This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Dynamic compensation of stray electric fields in an ion trap using machine learning and adaptive algorithm

Moji Ghadimi1 [email protected]    Alexander Zappacosta1    Jordan Scarabel1    Kenji Shimizu1    Erik W Streed1,2    Mirko Lobino1,3 1Center for Quantum Dynamics, Griffith University, Nathan QLD Australia
2Institute for Glycomics, Griffith University, Southport QLD Australia
3Queensland Micro Nanotechnology Centre, Nathan QLD Australia
Abstract

Surface ion traps are among the most promising technologies for scaling up quantum computing machines, but their complicated multi-electrode geometry can make some tasks, including compensation for stray electric fields, challenging both at the level of modeling and of practical implementation. Here we demonstrate the compensation of stray electric fields using a gradient descent algorithm and a machine learning technique, which trained a deep learning network. We show automated dynamical compensation tested against induced electric charging from UV laser light hitting the chip trap surface. The results show improvement in compensation using gradient descent and the machine learner over manual compensation. This improvement is inferred from an increase of the fluorescence rate of 78% and 96% respectively, for a trapped 171Yb+ ion driven by a laser tuned to -7.8 MHz of the 2S21/2{}_{1/2}\leftrightarrow^{2}P1/2 Doppler cooling transition at 369.5 nm.

preprint: APS/123-QED

I Introduction

Trapped ions are one of the most promising candidates for implementing scalable quantum information processing [1, 2, 3, 4]. Traditionally macroscopic rod or needle traps were used for electrically trapping laser cooled ions. These Paul traps have a deep potential well that enables the trapping of multiple ions in a chain, but their distant electrodes lack the fine spatial control of electric fields necessary to efficiently split and combine small ion crystals [5]. One of the proposed architectures for scaling up the number of trapped ions in a quantum processor is the use of microfabricated surface chip traps with multiple DC electrodes that are able to manipulate multiple ions individually [6, 7, 8, 9, 10].

In these chip traps the ions are confined to a potential formed by the node of an oscillating RF electric field in two dimensions while a DC field generated from multiple electrodes provide a finely spatially adjustable potential in the third dimension. It is well-known in Paul trapping that stray DC electric fields push the ion off the RF node and induce micro-motion that degrades or outright prevents effective confinement and laser-cooling [11]. The solution to this problem is to use the DC electrodes to generate an opposite DC electric field that compensates for these stray fields [9].

Micro-motion amplitude can be measured and dealt with in multiple ways [12, 13, 14, 9, 15]. A simple and quick proxy for small amounts of micromotion is to measure the ion’s fluorescence under laser cooling near the Doppler limit. In fact for low magnitudes of micro-motion, if a laser is red-detuned near the natural linewidth of the atom, the ion fluorescence rate increases as the micro-motion amplitude decreases (see Section II and Fig. 2). While stray electric fields can be readily compensated in simple fixed trap geometries [16], this task becomes more challenging for multi-electrode designs [17] and in proximity of surfaces which are vulnerable to laser induced charging [18].

Automation, optimization, and machine learning have been used to improve different manual tasks in atom and ion traps [19, 20, 21] and they are also useful tools for optimizing the individual electrode voltages which generate the trapping electric field. In this work we show how the voltages of an array of electrodes in a surface chip trap can be optimized using a fully automated, machine learning driven process to minimize micromotion and maximize florescence. The chip trap in our experiment has 44 segmented DC electrodes and 2 RF rails, as shown in Fig. 1. The trap incorporates nearly aberration free diffractive mirrors directly fabricated onto the central ground electrode for efficient fluorescence collection from the ion [22]. A Gradient descent (ADAM [23]) and a machine learning algorithm (MLOOP [20]) were compared on the basis of the versatility and time taken to find the optimal compensation and highest fluorescence rate of the ion. These methods were applied to a trapped 171Yb+ ion and compared to a manual optimization performed by an experienced operator.

With ADAM we were able to improve the fluorescence by 78% starting from voltages already optimized by manual adjustment, while MLOOP achieved a 96% improvement. Subsequently we tested the versatility and adaption of this procedure by deliberately charging the trap with UV light to drop the fluorescence rate by around 35% and compensating back to the optimal with ADAM. This was not tested on MLOOP as it was observed to be quite sensitive to the noise in the fluorescence signal and did not make a solid candidate for final optimization of the artificially charged trap.

II Experimental Setup

