## Abstract

Efficient cooling of trapped charged particles is essential to many fundamental physics experiments^{1,2}, to high-precision metrology^{3,4} and to quantum technology^{5,6}. Until now, sympathetic cooling has required close-range Coulomb interactions^{7,8}, but there has been a sustained desire to bring laser-cooling techniques to particles in macroscopically separated traps^{5,9,10}, extending quantum control techniques to previously inaccessible particles such as highly charged ions, molecular ions and antimatter. Here we demonstrate sympathetic cooling of a single proton using laser-cooled Be^{+} ions in spatially separated Penning traps. The traps are connected by a superconducting LC circuit that enables energy exchange over a distance of 9 cm. We also demonstrate the cooling of a resonant mode of a macroscopic LC circuit with laser-cooled ions and sympathetic cooling of an individually trapped proton, reaching temperatures far below the environmental temperature. Notably, as this technique uses only image–current interactions, it can be easily applied to an experiment with antiprotons^{1}, facilitating improved precision in matter–antimatter comparisons^{11} and dark matter searches^{12,13}.

## Main

Measurements of the charge-to-mass ratio and *g*-factor of the proton and antiproton, a prominent, stable particle–antiparticle system, are limited by cryogenic particle temperatures^{14,15,16}. However, with no electronic structure, (anti)protons (protons and antiprotons) are not amenable to standard laser cooling techniques. Moreover, ions that are easily laser cooled are not readily trapped in the same potential well as negatively charged antiprotons or antimatter molecular ions (for example, \({\bar{{\rm{H}}}}_{2}^{-}\))^{17}. Sympathetic laser cooling with negatively charged ions^{18,19,20} and with microscopically fabricated trapping potentials^{1,21,22} have been proposed. Another technique, proposed over 30 years ago, extended laser cooling to exotic systems by coupling via induced image currents in trap electrodes to ions with a well-suited cooling transition^{9}. Similarly, coupling laser addressable ions to systems with no optical structure is sought after in precision spectroscopy^{23,24}, mass measurements^{25}, quantum information^{10} and quantum engineering^{5}.

We demonstrate sympathetic cooling of a single proton, extending the image-current coupling technique with a superconducting LC circuit that resonantly enhances energy exchange between the proton and laser-cooled ions. We use a cryogenic multi-Penning-trap system to store a single proton in the proton trap and a cloud of Be^{+} ions in a beryllium trap, separated axially by around 9 cm (Fig. 1a). A homogeneous magnetic field *B* parallel to the electrode axis and an electric quadrupole potential at voltage *V*_{0} confines the particles and gives rise to circular magnetron and modified cyclotron motion in the radial plane, and harmonic axial motion, at frequencies *ν*_{−}, *ν*_{+} and *ν*_{z}, respectively^{26}. LC resonators with high quality factors (in our case *Q* ~ 15,000) are commonly used to detect image currents in trap electrodes^{27} and, as shown in Fig. 1a, we connect the resonator to both traps so that the two ion-trap systems are coupled via particle-induced image currents. The LC circuit with total capacitance *C*_{R} ≈ 36 pF and inductance *L*_{R} ≈ 3.0 mH has an equivalent parallel resistance, at resonance frequency *ν*_{0}, of *R*_{p} = 2*πν*_{0}*L*_{R}*Q*. The entire system is modelled by an equivalent circuit where the proton and Be^{+} ions are series LC circuits with capacitances and inductances *C*_{p}, *C*_{Be}, *L*_{p} and *L*_{Be} (ref. ^{28}) connected in parallel to the superconducting LC circuit (Fig. 1b). The Be^{+} ions are also damped by the cooling laser, represented as a variable virtual resistance *R*_{L} (ref. ^{9}). In contrast to the proposal in which energy is exchanged between the ion-trap systems via only a shared electrode^{9}, we use the LC-circuit resonator to couple the axial modes of the trapped particles. On resonance, the large inductance of the resonator coil compensates the electrode capacitance and enhances the ion-induced image current by the *Q*-value. With mechanically machined traps used for precision Penning-trap experiments, the ~10-mHz coupling rates expected from the non-resonant proposal^{11} require minute-scale cooling cycles and are limited by the loss of resonant coupling; for example, from voltage fluctuations of the trapping potential. For the parameters used in our resonant cooling demonstration, energy is exchanged between the proton and sympathetically laser-cooled resonator at a rate of 2.6 Hz (measured by the dip width on resonance) so that thermal equilibrium is reached within seconds and the axial frequencies of the two species are easily matched. Notably, by coupling the ion-trap systems via the resonator, the coupling does not rely on a shared electrode, so the energy exchange rate is not limited by the trap capacitance. Consequently, this cooling scheme can be realised over long distances and with several distributed ion traps.

The noise spectrum of this coupled system is shown in a fast Fourier transform (FFT) of the voltage signal of the resonator (Fig. 1c). The entire system is driven by a combination of the Johnson noise of the resonator and additional voltage noise from the cryogenic amplifier, resulting in an effective noise temperature *T*_{0} = 17.0(2.4) K (where the number in parentheses is the 1*σ* uncertainty). Here, the axial frequencies of the proton in the proton trap and the Be^{+} ion in the beryllium trap were set close to resonance with the LC circuit by adjusting the voltage of the axial potential in each trap, *ν*_{z} ∝ *V*_{0}^{1/2}. In the measured noise spectrum, the detector appears as a broad ~40 Hz (full-width at half-maximum (FWHM)) resonance while the proton and the Be^{+} ions short the parallel resistance of the resonator and appear as narrow dips with widths determined by the charge-to-mass ratios and the trap diameter^{28}.

