Quantum logic inspired techniques for spacetime-symmetry tests with (anti-)protons

Cosmological observations as well as theoretical approaches to physics beyond the Standard Model provide strong motivations for experimental tests of fundamental symmetries, such as CPT invariance. In this context, the availability of cold baryonic antimatter at CERN has opened an avenue for ultrahigh-precision comparisons of protons and antiprotons in Penning traps. This work discusses an experimental method inspired by quantum logic techniques that will improve particle localization and readout speed in such experiments. The method allows for sympathetic cooling of the (anti-)proton to its quantum-mechanical ground state as well as the readout of its spin alignment, replacing the commonly used continuous Stern-Gerlach effect. Both of these features are achieved through coupling to a laser-cooled `logic' ion co-trapped in a double-well potential. This technique will boost the measurement sampling rate and will thus provide results with lower statistical uncertainty, contributing to stringent searches for time dependent variations in the data. Such measurements ultimately yield extremely high sensitivities to CPT violating coefficients acting on baryons in the Standard-Model Extension, will allow the exploration of previously unmeasured types of symmetry violations, and will enable antimatter-based axion-like dark matter searches with improved mass resolution.

Cosmological observations as well as theoretical approaches to physics beyond the Standard Model provide strong motivations for experimental tests of fundamental symmetries, such as CPT invariance. In this context, the availability of cold baryonic antimatter at CERN has opened an avenue for ultrahigh-precision comparisons of protons and antiprotons in Penning traps. This work discusses an experimental method inspired by quantum logic techniques that will improve particle localization and readout speed in such experiments. The method allows for sympathetic cooling of the (anti-)proton to its quantum-mechanical ground state as well as the readout of its spin alignment, replacing the commonly used continuous Stern-Gerlach effect. Both of these features are achieved through coupling to a laser-cooled 'logic' ion co-trapped in a double-well potential. This technique will boost the measurement sampling rate and will thus provide results with lower statistical uncertainty, contributing to stringent searches for time dependent variations in the data. Such measurements ultimately yield extremely high sensitivities to CPT violating coefficients acting on baryons in the Standard-Model Extension, will allow the exploration of previously unmeasured types of symmetry violations, and will enable antimatter-based axion-like dark matter searches with improved mass resolution.

I. INTRODUCTION
The study of symmetries and their breakdown is closely intertwined with the development of the Standard Model (SM) and continues to be a cornerstone of experimental and theoretical explorations aiming to uncover more fundamental physics. While most discrete spacetime symmetries are known to be violated in weak interactions, there is no credible experimental evidence of departures from the combined transformation of charge conjugation (C), parity inversion (P) and time reversal (T) [1,2]. The theoretical underpinning of this result is provided by the CPT theorem [3], which states that a realistic Lorentzsymmetric quantum field theory must also be CPT invariant. This profound result entails an exact symmetry between the properties of matter-antimatter conjugates. However, the observable universe appears to exhibit an imbalance between baryons and antibaryons [4]. Although the SM contains the qualitative features for generating a baryon asymmetry, it struggles to explain the asymmetry quantitatively: the required amount of CP violation has thus far eluded solid experimental identification in various sectors of the SM [1,4], and the large Higgs mass can impede the necessary departures from thermal equilibrium in the early universe [1,5]. Plancksuppressed CPT breaking, on the other hand, may provide an alternative mechanism, as it can evade both of these conditions [6]. Moreover, a number of theoretical ideas involving physics beyond the SM can accommodate small deviations from CPT symmetry [7][8][9][10]. This situation has spurred the development of the Standard-Model Extension (SME) test framework for the theoretical description of such effects [11] and has led to intensified efforts to perform matter-antimatter comparisons in general and improved CPT tests with baryons in particular [2].
The rapid progress and interest in ultrahigh-precision studies with baryons is evident from a number of recent measurements involving the (anti-)proton. These include proton experiments performed by the BASE (Baryon-Antibaryon Symmetry Experiment) collaboration at Mainz [12,13] as well as measurements with both protons and antiprotons performed at CERN's Antiproton Decelerator (AD) by several collaborations. These efforts have yielded key measurements involving chargeto-mass ratios and magnetic moments of these particles.
One of these was performed by the TRAP collaboration at CERN's low-energy antiproton ring LEAR by cyclotron-frequency comparisons of antiprotons and protons [14], which was subsequently improved by using negatively-charged hydrogen ions [15]. Later, a more sensitive result for the proton-to-antiproton charge-tomass ratio (q/m) p /(q/m)p = 1 + 1(69) × 10 −12 was obtained from cyclotron-frequency comparisons of a single antiproton and a H − ion by using the BASE apparatus at CERN's AD [16]. BASE has also performed a 0.8 ppm [17] measurement of the antiproton magnetic moment improving on a previous measurement [18] from the ATRAP collaboration by a factor of six. Using single spin-flip resolution [19] and an advanced multi-trap arXiv:2106.06252v2 [hep-ph] 13 Jul 2021 method, BASE further improved this result with the most precise determination of the antiproton magnetic moment with a fractional precision of 1.5 ppb in units of the nuclear magneton [20]. This represents the first time that a measurement sensitivity for an antimatter property surpassed that of its matter conjugate. This work has also placed constraints on CPT-odd SME coefficients at the 10 −24 GeV level [20], which is comparable in sensitivity to the corresponding results for the electron [21,22], and outperforms limits extracted from muon g-factor measurements [23,24]. BASE-Mainz has subsequently continued to refine our knowledge of the proton magnetic moment at the sub-ppb level [13]. The BASE measurements have adapted techniques originally developed by Dehmelt for electrons and positrons [21,22,25], combined them with the doubletrap method [26], and pushed their limits for the extremely challenging case of the (anti-)proton. This has enabled the measurement of the proton and antiproton magnetic moment at a respective relative precision of around 30 and 3000 times better than the previous best values performed by the hydogen MASER experiment [27] and the ATRAP collaboration [18], respectively. More recently, this has been used to tighten constraints on a possible axion-antiproton coupling, improving by up to five orders of magnitude the astrophysicsbased previous measurements [28].
In parallel, a range of sophisticated laser-based techniques have been developed for the preparation, manipulation, and detection of atomic ions [29,30]. Early on, it was recognized that these techniques could also be used to support control and precision spectroscopy of other quantum systems [31], such as the (anti-)proton, ultimately leading to the experimental realization of quantum logic spectroscopy [32], the development of the single-ion aluminum clock [33,34], and more recently quantum logic spectroscopy of molecular ions [35,36] and even highly charged ions [37]. In Penning traps, applications to mass spectrometry have been proposed by using a single 40 Ca + ion as a sensor in a double Penning trap system [38]. More recently, the idea of using two Penning traps connected by a common electrode has been proposed to improve cooling in present g-factor experiments [39]. In previous work [40,41], we have already discussed concepts for coupling of atomic ions to single (anti-)protons using the double-well technique first suggested in Ref. [29] and experimentally demonstrated in radio-frequency (rf) traps [42,43]. As a first experimental step, we have reported on the prospects of coupling pairs of 9 Be + ions in a double-well potential and on the implementation of Doppler cooling of 9 Be + ions in a specifically constructed 5 T cryogenic Penning-trap setup [44]. Finally, we have already discussed the realization of elementary quantum logic operations with (anti-)protons in Penning traps [45].
In this work, we discuss in detail the experimental conditions of a scenario for implementing quantum logic spectroscopy of (anti-)protons. We also provide a few comments on the expected implications for constraints on SME coefficients.
The outline of this paper is as follows. Section II reviews the main techniques currently employed in Penning-trap experiments. A brief review of the interplay between CPT and Lorentz violation, the SME test framework, and the prediction of sidereal variations is presented in Sec. III. In Sec. IV, our experimental concept is presented in more detail. A time and error budget with projected constraints on various SME coefficients is presented in Sec. V. Finally, a short summary is contained in Sec. VI.

II. PRINCIPLES OF g -FACTOR MEASUREMENTS IN PENNING TRAPS
The Penning trap is an essential tool for high-precision measurements of fundamental properties of charged particles, such as charge-to-mass ratios or g-factors. In this device, a charged particle of interest is stored by superimposing an electric quadrupole field and a strong homogeneous magnetic field B in the axial direction. Figure 1 shows a sketch of a Penning trap with cylindrical shape where a quadrupole electric field is generated by applying a voltage V RE between the ring and endcap electrodes and a voltage V CE between the correction and endcap electrodes. By choosing a proper set of voltages as well as an optimized trap geometry, a harmonic potential shape along the trap axis V (z) ≈ V RE C 2 z 2 is obtained, where C 2 is the unique nonzero coefficient of a multipole expansion [46].
The trajectory of the ion in the trap can be described as a superposition of three motions, one in the axial direction and two in the radial plane, known as the modified cyclotron and magnetron motions. The axial motion has characteristic frequency where m and q are the mass and charge of the stored particle, respectively. In the radial plane, the characteristic frequencies are given by where ν + and ν − are the modified cyclotron and magnetron frequencies, respectively. The motional frequencies of the particle can be determined accurately through the image currents induced in the trap electrodes by the charge's motion [47,48]. They are related to the free cyclotron frequency which depends on the magnitude B = | B| of the magnetic field strength and the charge-to-mass ratio of the particle, through the invariance theorem [49] according to In contrast to the orbital oscillation frequencies of a trapped particle, the spin-precession frequency for particles with magnetic moments, so-called Larmor frequency ν L , cannot be determined directly with this image-current detection technique. Instead, the continuous Stern-Gerlach effect is used [50], where an inhomogeneous magnetic field with z component B z (z, ρ) = B 0 + B 2 (z 2 − ρ 2 /2), a so-called magnetic bottle, is superimposed on the trap center. This magnetic bottle couples the spin magnetic moment to the axial frequency and causes an observable shift of the axial frequency due to a spin transition [51]. Taking into account the net magnetic moment of the trapped particle, which is composed of its spin magnetic moment and the magnetic moment due to the two radial motions, the resulting shift in the axial frequency is given by ∆ν z,SF = 2πm where is the reduced Planck constant and n ± and m s are the quantum numbers associated with the modified cyclotron and magnetron motion and spin orientation, respectively. The Larmor frequency is then determined by measuring spin-flip probabilities for different externally applied drive frequencies through observation of the associated jumps in the axial frequency.
To observe clear jumps of the axial frequency due to a change in m s , spurious noise-driven transitions of n + (which scale with the value of n + directly) need to be suppressed so that the resulting frequency changes do not mask the effect of spin flips. Typically, this requires a thermal occupation of the cyclotron motion corresponding to a temperature of about 100 mK or below [13]. Cooling to such low temperatures is realized by thermalizing the cyclotron motion with a tank circuit at liquidhelium temperature for some time, and then monitoring the axial frequency in an alternating fashion. Only when the energy in the modified cyclotron oscillator is low enough, experimentation can proceed. Currently the best achieved median time range to prepare particles at modified cyclotron energies that provide spin state detection fidelities of about 85% is around 20 min [13,19]. This step consumes approximately 90% of the total measurement time budget, and thus represents the key target for substantial improvements through quantum logic spectroscopy by speeding up the above measurement procedure.
In a double-trap g-factor measurement the spin state of the particle is first initialized in the magnetic bottle trap. Afterwards, the particle is moved to the homogeneous precision trap where the cyclotron frequency is measured while a magnetic rf drive is applied to flip, once resonant, the particle's spin. In a next step the particle is moved back to the analysis trap and the spin state is analyzed. This method is used to scan over spin-flip drive frequencies ν rf yielding the spin-flip probability as a function of the frequency ratio ν rf /ν L from which the Larmor frequency can be extracted. The relation valid in the SM, then permits the determination of the gfactor in the homogeneous magnetic field of the precision trap.

III. SME SENSITIVITIES IN PENNING-TRAP EXPERIMENTS
For a conventional Lorentz-invariant charged Dirac fermion, such as the SM proton, this g-factor corresponds to a coupling that parametrizes the strength of the Pauli interaction [52] determined via the effective-Lagrangian contribution (g−2) m F µνψ σ µν ψ, where F µν is the electromagnetic field strength and ψ is a single Dirac fermion field. The CPT theorem guarantees that all models of this type are CPT symmetric, so that both a fermion and its associated antifermion necessarily possess the same g. But this reasoning also implies that CPT violation requires a framework different from the SM, in which the derivation of Eq. (6) is likely invalidated. Thus, for a clean interpretation and presentation of the protonantiproton ν L /ν c comparison, a more meaningful figure of merit might involve coefficients of a suitable test framework that allows for departures from CPT invariance.
If the test framework is to be a realistic field theory with conventional quantum mechanics, rigorous arguments can be made that CPT violation must be accompanied by deviations from Lorentz symmetry [53]. Such a framework can therefore be expected to predict a dependence of ν L /ν c on the orientation and velocity of the Penning trap. This expectation is borne out by the aforementioned SME, a framework that allows the general classification of such orientation and velocity dependence in all physical systems. The basic idea is that in addition to established physics, the SME Lagrangian contains small Lorentz-violating corrections parametrized by preferred directions in the form of fixed vectors and tensors [54]: , etc. The components of these vectors and tensors represent the coefficients of the SME, which precision experiments seek to constrain. About half of these coefficients also break CPT symmetry. The notation is chosen such that key physical properties become apparent. For example, the subscript p signifies that the coefficient is associated with the proton. We remark that just like m and g, these coefficients cannot be specified independently for a particle and its antiparticle. The subscript F signifies proportionality to F µν , which means in the present context that its physical effects will grow with the trapping B field, and the superscript in parentheses indicates mass dimension [54].
With these considerations, one generically expects ν L /ν c to be given by The CPT-violating corrections are encoded in the functions δf and δf , which can in general differ from one another; they can also contain additional physical quantities, such as the proton mass m and the trapping field B. We note again that both expressions depend on the same set of proton coefficients g, m, b µ p , g λµν p , b (5)λµν p etc. Explicit expressions for δf and δf can be inferred from the results in Ref. [54]. In the remainder of this section, we present a brief overview of these results.
We begin by considering an inertial Penning trap. The B field of such a trap will exhibit a fixed orientation with respect to the SME preferred directions b µ p , g λµν p , b (5)λµν p , . . . [54]. For this system, we may select coordinates such that B = (0, 0, B) points entirely along the x 3 axis. For leading-order SME effects, it can be argued that the electric field in the Penning trap can be neglected and that dominant SME contributions are given by corrections δE p n,s to the usual homogenous-B Landau levels [54].
With this setup, one finds for the proton [54]: Although each stack of Landau levels, determined by the spin state, is affected differently, the dominant SME shifts within a fixed stack are identical for each level, so that the transition frequencies ν p c and νp c are unaffected within the SME at leading order. For subdominant effects, see Ref. [55].
On the other hand, spin-flip transitions between Landau stacks do contain leading-order corrections: Here, the x 3 axis denotes the direction of the B * field in the antiproton Penning trap, which may differ in magnitude and direction from B in the proton trap. The tilde coefficients represent combinations of preferred-direction coefficients in the SME: and b * 3 The grouping of SME coefficients in these expression has been chosen so as to display their properties under rotations, as will become apparent below.
Since Lorentz-and CPT-violating directions represent 4-vectors and 4-tensors, their components can in principle be specified with respect to any basis. However, for practical reasons, such as for comparison between different experiments, measurements are quoted with respect to a standard coordinate system. The canonical choice for this system is the Sun-centered celestial equatorial frame (SCCEF) [56]. The above tilde componentsb 3 p , b 33 F,p ,b * 3 p , andb * 33 F,p must then be expressed in terms of SCCEF components, conventionally denoted with capital coordinates and indices X J = (X, Y, Z). For usual terrestrial Penning traps, the Lorentz transformation between the laboratory and the SCCEF can be decomposed into various contributions including the rotations and boosts arising from Earth's revolution about its axis and its orbital motion around the Sun. Below, we focus on the primary contribution only: the rotation transformation due to the laboratory's motion about the Earth's axis. Purely non-inertial effects due to the rotation of the laboratory frame, such as Coriolis-type corrections, are small and will be neglected. In this context, it is useful to introduce also standard local laboratory frames with a notation involving lower-case coordinates and indices x j = (x, y, z) with the corresponding unit vectors (x,ŷ,ẑ); the conventional choice of these axes hasẑ pointing straight up to the local zenith, andx andŷ are oriented horizontally and point in the local southern and eastern directions, respectively. In these coordinates, B may no longer point along the z direction, but rather alongn = B/B = (n 1 , n 2 , n 3 ). For example,n = (0, 0, 1) for B vertically upward, or n = −(cos ξ, sin ξ, 0) for horizontal B at an angle ξ west of north. The conversion between such laboratory frames and SCCEF coordinates is then given by Here, χ denotes the colatitude of the laboratory, ω ⊕ = 2π/(23 h 56 min) the Earth's sidereal angular frequency, and T ⊕ the local sidereal time [54].
With the above definitions, the Penning-trap frequency shifts (9) can be expressed with SCCEF components as follows: where starred or primed quantities refer to the antiproton laboratory, which may differ from the proton laboratory. With these results, explicit expressions for the corrections δf and δf in Eq. (7) can be obtained. Note that these expressions imply an oscillatory time dependence of the observable 2ν L /ν c in the form of sidereal variations [54]. Note also the dependence on the laboratory's location via both colatitude and local sidereal time, as well as the dependence on the B-field direction and magnitude.

IV. EXPERIMENTAL CONCEPT FOR QUANTUM LOGIC SPECTROSCOPY
From the discussion in Sec. II, it should be evident that g-factor measurements with (anti-)protons are extremely challenging because of the extraordinarily low noise levels required, the long cooling times, and the need to have outstanding control of systematic effects in the presence of considerable residual magnetic-field inhomogeneities. Measurements of sidereal variations in proton-antiproton ν L /ν c comparisons introduce an additional layer of complexity because of the amount of data required and the complication due to the interleaved long cooling times. Different sympathetic laser-cooling approaches have already the potential to reduce considerably the time required for cooling. Here, in addition, we discuss prerequisites and approaches for the implementation of quantum logic spectroscopy, which would allow one to probe all resonances essentially from the motional ground state of the particle, with cooling or spin-state detection happening on time scales of a second or less. The underlying principles were suggested in Refs. [29,31] and result in the following protocol, which is also illustrated in the bottom right part of Fig. 2. 1. Both particles are assumed to be cooled to the ground state and the 'logic' ion (blue) is prepared in the |↑ state as shown in Fig. 2. The spin state of the (anti-)proton (red) is unknown.
2. An rf pulse is applied to flip the (anti-)proton spin if and only if the (anti-)proton spin-state was |↓ previously, and to add one quantum of motion in the axial mode. This effectively implements a SWAP gate between the (anti-)proton's spin and motional (n p = 0 vs. n p = 1) qubit [45].
3. The motional state of the (anti-)proton is transferred to the 'logic' ion by a double-well potential technique [29,44]. This realizes a SWAP gate between the (anti-)proton's and the 'logic' ion's motional (n l = 0 vs. n l = 1) qubits.

4.
A laser pulse flips the 'logic' ion's spin state and removes a quantum of motion if and only if it was previously in the n l = 1 state. This realizes a SWAP gate between the 'logic' ion's motional qubit and its spin state.
5. The fluorescence signal induced by a resonant laser beam is used to detect the spin state of the 'logic' ion. If and only if the ion is |↑ at the end of the previous step, scattered photons are detected.
6. The 'logic' ion is re-initialized to |↑ by optical pumping, and the motion of both particles is reinitialized to the ground state, if necessary.
Repetition of steps 1 to 6, accompanied by variation of the frequency of a rf pulse to flip the (anti-)proton spin in step 1 yields a resonance curve from which ν L can be extracted. Implementing these ideas requires several ingredients: a. Spin-state discrimination for the 'logic' ion b. The ability to ground-state cool the 'logic' ion and carry out spin-motional state swap operations on the 'logic' ion c. A suitable trap stack for proton loading and antiproton injection d. A suitable coupling mechanism allowing the transfer of single quanta of motion between the 'logic' ion and the (anti-)proton Based on these requirements, a multi-trap method is favoured. Different points above impose mutually incompatible boundary conditions difficult to meet in a single trap. For example, probing the Larmor or cyclotron resonance would require a trap of relatively large scale (typically 5-10 mm), such that the systematics are well controlled [13]. Laser access, on the other hand, typically involves modifications to the trap geometry which slightly distort potentials and impose additional system-atic effects [57]. Furthermore, laser-induced charging may inhibit the refined control of potentials necessary for a precision trap [58]. Implementing a coupling mechanism between single (anti-)protons and the 'logic' ion at the single motional quantum level that works in a reasonably short amount of time typically requires a small length scale for the trap [41], so that both particles can be kept in close proximity, which is again at odds with the demands of a precision trap. Implementing the swap operation between the (anti-)proton's spin and motional states would typically also call for a small-scale trap and additional coupling fields [29], which might again conflict with the needs for a precision trap.
Looking at the list of requirements, the core functions can be realized in a trap stack, schematically illustrated in the upper left of Fig. 2, with four different zones for these functions. Items (a) and (b) can be carried out in a trap zone with laser-beam access and fluorescence detection. To allow for item (c), antiproton loading, one end of the trap stack needs to be open for the injection of antiprotons. Item (d) could be accomplished in a dedicated coupling trap with small length scales in order to reduce the distance between the 'logic' ion and the (anti-)proton. Item (e) could be performed in a sideband trap for spinmotional coupling with (anti-)protons, inspired by recent advances in rf or microwave driven quantum logic operations for trapped-ion qubits [45]. Item (f) could be realized using a precision trap as employed in present experiments. Below we discuss these different trap zones before highlighting how they come together to allow the implementation of quantum logic spectroscopy. These traps would be interconnected through physical transport of particles along the magnetic-field axis (item (g) above). In order to facilitate transport, traps should be built as cylindrical Penning traps. The physical location of both particles during each step of the measurement algorithm is indicated in the lower left part of Fig. 2.
We start by discussing the general structure of the apparatus. As in present g-factor measurements [13,20], the experiment should be operated at cryogenic temperatures to facilitate pre-cooling of (anti-)protons using cryogenic tank circuits, and it should feature a hermetically sealed secondary vacuum enclosure ('trap can') to support long storage times of antiprotons. In order to obtain magnetic fields as homogeneous as possible, radial optical access ports in the Penning-trap's superconducting magnet should be avoided, as they are likely to introduce additional inhomogeneities in the magnetic field [59]. As a consequence, all laser beams would have to be brought in along the axis of the superconducting magnet, and fluorescence detection would have to occur along that direction as well, and thus downstream from the antiproton source. The laser cooling and detection trap would thereby also have to be the last trap in the trap stack with the imaging optics downstream of it. In what follows, we consider a magnetic field of 5 T for higher cyclotron and Larmor frequencies and beryllium ions are assumed as 'logic' ions of choice.

a. Laser cooling and detection trap
The design of our current laser cooling and detection trap has been discussed in Ref. [44]. Beryllium ions are produced in the center of the laser trap by ablation of beryllium atoms from a solid target. A single trapped ion can be obtained from repeated splitting of the loaded sample, using suitable voltage ramps, and then discarding one part. This ion can be laser-cooled near the Doppler limit using two laser fields resonant with the | 2 S 1/2 , m J = +1/2 → | 2 P 3/2 , m J = +3/2 and | 2 S 1/2 , m J = −1/2 → | 2 P 3/2 , m J = +1/2 transitions in 9 Be + as is shown in Fig. 3. Details of the laser-system setup for this experimental arrangement can be found in Refs. [44,60]. Be + at a 5 T magnetic field. The transitions used for Doppler cooling and repumping as well as for the two-photon stimulated Raman process are shown. Laser detunings should be chosen for the mI = +3/2 nuclear spin state. Weak off-resonant transitions will eventually lead to the preparation of that nuclear spin state starting from mI = +3/2. Energy levels are not drawn to scale.

b. Ground-state cooling and spin-motional SWAP gate for 'logic' ion
Ground-state cooling and the implementation of the SWAP gate between the ion's spin and motional qubits will require the implementation of red and blue spinmotional sideband pulses [61]. In the case of 9 Be + , these will have to be carried out on the Zeeman qubit in the 9 Be + ground state, split by roughly 139 GHz at 5 T. In Ref. [62], we have demonstrated a pulsedlaser approach to bridge this big Zeeman splitting by using two laser beams as depicted in Fig. 3, and successfully implemented it at a low magnetic field of 22.3 mT. For sideband cooling, blue-sideband pulses on the | 2 S 1/2 , m J = +1/2 to | 2 S 1/2 , m J = −1/2 transition can be employed to remove motional energy from the ion, followed by repump pulses on the | 2 S 1/2 , m J = −1/2 to | 2 P 3/2 , m J = +1/2 transition to re-initialize the ion in the | 2 S 1/2 , m J = +1/2 state. In order for the laser beams to couple only to the axial direction for sideband cooling, they need to be irradiated with a wave-vector difference along the trap axis as illustrated in Fig. 2.

c. Proton loading and antiproton injection
With one end of the trap stack open, antiprotons can be injected and trapped using a catching and reservoir trap as in the BASE-CERN setup [41]. Proton loading will have to be carried out differently. In the BASE-CERN setup, the electron gun used for this operation is pointing in the opposite direction as the antiproton beam, a direction obstructed by the fluorescence detection system in our setup. We therefore plan to use an off-axis electron gun and an organic target for proton loading [63]. (Anti-)protons could be stored for use in a reservoir trap similar to that demonstrated at BASE-CERN [64].

d. Coulomb coupling using a double-well potential
A crucial feature of our experimental protocol is the coupling between a single (anti-)proton and a laser-cooled 9 Be + ion through the Coulomb interaction in free space. For this purpose, the two particles should be spatially separated and held in adjacent potential wells. We have already discussed this ion coupling mechanism in Ref. [41,44,65]. For two charges q 1 and q 2 trapped in two potential wells separated by a distance s 0 , if the radial component of the Coulomb force is neglected and a one-dimensional interaction is assumed, in the limit of weak coupling with both particles oscillating at the same axial frequency ν z ν 1 = ν 2 , the particles exchange energy at the time scale where Ω ex is the so-called coupling rate, 0 is the permittivity of free space, and m 1 and m 2 are the masses of particle 1 and 2, respectively. This coupling mechanism has been demonstrated for two 9 Be + ions in their motional ground states in a surface-electrode rf trap [42] and for 40 Ca + ions in Ref. [43]; it should also be feasible in Penning traps because both systems are analogous in the axial mode. Figure 4 shows the double-well potential shape along the axial direction for storing a 9 Be + ion and a single antiproton. For an axial frequency ν z = 4 MHz, a distance s 0 = 300 µm leads to an exchange time of t ex 3.7 ms [41]. The choice of 9 Be + as the 'logic' ion implies higher coupling rates due to its low mass compared to other readily laser-cooled ions. Also, it is desirable to have the ratio of the 'logic' ion mass to the proton mass as close to unity as possible in order to facilitate the engineering of a doublewell potential with two equal trap frequencies for both particles, which requires a potential-curvature ratio given by the mass ratio. Another advantage of 9 Be + is the relatively simple atomic structure without metastable states which would otherwise require a substantial amount of repumping fields as in the case of 40 Ca + . On the other hand, for ground-state cooling, 9 Be + requires stimulated Raman transitions, which can be challenging to implement given the Zeeman splittings in the Penning trap that typically exceed 100 GHz.
A further essential feature for achieving high coupling rates is a small distance s 0 between the two trapped particles, as is apparent from Eq. (14). The limitations in this respect are the trap construction together with anomalous motional heating rates [66], and the initial amplitude of radial motion of the (anti-)protons before they are brought into the coupling trap. To place the minima of the potential wells closer together in a cylindrical Penning trap, both thickness and diameter of the electrodes must be reduced. Microfabrication techniques can be employed for the construction of such traps [67], allowing distances s 0 on the order of a few hundred µm. However, this involves reducing the trap diameter, which would lead to an increase of the anomalous motional heating rate.
The cryogenic operation of the trap will suppress the anomalous motional heating. Recently, rates of a few quanta per hour have been reported for the motional heating of a single antiproton in a cryogenic Penning trap with a 3.6 mm inner diameter [68]. Assuming that magnetron amplitudes of (anti-)protons injected into the coupling trap could be on the order of a few tens of microns, it would seem reasonable to consider an ion-toelectrode distance in the coupling trap of no more than 400-500 µm.
Using such a double-well trap, sympathetic cooling of single (anti-)protons to the motional ground state could proceed as follows: 1. The 9 Be + ion is cooled to the motional ground state by using Doppler [69] and sideband cooling [70] techniques as discussed above. A single (anti-)proton is prepared in the precision trap and pre-cooled by a resistive cooling technique.
2. Both particles are adiabatically shuttled to the coupling trap. The axial motional modes of both particles are tuned into resonance for the time it takes to swap the motional states between the two particles.
3. The 'logic' ion is transported back to the laser cooling trap and re-cooled to the motional ground state using the procedure outlined in step 1.
This implements sympathetic cooling of the (anti-)proton to the motional ground state. The relevant time scales would be the time required for motional ground state cooling of the 9 Be + ions (a few milliseconds at most) and the motional exchange time of t ex = 3.7 ms, enabling ground state cooling of the axial degree of freedom of the (anti-)proton within milliseconds, significantly shorter than the time scales associated with resistive cooling. The motional state exchange could exhibit imperfections if the initially finite magnetron radius of the (anti-)proton (see above) leads to a different effective distance of the two particles. This effect should be small, because an (anti-)proton magnetron radius of a few tens of microns would not effectively change the relevant distance of the two particles of e.g. 300 µm by much. Other imperfections could arise from timing and resonancecondition imperfections or due to trap anharmonicities. 4. In the case of imperfect cooling, steps 2 and 3 could be performed several times.

5.
To cool the radial modes, a sideband coupling procedure between the radial and axial modes of the (anti-)proton must be carried out in the precision trap [71], followed by repetition of the above axialcooling protocol. . Potential shape along the z axis for storing a single 9 Be + ion (blue) and a single antiproton (red) at a distance of 300 µm. The coupling trap used to simulate this potential shape is composed of nine ring-shaped electrodes with an inner diameter of 800 µm and a thickness of 200 µm each. The distance between adjacent electrodes is 50 µm. The potential shape has been calculated for an axial frequency of 4 MHz for both particles in a 5 T Penning trap. The potential shape for the coupling between a 9 Be + ion and a single proton was reported in Ref. [41].

e. Sideband coupling between motional and spin states
To couple the motional and spin states of the (anti-)proton in the sideband trap and to carry out the SWAP operation between the (anti-)proton's spin state and its motional qubit, a rf pulse close to the (anti-)proton spin flip resonance ω SF = g p µ N B/ should be applied, where µ N is the nuclear magneton and g p is the (anti-)proton gyromagnetic ratio in units of µ N . For this purpose, the proposal of Ref. [31] envisioned the use of an oscillating-amplitude magnetic-field gradient to couple spin and motional degrees of freedom. This mechanism is conceptually similar to the oscillating near-field microwave quantum logic gates of Ref. [72]. However, with the particle at least a factor of ten further away from the conductor inducing the transition, without any current scaling, the gradient would be about two orders of magnitude smaller. In addition, compared to an electronic atomic magnetic moment, the magnetic moment of the (anti-)proton is lower by about three orders of magnitude, making effective Rabi rates about five orders of magnitude smaller.
This issue can be overcome by introducing a strong static magnetic-field gradient [73]. In a variation of the magnetic bottle approach, very large magnetic field gradients can be generated in the center of a Penning trap by introducing a ferromagnetic electrode as an endcap electrode, effectively trapping the particle in the wing of a magnetic bottle, where the magnetic-field curvature vanishes and the gradient is strongest. For traps with an inner diameter of about 1 mm, gradients as strong as 1000 T/m can be anticipated, much stronger than the ≈ 35 T/m that was used in Ref. [72], yet at considerably larger surface-to-electrode distance, thereby at least partly compensating for the smallness of the magnetic moment of the (anti-)proton. In Ref. [45], we have analyzed the behavior of an (anti-)proton in a large magnetic-field gradient superimposed on the Penning trap and found that as in Ref. [73], a homogeneous magnetic rf field at a spin-motional sideband frequency will now allow the spin and motional degrees of freedom to be coupled as required for the implementation of the SWAP gate.

f. Cyclotron-frequency measurement
In addition to the Larmor frequency, a g-factor determination based on Eq. 6 also requires a cyclotronfrequency measurement for the (anti-)proton. For this purpose, different approaches may be taken. One of these is based on the image-current technique in the precision trap. However, this method requires the particle to be in thermal equilibrium with the detection system [48], which may be a disadvantage for tests with the highest precision when residual magnetic-field inhomogeneities are present.
An alternative scheme that probes all transition out of the motional ground state at the level of a single quantum has already been discussed in Refs. [29,31]. It relies on the application of the same sympathetic ground-state cooling protocol outlined above and on an electric-field pulse to probe the motional resonances in the precision trap. A possible excitation as a result of the pulse can then be read out using the Coulomb-coupling technique and the 'logic' ion again. By varying the excitation-pulse frequency, the motional resonances can be recorded and the free cyclotron frequency determined.
A variant of this method would be carried out by us-ing a Ramsey-like experiment [74], where two electricfield pulses, separated by a free-precession time t f , are applied [75,76]. The first excitation pulse produces a coherent state by imposing a displacementD(α) = exp(αâ † −α * â ) on the motional mode of interest, where α is the displacement amplitude andâ andâ † represent the usual creation and annihilation operators of the trap harmonic oscillator. The expected state overlap is given by n D n = exp(−1/2|α| 2 )L n (|α| 2 ), where L n is the Laguerre polynomial [77]. After a free period t f , a relative phase φ = ∆×t f is accumulated in phase space, where ∆ is the detuning of the excitation. For φ = π(k + 1) with k an integer, the second pulse undoes the displacement and returns the system to |0 [76]. For any other outcome, the final state will have an excited state contribution different from |0 which can be detected using the axial Coulomb coupling technique, combined with π-pulses which can exchange the axial and cyclotron degrees of freedom [71].

g. Ion transport
Regarding ion transport, adiabatic transport is achieved for a harmonic trapping potential with a minimum that can be shuttled along the axial direction between two traps without any abrupt changes to the potentials. In the context of quantum-information processing with linear rf traps, this goal has already been achieved for particles in their motional ground state [78], even through complex trap geometries [79,80]. In principle, the axial degree of freedom in a Penning trap is purely electrostatic, as for the linear rf trap, so one would expect comparable results. With the methods available in the BASE experiment at CERN, transport heating has been characterized to 0(2) quanta per transport in the cyclotron mode [81]. Stimulated-Raman laser control would be able to place tighter bounds here, and would be expected to yield results similar to Paul traps with significantly less than one quantum of motion per round-trip, as observed in Paul traps. Note that this is a prerequisite to carry out the above protocols.

V. SENSITIVITIES AND SME IMPLICATION
By using the experimental protocol discussed in this work, a Larmor-frequency determination will consist of several experimental cycles. For each i-th cycle, the spin-flip probability is determined for a pulse excitation ν i L close to ν L . The principal steps for a single measurement cycle are summarized in Table I. The Larmor resonance is probed by using the Ramsey method with two π/2 pulses separated by a free-precession time t f . Each pulse must be a transverse magnetic rf field at an amplitude b rf and duration t p = πm/g p qb rf ≈ 3 ms for b rf ≈ 2 µT [82], which must be shorter than t f . Transport times will depend on the low-pass RC filter of the DC lines, whose cutoff frequency will be given by the magnetron frequency. After that, sideband coupling between the motional and spin states of the (anti-)proton can be produced in around 2 ms with low intrinsic error [45]. Coulomb coupling is discussed in detail in Sec. IV d, and finally, a SWAP gate between the Be ion's motional and internal qubits can be performed in several hundred µs [62], as can photon detection [42]. For t f ≈ 100 ms, the probe time is around 90 % of the total measurement time.
For a complete g-factor measurement, the axial and modified cyclotron frequencies will be determined by scanning external excitation frequencies ν i z and ν i + close to ν z and ν + , respectively. The cyclotron frequency can be calculated by using Eq. 4, where ν − = ν 2 z /2ν + [49]. As before, the Ramsey method can be used to measure these motional frequencies by using pulses with duration of a few tens of µs [76]. In addition, for the modified cyclotron frequency, a pulse excitation at ν + − ν z must be applied in order to couple the modified cyclotron and axial modes. The duration of this pulse can be inferred from Ref. [71]. It is important to note that the cooling of any of the ion motions can be performed at any moment of the complete experimental process without a significant increase in the measurement time.
Systematic shifts of the ion frequencies based on energy-and time-dependent variations in the precision trap [83], where the spectroscopy takes place, as well as the frequency stability of the ion transitions will be the fundamental sources of uncertainty. Since the (anti-)proton will be in a well-defined energy state close to the ground state, systematic corrections based on highorder electric-and magnetic-field imperfections as well as relativistic effects will contribute with fractional uncertainties at the sub-ppt level since they are related to the energy of the ion eigenstates [84,85]. In the case of the image charge effects, uncertainties due to this correction at the ppt level are expected [13]. However, our detection method would allow the use of larger trap sizes, which would reduce this error to levels around 1 ppt [86]. Timedependent variations of the electric and magnetic field can be reduced to the sub-ppt levels by using linear interpolations of several axial-frequency measurements [13] as well as several measurements of the qubit transition of the Be ion (see Sec. IV b) to interpolate the magneticfield variations.
Assuming that frequency uncertainties of the ion transitions at a level higher than 1 ppt are dominated by the frequency stability of the ion transitions, under some considerations [87], the stability for fractional frequency deviations y will be given by the quantum projection noise [88] where ν 0 is the transition frequency and τ the total measurement duration. In the case of the Larmor transition, for a magnetic field of 5 T, a stability of 2×10 −9 / τ /s is estimated. For the axial and modified cyclotron frequencies, the stabilities are 5×10 −7 / τ /s and 7×10 −9 / τ /s, respectively. This implies an expected stability at the 10 ppt level for the Larmor frequency after measurement times of around 12 hours. For the cyclotron frequency, the expected stability would increase by a factor of about 5.5 for the same measurement time, which should be divided equally between the axial and modified cyclotron frequencies. The magnetron frequency would contribute Table II. Constraints on variousb J p ,b * J p ,b JK F,p , andb * JK F,p SME coefficients. The second column contains estimates derived in Ref. [91] from two recent BASE measurements [13,20]. The third column lists constraints projected to be achievable with the methods discussed in this work.

Coefficient Estimated constraints Projected constraints
at the sub-ppt level by using the relation ν − = ν 2 z /2ν + . It should be noted, that the motional frequencies are not a two-level system and a factor 1/2|α| [76] in the frecuency stabilities is expected. Here we assume |α| ≈ 1 for the lowest sensitivity to field-inhomogeneity induced systematics. In addition, frequency stability degradation (see, e.g., Ref. [89]) due to dead time on the measurement cycle is presumed to be well below the quantum projection noise instability [87,90]. Therefore, uncertainties at the order of several tens of ppt in the (anti-)proton g-factor would be expected with a measurement duration on a sidereal-time scale.
These rapid measurements of the (anti-)proton frequencies would allow access to those SME coefficients discussed in Sec. III that are associated with sidereal variations, which remain experimentally unexplored in current g-factor proton experiments [54,91]. Following Ref. [54] and applying the techniques discussed in this work to a hypothetical experiment located at BASE CERN with a colatitude of 43.8 • and a magnetic-field strength of 5 T pointing to the local zenith, several bounds on SME coefficients given by for the second harmonic of sidereal variations, respectively, are obtained. Here, ∆b J p ≡ (b J p −b * J p )/2 and ∆b JK F,p ≡ (b JK F,p −b * JK F,p )/2 with J, K = X, Y, Z. In addition, a fractional precision of the antiproton and proton magnetic moments of 100 and 10 ppt are estimated for measurement of sidereal variations and constant effects, respectively. The resulting constraints on these tilde coefficients are summarized in Table II. Natural units with c = = 1 are used. The upper part of Table II shows the constant coefficients, which are measured with current techniques, and the lower part displays the previously unexplored coefficients in the SME's proton sector that could be accessed with the techniques presented in this work.

VI. SUMMARY
CPT violation is closely intertwined with a breakdown of Lorentz symmetry, so that sidereal modulations of physical observables in a terrestrial laboratory represent a key experimental signature for such effects. To exploit this insight for investigations in Penning traps, the temporal resolution for measurements should be at subsidereal time scales, a criterion difficult to meet in existing experimental setups involving the (anti-)proton. We have identified the complexity and duration of the currently used spin-readout techniques as the primary target for substantial improvements.
More specifically, we have argued that recent experimental advances in laser-based quantum logic spectroscopy have opened an avenue for shortening highfidelity (anti-)proton spin readout by three orders of magnitude. This idea can be realized via a co-trapped 9 Be + 'logic' ion: a SWAP operation essentially permits the conversion of the (anti-)proton's spin-state information to motional excitations of the Be + ion, which are readily detectable with conventional methods.
In this context, we have presented the key features of the experimental implementation of this idea. In particular, we have described the design of the required cryogenic multitrap setup consisting of interconnected laser, coupling, sideband, and precision traps, supplemented by an antiproton storage trap. We have further detailed the associated experimental protocol with particular focus on the physical underpinnings of the SWAP operation. Finally, we have provided estimates for time-and errorbudgets. These considerations have allowed us to project improvements by up to two orders of magnitude for existing constraints as well as sensitivity to previously unmeasured effects tripling the number of experimentally accessible CPT-violating SME coefficients at mass dimension up to six.

ACKNOWLEDGMENTS
We are grateful for discussions with D.J. Wineland, J.J. Bollinger, R.C. Thompson, and P.O. Schmidt. This work was supported by PTB, LUH, and DFG through the clusters of excellence QUEST and QuantumFrontiers as well as through the Collaborative Research Center 1227 (DQ-mat) and ERC StG "QLEDS." We also acknowledge financial support from the RIKEN Pioneering Project Funding and the MPG/RIKEN/PTB Center for Time, Constants and Fundamental Symmetries. R.L. acknowledges support by the Alexander von Humboldt Foundation.