Refer to caption
Figure 1: a: Transitions for laser cooling 171Yb+. The ion is laser cooled on the nearly closed 396.5 nm 2S1/2F=1 to 2P1/2F=0 transition. The small amount of off-resonant scattering into the 2S1/2 F=0 state is repumped with a 14.7 GHz sideband on the primary cooling laser. Decays from 2P1/2 can also populate the 2D3/2 level, which are repumped using a 935 nm laser. b: Schematics of the ion trap chip and the optimization devices. RF and DC electrodes are used to trap the ion and shuttle it above the diffracting mirror that collects photons from the ion and sends it to photo multiplayer tube (PMT). The PMT counts are processed by the optimiser and DC electrode voltages and cooling laser position are updated using a digital to analogue converter (DAC).

The ion trap used in these experiments and described in [24, 22], consists of a planar rectangular chip of 1400 μ\mum length and 800 μ\mum width. It has a set of aluminum electrodes patterned on top and is schematically shown in Fig. 1. The ion is trapped 60 microns above the surface, between the two RF electrode rails that extend for the full length of chip, and 44 DC electrodes on the sides. The DC electrodes are used to create a reconfigurable trapping potential along the length of the trap, compensating stray electric fields, and enabling ion crystals merging, separation, and shuttling.

An atomic oven underneath one end of the trap generates a beam of neutral Yb atoms which passes through a slit in the chip. The 171Yb is first excited by an isotopically selective 399 nm laser [25], and subsequently non-resonantly ionized by the a 369.5 nm laser. The ion is then Doppler cooled using the same 369.5 nm laser tuned nearly resonant with 2S1/2 F=12\leftrightarrow^{2}P1/2 F=0 transition. Occasionally off-resonant scattering from the 2P1/2 F=1 state will populate the dark 2S1/2 F=0 ground state, which is repumped by a small amount of light from a 14.7 GHz sideband added to the primary 369.5 nm cooling laser. There is also a probability of 0.5% for the atom to decay from the 2P1/2 into the meta-stable 2D3/2 state. In this case a 935 nm laser repumps the ions back into the cooling cycle. Figure 1 depicts the cooling and repumping transitions.

After initial trapping near the oven slit region, the ion is shuttled along the length of the trap by properly controlling the voltages of the array of DC electrodes with a 12-bit National Instruments PXI-6713 DAQ (output doubled to +/-20V with 10 mV resolution). For our experiments, the ion was shuttled above the surface of an integrated diffracting micro-mirror patterned on the ground electrode with a focal length of 60 µm, equal to the height of the ion above the trap, to improve florescence collection and coupling into a single mode fiber [22].

Refer to caption
Figure 2: Normalised florescence as a function of the dimensionless micromotion parameter β\beta (Eq. 1). Here δL\delta_{L}, the laser frequency detuning, is -7.8 MHz, γ\gamma, the natural linewidth is 20 MHz and Ω\Omega, the RF trap frequency, is 25.5 MHz. The 2P1/2 state population and therefore the fluorescence rate drops monotonically as β\beta increases as long as its value is not too large (smaller than 2.5). Our operating regime is far below this threshold.

For efficient Doppler cooling and subsequent quantum information processing experiments, it is important to minimize any micro-motion induced by stray electric fields. This is done by tuning the voltages of the DC electrodes so that they generate an opposing electric field that compensates for the stray fields. To generate a compensation field in the RF confinement plane xyxy, we use two waveforms (array of DC voltages for individual electrodes) which were calculated for generating an electric field of 100 V/m in the xx (horizontal) and yy (vertical) directions. Electrode voltages of 1-1 to +1+1 V for the xx direction and 2-2 to +2+2 V for the yy direction are needed to create electric field of 100 V/m at the location of the ion. These sets of voltages are multiplied by arbitrary weights wxw_{x} and wyw_{y} to compensate arbitrary fields in the RF trapping plane. In practice this process has imperfections and cannot fully compensate for the stray electric fields [17]. Instead, individual electrodes need to be tuned independently for optimal compensation. This optimum compensation is position dependent and there might be multiple solutions for a multi-electrode trap.

Two other DC waveforms are used and their relative weights tuned: a harmonic waveform, and RF plane trap axis rotation waveform. The former creates a harmonic potential to trap the ion in different locations along the chip length, while the latter rotates the trap’s secular confinement axes to lift the degeneracy and improve cooling [26].