We demonstrate that the proton, Be^{+} ions and resonator form a system of three coupled oscillators by measuring the noise spectrum at thermal equilibrium. We detune both ion species around three resonator linewidths away from the LC-circuit resonance frequency *ν*_{0} ≈ 479,000 Hz to observe coupling signatures via the FFT lineshape. In these measurements, we store the proton in the proton trap at constant axial frequency and gradually increase the axial frequency of a single Be^{+} ion in the beryllium trap. The resulting FFT spectra show observed particle frequencies at the dip positions (Fig. 2a, dark blue) and two of the normal modes of the coupled three-oscillator system at the maxima (Fig. 2a, red). Near the p-Be^{+} resonance, the axial motion of both particles is no longer determined by the trapping potential alone, and we observe the coupling signature in two of the normal modes of the three-oscillator system. This feature is consistent with the analytical solution derived from the impedance of the circuit model in Fig. 2b (see Methods) and appears when the three oscillators exchange energy. Using simulations (described in the Methods), we show the corresponding time-domain behaviour in the presence of the environmental noise (Fig. 2c). In the absence of environmental noise, the energy of each oscillator as a function of time is deterministic and can be found from the initial phases and energy exchange rates. With environmental noise included, energy is still exchanged and, as shown here when laser cooling is absent, the oscillator energies are determined by this equivalent noise temperature.

