Experimental sequence<\/p>\n
Each interrogation beginsonly after both ions have been subjected to Doppler cooling and prepared in the |g\u27e9 state with >99% probability using the conditional state preparation sequence described in ref.\u200953<\/a>. Both ions are interrogated with the HA\u2013HR sequence shown in Fig. 1d<\/a>, with the parameters \u03c4L\u00a0\u2248\u00a04\u2009ms, \u03c41\u00a0=\u00a0\u03c42\u00a0\u2248\u00a012\u2009ms, and a total Ramsey time given by TR\u00a0=\u00a03(T\u00a0+\u00a0\u03c41\u00a0+\u00a0\u03c42). Comparisons of the two independent frequency references are carried out by measuring the parity alternately for \\(\\phi =\\pm \\frac{{\\rm{\\pi }}}{2}\\) and, after N interrogation cycles, steering \\(\\varPi \\left(\\frac{{\\rm{\\pi }}}{2}\\right)-\\varPi \\left(-\\frac{{\\rm{\\pi }}}{2}\\right)\\) to zero by updating \u03b4f. For the comparisons at the same 0.1\u2009mT magnetic field, as is the case for all results shown in Fig. 2<\/a>, the microwave drive frequencies f1 and f2 are identical for both ions and fixed throughout. For comparisons at different magnetic fields, as when evaluating the quadratic Zeeman coefficient \u03b1z, the microwave frequencies are necessarily different to compensate quadratic Zeeman shifts. The frequency difference \u03b4f is set by an acousto-optic modulator (AOM), labelled AOM 2a in Fig. 1a<\/a>. Two auxiliary measurements are interleaved with the comparison servo: Rabi spectroscopy of the \\(| g\\rangle \\leftrightarrow | 8\\rangle \\) optical transition to ensure the 848\u2009nm clock laser is near resonance and measurement of the \\({|}^{3}{D}_{1},6,\\pm 1\\rangle \\) Zeeman splitting for each ion using microwave spectroscopy. The Zeeman splitting is used to infer the magnetic field, and feedback is applied to shim coils to compensate for any slow field drift, with further details in the\u00a0Supplementary Information<\/a>. HA\u2013HR interrogation time is typically 87% of the total duty cycle, including the auxiliary measurements.<\/p>\n Magnetic field stability and servo<\/p>\n The short-term stability of the magnetic field is shown in Extended Data Fig. 2<\/a> (blue points). This was evaluated from a separate experiment measuring a Zeeman splitting with fast servo attack time. Over the range of Ramsey times of interest in this work (about 1\u201350\u2009s), the magnetic field noise shows flicker instability of approximately 7\u2009nT. During comparison measurements, the magnetic field is steered to the set point of B0\u00a0=\u00a00.1\u2009mT by a compensation coil with a typical servo attack time of tser\u00a0\u2248\u00a080\u2009s using interleaved measurements of the \\({|}^{3}{D}_{1},6,\\pm 1\\rangle \\) Zeeman splitting. The orange points in Extended Data Fig. 2<\/a> show the field instability inferred from the in-loop measurements of the \\({|}^{3}{D}_{1},6,\\pm 1\\rangle \\) Zeeman splitting, which averages down \u221d(tser\/\u03c4) for reasons explained in the\u00a0Supplementary Information<\/a>. The true instability of the magnetic field is limited by both projection noise and deadtime \\(\\propto {({t}_{{\\rm{ser}}}\/\\tau )}^{1\/2}\\), which was verified in a test run using interleaved out-of-loop measurements of the \\({|}^{3}{D}_{1},8,\\pm 1\\rangle \\) Zeeman splitting (Extended Data Fig. 2<\/a>, green points). This out-of-loop stability is taken as an upper bound on the magnetic-field instability with compensation engaged. The magnetic field contributes uncertainty through the quadratic Zeeman shift of \\(\\frac{\\delta {\\nu }_{\\mathrm{QZ}}}{{\\nu }_{0}}=2{\\alpha }_{z}{B}_{0}\\delta B(\\tau )\\approx 2\\times 1{0}^{-19}\\times {(\\tau \/{\\rm{s}})}^{-1\/2}\\) beyond the servo attack time.<\/p>\n Detection and background gas collisions<\/p>\n To detect background collisions as much as possible and ensure the ions are sufficiently re-cooled before the next experiment cycle, we use the detection sequence shown in Extended Data Fig. 3a<\/a> at the end of every Ramsey experiment. The intervals di represent adaptive Bayesian state detection54<\/a> by 646\u2009nm fluorescence. If the ion is detected dark in the initial detection d0 after the Ramsey sequence, which is about 50% of events, then three attempts are made to shelve on the 848\u2009nm clock transition and detect. The variable 848\u2009nm clock pulse areas are tailored to maximize population transfer even for higher thermally occupied vibrational n states. The probability of at least one successful shelving is >99.5% for the expected thermal distribution accounting for ion heating during the Ramsey dark time TR. If the ion is not detected bright on any of d0…3, then it is assumed that a collision has occurred with sufficient energy transfer to significantly reduce either (1) the coupling on the 848\u2009nm transition or (2) the 646\u2009nm fluorescence rate. During clock comparison servo operation, if a collision is detected in this way on either Lu-1 or Lu-2, that interrogation cycle is considered invalid and repeated.<\/p>\n If a collision is detected (d0…3 all dark), the additional ending sequence shown in Extended Data Fig. 3a<\/a> is applied. First, a repump pulse (350\u2009nm, 622\u2009nm and 895\u2009nm) and another detection d4 is attempted. If bright, this is attributed to a lower-energy-transfer collision of type 1: sufficient to reduce efficiency of shelving on the clock transition but not detection. If dark, then the collision is of type 2: energetic enough to disrupt detection. In either case, a cycle of cooling and fluorescence detection with a high threshold, bi, is repeated until the ion is confirmed bright to ensure the ion is effectively cooled for the next experiment cycle. Extended Data Fig. 3b<\/a> shows the probability of outcomes (1) and (2) as a function of TR from which we extract rates \u0393(1) and \u0393(2). To estimate the detectable collision rate \u0393, we assume \\({\\varGamma }^{(1)}=\\frac{1-{P}_{{\\rm{d}}}}{2}\\varGamma \\) and \u0393(2)\u00a0=\u00a0Pd\u0393, where Pd is the probability a collision interferes with bright detection, resulting in outcome (2). Of the lower energy transfer collisions at rate (1\u00a0\u2212\u00a0Pd)\u0393, we assume half are not detected because the ion was in 3D1 at the end of the Ramsey experiment and half are detected as outcome (1). We infer the collision rates \u0393\u00a0=\u00a01.9\u00a0\u00d7\u00a010\u22123\u2009s\u22121 for Lu-1 and 5.9\u00a0\u00d7\u00a010\u22123\u2009s\u22121 for Lu-2.<\/p>\n After the comparison measurement campaign, additional experiments were performed to estimate the collision rate without the complications introduced by Ramsey spectroscopy. The ions were prepared in the 1S0 ground state, and after a delay of 5\u2009s, an attempt was made to reshelve and detect using a similar sequence as shown in Extended Data Fig. 3a<\/a>. From the probability that all of multiple shelving attempts on the clock transition fail, we directly infer the detectable collision rates \u0393\u00a0=\u00a02.9(3)\u00a0\u00d7\u00a010\u22123\u2009s\u22121 for Lu-1 and 6.5(6)\u00a0\u00d7\u00a010\u22123\u2009s\u22121 for Lu-2, in reasonable agreement with the rates inferred from the measurement campaign data.<\/p>\n Our extensive analysis of collision shifts in ion-based optical clocks is given in ref.\u200924<\/a>. The collision shift is evaluated with respect to a Langevin collision rate \u0393L, defined as the rate of collisions below a critical impact parameter that result in an inward-spiralling trajectory. We calculate the overestimate of the Langevin collision rate using eq.\u200931 of ref.\u200924<\/a> and with the cutoff velocity given by the 50% threshold of the Ramsey suppression factor as defined in section 2 of ref.\u200924<\/a>. Using \u03c9\u00a0\u2248\u00a01\u2009MHz for Lu-1 and \u03c9\u00a0\u2248\u00a0500\u2009kHz for Lu-2 yields \u0393L\u00a0=\u00a01.3\u00a0\u00d7\u00a010\u22123\u2009s\u22121 and 2.1\u00a0\u00d7\u00a010\u22123\u2009s\u22121 for the respective Langevin rates, corresponding to 3.5\u2009nPa and 5.8\u2009nPa background pressures of molecular hydrogen at 300\u2009K and consistent with pressure gauge readings. An inverted magnetron gauge on Lu-1 reads 4\u2009nPa, and the ion pumps on both chambers read \u2018low pressure\u2019 (<13\u2009nPa) at the limit of the sensitivity of the ion pump controllers. With \u03ba\u00a0=\u00a01 in eq.\u200960 of ref.\u200924<\/a>, we estimate collision shift bounds of 3.9\u00a0\u00d7\u00a010\u221220 and 3.2\u00a0\u00d7\u00a010\u221220 for Lu-1 and Lu-2, respectively.<\/p>\n Blackbody radiation<\/p>\n The differential dynamic polarizability, \u0394\u03b1(\u03c9), for the 848-nm clock transition in 176Lu+ has been well characterized39<\/a>,40<\/a>. The BBR shift is given by <\/p>\n $$\\frac{\\delta {\\nu }_{\\mathrm{BBR}}}{{\\nu }_{0}}=-4.90\\times 1{0}^{-19}{\\bar{T}}^{4}(1+1.77{\\bar{T}}^{2}),$$<\/p>\n \n (2)\n <\/p>\n where \\(\\overline{T}\\equiv T\/{T}_{0}\\) and T0\u00a0=\u00a0300\u2009K. Over the practical temperature range of 270\u2013330\u2009K, the uncertainty contribution from \u0394\u03b1(\u03c9) is well approximated by \\(9.8\\times 1{0}^{-20}{\\bar{T}}^{4}\\) (ref.\u200940<\/a>). A detailed temperature assessment is given in the\u00a0Supplementary Information<\/a> in which we bound the minimum and maximum temperatures to [299.8,\u00a0303.2]\u2009K for Lu-1 and [299.6,\u00a0301.9]\u2009K for Lu-2.<\/p>\n Trap secular frequencies<\/p>\n The Lu-1 trap used in previous work16<\/a> had an unusually high heating rate, which has been resolved after replacing the ion trap with one of the same design53<\/a>. Moreover, the helical resonators for both Lu-1 and Lu-2 were changed to reduce the RF-drive frequency \u03a9RF and obtain higher trap confinement without increased RF power. For 0.25\u2009W of RF power at \u03a9RF\u00a0=\u00a02\u03c0\u00a0\u00d7\u00a09.4\u2009MHz and 11.2\u2009MHz, we obtain secular trapping frequencies of (207, 1,063, 1,134)\u2009kHz and (138, 492, 543)\u2009kHz for Lu-1 and Lu-2, respectively, in which the weakest confinement corresponds to the axial direction. The combined effect of reduced heating rates and the smaller Lamb\u2013Dicke parameters resulting from the increased radial confinement means ion heating is no longer a significant limitation to the interrogation time. The simulated contrast loss shown in Fig. 2a<\/a> is evaluated by averaging the final optical pulses over a thermal distribution of radial mode occupation accounting for the heating during the Ramsey time.<\/p>\n Thermal second-order Doppler<\/p>\n Thermal motion gives rise to an SODS given by <\/p>\n $$\\frac{\\delta {\\nu }_{\\mathrm{SODS}}}{{\\nu }_{0}}=-\\frac{\\langle {v}^{2}\\rangle }{2{c}^{2}}=-\\frac{{k}_{{\\rm{B}}}}{2m{c}^{2}}\\sum _{i}({\\kappa }_{i}{T}_{i}),$$<\/p>\n \n (3)\n <\/p>\n where Ti are the average temperatures during the Ramsey interrogation for the axial (i\u00a0=\u00a01) and radial (i\u00a0=\u00a02,\u00a03) principal axes, respectively, and kB is Boltzmann\u2019s constant. The geometry factor \u03ba2\u00a0=\u00a0\u03ba3\u00a0\u2248\u00a02 accounts for equal contributions from the secular motion and intrinsic micromotion for the two radial principal axes, which is valid when the confinement is predominantly ponderomotive55<\/a>. Axial confinement is predominantly due to the static potential, and so only secular motion contributes (\u03ba1\u00a0\u2248\u00a01).<\/p>\n The initial temperatures and heating rates are determined by measuring the temperature immediately after Doppler cooling and after a fixed delay time. The axial temperature is measured by fitting the thermal dephasing of Rabi flopping on the 804\u2009nm E2 transition. Although the 804\u2009nm laser wave vector is at 45\u00b0 with respect to the ion trap axis and has a projection onto all principal axes, the radial modes have sufficiently low thermal occupation \\(\\overline{n}\\) and a small Lamb\u2013Dicke parameter \u03b7 that they contribute negligibly to the thermal dephasing. From the weighted mean of three measurements over the course of the campaign, we find the initial temperatures T1,0\u00a0=\u00a0190.5(5.9)\u2009\u03bcK for Lu-1 and 189.6(4.7)\u2009\u03bcK for Lu-2 and heating rates \\(\\frac{{\\rm{d}}{T}_{1}}{{\\rm{d}}t}=73(14)\\,\\mathrm{\\mu K}\\,{{\\rm{s}}}^{-1}\\) for Lu-1 and 119(17)\u2009\u03bcK\u2009s\u22121 for Lu-2 (Extended Data Fig. 5<\/a>). The average temperature during the Ramsey experiment for each mode is given by \\({T}_{i}={T}_{i,0}+\\frac{{\\rm{d}}{T}_{i}}{{\\rm{d}}t}\\frac{{T}_{{\\rm{R}}}}{2}\\).<\/p>\n The radial temperatures are measured by spectroscopy on the secular-motion sidebands using the 804\u2009nm E2 transition. The average thermal occupation \\(\\overline{n}\\) of one of the radial modes is extracted from the ratio of population transferred on the red and blue sidebands56<\/a>: <\/p>\n $$\\frac{{P}_{\\mathrm{red}}(t)}{{P}_{\\mathrm{blue}}(t)}=\\frac{\\bar{n}}{\\bar{n}+1}.$$<\/p>\n \n (4)\n <\/p>\n For both traps, the initial \\(\\overline{n}\\) measurements on both radial modes are consistent with the temperature 150(30)\u2009\u03bcK, which is also consistent with the initial axial temperatures to within statistical uncertainty and corresponds to approximately three times the Doppler cooling limit. The sideband ratio method is most sensitive for low \\(\\overline{n}\\), so for measuring the radial heating rates, we first apply Zeeman-degenerate Raman sideband cooling57<\/a> to prepare in the motional ground state and then measure the sideband ratio after some delay. Extended Data Fig. 6<\/a> shows the results of radial heating measurements in both chambers, which yield heating rates \\(\\left[\\frac{{\\rm{d}}{T}_{2}}{{\\rm{d}}t},\\frac{{\\rm{d}}{T}_{3}}{{\\rm{d}}t}\\right]=[110(15),29.1(2.9)]\\,\\mathrm{\\mu K}\\,{{\\rm{s}}}^{-1}\\) for Lu-1 and [85.5(5.5),\u00a070.4(5.5)]\u2009\u03bcK\u2009s\u22121 for Lu-2.<\/p>\n Equation (3<\/a>) is the SODS for a thermal state in the high-temperature limit, \\({\\bar{n}}_{i}\\approx \\frac{{k}_{{\\rm{B}}}{T}_{i}}{\\hbar {\\omega }_{i}}\\), and does not account for the zero-point fluctuations for the quantum ground state for both secular and micromotion58<\/a>,59<\/a>, which contribute additional SODS of \\(-\\sum _{i}{\\kappa }_{i}\\left(\\frac{\\hbar {\\omega }_{i}}{4m{c}^{2}}\\right)\\). This term contributes only a 6% correction to the total SODS for Lu-1 but has been included for completeness.<\/p>\n We estimate the total SODS for a Ramsey time TR as follows: <\/p>\n $$\\begin{array}{c}(\\mathrm{Lu}\\text{-1})\\,\\frac{\\delta {\\nu }_{\\mathrm{SODS}}}{{\\nu }_{0}}=-2.4(2)\\times 1{0}^{-19}+{T}_{{\\rm{R}}}[-4.6(4)\\times 1{0}^{-20}\\,{{\\rm{s}}}^{-1}],\\\\ (\\mathrm{Lu}\\text{-2})\\frac{\\delta {\\nu }_{\\mathrm{SODS}}}{{\\nu }_{0}}=-2.4(2)\\times 1{0}^{-19}+{T}_{{\\rm{R}}}[-5.7(3)\\times 1{0}^{-20}\\,{{\\rm{s}}}^{-1}].\\end{array}$$<\/p>\n Excess micromotion<\/p>\n The micromotion shift has two components: an SODS due to oscillatory motion and the a.c. Stark shift due to the RF electric field. Owing to the very low differential polarizability of the clock transition, the a.c. Stark contribution is more than two orders of magnitude smaller than the SODS component and therefore negligible. The micromotion amplitude is evaluated using phase-modulated sideband spectroscopy on the 804-nm E2 clock transition22<\/a>. The excess micromotion (EMM) SODS is given by <\/p>\n $$\\frac{\\delta {\\nu }_{\\mathrm{EMM}}}{{\\nu }_{0}}=-{\\left(\\frac{{\\varOmega }_{\\mathrm{RF}}}{2c{k}_{804}}\\right)}^{2}|\\beta {|}^{2},$$<\/p>\n \n (5)\n <\/p>\n where k804 is the wavenumber of the 804-nm clock transition, and \u03b2 is the modulation depth. As phase-modulated sideband spectroscopy is sensitive to the phase of the micromotion, we distinguish between two quadrature components \u03b2\u00a0=\u00a0\u03b2m\u00a0+\u00a0i\u03b2p, where \u03b2m is the excess micromotion due to displacement by stray fields, and \u03b2p is due to a phase shift between RF electrodes, which we are not able to compensate. Intrinsic micromotion is accounted for separately with the thermal SODS.<\/p>\n Following the procedure described in ref.\u200922<\/a>, we measure the modulation depth in three orthogonal directions to evaluate \\({\\beta }_{{\\rm{m}}}^{2}=\\sum _{i}{\\beta }_{{\\rm{m}},i}^{2}\\). The results of all measurements of \u03b2m,i taken for both ion traps are shown in Extended Data Fig. 4a,b<\/a>, in which the open circles are measured before compensating for the d.c. stray field and closed circles immediately after. To estimate the EMM shift for each comparison measurement, we linearly interpolate the \u03b2m,i measurements, assuming linear growth of the d.c. stray field between measurements. Extended Data Fig. 4c<\/a> shows the EMM shifts evaluated from the time average of \\({\\beta }_{{\\rm{m}}}^{2}\\) over each comparison interval. We estimate an average relative shift of \u22121.6(1.3)\u00a0\u00d7\u00a010\u221220 for Lu-1 and \u22121.4(0.4)\u00a0\u00d7\u00a010\u221220 for Lu-2.<\/p>\n As reported in ref.\u200922<\/a>, \u03b2p is negligibly small for Lu-1. For Lu-2, the three orthogonal components of \u03b2p are measured to be [2.22(13),\u00a01.43(12),\u00a01.92(13)]\u00a0\u00d7\u00a010\u22122, corresponding to an EMM shift of \u22122.42(19)\u00a0\u00d7\u00a010\u221219. \u03b2p was remeasured at the end of the campaign and found to be in statistical agreement with the original measurement.<\/p>\n a.c. Zeeman, RF<\/p>\n Time-varying magnetic fields give rise to an a.c. Zeeman shift, with the dominant contributions coming from currents in the electrodes driven by the RF-trapping potential. The contribution depends on the component of the RF magnetic field perpendicular to the applied d.c. field60<\/a>. The clock shift given by <\/p>\n $$\\delta {\\nu }_{\\mathrm{RF}}={\\alpha }_{\\perp }\\langle {B}_{\\perp }^{2}\\rangle ,$$<\/p>\n \n (6)\n <\/p>\n where \\(\\sqrt{\\langle {B}_{\\perp }^{2}\\rangle }\\) is the root-mean-square RF magnetic field amplitude perpendicular to the d.c. field and \u03b1\u22a5\u00a0=\u00a00.20\u2009mHz\u2009\u03bcT\u22122 is the sensitivity coefficient after 3D1 hyperfine averaging60<\/a>.<\/p>\n Although in previous work B\u22a5 was measured using an Autler\u2013Townes splitting on the Ba+ clock transition41<\/a>, here we use an Autler\u2013Townes splitting on the Lu+3D2 |6,\u20090\u27e9 to |5,\u20090\u27e9 microwave transition. The 3D2 F\u00a0=\u00a05 Land\u00e9 g-factor, g5\u00a0=\u00a0\u22120.38575750(19) (ref.\u200953<\/a>), is the largest of all 3D1 and 3D2 F manifolds, and for \u03a9RF\u00a0\u2248\u00a02\u03c0\u00a0\u00d7\u00a010\u2009MHz, the linear Zeeman splitting can be brought into resonance with \u03a9RF at experimentally accessible magnetic fields of B\u00a0\u2272\u00a02\u2009mT.<\/p>\n The experimental measurement procedure is as follows. Neodymium permanent magnets are used to bias the magnetic field to the required 1.7\u2009mT for Lu-1 and 2.0\u2009mT for Lu-2. Three pairs of coils are used to fine-tune the magnetic field amplitude and optimize the orientation for the \\({|}^{3}{D}_{1},7,0\\rangle \\) state preparation by aligning to a \u03c0-polarized 646-nm laser. This ensures the d.c. magnetic field is aligned in the same direction as during the comparison measurements. From \\({|}^{3}{D}_{1},7,0\\rangle \\), we transfer sequentially to the \\({|}^{3}{D}_{1},8,0\\rangle \\) state using a microwave \u03c0 pulse and then the \\({|}^{3}{D}_{2},6,0\\rangle \\) state using an 804-nm\u2013848-nm Raman pair. We observe the Autler\u2013Townes splitting by interrogating the \\({|}^{3}{D}_{2}\\,| 6,0\\rangle \\) to |5,\u20090\u27e9 microwave transition with a \u03c0 pulse for a range of detunings, after which the remaining \\({|}^{3}{D}_{2},6,0\\rangle \\) population is reshelved to \\({|}^{3}{D}_{1}\\rangle \\) by Raman transfer for detection.<\/p>\n Following the treatment in ref.\u200960<\/a>, the effective Hamiltonian in the rotating wave approximation describing the microwave interrogation is <\/p>\n $$H=\\frac{\\hbar }{2}\\left(\\begin{array}{cccc}-2{\\varDelta }_{{\\rm{M}}} & {\\varOmega }_{M} & 0 & 0\\\\ {\\varOmega }_{M} & 0 & {\\varOmega }_{B} & {\\varOmega }_{B}\\\\ 0 & {\\varOmega }_{B} & 2{\\varDelta }_{1} & 0\\\\ 0 & {\\varOmega }_{B} & 0 & -2{\\varDelta }_{-1}\\end{array}\\right),$$<\/p>\n \n (7)\n <\/p>\n where \u03a9M is the microwave coupling, \u0394M is the microwave detuning, \\({\\varOmega }_{B}=\\frac{1}{\\hbar }\\frac{\\sqrt{15}}{2}{g}_{5}{\\mu }_{B}{B}_{\\perp }\\) is the coupling between \\({|}^{3}{D}_{2},5,0\\rangle \\) and |5,\u2009\u00b11\u27e9 states due to the transverse oscillating magnetic field at frequency \u03a9RF, and <\/p>\n $${\\varDelta }_{\\pm 1}={\\varOmega }_{{\\rm{R}}{\\rm{F}}}\\pm {\\omega }_{\\pm 1}={\\varOmega }_{{\\rm{R}}{\\rm{F}}}+{g}_{5}{\\mu }_{B}B\/\\hbar \\pm 2\\pi \\times \\varDelta \\alpha {B}^{2}$$<\/p>\n \n (8)\n <\/p>\n is the detuning of the RF field from the |5,\u20090\u27e9 to |5,\u2009\u00b11\u27e9 Zeeman splittings, \u03c9\u00b11. Degeneracy of \u03c9\u00b11 is lifted by the differential quadratic Zeeman shift for which \u0394\u03b1\u00a0\u2248\u00a01.23\u2009kHz\u2009mT\u22122.<\/p>\n Extended Data Fig. 7a<\/a> shows the simulated spectrum of the Autler\u2013Townes splittings for the operating parameters of Lu-1 and \u03a9B\u00a0=2\u03c0\u00a0\u00d7 3\u2009kHz. For both Lu-1 and Lu-2, several spectra were taken at a fixed magnetic field in the vicinity of the Autler\u2013Townes splitting, of which three from Lu-1 are shown in Extended Data Fig. 7b\u2013d<\/a>. We fit the measured data to simulated spectra to obtain \u03a9B. Taking the mean of the \u03a9B fit values, with the standard deviation of fit values as the uncertainty, we find \\(\\sqrt{\\langle {B}_{\\perp }^{2}\\rangle }\\) to be 0.220(5)\u2009\u03bcT and 0.303(4)\u2009\u03bcT for Lu-1 and Lu-2, respectively. The RF-drive voltage to the trap is monitored over the campaign and for both Lu-1 and Lu-2 varied by less than 0.5%.<\/p>\n a.c. Zeeman, microwave<\/p>\n When applying the microwave fields during the Ramsey sequence for hyperfine averaging, there is a probe-induced shift because of the \u03c3\u00b1 polarization components off-resonantly coupling to m\u00a0=\u00a0\u00b11 Zeeman states. Evaluation of this shift is discussed in detail in the supplementary material of ref.\u200916<\/a>. The total microwave a.c. Zeeman shift is given by <\/p>\n $$\\delta {\\nu }_{\\mu }=\\frac{({\\varDelta }_{1,7}+{\\varDelta }_{1,8}){\\tau }_{1}+({\\varDelta }_{2,6}+{\\varDelta }_{2,7}){\\tau }_{2}}{{T}_{{\\rm{R}}}},$$<\/p>\n \n (9)\n <\/p>\n where \u0394k,F is the shift of |F,\u2009m\u2009=\u20090\u27e9 when microwave coupling \u03a9k is on. The microwave field polarizations are characterized by measuring the \u03c0 times \u03c4kq for the field k and polarization q, at fixed microwave power, from which the shifts \u0394k,F are evaluated as<\/p>\n $$\\begin{array}{c}{\\varDelta }_{1,8}=-\\frac{7}{9}\\frac{{\\varOmega }_{1}^{2}}{4{\\omega }_{7}}\\left[{\\left(\\frac{{\\tau }_{10}}{{\\tau }_{1+}}\\right)}^{2}-{\\left(\\frac{{\\tau }_{10}}{{\\tau }_{1-}}\\right)}^{2}\\right],\\\\ {\\varDelta }_{1,7}=\\frac{{\\varOmega }_{1}^{2}}{4{\\omega }_{8}}\\left[{\\left(\\frac{{\\tau }_{10}}{{\\tau }_{1+}}\\right)}^{2}-{\\left(\\frac{{\\tau }_{10}}{{\\tau }_{1-}}\\right)}^{2}\\right],\\\\ {\\varDelta }_{2,7}=\\frac{{\\varOmega }_{2}^{2}}{4{\\omega }_{6}}\\left[{\\left(\\frac{{\\tau }_{20}}{{\\tau }_{2+}}\\right)}^{2}-{\\left(\\frac{{\\tau }_{20}}{{\\tau }_{2-}}\\right)}^{2}\\right],\\\\ {\\varDelta }_{2,6}=\\frac{4}{3}\\frac{{\\varOmega }_{2}^{2}}{4{\\omega }_{7}}\\left[{\\left(\\frac{{\\tau }_{20}}{{\\tau }_{2+}}\\right)}^{2}-{\\left(\\frac{{\\tau }_{20}}{{\\tau }_{2-}}\\right)}^{2}\\right],\\end{array}$$<\/p>\n where \u03c9F\u00a0>\u00a00 is the Zeeman splitting for the hyperfine level F.<\/p>\n This shift is suppressed by ensuring \u03c4k\u00a0\u226a\u00a0TR, balancing the circular polarization components (\u03c4k+\u00a0\u2248\u00a0\u03c4k\u2212), and maximizing the \u03c0 coupling (\u03c4k0\u00a0\u226a\u00a0\u03c4k\u00b1). All microwave horns are mounted on rotation mounts and are tuned, through a combination of translation and rotation, to the condition \u03c4k+\u00a0\u2248\u00a0\u03c4k\u2212 as much as possible at the start of the campaign. All \u03c4kq were characterized at three points in the campaign, and the evaluated clock shifts were consistently below 1\u00a0\u00d7\u00a010\u221220.<\/p>\n Quadrupole shift<\/p>\n The hyperfine-averaged residual quadrupole moment has been reported as \\(\\widetilde{\\varTheta }=-2.54(0.25)\\times 1{0}^{-4}\\,e{a}_{0}^{2}\\) (ref.\u200942<\/a>), which has been independently confirmed by an alternate method23<\/a>. The residual quadrupole shift, \\(\\delta {\\widetilde{\\nu }}_{Q}\\), is evaluated by measuring the quadrupole shift, \u03b4\u03bdQ,1, on the 3D1 |7,\u20090\u27e9 to |8,\u20090\u27e9 microwave clock transition, which are related by <\/p>\n $$\\delta {\\widetilde{\\nu }}_{Q}=\\frac{\\widetilde{\\varTheta }}{\\frac{8}{5}\\varTheta {(}^{3}{D}_{1})}\\delta {\\nu }_{Q,1},$$<\/p>\n \n (10)\n <\/p>\n where \\(\\varTheta {(}^{3}{D}_{1})=0.63862(74)\\,e{a}_{0}^{2}\\) (ref.\u200961<\/a>). We measure the microwave transition frequencies for both traps by microwave Ramsey spectroscopy with a 10\u2009s interrogation time. The quadrupole shifts \u03b4\u03bdQ,1 are inferred using the unperturbed microwave frequencies23<\/a> and accounting for the second-order Zeeman shifts. The residual quadrupole shifts on Lu-1 and Lu-2 are evaluated to be 3.3(3)\u00a0\u00d7\u00a010\u221220 and 1.35(14)\u00a0\u00d7\u00a010\u221219, respectively. Given the accuracy of the measured unperturbed microwave frequencies23<\/a>, the quadrupole shift may be easily suppressed further by setting the magnetic field angle so as to null the microwave quadrupole shifts as much as is required.<\/p>\n Microwave coupling errors<\/p>\n Incorrect microwave pulse areas due to uncertainty in the respective couplings give rise to timing errors in HA\u2013HR spectroscopy. As derived in the supplementary material of ref.\u200916<\/a>, the shifts due to microwave coupling errors are bound by <\/p>\n $$\\delta {\\nu }_{k}\\approx -\\frac{1}{{T}_{{\\rm{R}}}}\\left({\\left(\\frac{{\\varDelta }_{k}{\\tau }_{k}}{\\pi }\\right)}^{2}+{\\left(\\frac{\\pi {q}_{k}}{2}\\right)}^{2}\\right),$$<\/p>\n \n (11)\n <\/p>\n where \u0394k is the microwave detuning, \u03c4k is the microwave \u03c0 pulse duration, and qk is the fractional error in the microwave coupling for the k\u00a0=\u00a01,\u00a02 microwave transitions.<\/p>\n These shifts are suppressed for Ramsey times TR\u00a0\u226b\u00a0\u03c4L,\u00a0\u03c4k, as is the case for the comparison experiments reported here. The magnetic field was actively steered to a fixed value in both chambers so the microwave detunings \u0394k are stable and known precisely from the quadrupole shift assessment. The microwave couplings were measured at the beginning, middle and end of the campaign with about 0.2% measurement precision and found to deviate by less than 0.5%, with the exception of \u03a92 on Lu-2, which drifted by 1.5% from the beginning to the end of the campaign. Extended Data Fig. 8<\/a>, for example, shows this shift evaluated for the worst case of q2\u00a0=\u00a00.015 as a function of the detuning \u03942. Although the \u0394k are precisely known at the <10\u2009mHz level, which, in principle, allows for a more accurate estimation of this shift, we take a conservative approach and use the full width of the bounding envelope, given by the second term in equation (11<\/a>), as the uncertainty.<\/p>\n a.c. Stark<\/p>\n For an optical \u03c0 time of \\({\\tau }_{{\\rm{L}}}=\\frac{{\\rm{\\pi }}}{{\\varOmega }_{{\\rm{L}}}}=4\\,\\mathrm{ms}\\), the a.c. Stark shift on the |g\u27e9 to |8\u27e9 optical transition when interrogating with the 848-nm laser is approximately \u0394S\u00a0=\u00a02\u03c0\u00a0\u00d7\u00a025\u2009Hz. In the absence of additional effects, both the HR37<\/a> sequence, Extended Data Fig. 9a<\/a>, which was used in previous work16<\/a>, and the time-reversed HR sequence, Extended Data Fig. 9b<\/a>, which is used in this work, equivalently suppress the resulting clock shift to \\(\\frac{2}{{\\rm{\\pi }}{T}_{{\\rm{R}}}}{\\left(\\frac{\\varDelta }{{\\varOmega }_{{\\rm{L}}}}\\right)}^{3}\\), in Hz, where \u0394\u00a0=\u00a0\u0394SP\u00a0\u2212\u00a0\u0394S is the error in the frequency step, \u0394SP, applied during the optical interrogation pulses. To set the laser frequency step, \u0394SP, the a.c. Stark shift was measured to better than 1% for both Lu-1 and Lu-2 at the start of the campaign by interleaved self-comparison of Rabi and HR spectroscopy. The value of the \u0394SP was fixed throughout the campaign, and we bound the uncertainty in \u0394 to 2% of |\u0394S| for Lu-1 and 5% of |\u0394S| for Lu-2 based on the maximum observed variation in the optical couplings, which were measured at the beginning, middle and end of the campaign with measurement precision of about 0.2%.<\/p>\n When including the effects of ion heating, HR schemes present a weak linear dependence on \u0394 (ref. 62<\/a>). This linear dependence is evaluated by simulation with the results for the typical operating parameters of Lu-2 shown in Extended Data Fig. 9c<\/a>. It is noted that the linear slope has the opposite sign for HR and HR reverse schemes, and that a hybrid error signal constructed by averaging the two schemes cancels the thermal effect, restoring the cubic dependence on \u0394. The idea of combining sequences to generate a more robust error signal is similar to modified HR63<\/a> and generalized HR64<\/a>, but unlike those schemes for which ion heating results in either an offset or increased sensitivity compared with the basic HR scheme62<\/a>, this simple hybrid combination of HR and HR reverse heavily suppresses the linear dependence on ion heating, at least in the regime \u03c4L\u00a0\u226a\u00a0TR, in which heating during the optical pulse is negligible compared with during the dark time TR.<\/p>\n For the experiments reported here, the shift was already sufficiently small that only HR reverse was used. From the measured radial heating rates, uncertainty in \u0394 and the linear dependence determined from simulation, we estimate a total uncertainty from the a.c. Stark effect to be 3.8\u00a0\u00d7\u00a010\u221221 for Lu-1 and 2.5\u00a0\u00d7\u00a010\u221220 for Lu-2 for TR\u00a0=\u00a05\u2009s.<\/p>\n a.c. Quadrupole, RF<\/p>\n The oscillating RF quadrupole field couples off-resonantly to hyperfine transitions giving rise to an a.c. quadrupole shift that is not cancelled by hyperfine averaging65<\/a>. We estimate this shift by assuming an ideal linear Paul trap with RF potential of the form \u03a6(x,\u00a0y,\u00a0z)\u00a0=\u00a0\u03f5(x2\u00a0\u2212\u00a0y2) in the principal-axis frame, where the electric field gradient \\({\\epsilon }=\\frac{m{\\varOmega }_{\\mathrm{RF}}{\\omega }_{{\\rm{r}}}}{e\\sqrt{2}}\\) is determined by the radial pseudo-potential confinement frequency \u03c9r\u00a0\u2248\u00a02\u03c0\u00a0\u00d7\u00a01,100\u2009kHz for Lu-1 and 510\u2009kHz for Lu-2. The magnetic field in both chambers is aligned to approximately 33(3)\u00b0 with respect to the ion trap axis (z). Under these conditions, we evaluate65<\/a> a relative clock shift of the HA frequency of \u22122.0\u00a0\u00d7\u00a010\u221221 for Lu-1 and \u22125.7\u00a0\u00d7\u00a010\u221222 for Lu-2.<\/p>\n Differential AOM chirp<\/p>\n As shown schematically in Fig. 1a<\/a>, the differential path length is actively stabilized with reference to retroreflecting mirrors below the respective experimental chambers that are near the table surface and approximately 20\u2009cm below the trapped ions. The differential phase stability is characterized by out-of-loop measurement of the optical phase with the clock beams directed to a common beam splitter instead of the ions, requiring approximately 2\u2009m of additional unstabilized optical path length. When simultaneously switching on AOMs 1a and 2a, as for the Ramsey pulses in the comparison interrogation sequence, a differential phase chirp is induced by the lock circuitry as shown in Extended Data Fig. 10a<\/a>. This is well modelled by a damped harmonic oscillation with the fit parameters shown in Extended Data Fig. 10a<\/a>. By straightforward extension of the analysis given in ref.\u200966<\/a>, we evaluate a \u22124.8(8)\u00a0\u00d7\u00a010\u221221 systematic shift to the difference frequency between Lu-1 and Lu-2 for TR\u00a0=\u00a05\u2009s.<\/p>\n First-order Doppler<\/p>\n We consider two sources of first-order Doppler shift (FODS): drift of the ion with respect to the trap along the clock probe (vertical) direction and differential drifts of the optical phase.<\/p>\n The measured micromotion modulation depth \u03b2m,3 can be directly related to the displacement of the ion from the trap RF null along the clock interrogation direction22<\/a>. For small modulation depth, the relationship is linear with scale factors of approximately 0.77\u2009\u03bcm\/\u03b2m,3 for Lu-1 and 2.0\u2009\u03bcm\/\u03b2m,3 for Lu-2. Extended Data Fig. 4d<\/a> shows the evaluated FODS shift for each comparison measurement, where the average velocity is estimated from the inferred displacement between the EMM measurements. The weighted FODS for the entire campaign is estimated to be 2.4(3.5)\u00a0\u00d7\u00a010\u221223 for Lu-1, \u22120.3(2.0)\u00a0\u00d7\u00a010\u221222 for Lu-2 and 0.5(2.0)\u00a0\u00d7\u00a010\u221222 for the difference.<\/p>\n Using the same measurement setup as for AOM chirp, we measure the out-of-loop differential phase for long durations for which we observed an instability asymptote of 3.0\u00a0\u00d7\u00a010\u221216\u2009(\u03c4\/s)\u22121 as shown in Extended Data Fig. 10b<\/a> (blue). For typical servo operation at 87% interrogation duty cycle, we evaluate a contribution of 2.0\u00a0\u00d7\u00a010\u221217\u2009(\u03c4\/s)\u22121\/2 (orange). This is well below the comparison instability on all timescales. When operating as a single clock, the path length is actively stabilized from the laser to reference mirror instead of differentially between chambers, but can be expected to contribute at the same level or lower because of the shorter unstabilized path.<\/p>\n RF synthesis<\/p>\n The RF synthesizers used for the AOMs are based on the AD9912 chip and have approximately 7\u2009\u03bcHz resolution. This contributes about 1\u00a0\u00d7\u00a010\u221220 fractional uncertainty, half of the minimum step size, to the comparison servo. The synthesized microwave frequencies, f1 and f2, were identical and generated from the same sources for Lu-1 and Lu-2, set to the nearest 1\u2009mHz. All synthesizers are referenced to a common hydrogen maser, which is calibrated to \u22722\u00a0\u00d7\u00a010\u221215 at the start of the measurement campaign67<\/a>. For the comparison measurement, any error in maser accuracy is common mode and contributes no uncertainty to the difference.<\/p>\n However, we note that when realizing a Lu+ standard using HA\u2013HR spectroscopy as \\({\\nu }_{0}={f}_{{\\rm{L}}}+\\frac{1}{3}(2{f}_{1}+{f}_{2})\\), the accuracy of the synthesized microwaves fk is also not a limitation to optical clock accuracy because the servo action steers the laser to compensate any errors in the microwave fields such that this specific linear combination of the fields is the HA frequency. The HA\u2013HR servo is sensitive to the instability of both the laser and microwave fields at the same level, in Hz. In practice, this places a more relaxed requirement on the microwaves as compared with the laser; for example, for a laser with fractional instability \\(\\frac{\\delta {f}_{{\\rm{L}}}}{{f}_{{\\rm{L}}}}\\approx 1{0}^{-16}\\), the microwaves need only a fractional stability better than \\(\\frac{\\delta {f}_{k}}{{f}_{k}}\\lesssim \\frac{\\delta {f}_{{\\rm{L}}}}{{f}_{k}} \\sim 1{0}^{-12}\\) to not limit the clock instability.<\/p>\n Gravitational redshift<\/p>\n The differential redshift between the ions is given by \\(\\frac{g\\delta h}{{c}^{2}}\\), where g\u00a0\u2248\u00a09.776\u2009m\u2009s\u22122 is the local gravitational acceleration and \u03b4h\u00a0=\u00a0h1\u00a0\u2212\u00a0h2 is the height difference of the ions. We measure the heights of the ions relative to a laser levelling assembly fixed to the optical table, which is precisely aligned to each ion by the same imaging optics used for state detection. We determine the ion heights h1\u00a0=\u00a07.44(10)\u2009mm and h2\u00a0=\u00a011.41(10)\u2009mm. Measurements of the table levelling by a precision spirit level with 0.02\u2009mm\u2009m\u22121 resolution are limited by the table surface flatness but bound the table tilt to <0.2\u2009mm\u2009m\u22121. The ion traps are separated by 1.2\u2009m, and table levelling contributes the dominant uncertainty of 240\u2009\u03bcm. The height difference of the ions is evaluated to be \u03b4h\u00a0=\u00a0\u22123.97(27)\u2009mm, corresponding to a \u22124.31(29)\u00a0\u00d7\u00a010\u221219 differential shift.<\/p>\n Clock stability<\/p>\n The stability of a single-ion optical clock is quantum-projection-noise limited at the longest achievable interrogation time, which is set by either the optical coherence of the local oscillator or the atomic coherence, whichever is shorter. State-of-the-art cavity-stabilized lasers reach fractional instabilities in the low 10\u221217 at 1\u2009s (refs.\u200928<\/a>,33<\/a>) and have enabled demonstrated single-ion clock instabilities as low as 3.5\u00a0\u00d7\u00a010\u221216\u2009(\u03c4\/s)\u22121\/2 (ref.\u200912<\/a>). Although our laboratory lacks a local oscillator of this quality with which to challenge ion clock stability at this level directly, our comparison does probe the limits set by atomic coherence, allowing us to infer the technical limits to interrogation time independent of any clock laser. Comparison of two references by correlation spectroscopy carries a factor of \u221a2 higher instability than a comparison of independent clocks, and a factor of 2 higher than a single standalone clock, for equal interrogation time. Insofar as the present comparison is limited by uncorrelated magnetic field noise, these systems, when combined with state-of-the-art laser stabilization, would support a single-clock instability of about 2.4\u00a0\u00d7\u00a010\u221216\u2009(\u03c4\/s)\u22121\/2. Because the 3D1 lifetime imposes no practical limit, substantially longer interrogation times remain possible, contingent on overcoming magnetic field noise, ion trap heating and background gas collisions. Reducing the d.c. magnetic field lowers the linear sensitivity to field noise at the cost of smaller Zeeman splittings; nevertheless, operation at 25\u2009\u03bcT is achievable with modest changes to our laser configuration. By our assessment, this would permit interrogation times up to 30\u2009s in the present system. Even longer coherence times can be anticipated with the use of magnetic shielding and cryogenics, albeit at the cost of increased experimental complexity.<\/p>\n","protected":false},"excerpt":{"rendered":"Experimental sequence Each interrogation beginsonly after both ions have been subjected to Doppler cooling and prepared in the…\n","protected":false},"author":2,"featured_media":634062,"comment_status":"","ping_status":"","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[24],"tags":[13734,2026,61,60,2027,33193,248,82],"class_list":["post-634061","post","type-post","status-publish","format-standard","has-post-thumbnail","category-physics","tag-atomic-and-molecular-physics","tag-humanities-and-social-sciences","tag-ie","tag-ireland","tag-multidisciplinary","tag-optical-spectroscopy","tag-physics","tag-science"],"_links":{"self":[{"href":"https:\/\/www.newsbeep.com\/ie\/wp-json\/wp\/v2\/posts\/634061","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.newsbeep.com\/ie\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.newsbeep.com\/ie\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.newsbeep.com\/ie\/wp-json\/wp\/v2\/users\/2"}],"replies":[{"embeddable":true,"href":"https:\/\/www.newsbeep.com\/ie\/wp-json\/wp\/v2\/comments?post=634061"}],"version-history":[{"count":0,"href":"https:\/\/www.newsbeep.com\/ie\/wp-json\/wp\/v2\/posts\/634061\/revisions"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/www.newsbeep.com\/ie\/wp-json\/wp\/v2\/media\/634062"}],"wp:attachment":[{"href":"https:\/\/www.newsbeep.com\/ie\/wp-json\/wp\/v2\/media?parent=634061"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.newsbeep.com\/ie\/wp-json\/wp\/v2\/categories?post=634061"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.newsbeep.com\/ie\/wp-json\/wp\/v2\/tags?post=634061"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}