We infer the magnitude of the stray electric field from the fluorescence rate of the ion when the cooling laser is red-detuned near half its natural linewidth for optimal cooling. When driving an ideal cold two-level atom with this laser, the fluorescence is proportional to the population of excited state, PeP_{e}, that can be calculated from the optical Bloch equation [27]. The relationship between PeP_{e} and the strength of the stray electric field is given by [15]:

Pe(β)=Cm=Jm2(β)(δLγ+mΩγ)2+14P_{e}(\beta)=C\sum_{m=-\infty}^{\infty}{\frac{J_{m}^{2}(\beta)}{(\frac{\delta_{L}}{\gamma}+m\frac{\Omega}{\gamma})^{2}+\frac{1}{4}}} (1)

where JmJ_{m} is the mth order Bessel function of the first kind, δL=ωatomωlaser\delta_{L}=\omega_{atom}-\omega_{laser} is the laser frequency detuning, γ\gamma is the natural linewidth and Ω\Omega is the RF trap frequency. β\beta is a dimensionless measure of the ion’s coherently driven motion from the stray electric fields and RF phase imbalances (AC micromotion) [15]. CC depends on the strength of the coupling between the levels that is a constant here since we keep laser intensity and direction fixed. Figure 2 shows a graph of Pe(β)/Pe(0)P_{e}(\beta)/P_{e}(0) as a function of β\beta. The graph shows that at our detuning (-7.8 MHz) the population decreases monotonically as β\beta increases if β\beta is less than 2.5. For larger β\beta micromotion induced local maxima arise that prevents inferring the magnitude of micromotion from fluorescence rate. Since our operating regime is below this threshold, we can use the change in fluorescence rate to detect the change in magnitude of the stray electric field.

III Results

A gradient descent algorithm (ADAM) and a deep learning network (MLOOP) were tested for compensating stray fields in different working regimes. The source code used for the experiments is available in [28]. The software controlled the voltages using the PXI-6713 DAQ and read the fluorescence counts from a photo-mutliplier-tube (PMT) through a time tagging counter (IDQ id800). All software was written in python and interfaced with the DAQ hardware using the library NI-DAQmx Python. A total of 44 DC electrodes and the horizontal position of the cooling laser were tuned by the program, resulting in a total of 45 input parameters.

III.1 Gradient Descent optimizer

The first compensation test was performed by ADAM gradient descent algorithm. This is a first order optimizer that uses the biased first and second order moments of the gradient to update the inputs of an objective function, and was chosen for its fast convergence, versatility in multiple dimensions and tolerance to noise[23]. Our goal was to maximize the fluorescence of the ion which was described by a function f(α)f(\vec{\alpha}{\,}), where α=(α1,α2,α3,α45)\vec{\alpha}{\,}=(\alpha_{1},\alpha_{2},\alpha_{3},…\alpha_{45}) represents the array of parameters to be optimized. To find the optimal α\vec{\alpha}{\,}, the algorithm needs to know the values of the partial derivatives for all input parameters. Because we do not have an analytic expression for f(α)f(\vec{\alpha}{\,}), the values of its derivatives were estimated from experimental measurements by sequentially changing each input αi\alpha_{i}, and reading the associated change in fluorescence ff. This data were used as inputs to ADAM for finding the optimal α\vec{\alpha}{\,} which maximized ff.

Refer to caption
Figure 3: The photon count improvement using ADAM starting from manually optimized waveform weights. This graph demonstrates a 78±1%78\pm 1\% improvement in ion fluorescence. The error bars are a combination of the readout error of the PMT and peak to peak variations in the photon count whilst not optimizing. The background photon count here is 1184±\pm34 counts/s.

Before running the automated compensation, we manually adjusted the 4 weights of the waveforms used for compensation described in the previous section. We also tried to run ADAM to optimize these 4 parameters but the increase in fluorescence was limited to 6%. After manual compensation, we ran ADAM on all 45 inputs with the algorithm parameters given in the source code [28]. Each iteration took 12 s, where 9.8 s were the photon readout (0.1 s×2 readouts per parameter plus 2×0.1 s readouts at the beginning and end of the iteration), and the rest of the time was the gradient computation. If the photon count reduced by more than 40% of its initial value, the algorithm terminated and applied the previously found optimum. This acted as the safety net for the program, ensuring the ion was not lost while optimizing the 45 inputs. We need this safety net because if the ion is heated past the capture range for the used cooling detuning, it will be ejected from the trap.