We further demonstrate that the temperature of the proton can be modified by coupling to a cloud of excited Be^{+} ions, here consisting of around 15 ions. To this end, we apply a parametric rf drive at 2*ν*_{0} in the beryllium trap, which excites the Be^{+} ions if *ν*_{z,Be} = *ν*_{0} but, as confirmed by background measurements (see Methods), has no direct effect on the proton in the proton trap. By bringing the proton into resonance with the weakly excited Be^{+} ions the Be^{+} ions appear as a broad, shallow dip and the sympathetically excited proton appears as a narrow peak (see Methods). To quantify the energy transferred to the proton, we measure the axial frequency of the proton before and after coupling to the excited Be^{+} ions. Coupling the excited axial mode to the cyclotron mode with a sideband drive at *ν*_{+}− *ν*_{z} transfers the energy of the axial mode to the cyclotron mode with resulting energy *E*_{+} = (*ν*_{+}/*ν*_{z})*E*_{z} (ref. ^{29}). Similar to the continuous Stern–Gerlach effect^{30}, the quadratic component of the magnetic field in the proton trap, *B*_{2} = −0.39(11) T m^{−}^{2}, interacts with the magnetic moment of the modified cyclotron mode at energy *E*_{+,p}, producing the axial frequency shift Δ*ν*_{z} ∝ *B*_{2}Δ*E*_{+,p} (see Methods), which we measure to determine the change in axial energy of the proton. We show the evolution of the standard deviation of the change in the proton energy *E*_{+,p} while the excited Be^{+} ions are tuned to resonance with the proton (Fig. 2d, orange) and when detuned from the resonator (Fig. 2d, blue). During this experiment we interleave on-resonance and off-resonance measurements and see a clear increase in proton energy that arises from the remote, resonator-mediated coupling to the excited ions. In contrast to the off-resonant points, the resonance points exhibit scatter that is nearly three orders of magnitude larger, with a statistical significance of more than 20*σ*. The excitation drive remains on during both of the interleaved measurements to ensure that the increased scatter is due only to the ion–proton coupling. In addition, we constrain variations in the resonator temperature due to off-resonant coupling to the excitation drive (see Methods). For comparison, the scatter when the drive is off is also shown in Fig. 2d (green points).

Our demonstration of sympathetic cooling employs similar axial frequency shift measurements in the presence of a continuously laser-cooled Be^{+} ion cloud. The Be^{+} ions are cooled with the closed ^{2}S_{1/2} → ^{2}P_{3/2} transition and tuned to resonance with the superconducting circuit and the proton (Fig. 3a). The cooling laser damps the axial motion, increasing the equivalent resistance *R*_{L} (Fig. 1b) and reducing the signal of the broad Be^{+} dip. The laser-cooled ions reduce the effective noise temperature in the entire circuit and lower the temperature in a narrow frequency range. Using the narrow proton dip as a temperature sensor for the cooled common mode of the system, we determine the temperature reduction experimentally with well understood energy-dependent shifts of the axial dip and develop further insight into the cooling using time-domain simulations.

A symmetric, cylindrical Penning trap provides a high degree of control over the trapping potential. We use a deliberately introduced trap anharmonicity in the proton trap that shifts the axial frequency by

Here, *C*_{n}(TR) are the coefficients of the expansion of the local trapping potential along the trap axis that depend on the ratio of voltages applied to the central ring electrode (*V*_{0}), and the two nearest correction electrodes (*V*_{CE})^{26,31,32}, referred to as the tuning ratio, TR = *V*_{CE}/*V*_{0}. When the laser-cooled Be^{+} ions are tuned to resonance, the noise energy of the common mode of the proton, resonator and Be^{+} ions is reduced from the noise temperature of the environment resulting in an axial frequency shift,

where *ν*_{z,1}(TR,*T*_{0}) = *ν*_{z} + *δν*_{z}(TR,*T*_{0}) is the axial frequency measured at *T*_{0} when the Be^{+} ions are detuned and *ν*_{z,2}(TR,*T*_{p}) = *ν*_{z} + *δν*_{z}(TR,*T*_{p}) is the axial frequency measured when laser-cooled ions are in resonance and reduce the temperature to *T*_{p}. The trap anharmonicity is characterized by an offset from the ideal tuning ratio ΔTR = TR − TR(*C*_{4} = 0) and a constant determined from the trap geometry *κ* = 45.4 Hz K^{−}^{1} that we crosscheck with additional measurements that use electronic feedback to change the temperature of the resonator. We measure Δ*ν*_{z} as a function of ΔTR and the measured slope *s* determines the change in temperature, Δ*T* = *T*_{0} − *T*_{p} = −*s*/*κ*. The results of an example measurement are shown in blue in Fig. 3b. With ten Be^{+} ions in resonance, we measure a slope *s* = −350(14) Hz and in a background measurement with the Be^{+} ions detuned (Fig. 3b, orange), obtain a slope *s* = 4(13) Hz. This corresponds to a temperature reduction of Δ*T* = 7.7(0.3) K and demonstrates sympathetic laser cooling of a single trapped proton. With a significance of more than 20 standard deviations, this is also a demonstration of remote, image-current mediated sympathetic cooling, applicable to any charged particle without convenient cooling transitions.

The temperature of the proton is determined by the noise power dissipated by the laser-cooled ions. In the circuit representation, increasing the damping of the laser cooling *γ*_{L} increases *R*_{L} and has the effect of lowering the coupling rate of the Be^{+} ions to the resonator \({\mathop{\gamma }\limits^{ \sim }}_{{\rm{Be}}}\), which, in the absence of laser cooling is given by the dip width *γ*_{Be} ∝ *N*_{Be}. For a given number of laser-cooled ions *N*_{Be}, *γ*_{L} must be optimized and in the limiting case when *γ*_{L} ≪ *γ*_{Be} the Be^{+} ions are driven by the resonator and the dip signal is unchanged. Likewise, when *γ*_{L} ≫ *γ*_{Be}, the Be^{+} ions are decoupled from the resonator and the dip signal vanishes. In both limiting cases, the temperature of the resonator and the proton remain unchanged.

However, increasing *N*_{Be} increases \({\mathop{\gamma }\limits^{ \sim }}_{{\rm{Be}}}\) and laser cooling reduces the temperature of the resonator and the proton, even at large *γ*_{L}. To lower the temperature and to investigate the scaling of *T*_{p}, we performed a series of further measurements with varying *N*_{Be} and laser detuning *δ* (Fig. 4). We additionally analysed the temperature scaling by comparing to a temperature model in which the common mode temperature *T*_{CM} of the equivalent circuit arises from competing dissipation sources; the noise temperature of the environment, *T*_{0}, at a coupling rate given by the width of the LC resonance *γ*_{D}, and to the Be^{+} ions at temperature *T*_{Be}. As a result, the system comes to thermal equilibrium at

When \({\mathop{\gamma }\limits^{ \sim }}_{{\rm{B}}{\rm{e}}}\gg {\gamma }_{{\rm{D}}}\), the proton temperature is approximated as

reproducing the 1/*N*_{Be} scaling, by *γ*_{Be} ∝ *N*_{Be}, appearing in the non-resonant proposal^{9,11} and related proposals in the context of trapped ion quantum information^{10,33,34}. In these measurements, the laser detuning *δ* can be viewed as a tuning parameter that changes the *γ*_{L} and subsequently \({\mathop{\gamma }\limits^{ \sim }}_{{\rm{Be}}}\) (see Methods). As a result, the lowest proton temperatures are not found by minimizing the Be^{+} temperature, which would correspond to lower laser detunings, but by maximizing the coupling of the ions to the detector, corresponding to larger \({\mathop{\gamma }\limits^{ \sim }}_{{\rm{Be}}}\). For the experimental parameters used here, *δ* = −90 MHz is the largest laser detuning at which the proton dip is still visible. The largest ion cloud *γ*_{Be} = 164(5) Hz and largest detuning from the centre of the cooling transition *δ* = −90 MHz (Fig. 4) is representative of the lowest temperatures observed in our measurements. We achieve a temperature reduction of

and using the environment temperature *T*_{0}, we obtain

with uncertainty dominated by the one of *T*_{0}. This measurement demonstrates a temperature reduction of 85%.

Lower temperatures can be achieved by lowering the noise temperature of the amplifier, *T*_{0}, increasing the *Q*-value of the resonator, or by operating with smaller traps that increase *γ*_{Be} quadratically with lower radius. In addition, by performing these demonstration measurements fully on-resonance for maximal coupling rates, the balance of heating by the resonator to cooling by the Be^{+} ions is maximally inefficient, and future cooling work will be done off-resonantly to balance the coupling rate and the temperature limit with engineered cooling sequences^{35}.

In the context of our experimental goals, this technique can be readily applied to sympathetically laser cool protons and antiprotons in the same large macroscopic traps that enable precision measurements of the charge-to-mass ratio and *g*-factor^{1,11}. In addition, while we measure the axial temperature, sideband coupling^{29} or axialization^{36} can be used to cool the radial motion of the antiproton. In measurements of nuclear magnetic moments, this will enable nearly 100% spin-flip fidelity^{11,15,16,37}, and can reduce the dominant systematic effect proportional to the particle temperature in the highest precision mass measurements^{38,39,40}. In addition, this technique can be used to cool other exotic systems such as highly charged^{23,24} or molecular ions^{17,41} and the sympathetically cooled resonator can enhance the sensitivity of dark matter searches^{13,41,42}. Ultimately, this demonstration realises a long-sought experimental technique that will enable precision experiments of any charged species at lower temperatures.

## Methods

### Equations of motion for the coupled ion-trap systems

The axial motions of the trapped proton and Be^{+} ions are described by a harmonic oscillator driven by the oscillating voltage on the trap electrodes connected to the resonant LC circuit, *V*_{LC}, and in the case of the Be^{+} ion(s) by an additional photon scattering force from the cooling laser, *F*_{L}:

*V*_{LC} is composed of voltage noise from the environment, *V*_{noise}, and the voltage arising from image currents induced by the proton and Be^{+} ions, *I*_{p} and *I*_{Be}, respectively. On resonance *ω*_{R} = *ω*_{z,p} = *ω*_{z,Be}, the impedance of the LC circuit is given by its equivalent parallel resistance, *R*_{p}, and

*I*_{noise} is the noise current from the environment, and *I*_{p} and *I*_{Be} are the induced image currents of the proton and the Be^{+} ions, respectively:

where *D*_{p,Be} are the trap-dependent effective electrode distances^{28}.

The equations of motion and the equation for the voltage in the LC circuit form a set of coupled stochastic differential equations without closed analytical solutions available. As a result, we analyse the frequency response of the system by calculating the impedance of the equivalent circuit in Fig. 1b and estimate the energy of the proton by calculating the temperatures of each component based on their energy exchange rates. Finally, we numerically integrate the differential equations in simulations that allow the comparison of FFT spectra and the visualization of the time-domain behaviour in the system.

### Impedance analysis of the equivalent circuit

The FFT spectrum in Fig. 1c results from the noise on the image-current detector \({u}_{n}^{2}=4{k}_{{\rm{B}}}{T}_{0}{\rm{R}}{\rm{e}}[Z(\omega )]\Delta f\) at effective noise temperature *T*_{0}, FFT bandwidth Δ*f*, and the impedance *Z*(*ω*) of the circuit in Fig. 1b. The lineshape of resistively cooled particles stored in a single trap based on the impedance of the equivalent circuit is well understood^{28,43}. Here, we evaluate the impedance for two independently biased ion-trap systems as:

where *k*_{L} = *R*_{L}/*R*_{p} allows for additional damping in one of the traps. The lineshapes of the individual components arise from *δ*_{i}(*ω*) = 2(*ω*−*ω*_{i})/*γ*_{i}, which are parameters proportional to the ratio of the frequency detuning (*ω*−*ω*_{i}) to the oscillator linewidth, *γ*_{i}. The index *i*∈{R, Be, p} relates to the resonator, the Be^{+} ions and the proton, respectively. In the absence of additional damping *k*_{L} = 0, the impedance simplifies to

which describes the lineshape of the data shown in Fig. 1c. Similarly, the heat maps in Fig. 2a and Fig. 2b compare the FFT spectra from experiment to ones calculated with *Z*(*ω*), and show consistent behaviour.

With laser cooling included, *R*_{L} > 0 and the dip feature of the Be^{+} ions is modified as shown in Fig. 3a. The corresponding impedance is calculated for varying *R*_{L} in Extended Data Fig. 1. In both cases, regardless of the value of *R*_{L}, the proton shorts the noise of the LC circuit on resonance. The Be^{+} ions decouple from the LC circuit as *R*_{L} reduces the fraction of noise power dissipated in the series LC circuit of the Be^{+} ions—ultimately leading to a vanishing dip signal. This decoupling effect is well known from other coupled oscillator systems^{44} and motivates the reduced coupling of the Be^{+} ions to the LC circuit, \({\mathop{\gamma }\limits^{ \sim }}_{{\rm{Be}}} < {\gamma }_{{\rm{Be}}}\).

### Temperature model

The temperature model presented here is described in ref. ^{45} and assumes that each component of the three-oscillator system consisting of the trapped proton, the trapped Be^{+} ion(s) and the resonator comes to thermal equilibrium with the rest of the system at temperatures defined by the energy exchange rates in the system.

The Be^{+} ions are damped by the resonator as well as the cooling laser, and the power transmitted by the Be^{+} ions is then written

where \({\langle \frac{d{E}_{{\rm{Be}}}}{dt}\rangle }_{{\rm{laser}}}\) is the the power dissipated by scattered photons. An identical analysis applies to the resonator, which is coupled to the environment with a coupling rate *γ*_{D} given by the width of the resonance, or the *Q*-value, and to the Be^{+} ions with a coupling rate \({\mathop{\gamma }\limits^{ \sim }}_{{\rm{Be}}}\). These relations produce the system of equations shown in the main text,

The power dissipated by the resonator while the Be^{+} ions are laser cooled can be written as

and in combination with the power dissipated by the resonator in the absence of laser cooling,

allows the reduced coupling rate to be written as

where *k* is defined by the ratio

Although *k* can, in principle, be extracted from the FFT spectrum, the extraction of individual *k* values is imprecise and *k* and *T*_{Be} are treated as constant fit parameters in Fig. 4. A more accurate determination of \({\langle \frac{d{E}_{{\rm{Be}}}}{dt}\rangle }_{{\rm{laser}}}\) can be performed by measuring the photon scattering rate^{9} and is planned for future measurements.

### Simulations and time-domain behaviour

We access the time-domain behaviour of the proton–ion–resonator system through simulations, which are performed by numerically integrating equation (9) and equation (10). By replacing *V*_{LC} = *L*_{R}\({\dot{I}}_{L}\), where *I*_{L} is the current flowing through the inductance *L*_{R}, these equations can be rewritten as

The integration is performed in time steps of Δ*t* = 1 ns for most simulations and the equivalent thermal noise \(\langle {I}_{{\rm{n}}{\rm{o}}{\rm{i}}{\rm{s}}{\rm{e}}}^{2}\rangle =4{k}_{B}{T}_{0}\Delta f/{R}_{{\rm{p}}}\) is computed in each step *n* as

where, owing to the discrete time steps, the noise bandwidth is defined as Δ*f* = 1/(2Δ*t*). *G*_{n}(*μ* = 0,*σ* = 1) is a Gaussian distribution with mean *μ* = 0 and standard deviation *σ* = 1 which is sampled every step, conserving the standard deviation of the noise while fulfilling the criterion that two subsequent values must be uncorrelated.

We implement laser cooling in the simulations by assigning a photon absorption probability to an ion in its electronic ground state at each time step. The laser, at wavelength *λ*, with wave vector *k*_{L} in the axial direction, is detuned from the centre of the transition frequency at *f*_{0} by a detuning *δ* and we assume that the linewidth of the laser is negligible compared to the transition linewidth Γ. The discrete photon absorption probability depends on the velocity of the ion due to the Doppler effect and, in the low saturation limit, can be written as

where the saturation intensity is \({I}_{{\rm{s}}{\rm{a}}{\rm{t}}}=\frac{2{\pi }^{2}\hslash \Gamma }{2{\lambda }^{3}}\) and *I* is a free parameter that is tuned to match the laser intensity in the experiment. Upon absorption of a photon the ion transitions to the excited electronic state and receives a momentum kick −*ħk*_{L}. The ion decays to the ground state via spontaneous emission with probability ΓΔ*t* and receives a momentum kick in the axial direction of *ħk*cos*θ* where the angle *θ* accounts for radial momentum of the emitted photon. The ion can also decay via stimulated emission, in which case the momentum kick is + *ħk*_{L}.

Data preparation and analysis are performed in R^{46}, while the intensive part of the calculation is performed using C++ via the Rcpp-package^{47}. We use a fourth order symplectic integrator^{48} to calculate the particle trajectories and the voltage across the RLC circuit to ensure that energy is, on average, conserved for numeric integration with more than 10^{10} steps.

In simulations of the p-Be^{+} resonator system, we apply the conditions of the experiments described in the main text to reproduce the frequency-domain behaviour in Extended Data Fig. 2, with Figs. 2a, b corresponding to the experimental results shown in the inset of Figs. 2d and 3a), respectively. The evolution of the oscillator energies with ten Be^{+} ions, *N*_{Be} = 10, is shown in Extended Data Fig. 2c. Starting from *t* = 7.5 s, a parametric drive is applied, resulting in a significant increase in the energy of the proton and the Be^{+} ions from an initial temperature of 17 K. Similarly, Extended Data Fig. 2d shows the energy exchange between a single proton, 80 Be^{+} ions, and the resonator all on resonance, where the cooling laser is applied from *t* = 10 s, resulting in rapid cooling of the Be^{+} ions and a temperature reduction of the proton.

### Axial frequency shifts

A particle in a Penning trap is subjected to shifts in the mode frequencies due to the inhomogeneity of the magnetic field and the anharmonic contributions to trapping potential^{26,32}. The magnetic field in the trap centre can be written with the lowest order corrections as

where a quadratic gradient *B*_{2} shifts the axial frequency as a function of the radial energy by

We use this effect to demonstrate the energy exchange between the heated Be^{+} ions and the proton in Fig. 2c. Here, the proton axial mode and modified cyclotron mode are sideband-coupled with a quadrupolar rf drive, so that after the sideband coupling, the proton cyclotron energy freezes out at an energy *E*_{+} = (*ν*_{+}/*ν*_{z})*E*_{z}, where *E*_{z} is the axial energy while coupling the axial mode to the excited Be^{+} ions.

timized to have a homogeneous magnetic field and is unsuited for energy measurements using equation (25) at low energy, with a temperature resolution of <0.1 mHz K^{−1}. In the the sympathetic cooling measurements presented here, we instead used the trapping potential anharmoncity that we introduced in the proton trap to determine the temperature of the trapped particle. The trapping potential can be expanded in terms of *C*_{n} coefficients^{26,31,32}, and the higher-order terms *C*_{2n}, *n* ≥ 2 shift the trap frequencies *ν*_{i} by

where *E*_{j} is the energy of a trap mode. The coefficients *C*_{n} can be written in terms of a ‘tuning ratio’, TR, defined by the ratio of the voltage applied to the central ring electrode to the voltage applied to a correction electrode, as

*D*_{n} can be calculated from the trap geometry, and the axial frequency shift due to the leading energy-dependent trap anharmonicity *C*_{4} can be written as

where ΔTR is the offset in applied tuning ratio from the ideal tuning ratio at which *C*_{4} = 0. Ultimately, the axial frequency shift as a function of TR and the axial energy *E*_{z} can be expressed as

where, for a proton stored in the proton trap, *κ*_{D4} = 45.4ΔTR Hz K^{−1}. This effect is used to determine the change of the proton axial temperature while the resonator is cooled with the laser-cooled Be^{+} ions, and is the underlying method for the data shown in Figs. 3 and 4.

Temperature measurements using this method are limited by the determination of *T*_{0} to the ~ 1 K level. We have previously performed higher precision temperature measurements using a dedicated, spatially distant trap with a ferromagnetic ring electrode that uses the shift of equation (25) to obtain a cyclotron energy resolution of up to 80 Hz K^{−1} (refs. ^{15,49}) and have developed a similar trap to reach 10 mK temperature resolution in future cooling measurements.

We also note that equation (28) and equation (29) cause the axial dip to spread out during an FFT averaging window and decrease the dip signal-to-noise ratio. This is reflected in the increased uncertainties at larger ΔTR that can be seen in Fig. 3b.

### Parasitic drive heating

The demonstration of remote energy exchange was performed by exciting a small cloud of Be^{+} with an rf drive at twice the resonance frequency, 2*ν*_{0}, with results presented in Fig. 2d. To confirm that the proton is excited only by the resonantly coupled ions, we performed a series of background measurements. These control measurements show that the proton is excited only when resonant with the Be^{+} ions and that the proton is unaffected by the excitation drive when the Be^{+} ions are detuned (Extended Data Fig. 3).

Although the measurements presented in Fig. 2d are performed by interleaving the on- and off-resonant configurations, we additionally analysed the temperature of the proton in the presence of the drive. During these measurements, we transferred the axial energy to the modified cyclotron mode and measured the resulting energy dependent axial frequency shift, described in the main text. We measured an axial temperature of

in the presence of the drive, and an axial temperature of

in the absence of the drive, where the error comes from the fit uncertainty of the frequency scatter distribution. As a result, we constrain the possible increase in axial temperature due to the excitation drive to no more than a factor of two. From the spectra shown in Extended Data Fig. 3, we see that the signal-to-noise ratio of the proton dip is unaffected by the drive and conclude that an increase in axial temperature would come not from direct coupling of the excitation drive to the proton but via an increase in the equivalent noise temperature of the resonator. We further note that the resonator temperature of 17.0 K given in the main text comes from the weighted mean of several temperature measurements performed with several methods. Importantly, residual heating due to the drive is far lower than the energy scatter shown in Fig. 2d; approximately 40,000 K k_{B}^{−1}.

## Data availability

The datasets generated and/or analysed during this study are available from the corresponding authors on request. Source data are provided with this paper.

## Code availability

The code used during this study is available from the corresponding authors on reasonable request.

## References

- 1.
Smorra, C. et al. BASE–the baryon antibaryon symmetry experiment.

*Eur. Phys. J. Spec. Top.***224**, 3055–3108 (2015). - 2.
Sturm, S. et al. The ALPHATRAP experiment.

*Eur. Phys. J. Spec. Top.***227**, 1425–1491 (2019). - 3.
Huntemann, N., Sanner, C., Lipphardt, B., Tamm C.H.R. & Peik, E. Single-ion atomic clock with 3 × 10

^{−18}systematic uncertainty.*Phys. Rev. Lett.***116**, 063001 (2016). - 4.
Brewer, S. M. et al.

^{27}Al^{+}quantum-logic clock with a systematic uncertainty below 10^{−18}.*Phys. Rev. Lett.***123**, 033201 (2019). - 5.
Kotler, S., Simmonds, R. W., Leibfried, D. & Wineland, D. J. Hybrid quantum systems with trapped charged particles.

*Phys. Rev. A***95**, 022327 (2017). - 6.
Bruzewicz, C. D., Chiaverini, J., McConnell, R. & Sage, J. M. Trapped-ion quantum computing: progress and challenges.

*Appl. Phys. Rev.***6**, 021314 (2019). - 7.
Brown, K. R. et al. Coupled quantized mechanical oscillators.

*Nature***471**, 196–199 (2011). - 8.
Harlander, M., Lechner, R., Brownnutt, M., Blatt, R. & Hänsel, W. Trapped-ion antennae for the transmission of quantum information.

*Nature***471**, 200–203 (2011). - 9.
Heinzen, D. J. & Wineland, D. J. Quantum-limited cooling and detection of radio-frequency oscillations by laser-cooled ions.

*Phys. Rev. A*,**42**, 2977 (1990). - 10.
Daniilidis, N., Lee, T., Clark, R., Narayanan, S. & Haffner, H. Wiring up trapped ions to study aspects of quantum information.

*J. Phys. B***42**, 154012 (2009). - 11.
Bohman, M. et al. Sympathetic cooling of protons and antiprotons with a common end-cap Penning trap.

*J. Mod. Opt.***65**, 568–576 (2018). - 12.
Smorra, C. et al. Direct limits on the interaction of antiprotons with axionlike dark matter.

*Nature***575**, 310–314 (2019). - 13.
Devlin, J. A. et al. Constraints on the coupling between axionlike dark matter and photons using an antiproton superconducting tuned detection circuit in a cryogenic Penning trap.

*Phys. Rev. Lett.***126**, 041301 (2021). - 14.
Ulmer, S. et al. High-precision comparison of the antiproton-to-proton charge-to-mass ratio.

*Nature***524**, 196–199 (2015). - 15.
Schneider, G. et al. Double-trap measurement of the proton magnetic moment at 0.3 parts per billion precision.

*Science***358**, 1081–1084 (2017). - 16.
Smorra, C. et al. A parts-per-billion measurement of the antiproton magnetic moment.

*Nature***550**, 371–374 (2017). - 17.
Myers, E. G. CPT tests with the antihydrogen molecular ion.

*Phys. Rev. A***98**, 010101 (2018). - 18.
Yzombard, P., Hamamda, M., Gerber, S., Doser, M. & Comparat, D. Laser cooling of molecular anions.

*Phys. Rev. Lett.***114**, 213001 (2015). - 19.
Tang, R. et al. Candidate for laser cooling of a negative ion: high-resolution photoelectron imaging of Th

^{−}.*Phys. Rev. Lett.***123**, 203002 (2019). - 20.
Cerchiari, G., Yzombard, P. & Kellerbauer, A. Laser-assisted evaporative cooling of anions.

*Phys. Rev. Lett.***123**, 103201 (2019). - 21.
Wineland, D. J. et al. Experimental issues in coherent quantum-state manipulation of trapped atomic ions.

*J. Res. Natl Inst. Stand. Technol.***103**, 259–328 (1998). - 22.
Niemann, M., Paschke, A. -G., Dubielzig, T., Ulmer, S. & Ospelkaus, C. CPT test with (anti)proton magnetic moments based on quantum logic cooling and readout. In

*Proc. Sixth Meeting on CPT and Lorentz Symmetry*(Ed. Kostelecký, V. A.) 41–44 (World Scientific, 2014); https://www.worldscientific.com/doi/abs/10.1142/9789814566438_0011 - 23.
Egl, A. et al. Application of the continuous Stern-Gerlach effect for laser spectroscopy of the 40Ar

^{13+}fine structure in a Penning trap.*Phys. Rev. Lett.***123**, 123001 (2019). - 24.
Micke, P. et al. Coherent laser spectroscopy of highly charged ions using quantum logic.

*Nature***578**,60–65 (2020). - 25.
Gutiérrez, M. J. et al. Dynamics of an unbalanced two-ion crystal in a Penning trap for application in optical mass spectrometry.

*Phys. Rev. A***100**, 063415 (2019). - 26.
Brown, L. S. & Gabrielse, G. Geonium theory: physics of a single electron or ion in a Penning trap.

*Rev. Mod. Phys.***58**, 233–311 (1986). - 27.
Nagahama, H. et al. Highly sensitive superconducting circuits at 700 kHz with tunable quality factors for image-current detection of single trapped antiprotons.

*Rev. Sci. Instrum.***87**, 113305 (2016). - 28.
Wineland, D. J. & Dehmelt, H. G. Principles of the stored ion calorimeter.

*J. Appl. Phys.***46**, 919–930 (1975). - 29.
Cornell, E. A., Weisskoff, R. M., Boyce, K. R. & Pritchard, D. E. Mode coupling in a Penning trap: π pulses and a classical avoided crossing.

*Phys. Rev. A***41**, 312– 315 (1990). - 30.
Dehmelt, H. Continuous Stern-Gerlach effect: principle and idealized apparatus.

*Proc. Natl Acad. Sci.***83**, 2291–2294 (1986). - 31.
Gabrielse, G., Haarsma, L. & Rolston, S. L. Open-endcap Penning traps for high precision experiments.

*Int. J. Mass Spectrom. Ion Processes***88**, 319–332 (1989). - 32.
Ketter, J., Eronen, T., Höcker, M., Streubel, S. & Blaum, K. First-order perturbative calculation of the frequency-shifts caused by static cylindrically-symmetric electric and magnetic imperfections of a Penning trap.

*Int. J. Mass spectrom.***358**, 1–16 (2014). - 33.
Daniilidis, N. & Häffner, H. Quantum interfaces between atomic and solid-state systems.

*Annu. Rev. Condens. Matter Phys.***4**, 83–112 (2013). - 34.
Möller, S. A.

*How to Couple Trapped Ions to Superconducting Resonators: Towards Hybrid Quantum Devices*. PhD thesis, Univ. California (2016). - 35.
Tu, B. et al. Tank-circuit assisted coupling method for sympathetic laser cooling.

*Adv. Quantum Technol.***4**, 2100029 (2021). - 36.
Powell, H. F., Segal, D. M. & Thompson, R. C. Axialization of laser cooled magnesium ions in a Penning trap.

*Phys. Rev. Lett.***89**, 093003 (2002). - 37.
Schneider, A. et al. A novel Penning-trap design for the high-precision measurement of the

_{3}He^{2+}nuclear magnetic moment.*Ann. Phys.***531**, 1800485 (2019). - 38.
Heiße, F. et al. High-precision measurement of the proton’s atomic mass.

*Phys. Rev. Lett.***119**, 033001 (2017). - 39.
Smith, J. A., Hamzeloui, S., Fink, D. J. & Myers, E. G. Rotational energy as mass in H3

^{+}and lower limits on the atomic masses of d and^{3}He.*Phys. Rev. Lett.***120**, 143002 (2018). - 40.
Rau, S. et al. Penning trap mass measurements of the deuteron and the HD+ molecular ion.

*Nature***585**, 43–47 (2020). - 41.
Cairncross, W. B. et al. Precision measurement of the electron’s electric dipole moment using trapped molecular ions.

*Phys. Rev. Lett.***119**, 153001 (2017). - 42.
Roussy, T. S., Palken, D. A., Cairncross, W. B. et al. Experimental constraint on axion-like particle coupling over seven orders of magnitude in mass. Preprint at https://arxiv.org/abs/2006.15787 (2020).

- 43.
Feng, X., Charlton, M., Holzscheiter, M., Lewis, R.A. & Yamazaki, Y. Tank circuit model applied to particles in a Penning trap.

*J. Appl. Phys.***79**, 8–13 (1996). - 44.
Zeuthen, E.

*Extending Cavity-Mechanical Cooling via a Hot LC Electrical Circuit Framework for Electromechanical Coupling Calculations*. MSc thesis, Niels Bohr Institute, Univ. Copenhagen (2012). - 45.
Bohman, M.

*Sympathetic Cooling of a Proton with Resonant Image Current Coupling*. PhD thesis, Ruperto-Carola Univ. Heidelberg (2020). - 46.
R Core Team.

*R: a Language and Environment for Statistical Computing.*https://www.R-project.org/ (R Foundation for Statistical Computing, 2020). - 47.
Eddelbuettel, D. & François, R. Rcpp: seamless R and C++ integration.

*J. Stat. Softw.***40**, 1–18 (2011). - 48.
Yoshida, H. Construction of higher order symplectic integrators.

*Phys. Lett. A***150**, 262–268 (1990). - 49.
Borchert, M. J. et al. Measurement of ultralow heating rates of a single antiproton in a cryogenic Penning trap.

*Phys. Rev. Lett.***122**, 043201 (2019).

## Acknowledgements

This study comprises parts of the PhD thesis work of M.B. We acknowledge the contributions of G. L. Schneider and N. Schön to the design and construction of the Penning-trap system. We thank S. Sturm for helpful discussions regarding the cooling method presented here and acknowledge similar developments toward cooling highly charged ions in the ALPHATRAP collaboration. We acknowledge financial support from the RIKEN Chief Scientist Program, RIKEN Pioneering Project Funding, the RIKEN JRA Program, the Max-Planck Society, the Helmholtz-Gemeinschaft, the DFG through SFB 1227 ‘DQ-mat’, the European Union (Marie Skłodowska-Curie grant agreement number 721559), the European Research Council under the European Union’s Horizon 2020 research and innovation programme (Grant agreement numbers 832848-FunI and 852818-STEP) and the Max-Planck-RIKEN-PTB Center for Time, Constants and Fundamental Symmetries.

## Funding

Open access funding provided by Max Planck Society.

## Author information

### Affiliations

### Consortia

### Contributions

M.B., A.M., S.U., J.W. and M.W. designed the experimental apparatus. M.B., A.M. and M.W. assembled the trap and laser systems. M.B., V.G., C.S. and M.W. contributed to the experiment run. M.B. and C.S. implemented the methods, and recorded and evaluated the experimental data. C.W. developed the simulation code and provided the simulation results. M.B., C.S., K.B. and S.U. prepared the manuscript, which was discussed and approved by all authors.

### Corresponding author

## Ethics declarations

### Competing interests

The authors declare no competing interests.

## Additional information

**Peer review information** *Nature* thanks Manas Mukherjee, Richard Thompson and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.

**Publisher’s note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Extended data figures and tables

### Extended Data Fig. 1 Calculated impedance of the equivalent circuit.

The normalised impedance \({\rm{R}}{\rm{e}}[Z(\omega )]/{R}_{{\rm{p}}}\) is plotted for different values of *R*_{L}, where no damping corresponds to *R*_{L} = 0, and no Be^{+} ions to *R*_{L} → ∞. Here, *γ*_{R} = 2*γ*_{Be} = 20*γ*_{p}, corresponding in experiment to about 30 Be^{+} ions and a single proton.

### Extended Data Fig. 2 Simulation results.

**a)** A computed FFT spectrum is shown simulating the experimental conditions of Fig. 2d) in the main text. **b)** Representative time domain behaviour for these measurements is shown where the excitation drive is applied at time at time *t* = 10s. **c)** A computed FFT spectrum is shown simulating the experimental conditions of Fig. 3a) in the main text. **d)** Representative time domain behaviour for these measurements is shown where the cooling laser is applied at time *t* = 10s.

### Extended Data Fig. 3 Excitation drive background measurements.

**a)** An FFT spectrum while the excitation drive is off, the proton is on resonance with the resonator and the Be^{+} ions are off resonance. **b)** An FFT spectrum while the excitation drive is on, the proton is off resonance with the resonator and the Be^{+} ions are on resonance. **c)** An FFT spectrum while the excitation drive is on, the proton is on resonance with the resonator and the Be^{+} ions are on resonance. **d)** An FFT spectrum while the excitation drive is on, the proton is on resonance with the resonator and the Be^{+} ions are off resonance.

## Rights and permissions

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

## About this article

### Cite this article

Bohman, M., Grunhofer, V., Smorra, C. *et al.* Sympathetic cooling of a trapped proton mediated by an LC circuit.
*Nature* **596, **514–518 (2021). https://doi.org/10.1038/s41586-021-03784-w

Received:

Accepted:

Published:

Issue Date:

DOI: https://doi.org/10.1038/s41586-021-03784-w

## Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.