A Single-Ion Reservoir as a High-Sensitive Sensor of Electric Signals

A single-ion reservoir has been tested, and characterized in order to be used as a highly sensitive optical detector of electric signals arriving at the trapping electrodes. Our system consists of a single laser-cooled 40Ca+ ion stored in a Paul trap with rotational symmetry. The performance is observed through the axial motion of the ion, which is equivalent to an underdamped and forced oscillator. Thus, the results can be projected also to Penning traps. We have found that, for an ion oscillator temperature T axial ≲ 10 mK in the forced-frequency range ω z = 2π × (80,200 kHz), the reservoir is sensitive to a time-varying electric field equivalent to an electric force of 5.3(2) neV/μm, for a measured quality factor Q = 3875(45), and a decay time constant γ z = 88(2) s−1. This method can be applied to measure optically the strength of an oscillating field or induced (driven) charge in this frequency range within times of tens of milliseconds. Furthermore the ion reservoir has been proven to be sensitive to electrostatic forces by measuring the ion displacement. Since the heating rate is below 0.3 μeV/s, this reservoir might be used as optical detector for any ion or bunch of charged particles stored in an adjacent trap.


Results
The experiments have been carried out using the open-ring Paul trap described in ref. 22, but now driving the trapping field with a radiofrequency ω RF = 2π × 1.47 MHz and applying a compensation DC voltage (U DC ) to some of the trap electrodes. The 40 Ca + ion is created by photoionization at a background pressure of ≈10 −10 mbar. The level structure of 40 Ca + shown in Fig. 1 contains two electric dipole-allowed transitions, which will contribute to the dynamics. Three laser beams have been used for this experiment as shown schematically in the right-hand side of Fig. 1, with tunable wavelengths around 397 nm (B1 and B2), and 866 nm (R1). B1 and B2 are tuned to drive the cooling transition 4 s 2 S 1/2 → 4 s 2 P 1/2 in the radial and axial directions, respectively. Moreover, R1, which is only applied in the radial direction, is needed for the transition 3d 2 D 3/2 → 4 s 2 P 1/2 to pump the metastable state 3d 2 D 3/2 , unavoidably populated with a probability of 7%. The natural linewidths of the 3d 2 D 3/2 → 4 s 2 P 1/2 and 4 s 2 S 1/2 → 4 s 2 P 1/2 transitions are Γ R = 2π × 1.35 MHz, and Γ B = 2π × 21.58 MHz, respectively. The fluorescence light from the cooling transition (397 nm) is collimated using a commercial system from Thorlabs (MVL12 × 12Z), which provides a magnification of .
6 75 (5) in the focal plane of an EMCCD (Electron Multiplier Charged Couple Device) from Andor (IXON3). The EMCCD sensor is made of 512 × 512 pixels. The size of each pixel is 16 μm × 16 μm. An interference filter is located in front of the EMCCD to collect only photons with wavelength around 397 nm. The effective pixel size has been obtained through the ratio between the minimum distance, due to Coulomb interaction between two trapped ions πε ω ∆ = z e m ( /2 ) z 2 0 2 1/3 23 , and the pixel distance between the centers of the distributions of the ions projected on the axial direction. The radiofrequency voltage (V RF ) has been varied from ≈520 to 1430 V pp . The trap frequencies are given by where q u = f(V RF , ω RF , u 0 ) and a u = f(U DC , ω RF , u 0 ) with u 0 representing the characteristic dimensions of the trap 24 . The trap frequencies in the axial direction ω z were measured to be in the range ≈2π × 60-250 kHz, and the corresponding trap frequencies in the radial direction ω r are a factor of two smaller provided no DC potential is applied. q z in the axial direction varies from 0.18 to 0.49 for different values of V RF and a z = −0.0088 for the measurements carried out with ω π ×  2 108 z kHz(q z ≃ 0.25). This frequency is close to the value originally proposed 20 and within a range where the adiabatic approximation is valid 24 . For the sensitivity measurements to electrostatic potentials, a z = −0.0184 and ω π ×  2 80 z kHz for the same q z . Our trapped ion can be considered, to a good approximation, a harmonic oscillator in all three directions. The effect of the described laser configuration will be modeled as a damping force acting on the direction of the cooling lasers, which is the standard treatment in the theory of Doppler cooling 8,25,26 . As a consequence, the dynamics of the ion in the axial direction is given by a damped harmonic oscillator, with equation of motion where ρ z and γ z are, respectively, the position of the ion and the damping coefficient in the axial direction. A term given by αρ z 3 , might be added to the left part of Eq. (2) to account for possible deviations of the quadrupole potential due to the DC voltages applied for compensation.
When the ion reaches a temperature given by the Doppler limit, the cooling effect of the lasers stops and the ion reaches a steady state. In a three-level system, this limit temperature will depend on the intensities and frequencies of the lasers driving the cooling and pumping transitions. In this steady state, the ion is sensitive to external oscillating electric fields acting on the electrodes of the trap. This signal could be originated by one or Figure 1. Atomic levels and lasers configuration to perform Doppler cooling on 40 Ca + . 40 Ca + atomic levels for Doppler cooling and schematic representation of the interaction between a trapped ion and laser beams in the radial and axial potential wells of the open-ring trap 22 . The double-ended thick arrow depicts the application of an external dipole field with amplitude V dip equivalent to a force F e /m, where m is the mass of the ion. more ions in a second trap. However, in all experiments carried out here, this electric field is artificially generated. This will test the sensitivity of our oscillator, which ultimately would act as a weighing device to measure the mass of one or many ions in the other trap 21 .
In Fig. 2a, we show the distribution of photons emitted by the ion as it oscillates in the axial direction while in its steady state. We find that the width of the photon distribution is related to the variance of the projected Gaussian-like distribution, and in this way, to the amplitude of the oscillation ρ z,max . To explain this, we assume that the photons emitted by the ion at each point in the axial line are collected at the CCD camera with a Lorentzian distribution. The width of this distribution is obtained from a χ 2 fit to the profiles of the fluorescence measurements for different oscillation amplitudes when applying a driving force as shown in Fig. 2b for ρ z,max = 21 μm. The function for the χ 2 fit is built by the convolution of the probability distribution of the oscillator in the axial direction P(ρ z ), and the Lorentzian distribution, considering on one hand ρ z,max , and on the other hand the Lorentzian width, as free parameters. This demonstrates that we can optically detect the effect of a time-varying electric signal applied to the electrodes of the trap while we obtain the optical response of our system. The uncertainties in ρ z,max are the quadratic sum of the statistical uncertainty (from one measurement), and the systematic uncertainties, arising from the election of the position of the central pixel, from the resolution and magnification of the imaging system, and from the step size in ρ z,max . We will now systematically use this detection method to calibrate the sensitivity of our device. Additionally, we can associate an axial temperature to the ion in equilibrium according to ω ρ ω σ 8 . In this manner, we measure an interval of Doppler limit temperatures for the trap frequencies tested in this work from 8.1(3.8) to 10.7(2.9) mK. The same can be done for the radial direction, provided one knows the radial frequency. In order to characterize the sensitivity of our single-ion reservoir, we measure its response to dipolar fields of different amplitudes and frequencies, as shown in Fig. 3. In the steady state, the amplitude of the driven-damped harmonic oscillator is given by  where F e and ω dip are, respectively, the amplitude and the frequency of the harmonic driving force. There is no analytical expression for ρ z,max if one introduces the anharmonicity term αρ z 3 . In such case however, it is possible to perform a numerical fit. Fitting the experimental data (V 2 in Fig. 3) with Eq. (3), we obtain F e,2 = 5.3(2) neV/μm, and a damping coefficient γ z,2 = 88(2) Hz, to which a quality factor = = Q 3875(45) w g 2 z z can be associated. Introducing a term to account for an anharmonicity and performing a numerical fit, one obtains values consistent with the latter. In case the amplitude of the external field is V 3 , we obtained F e,3 = 3.6(1) neV/μm. The ratio of the electric field the ion sees F e,3 /F e,2 = 0.68(5) equals the ratio of the amplitudes applied to the electrodes of V 3 /V 2 = 0.7, thus by reducing this amplitude, it is possible to infer an electric field on the ion of 0.1F e,2 , when a radiofrequency field with an amplitude V 4 = 0.1V 2 is applied. The resonance curve is shown for a frequency of 2π × 202 kHz in the lower-right size of Fig. 3. The signal-to-noise ratio is ~2.
The application of this ion reservoir to measure ion-motional frequencies is subject to the possibility of observing fluorescence signals in short periods of time. Figure 4a shows the image collected for different delay times after the ion has been excited. Each image is collected during 20 ms, which gives a good signal-to-noise ratio  (24) Hz, respectively, from the analytical fit. Middle right: sensitivity of the ion reservoir to electric forces in neV/μm. * shows the field generated by a single-charged ion oscillating in a second trap, and ** the field generated by e.g. 50 antiprotons or 50 single-charged superheavy-element ions. The lower right side of the figure, shows the variance of the ion distribution as a function of (ω dip − ω z )/2π for an amplitude V 4 = 0.1V 2 .
The cyan-colored area shows the noise level, which is defined by the average variances of the distributions off resonance. Since the invariance for each data point is obtained directly from the Gaussian fit to the projection of the fluorescence into the axial direction, the uncertainty of each data point is the quadratic sum of the statistical and systematic uncertainties, the first provided by the fit and the second by the observed variations.
Scientific RepoRts | 7: 8336 | DOI:10.1038/s41598-017-08782-5 without affecting the time resolution. The observed decay rates applying different amplitudes of the driving field are in agreement with the ones computed from the fitting of Fig. 3 (γ z,2 ).
We have demonstrated that the resonance of the ion is clearly visible even in cases when the external field amplitude is the lower available in our laboratory (V 4 in Fig. 3). This field should be equivalent to a force of 0.53 neV/μm, in contrast to the force that a single ion stored in an adjacent trap would originate, which is of F e = 0.01 neV/μm, considering an electronic charge state of 1 + and an oscillation amplitude ρ z,max = 100 μm in the micro-trap described in ref. 27. The minimum temperature reached in these measurements corresponds to a mean number of ~1900 phonons in the axial direction (E phonon = 4.5 neV for ω z = 2π × 108 kHz). Using the scheme based on fluorescence photons depicted in Fig. 4, one needs to couple the motion of the 40 Ca + ion with the motion of the ion to be studied and observe the evolution of the system in the absence of lasers, i.e. without cooling 20 . The left part of Fig. 4b shows the image in the sensor and the axial projection when an excitation has been previously applied, in the absence of lasers, during 50 ms. From the Gaussian fit to the axial distribution, and considering γ z = 88(2) Hz, one obtains, after excitation, an initial ion oscillation amplitude ρ z,max = 15.6 μm (T z = 230 mK). This value, compared with the result when no excitation is applied, while there is no interaction with the laser beams, yields a rate for the increase in energy of the order of 10 μeV/s, much larger than the measured heating rate, below 0.3 μeV/s. This upper limit for the heating rate has been obtained by measuring the response of the laser-cooled trapped ion to different amplitudes of applied white noise (V noise ) and fitting linearly the axial temperature with the function given by where ζ · V 2 is the heating rate in units of J s −1 , which is multiplied by the square of the amplitude of the noise signal and K is a constant from which we recover the Doppler limit, in the ideal case where V noise = 0. Then, the square of the variance σ ζ ω . The upper limit obtained is just the value corresponding to the smallest amplitude applied.
For completeness, we have also studied the response of the system to DC fields. The electrostatic force in the axial direction is given by ω = F m z z z 2 . Figure 5 shows the electrostatic force as a function of the ion displacement, which is defined as the center of the Gaussian distributions resulting from the projections of the fluorescence images in the axial direction (see the insets). From the fit, one can get F z /z = 0.10455(0) neV/μm/nm. The center of the distributions is determined with an uncertainty of 67 nm, which defines the uncertainty in the measurement of the electrostatic force of 7 neV/μm. The minimum force measured from the ion displacement in the axial direction shown in Fig. 5, marked with (a) is 42(7) neV/μm. Smaller displacements have been obtained in the radial direction but the radial frequency was not measured.

Discussion
Summarizing, this article experimentally realizes an atomic ion reservoir composed of a trapped ion as a high-sensitive measurement device. The Doppler-cooled ion will allow one to detect weak electric fields as well as the motional frequencies of charged particles in adjacent traps. Our experiment paves the way towards more sensitive mass and frequency spectroscopy with charged particles. This method can be extended to a broader frequency range, where ground-state cooling on the same ion can be applied. This might allow increasing the sensitivity of the optical method in the DC and radiofrequency regimes, as well as controlling and measuring the phase of the axial motion. Figure 5. Electrostatic force sensing. Electrostatic force exerted on the laser-cooled ion as a function of the ion displacement. The potential applied to a specific electrode was varied from 0.5 V in (a) to 4 V in (b). The axial frequency was 80 kHz. The inset shows the images collected for the cases (a) and (b). Each of the images was taken during 35 seconds. The red line is a linear fit to the data points (black-solid circles), and the blue area delimit the uncertainties. The white-solid circles are obtained after the fit.