Refer to caption
Figure 4: Voltage deviation from the original starting point during optimization with ADAM. a: uncharged trap (section III.1). b: During UV charging (section III.4). Top graphs show odd electrode numbers corresponding to top DC electrodes in Fig. 1 and the bottom graphs show the even electrode numbers. The values were determined by subtracting the voltage at each iteration by the starting voltage ΔV=VnV0\Delta V=V_{n}-V_{0}. Changes can be seen in almost all the electrodes of the trap.

In our implementation of the algorithm we removed the reduction in the step size of the optimization algorithm as iterations progressed. This step reduction, which is present in the standard version of ADAM, is not ideal when stray fields change with time since the optimal values of the voltages also drift in time. The removal caused some fluctuations in the photon readout near the optimal settings. Adding to these fluctuations, other sources of noise, such as wavemeter laser locking [29], and mechanical drift in the trap environment, resulted in daily photon count variations of around 5%. Despite these fluctuations, and the fact that stray fields change every day, the algorithm demonstrated an increase in fluorescence collection up to 78±1% (Fig. 3) when starting from a manually optimized configuration.

The ADAM algorithm was fast and reliable (the ion was never lost during optimization), even in extremely volatile conditions like having time-dependent charging and stray electric field buildup. Figure 4a shows a colourmap of the voltages and laser position adjustments, where most of the improvement came from adding the same voltage to all DC electrodes indicating the ion was not at optimal height. The volatility of the ion-trap environment causes the fluorescence rate to oscillate around the optimal point. To get the best value, instead of using the values of the final iteration, the software saved all voltage combinations and applied the setting with the highest photon count after all iterations were finished. Despite picking the best value it can be seen in Fig. 3 that the fluorescence for some iterations during the optimization are higher than the final point selected by the software. This is because when the settings are changed, the ion fluorescence rate may transiently increase and subsequently stabilize to a slightly lower value for the same voltage settings.

III.2 Deep Learning Network for the Ion trap

Refer to caption
Figure 5: MLOOP deep learning network. Differential Evolution explores the input space (blue points) and the neural network creates a model of the data and predicts an optimum (red points). Maximum photon count of the neural network points is 96±\pm1% higher than manual optimization. Differential Evolution continues to explore the input space and has varied photon counts. The beginning point for the process (found by manually adjusting the 4 waveform weights) was at 33700 counts/s and the highest photon count found by the neural network was at 66200 counts/s.

The second algorithm tested was a deep learning network using the python based optimization and experimental control package MLOOP[20]. MLOOP uses Differential Evolution [30] for exploring and sampling data. The blue points in Fig. 5 corresponds to these samples and it can be seen that even at the end of optimization, they can have non-optimum fluorescence rates.

MLOOP also trains a neural network using the data collected by Differential Evolution and creates an approximate model of the experimental system. It then uses this model to predict an optimum point. The red points in Fig. 5 shows the optimum points predicted by the neural network model. It can be seen that this section starts later than Differential Evolution, as it requires some data for initial neural network training, and gradually finds the optimum and stays near it. For training of the neural network, the inbuilt ADAM optimizer is used to minimize the cost function.

The sampling in MLOOP does not require a gradient calculation which greatly improves the sampling time. Even though the sampling is fast, training the network to find an optimal point requires a minimum of 100 samples and that makes MLOOP slower than ADAM. With our settings for MLOOP, each iteration took 0.7 s on average and therefore 700 s was needed to take 1000 samples shown in Fig. 5.

In our test the neural network in MLOOP had 5 layers with 45 nodes each, all with Gaussian error correction. The neural network structure was manually optimized and tested on a 45-dimensional positive definite quadratic function before being used for the experiment.

Once the ion was trapped, positioned above the integrated mirror [22], and photon counts were read, the program started sampling 100 different voltage combinations around its initial point. Then, the network started training on the initial data and making predictions for the voltages that maximise fluorescence.

Since the ion trap setup is very sensitive to changes in the electric field, the voltages were set to move a maximum of 1% of their previous value in each iteration to reduce the chance of losing the ion. As a step size value could not be explicitly defined, this percentage was chosen to make the changes similar to the step size used for ADAM.

A small percentage of our initial trials with the maximum change of a few percent (instead of 1%) led to an unstable ion during the parameter search sequence. This is because MLOOP is a global optimizer and can set the voltages to values far from the stable starting point. Since the ion trap is a complicated system that can only be modelled for a specific range of configurations, moving away from these settings can lead to unpredictable and usually unstable behavior.

