In this paper we consider two problems in diagnostics of trapped ion crystals in which an analysis of the ions' collective oscillatory motion yield potentially useful results. When one of the ions in a linear crystal undergoes a collision, observation of the subsequent motion allows one to deduce the identity of which ion sustained the collision. When a linear ion crystal is formed with a dark impurity ion, analysis of the ions' motion can identify the mass (and thus give an important clue to the species) of the impurity.
In this paper we consider two problems in diagnostics of trapped ion crystals in which an analysis of the ions' collective oscillatory motion yield potentially useful results. When one of the ions in a linear crystal undergoes a collision, observation of the subsequent motion allows one to deduce the identity of which ion sustained the collision. When a linear ion crystal is formed with a dark impurity ion, analysis of the ions' motion can identify the mass (and thus give an important clue to the species) of the impurity.
In this paper we consider two problems in diagnostics of trapped ion crystals in which an analysis of the ions' collective oscillatory motion yield potentially useful results. When one of the ions in a linear crystal undergoes a collision, observation of the subsequent motion allows one to deduce the identity of which ion sustained the collision. When a linear ion crystal is formed with a dark impurity ion, analysis of the ions' motion can identify the mass (and thus give an important clue to the species) of the impurity.
We have developed a novel scheme to measure the secular motion of trapped ions. Employing pulsed excitation and analysis of the fluorescence of laser-cooled ions, we have measured the centre-of-mass mode frequency of single as well as entire ion crystals with a frequency precision better than 5 × 10 −4 within an interrogation time on the order of seconds, limited only by the fluorescence collection efficiency and the background noise. We have used this method to measure the mass of ions and observed charge exchange collisions between trapped calcium isotopes.
Coulomb crystals have applications in many areas such as quantum computing and simulation, quantum logic spectroscopy, nonlinear dynamics and chaos, phase transitions, and chemical reaction process. The structure of the Coulomb crystal and the trajectory of each trapped ion are typically determined by the parameters of the trap and the ion species. However, dark ions are often inevitable in experiment, which introduces uncertainty into the desired crystal structures and ion trajectories. Few researches have been conducted to investigate the configuration change of the crystal in the presence of dark ions and the influence of a dark ion on its surrounding ion trajectories in a multi-ion system. In this work, we utilize the molecular dynamics simulation software LAMMPS and the (py)LIon package (modified to adapt the semi-classical theory of laser cooling) for simulating the three-dimensional ion trajectories of Coulomb crystals. The formation process of 40Ca+ Coulomb crystal in a linear trap is simulated. With the micromotion and secular motion trajectories of each ion, we calculate the temperature of Coulomb crystal and the average velocity of specific ions. It is observed that the crystal structure exhibits obvious layering phenomenon when the trapped ions yield a large difference in their charge-to-mass ratio (CMR), however, layering is not obvious with a small difference in the CMR. In addition, we simulate and compare the Coulomb crystal structure formed by pure 40Ca+ ions with that formed by 40Ca+ ions mixed with a small number of dark ions including isotopic ions (44Ca+) and impurity ions (CaH+). Three different cases are investigated, namely the one-dimensional ion string, two-dimensional planar structure and three-dimensional helical structure. The results show that the ions in the neighborhood of a dark ion exhibit around micron-order position change compared with their positions before the dark ion is formed. Such a change can be measured in experiment through microscopic imaging, thereby providing a way to identify the formation of dark ions in Column crystals with a large ion number.
Coulomb crystal yields a wide range of applications including quantum computing and simulation, phase transitions, quantum logic spectroscopy, nonlinear dynamics and chaos, chemical reaction process, etc. The structure of the Coulomb crystal and the trajectory of each trapped ion are typically determined by the parameters of the trap and the ion species. However, dark ions are often inevitable in experiments, which introduces uncertainty to the desired crystal structures and ion trajectories. Few research has been conducted to investigate the configuration change of the crystal in the presence of dark ions and the influence of a dark ion on its surrounding ion trajectories in a multi-ion system. In this paper, we utilize the molecular dynamics simulation software LAMMPS and the (py)Lion package (modified to adapt the semi-classical theory of laser cooling) for the simulation of the 3D ion trajectories of Coulomb crystals. The formation process of 40Ca+ Coulomb crystals in a linear trap is simulated. With the micromotion and secular motion trajectories of each ion, we calculate the temperature of Coulomb crystal and the average velocity of specific ions. It’s observed that the crystal structure exhibits obvious layering phenomenon when the trapped ions yield a large difference in their charge-to-mass ratio (CMR), however, layering is not obvious with a small difference in the CMR. In addition, we simulate and compare the Coulomb crystal structure formed by pure 40Ca+ ions and that formed by 40Ca+ ions mixed with a small number of dark ions including isotopic ions (44Ca+) and impurity ions (CaH+). Three different cases are investigated, namely the 1D ion string, 2D plane structure and 3D helical structure. Results show that ions in the neighborhood of a dark ion exhibit around micron-order position changes compared to their positions before the dark ion is formed. Such a change can be measured in experiment through microscopic imaging, providing a way to identify dark ions in Column crystals with a large number of ions.
We examine the time-dependent dynamics of ion crystals in radiofrequency traps. The problem of stable trapping of general three-dimensional crystals is considered and the validity of the pseudopotential approximation is discussed. We analytically derive the micromotion amplitude of the ions, rigorously proving well-known experimental observations. We use a recently proposed method to find the modes that diagonalize the linearized time-dependent dynamical problem. This allows one to obtain explicitly the ('Floquet–Lyapunov') transformation to coordinates of decoupled linear oscillators. We demonstrate the utility of the method by analyzing the modes of a small 'peculiar' crystal in a linear Paul trap. The calculations can be readily generalized to multispecies ion crystals in general multipole traps, and time-dependent quantum wavefunctions of ion oscillations in such traps can be obtained.
We firstly discuss classical stability for a dynamical system of two ions levitated in a 3D Radio-Frequency (RF) trap, assimilated with two coupled oscillators. We obtain the solutions of the coupled system of equations that characterizes the associated dynamics. In addition, we supply the modes of oscillation and demonstrate the weak coupling condition is inappropriate in practice, while for collective modes of motion (and strong coupling) only a peak of the mass can be detected. Phase portraits and power spectra are employed to illustrate how the trajectory executes quasiperiodic motion on the surface of torus, namely a Kolmogorov-Arnold-Moser (KAM) torus. In an attempt to better describe dynamical stability of the system, we introduce a model that characterizes dynamical stability and the critical points based on the Hessian matrix approach. The model is then applied to investigate quantum dynamics for many-body systems consisting of identical ions, levitated in 2D and 3D ion traps. Finally, the same model is applied to the case of a combined 3D Quadrupole Ion Trap (QIT) with axial symmetry, for which we obtain the associated Hamilton function. The ion distribution can be described by means of numerical modeling, based on the Hamilton function we assign to the system. The approach we introduce is effective to infer the parameters of distinct types of traps by applying a unitary and coherent method, and especially for identifying equilibrium configurations, of large interest for ion crystals or quantum logic.