MLOOP also has an in-built mechanism that handles function noise using a predefined expected uncertainty. We set this uncertainty to the peak-to-peak noise of the photon readout when no optimization was running.

Since MLOOP is a global optimizer it was able to find optimum points different from the points found by ADAM. For trials where low numbers of initial training data points were used, these configurations proved to be unstable and in most cases resulted in the loss of the ion. Unstable states were also observed occasionally if the optimizer was run for too long. With moderate-size training sets, MLOOP was able to find voltage settings with fluorescence rates similar or higher than optimum points found by ADAM as shown in Fig. 5.

Considering the long duration of the MLOOP iteration sequence and the possibility of finding unstable settings in volatile conditions, the test of optimization with induced changing stray fields (section III.4) was only performed with the ADAM optimizer as the gradient based search method proved to be more robust against fluctuations in the ion environment.

III.3 Ion Properties before and after optimization

Refer to caption
Figure 6: Fluorescence vs laser frequency detuning from the resonance for inital setting and after different optimizations. It can be seen from (a) that the experimental values are very close to the theoretical Lorentzian fit [31, 32, 33]. This shows the heating is low before and after optimization and therefore the change in fluorescence can be used to infer the change in heating. Deviation from the theory near the resonance shown in (b) is a sign of small heating instability.
PsatP_{sat} (μW\mu W) η\eta (%) Bkg¯(1/s)\bar{Bkg}(1/s)
Initial 4.11(1) 0.971(2) 1184(34)
ADAM-45 1.86(1) 0.938(1) 1227(35)
MLOOP-45 1.86(1) 0.989(1) 1291(35)
Table 1: Saturation power (PsatP_{sat}), detection efficiency (η\eta) and average background count rate (BkgBkg) before and after optimization with different methods. PsatP_{sat} and η\eta values are calculated by fitting the theoretical formula for Fluorescence Rate vs Laser Power to the experimental data [31, 32, 33]. A clear drop in saturation power is observed after ADAM and MLOOP optimization of individual electrodes indicating a reduction of the micromotion.

To test the effectiveness of the protocols, the saturation power, PsatP_{sat}, was measured before and after the optimization process. The PsatP_{sat} is the laser power at which the fluorescence rate of a two-level system is half the fluorescence at infinite laser power. We also measured the overall detection efficiency η\eta, the fraction of emitted photons which resulted in detection events. Table. 1 shows PsatP_{sat} decreased (ion photon absorption was improved) using both ADAM and MLOOP. The detection efficiency was approximately the same for all runs as expected.

Another test was done by measuring fluorescence vs laser detuning before and after optimization. Figure 6 shows that the measured values follows the expected Lorentzian profile [31, 32, 33] and associated linewidth before and after optimization. This indicates that the initial micromotion magnitude β\beta was sufficiently small for fluorescence to be a good optimization proxy. Clear increase in florescence can be seen after optimizing 44 electrodes individually both with ADAM and MLOOP. The fit residual curve (difference between the experimental values and the theoretical fit) shows that optimizing individual electrodes, resulted in slight increase in heating instability near the resonance.

III.4 Testing under poor trap conditions

To test the live performance of the optimization protocol in a non-ideal situation, we deliberately charged the trap by shining 369.5 nm UV laser light onto the chip for 70 minutes. The power of the laser was 200±15μW200\pm 15\mu W and the Gaussian diameter of the focus was 120±10μm120\pm 10\mu m. This process ejects electrons due to the photo-electric effect [34] and produces irregular and potentially unpredictable slow time varying electric fields within the trap. The process charged the trap significantly and made a noticeable reduction to the photon count. The ADAM algorithm was then tested both during charging and after charging was stopped. In both cases an improvement of fluorescence rate was observed.

Refer to caption
Figure 7: Real time compensation with ADAM of laser charging induced stray electric field. The ion was optimized using ADAM (left blue points) then the photon count was noted whilst charging for 70 minutes (red points) then re-optimized (right blue points). Initial improvement from manually optimized settings was 27%. The second optimization improved the fluorescence by 58% from the charged conditions and returned it back to the optimum value of the first optimization within the error.
Refer to caption
Figure 8: The trap was charged by hitting a UV laser to destabilize the ion and individual electrodes optimized using ADAM simultaneously for 70 minutes. The photon count fluctuates as a result of combination of fluctuations of power of the cooling laser power, algorithm search and charging irregularities. The optimizer keeps the fluorescence at the photon count similar to the case of optimizing after the charging is stopped (third section of Fig.7).

The first experiment was performed to test the optimizing process after charging. In this test, starting with the optimal manual setting, ADAM individual electrode optimizer was able to obtain 27% improvement in the fluorescence rate (blue points on the left side of Fig. 7). Then charging was induced onto the trap for 70 min and a clear decrease in photon count was seen that went even lower than the initial value (red points in Fig. 7). At this point charging was stopped and ADAM was run again and fluorescence rate returned back to the previous optimum, within the error, in approximately 12 min. During the second optimization, the fluorescence goes higher than the stable final value for some iterations before the final. This is because of the same effect explained in section III.1 that the fluorescence might spike right after a change but go down slightly after stabilizing. Looking at the changes of individual electrodes, shown in Fig. 4b, we see that the main electrodes adjusted were those around the ion and some throughout the trap. The change in the laser horizontal position was negligible.

Another experiment was done by running ADAM during continuous charging for real-time compensation. Since we induce charging via laser scattering from the trap, the collected photons are both from the ion and the scattered laser and fluctuations in the intensity of scattered light confuses the optimizer. Despite that the optimizer did not lose the ion nor needed to abort the process. Figure 8 shows that the fluorescence rate, even after a 70-minute charging session, remained near the optimum value. After stopping the charging, the ion remained trapped for more than 8 hours and was intentionally removed from the trap after this time.

IV Conclusion

Comparing the results of the gradient descent method and MLOOP shows that ADAM, the local optimizer, was a better performer in volatile conditions. This is because unlike MLOOP, ADAM does not need any training and can react faster to changes. MLOOP as a global optimizer was able to find settings with higher fluorescence but occasionally optimums were unstable if the defined parameter search range for the settings were large or if the machine learner ran for too long. The main drawback of MLOOP was that it was much slower than ADAM. In both cases photon collection times were the limiting factor.

Optimizing the individual electrodes improved the fluorescence by 78% for ADAM and 96% for MLOOP. ADAM showed its reliability in the charging test, returning the ion to the optimal fluorescence in 12 minutes after 70 minutes of purposely charging the trap with a UV laser. Improvement amounts depended greatly on the overall stability of the cooling laser and efficiency of the photon counting system.

To improve the speed of the gradient based optimizer, shorter readout times can be used although this increases variability in the value of the photon count. It is also possible to try local optimizers like SPSA [35] that approximate the gradient using only 2 readouts instead of 90. This improves the time needed for each iteration but because using this approximate gradient reduces the improvement in each iteration, the overall speed comparison to reach the same optimum can only be made when both are tried in practice. One possible problem with SPSA is that it cannot handle noise like ADAM and multiple readouts and averaging are needed in noisy conditions. If high amount of averaging is required, the algorithm may become very slow.

IV.1 Acknowledgement

This research was financially supported by the Griffith University Research Infrastructure Programme, the Griffith Sciences equipment scheme, and Australian Research Council Linkage (LP180100096) Projects, ML is supported by Australian Research Council Future Fellowship (FT180100055), AZ was supported by Centre for Quantum Dynamics, Griffith University, Summer Scholarship, JS, and KS were supported by the Australian Government Research Training Program Scholarship.

References

  • Debnath et al. [2016] S. Debnath, N. M. Linke, C. Figgatt, K. A. Landsman, K. Wright, and C. Monroe, Demonstration of a small programmable quantum computer with atomic qubits, Nature 536, 63 (2016).
  • Lekitsch et al. [2017] B. Lekitsch, S. Weidt, A. G. Fowler, K. Mølmer, S. J. Devitt, C. Wunderlich, and W. K. Hensinger, Blueprint for a microwave trapped ion quantum computer, Science Advances 3, e1601540 (2017).
  • Monroe and Kim [2013] C. Monroe and J. Kim, Scaling the Ion Trap Quantum Processor, Science 339, 1164 (2013).
  • Blinov et al. [2004] B. B. Blinov, D. Leibfried, C. Monroe, and D. J. Wineland, Quantum Computing with Trapped Ion Hyperfine Qubits, Quantum Information Processing 3, 45 (2004).
  • Wineland et al. [1998] D. J. Wineland, C. Monroe, W. M. Itano, D. Leibfried, B. E. King, and D. M. Meekhof, Experimental Issues in Coherent Quantum-State Manipulation of Trapped Atomic Ions, Journal of Research of the National Institute of Standards and Technology 103, 259 (1998).
  • Kielpinski et al. [2002] D. Kielpinski, C. Monroe, and D. J. Wineland, Architecture for a large-scale ion-trap quantum computer, Nature 417, 709 (2002).
  • Chiaverini et al. [2005] J. Chiaverini, R. B. Blakestad, J. Britton, J. D. Jost, C. Langer, D. Leibfried, R. Ozeri, and D. J. Wineland, Surface-electrode architecture for ion-trap quantum information processing, Quantum Information & Computation 5, 419 (2005).
  • Seidelin et al. [2006] S. Seidelin, J. Chiaverini, R. Reichle, J. J. Bollinger, D. Leibfried, J. Britton, J. H. Wesenberg, R. B. Blakestad, R. J. Epstein, D. B. Hume, W. M. Itano, J. D. Jost, C. Langer, R. Ozeri, N. Shiga, and D. J. Wineland, Microfabricated Surface-Electrode Ion Trap for Scalable Quantum Information Processing, Physical Review Letters 96, 253003 (2006).
  • Doret et al. [2012] S. C. Doret, J. M. Amini, K. Wright, C. Volin, T. Killian, A. Ozakin, D. Denison, H. Hayden, C.-S. Pai, R. E. Slusher, and A. W. Harter, Controlling trapping potentials and stray electric fields in a microfabricated ion trap through design and compensation, New Journal of Physics 14, 073012 (2012).
  • Amini et al. [2010] J. M. Amini, H. Uys, J. H. Wesenberg, S. Seidelin, J. Britton, J. J. Bollinger, D. Leibfried, C. Ospelkaus, A. P. VanDevender, and D. J. Wineland, Toward scalable ion traps for quantum information processing, New Journal of Physics 12, 033031 (2010).
  • Kielpinski [2001] D. Kielpinski, Entanglement and decoherence in a trapped-ion quantum registerPh.D. thesis (2001).
  • Tanaka et al. [2012] U. Tanaka, K. Masuda, Y. Akimoto, K. Koda, Y. Ibaraki, and S. Urabe, Micromotion compensation in a surface electrode trap by parametric excitation of trapped ions, Applied Physics B 107, 907 (2012).
  • Allcock et al. [2012] D. T. C. Allcock, T. P. Harty, H. A. Janacek, N. M. Linke, C. J. Ballance, A. M. Steane, D. M. Lucas, R. L. Jarecki, S. D. Habermehl, M. G. Blain, D. Stick, and D. L. Moehring, Heating rate and electrode charging measurements in a scalable, microfabricated, surface-electrode ion trap, Applied Physics B 107, 913 (2012).
  • Narayanan et al. [2011] S. Narayanan, N. Daniilidis, S. A. Möller, R. Clark, F. Ziesel, K. Singer, F. Schmidt-Kaler, and H. Häffner, Electric field compensation and sensing with a single ion in a planar trap, Journal of Applied Physics 110, 114909 (2011).
  • Berkeland et al. [1998] D. J. Berkeland, J. D. Miller, J. C. Bergquist, W. M. Itano, and D. J. Wineland, Minimization of ion micromotion in a Paul trap, Journal of Applied Physics 83, 5025 (1998).
  • Noek et al. [2013] R. Noek, T. Kim, E. Mount, S.-Y. Baek, P. Maunz, and J. Kim, Trapping and cooling of 174Yb+ ions in a microfabricated surface trap, Journal of the Korean Physical Society 63, 907 (2013).
  • Allcock et al. [2010] D. T. C. Allcock, J. A. Sherman, D. N. Stacey, A. H. Burrell, M. J. Curtis, G. Imreh, N. M. Linke, D. J. Szwer, S. C. Webster, A. M. Steane, and D. M. Lucas, Implementation of a symmetric surface-electrode ion trap with field compensation using a modulated Raman effect, New Journal of Physics 12, 053026 (2010).
  • Harlander et al. [2010] M. Harlander, M. Brownnutt, W. Hänsel, and R. Blatt, Trapped-ion probing of light-induced charging effects on dielectrics, New Journal of Physics 12, 093035 (2010), arXiv: 1004.4842.
  • Seif et al. [2018] A. Seif, K. A. Landsman, N. M. Linke, C. Figgatt, C. Monroe, and M. Hafezi, Machine learning assisted readout of trapped-ion qubits, Journal of Physics B: Atomic, Molecular and Optical Physics 51, 174006 (2018).
  • Wigley et al. [2016] P. B. Wigley, P. J. Everitt, A. van den Hengel, J. W. Bastian, M. A. Sooriyabandara, G. D. McDonald, K. S. Hardman, C. D. Quinlivan, P. Manju, C. C. N. Kuhn, I. R. Petersen, A. N. Luiten, J. J. Hope, N. P. Robins, and M. R. Hush, Fast machine-learning online optimization of ultra-cold-atom experiments, Scientific Reports 6, 25890 (2016).
  • Tranter et al. [2018] A. D. Tranter, H. J. Slatyer, M. R. Hush, A. C. Leung, J. L. Everett, K. V. Paul, P. Vernaz-Gris, P. K. Lam, B. C. Buchler, and G. T. Campbell, Multiparameter optimisation of a magneto-optical trap using deep learning, Nature Communications 9, 4360 (2018).
  • Ghadimi et al. [2017] M. Ghadimi, V. Blūms, B. G. Norton, P. M. Fisher, S. C. Connell, J. M. Amini, C. Volin, H. Hayden, C.-S. Pai, D. Kielpinski, M. Lobino, and E. W. Streed, Scalable ion–photon quantum interface based on integrated diffractive mirrors, npj Quantum Information 3, 1 (2017).
  • Kingma and Ba [2017] D. P. Kingma and J. Ba, Adam: A Method for Stochastic Optimization, arXiv:1412.6980 [cs]  (2017), arXiv: 1412.6980.
  • Shappert et al. [2013] C. M. Shappert, J. T. Merrill, K. R. Brown, J. M. Amini, C. Volin, S. C. Doret, H. Hayden, C.-S. Pai, K. R. Brown, and A. W. Harter, Spatially uniform single-qubit gate operations with near-field microwaves and composite pulse compensation, New Journal of Physics 15, 083053 (2013).
  • Blūms et al. [2020] V. Blūms, J. Scarabel, K. Shimizu, M. Ghadimi, S. C. Connell, S. Händel, B. G. Norton, E. M. Bridge, D. Kielpinski, M. Lobino, and E. W. Streed, Laser stabilization to neutral Yb in a discharge with polarization-enhanced frequency modulation spectroscopy, Review of Scientific Instruments 91, 123002 (2020).
  • Norton et al. [2011] B. G. Norton, E. W. Streed, M. J. Petrasiunas, A. Jechow, and D. Kielpinski, Millikelvin spatial thermometry of trapped ions, New Journal of Physics 13, 113022 (2011).
  • Budker et al. [2004] D. Budker, D. Kimball, D. F. Kimball, and D. P. DeMille, Atomic Physics: An Exploration Through Problems and Solutions (Oxford University Press, 2004).
  • Zappacosta [2021] A. Zappacosta, Auto ion compensation, github.com/alexzappa11/Auto_Ion_Compensation (2021).
  • Ghadimi et al. [2020] M. Ghadimi, E. M. Bridge, J. Scarabel, S. Connell, K. Shimizu, E. Streed, E. Streed, M. Lobino, and M. Lobino, Multichannel optomechanical switch and locking system for wavemeters, Applied Optics 59, 5136 (2020).
  • Price et al. [2005] K. Price, R. M. Storn, and J. A. Lampinen, Differential Evolution: A Practical Approach to Global Optimization, Natural Computing Series (Springer-Verlag, Berlin Heidelberg, 2005).
  • Metcalf et al. [1999] H. J. Metcalf, P. v. d. Straten, and P. v. d. Straten, Laser Cooling and Trapping, Graduate Texts in Contemporary Physics (Springer-Verlag, New York, 1999).
  • Ejtemaee et al. [2010] S. Ejtemaee, R. Thomas, and P. C. Haljan, Optimization of yb171+ fluorescence and hyperfine-qubit detection, Physical Review A 82, 063419 (2010).
  • Berkeland and Boshier [2002] D. J. Berkeland and M. G. Boshier, Destabilization of dark states and optical spectroscopy in Zeeman-degenerate atomic systems, Physical Review A 65, 033413 (2002).
  • Wang et al. [2011] S. X. Wang, G. H. Low, N. S. Lachenmyer, Y. Ge, P. F. Herskind, and I. L. Chuang, Laser-induced charging of microfabricated ion traps, Journal of Applied Physics 110, 104901 (2011), arXiv: 1108.0092.
  • Spall [1998] J. C. Spall, Implementation of the simultaneous perturbation algorithm for stochastic optimization, IEEE Transactions on Aerospace and Electronic Systems 34, 817 (1998).