## Impact of ion motion on atom-ion confinement-induced resonances in hybrid traps

Vladimir S.Melezhik, 1, 2, ∗ Zbigniew Idziaszek, 3, † and Antonio Negretti 4, ‡

1 Bogoliubov Laboratory of Theoretical Physics, Joint Institute for Nuclear Research, Dubna, Moscow Region 141980, Russian Federation 2 Dubna State University, 19 Universitetskaya Street, Moscow Region 141982, Russian Federation 3 Faculty of Physics, University of Warsaw, Pasteura 5, 02-093 Warsaw, Poland 4 Zentrum f¨r Optische Quantentechnologien, u Fachbereich Physik, and Hamburg Center for Ultrafast Imaging, Universit¨t a Hamburg, Luruper Chaussee 149, 22761 Hamburg, Germany (Dated: October 28, 2019)

We investigate confinement-induced resonances in atom-ion quantum mixtures confined in hybrid traps for small atom-ion mass ratios. Specifically, we consider an ion confined in a time-dependent radio-frequency Paul trap with linear geometry, while the atom is constrained to move into a quasione-dimensional optical waveguide within the ion trap. We evaluate the impact of the ion intrinsic micromotion on the resonance position. Thus, we solve the atom-ion dynamics semiclassically, namely the atom dynamics is governed by the three-dimensional time-dependent Schr¨dinger equao tion, whereas the ion motion is described by the classical Hamilton equations. We find that the energy of the ion provided by the oscillating radiofrequency fields can affect the resonance position substantially. Notwithstanding, the peculiar phenomenology of those resonances regarding perfect transmission and reflection is still observable. These findings indicate that the intrinsic micromotion of the ion is not detrimental for the occurrence of the resonance and that its position can be controlled by the radiofrequency fields. This provides an additional mean for tuning atom-ion interactions in low spatial dimensions. The study represents an important advancement in the scattering physics of compound atomic quantum systems in time-dependent traps.

PACS numbers: 32.60.+i,33.55.Be,32.10.Dk,33.80.Ps

## I. INTRODUCTION

Compound atom-ion systems afford a new platform to study quantum physics in which multi-energy and multilength scales are involved. In particular, atom-ion systems allow investigating condensed-matter systems more closely. For instance, an important component of a solidstate system is the electron-phonon coupling, which is mimicked naturally in an atom-ion system. Indeed, the atom-ion interaction has the effect that the passage of an atom in the proximity of an ion crystal influences the state of the ions via the exchange of phonons as in a real solid-state system. This feature is absent in ultracold atoms trapped in an optical lattice, where there is no back-action of the atoms on the lattice. Furthermore, with such compound system it is also possible to study the formation of mesoscopic molecular ions [1, 2] and charge transport [3, 4], to mention a few examples (for a detailed overview, see Ref. [5]). The realisation of compound atom-ion systems in the laboratory, however, is quite challenging, since it combines different trap technologies such as radiofrequency traps for ions and optical dipole traps for atoms and this limits fundamentally the attainable temperatures. Albeit the s-wave collisional

∗ melezhik@theor.jinr.ru

† Zbigniew.Idziaszek@fuw.edu.pl

‡ anegrett@physnet.uni-hamburg.de

quantum regime (i.e. the s-wave regime) has been not yet reached in current atom-ion laboratories, a significant experimental effort has been put forward in very recent years. Specifically, it has been shown that submicrokelvin temperatures can be attained when ionising a Rydberg atom in a Bose-Einstein condensate [6] and that the required collisional energies to enter the s-wave limit in radiofrequency traps are attained when choosing a small atom-ion mass ratio [7-11], which is within experimental reach [12]. Furthermore, proof-of-principle experiments have demonstrated laser-controlled atom-ion interactions [13-15], therefore opening new possibilities for controlling interactions and developing light-matter interfaces for quantum information processing [16, 17]. We underscore that the attainment of the s-wave regime is crucial for the observation of atom-ion Feshbach resonances [18] as well as for quantum technological applications with atom-ion systems such as quantum gates [19], quantum simulation of the electron-phonon coupling [20], and for reaching the strong-coupling polaron regime [21].

Confinement-induced resonances (CIRs) [22-27] have been pivotal in entering the regime of strongly correlated atomic matter [28]. A CIR occurs in an atomic trap when the atom-atom scattering length in free space, a s , becomes comparable to the transversal width, a ⊥ , of the trap: a ⊥ /a s → 1.4603 [22]. At that ratio the bound state in the closed channel, that is, the first excited state of the transverse confinement, coincides with the collision threshold of the entrance channel. In this case, the scattering amplitude f ( a ⊥ , a s ) of the atom on atom

tends to -1, and therefore the transmission T = 1 + | f | 2 approaches zero [22-25]. Thus, by varying a ⊥ and a s near the CIR, it is possible to control the effective atomatom interaction of the confined atomic system. As for the neutral atomic counterpart, atom-ion or atom-dipole CIRs can be an additional 'knob' for manipulating the mixture's interaction in low spatial dimensions. For instance, CIRs can be utilised to steer the particle flow in Josephson junctions [9, 29, 30], to tune interactions in a Tomonaga-Luttinger liquid in which a linear ion crystal is immersed [31], for precise magnetometry [32, 33], and quantum simulation [20, 34-38]. This is especially important in the atom-ion setting when state-dependent atom-ion interactions are needed to perform particular quantum information processing tasks. Indeed, experiments [39-41] have shown that spin-exchange collisions can occur after a few Langevin collisions because of spinorbit couplings [42]. This effect has a minor impact in a quasi-one dimensional geometry.

A first theoretical study investigating the possibility of realising CIRs in an atom-ion system was done in Ref. [43] for the case of time-independent atom and ion traps with identical frequencies. In Ref. [25] the CIRs in atom-ion hybrid systems were predicted and the conditions for the atom-ion CIR appearance were obtained in the static ion approximation, where the ion is pinned in a precise position in space and cannot move. Interestingly, it has been found an isotope-like effect of the resonance position. While in the neutral setting the resonance condition is attained when a ⊥ /a s → 1.4603, as previously discussed, in atom-ion systems, when the effective spatial range, R ∗ , of the polarisation potential, -C /r 4 4 , is comparable to a ⊥ (or even larger), the position of the CIR strongly relies on the atom-ion mass ratio [25]. This effect is much harder to observe in neutral atomic systems, as it requires large trap frequencies (hundreds of kHz) in order to produce small trap widths a ⊥ , while atom-ion systems require frequencies ≤ 100 kHz [25].

In view of the foregoing, it becomes relevant to go beyond the static approximation in the problem of atomion CIR. So far, however, atom-ion collisions in radiofrequency traps have been treated either purely classically [7] or within the Markovian quantum master equation formalism [44]. Both studies indicated that the most favourable atom-ion species for reaching the quantum regime is Li/Yb + . In the study of Ref. [44], however, the atom-ion interaction has been treated perturbatively, i.e. in the Born approximation, which omits, for instance, the resonance effects in the atom-ion scattering. In this work we go one step further in the quantum mechanical treatment of atom-ion collisions in Paul traps and extend the previous work of two of us [25]. Indeed, we treat the atom dynamics fully quantum mechanically with the time-dependent Schr¨dinger o equation, whereas the ion motion is described classically via the Hamilton equations. Such a semiclassical treatment is well justified when the ion is much heavier than the colliding atom. As we explain in detail below, the equations of motion of the atom and the ion are coupled via the atom-ion interaction. In the atom Schr¨dinger equation the ion o position is treated as a time-dependent parameter, while the interaction in the ion Hamilton equations is averaged quantum mechanically over the instantaneous quantum state of the atom.

Towards this aim, the quantum-semiclassical computational method [45-48] specifically designed for particle collisions such as the problem of ionisation of the helium ion colliding with protons [47] and antiprotons [48] has been employed and extended to the time-dependent domain, as our radiofrequency ionic confinement requires. Moreover, our analyses focus on the specific Li/Yb + atom-ion pair, since it is the most promising atomic pair to reach the s-wave regime in Paul traps and it is currently under intense experimental investigations [12, 40, 49]. We note, however, that with regard to the mass ratio a good atom-ion pair is also Li/Ca + [5052]. Our analysis shows that the intrinsic micromotion of the ion, which is unavoidable in Paul traps (differently from optical trapping of ions [53]), is not detrimental for the occurrence of an atom-ion CIR. We find that the CIR position strongly relies on the ion kinetic energy. This implies that the atom-ion interaction can be controlled not only by the short-range atom-ion physics or by means of magnetic Fano-Feshbach resonances as well as the width of the atomic waveguide, but also by the ion kinetic energy, which can be manipulated via the external driving, e.g., by changing the radiofrequency of the Paul trap.

The paper is organised as follow. In Sec. II we introduce our microscopic atom-ion Hamiltonian and the corresponding equations of motion, namely the coupled Schr¨dinger-Hamilton equations. Beside this, we discuss o how to determine the scattering amplitude in the timedependent scenario. In Sec. III we present our findings and discuss the physical implications. In Sec. IV we draw our conclusions and provide an outlook for future work.

## II. PROBLEM FORMULATION AND METHODOLOGY

In this section we describe theoretically the atom-ion system including the equations of motion for the atom and the ion, how do we face the scattering problem in time-dependent radiofrequency traps, and what are the scattering quantities of interest in order to assess the CIR position in reliance of the ion kinetic energy and atomic waveguide width. A schematic view of the system under investigation is given in Fig. 1.

## A. The atom quantum Hamiltonian

The atom is described by the Hamiltonian

FIG. 1: (color online) Pictorial illustration of the system under investigation in this study. The light- and dark-grey electrodes (the big bars in the figure) of the Paul trap generate the time-dependent electric fields needed to confine the ion (red sphere) transversally, whereas longitudinally a static voltage is applied to ensure confinement (not shown). The atom (green sphere) is injected from the right to the left into a waveguide (light-blue), whose centre hosts the ion. The waveguide is orientated along the longitudinal axis, z , of the linear Paul trap. In the transverse directions, x, y , the confining potential both for the atom and the ion is strong.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01151v2-figure-f316c1fa-b04a-4ded-b35a-1f678f56385e.jpg)

$$\hat { H } _ { a } ( \hat { r } _ { a }, t ; \mathbf r _ { i } ) = - \frac { \hslash ^ { 2 } } { 2 m _ { a } } \nabla _ { a } ^ { 2 } & + \frac { m _ { a } \omega _ { \perp } ^ { 2 } } { 2 } ( \hat { x } _ { a } ^ { 2 } + \hat { y } _ { a } ^ { 2 } ) \\ & + V _ { a i } ( | \hat { \mathbf r } _ { a } - \mathbf r _ { i } ( t ) | ), \quad ( 1 ) \quad \text{ Here} \\ & \quad \underset { \dots } { \omega _ { i } }$$

where ˆ r a ≡ (ˆ x , y a ˆ a , z ˆ ) and a r i ≡ ( x , y , z i i i ) are the atom position operator and ion position vector, respectively, m a is the atom mass, and ∇ 2 a ≡ ∂ x 2 a + ∂ y 2 a + ∂ z 2 a is the laplacian operator ( ∂ x a denotes the partial derivative, for instance, in the x a direction). The second term in Eq. (1) describes the harmonic potential of frequency ω ⊥ due to the waveguide, whereas along the longitudinal direction, z a , the atom does not experience any confinement. The transverse potential is assumed to be so tight to render the atom motion quasi-one-dimensional (quasi-1D). We note that in the Hamiltonian (1) we explicitly introduced the parametric dependence on the ion position.

The last term in Eq. (1), V ai ( ˆ | r a -r i ( ) ), t | describes the interaction between the atom and ion, whose asymptotic behaviour at large atom-ion separations r t ( ) ≡ | ˆ r a -r i ( ) t | →∞ has the form

$$V _ { a i } ( r ( t ) ) \simeq - \frac { C _ { 4 } } { r ( t ) ^ { 4 } } \,, \quad \quad ( 2 )$$

i.e. it behaves as the polarisation potential. The dispersion coefficient is given by C 4 = αe 2 2 1 4 πglyph[epsilon1] 0 (in SI units)

with α being the static atom polarisability, e the electron charge, and glyph[epsilon1] 0 the vacuum permittivity. At shortrange distances, when the two electronic clouds do overlap, the potential is repulsive, albeit its form is generally unknown. For the sake of numerical convenience, we use the regularised potential [44]

$$V _ { a i } ( r ( t ) ) = - \frac { [ r ^ { 2 } ( t ) - c ^ { 2 } ] } { [ r ^ { 2 } ( t ) + c ^ { 2 } ] } \frac { C _ { 4 } } { [ r ^ { 2 } ( t ) + b ^ { 2 } ] ^ { 2 } } \,. \quad ( 3 )$$

This potential reproduces at large distances the polarisation potential (2), whereas at r = 0 it assumes a large, but finite, numerical value, contrarily to the typically employed C 12 /r 12 potential. Variation of the two parameters b and c permits to tune the atom-ion interaction for any value -∞ &lt; a s &lt; + ∞ of the atom-ion scattering length in free space in the zero-energy limit. Wenote that the time-dependence of the Hamiltonian (1) enters via the ion trajectory r i ( ) t in the atom-ion interaction (3). The ion trajectory is obtained by simultaneous solving the classical equations of motion for a charge particle in a Paul trap, as explained in the next section.

## B. The ion classical Hamiltonian

The ion is assumed to be confined in a linear Paul trap, whose electric fields read as [54]

$$\dot { \mathbf E } _ { s } & = \frac { m _ { i } } { 2 | e | } \omega _ { i } ^ { 2 } \left ( x _ { i }, y _ { i }, - 2 z _ { i } \right ), \\ \mathbf E _ { r f } & = \frac { m _ { i } \Omega _ { r f } ^ { 2 } q } { 2 | e | } \cos ( \Omega _ { r f } t ) \left ( x _ { i }, - y _ { i }, 0 \right ) \,. \quad \ \ ( 4 )$$

Here, m i is the ion mass, Ω rf is the radiofrequency (rf), ω i = Ω rf √ a/ 2, q and a are dimensionless geometric parameters (i.e. q z = 0, q y = -q x ≡ q , -a / z 2 = a x = a y ≡ a , and a glyph[lessmuch] q 2 &lt; 1). Hereafter, we assume that the axis of the waveguide in which is travelling the colliding atom is precisely the z -axis of the Paul trap (see also Fig. 1). The corresponding non-conservative potential is given by

$$U ( { \mathbf r } _ { i }, t ) & = \frac { m _ { i } \omega _ { i } ^ { 2 } } { 2 } \left ( z _ { i } ^ { 2 } - \frac { x _ { i } ^ { 2 } + y _ { i } ^ { 2 } } { 2 } \right ) \\ & + \frac { m _ { i } \Omega _ { r f } ^ { 2 } } { 2 } q \cos ( \Omega _ { r f } t ) \left ( \frac { y _ { i } ^ { 2 } } { 2 } - \frac { x _ { i } ^ { 2 } } { 2 } \right ) \,. \quad ( 5 )$$

Hence, the classical Hamiltonian describing an ion in a Paul trap is given by

$$H _ { i } ^ { t r a p } ( { \mathbf p } _ { i }, { \mathbf r } _ { i }, t ) = \frac { { \mathbf p } _ { i } ^ { 2 } } { 2 m _ { i } } + U ( { \mathbf r } _ { i }, t ) \,. \quad \quad ( 6 )$$

## C. The atom-ion equations of motion

When the atom is confined in the waveguide within the Paul trap, the ion will experience its presence via

the atom-ion interaction modelled by Eq. (3). The full classical ion Hamiltonian is therefore given by

$$H _ { i } ( \mathbf p _ { i }, \mathbf r _ { i }, t ; \mathbf r _ { a } ) = H _ { i } ^ { t r a p } ( \mathbf p _ { i }, \mathbf r _ { i }, t ) \\ + \langle V _ { a i } ( | \hat { \mathbf r } _ { a } - \mathbf r _ { i } ( t ) | ) \rangle \,, \quad \quad ( 7 ) \quad \text{Pau} \\ \Lambda _ { \dots }$$

where

$$\langle V _ { a i } ( | \hat { \mathbf r } _ { a } - \mathbf r _ { i } ( t ) | ) \rangle = \int _ { \begin{array} { c } \text{out} \\ \text{$1$} \end{array} } \left \langle V _ { a i } ( | \hat { \mathbf r } _ { a } - \mathbf r _ { i } ( t ) | ) | \Psi ( \mathbf r _ { a }, t ; \mathbf r _ { i } ) \right \rangle \quad ( 8 ) \quad \ \ \ \ \ t = ( \begin{array} { c } \text{out} \\ \text{$1$} \end{array} \right \rangle \end{array} \right \rangle$$

is the quantum mechanical average of the atom-ion interaction. As for the atom case, also for the ion Hamiltonian we explicitly emphasise the parametric dependence on the atom position.

The atom wavefunction Ψ( r A , t ) is governed by the time-dependent Schr¨dinger equation o

$$i \hbar { \frac { \partial } { \partial t } } \Psi ( { \mathbf r } _ { a }, t ; { \mathbf r } _ { i } ) = \hat { H } _ { a } ( \hat { \mathbf r } _ { a }, t ; { \mathbf r } _ { i } ) \Psi ( { \mathbf r } _ { a }, t ; { \mathbf r } _ { i } ) \,, \quad ( 9 )$$

where ˆ H a (ˆ r a , t ; r i ) is defined by Eq. (1). Hence, the dynamics of the ion is governed by the Hamilton equations

$$\frac { d } { d t } \mathbf p _ { i } & = - \frac { \partial } { \partial \mathbf r _ { i } } H _ { i } ( \mathbf p _ { i }, \mathbf r _ { i }, t ; \mathbf r _ { a } ) \,, \quad \text{ with } \quad \text{ and } \quad \text{agr} \\ \frac { d } { d t } \mathbf r _ { i } & = \frac { \partial } { \partial \mathbf p _ { i } } H _ { i } ( \mathbf p _ { i }, \mathbf r _ { i }, t ; \mathbf r _ { a } ) \,. \quad \text{ (10)}$$

This set of classical equations together with the atom Schr¨dinger equation (9) forms the complete set of dyo namical equations for describing the confined atom-ion collision in hybrid traps. In the present study we consider collisions of a light atom with a much heavier ion in the range of very low atomic colliding energies E coll (ultracold atoms), where the relation p a = √ 2 m E a coll glyph[lessmuch] p i for their momentums is satisfied. In addition, we require that E i = p / 2 i (2 m i ) glyph[greatermuch] glyph[planckover2pi1] ω i , which further justifies the application of the classical description for the ion. Additionally, because of the separation of energy (i.e. time) scales among the ion and atom dynamics, in the Hamilton equations (10) we neglect the functional derivatives owed to the parametric reliance of the atom wavefunction on the ion position, that is,

$$\frac { \partial } { \partial \xi _ { i } } \Psi ( { \mathbf r } _ { a }, t ; { \mathbf r } _ { i } ) \equiv 0 \quad \xi _ { i } = x _ { i }, \, y _ { i }, \, z _ { i }. \quad ( 1 1 )$$

In order to integrate simultaneously the equations (9) and (10), we need proper initial conditions with physical significance. At the beginning of the collisional process, the atom and the ion are assumed to be far away from each other such that they do not interact ( V ai = 0). In particular, the atom is initially in the ground state of the atomic trap with the low longitudinal colliding energy, that is, E coll glyph[lessmuch] 2 glyph[planckover2pi1] ω ⊥ , whereas the ion performs fast (with respect of atom motion) oscillations in the Paul trap with mean transversal ¯ E ⊥ and longitudinal ¯ E ‖ energies. Since the atom approaches the region of interaction with the ion very slowly ( E coll / glyph[planckover2pi1] glyph[lessmuch] ω ⊥ glyph[lessmuch]

ω , i Ω rf ), the initial position of the ion does not influence the scattering process itself, which depends only on ¯ E ⊥ and ¯ E ‖ . Specifically, the classical solution of the ion equations of motion (Mathieu equation) in the Paul trap (without the atom) are well approximated by: A j cos( ω t i + φ j )[1+ q j cos(Ω rf t ) / 2] ∀ j = x, y, z [55]. The associated kinetic energy depends on the amplitude A j , but not on the phase φ j [68]. Therefore, we choose, without loss of generality, the ion position at the initial time t = 0 in the trap centre with transversal energy, E ⊥ , and longitudinal energy, E ‖ . This can be summarised with the following set of initial conditions:

$$\mathbf r _ { i } ( t & = 0 ) & = ( 0, 0, 0 ), \\ p _ { i, x } ( t & = 0 ) & = \sqrt { 2 m _ { i } E _ { \perp } }, \\ p _ { i, y } ( t & = 0 ) & = 0, \\ p _ { i, z } ( t & = 0 ) & = \sqrt { 2 m _ { i } E _ { \| } }.$$

These initial conditions set the mean values of the ion transversal and longitudinal energies as ¯ E ⊥ = 1 64 . E ⊥ (calculated numerically for an ytterbium ion in a trap with Ω rf = 2 π × 2MHz, ω i = 2 π × 63kHz, a = 0 002 . and q = 0 08) and . ¯ E ‖ = E / ‖ 2, which is in qualitative agreement with the estimate

$$\bar { E } _ { \perp } = \frac { E _ { \perp } } { 2 } \left [ 1 + \left ( \frac { q \Omega _ { r f } } { 2 \omega _ { i } } \right ) ^ { 2 } \right ] \simeq 1. 3 E _ { \perp } \quad ( 1 3 )$$

from the first-order solution of the Mathieu equation [55, 56]. We also note that such a choice for the initial transversal momentum orientation, that is, along the x -axis, is not relevant for the scattering problem we are interested in. Indeed, as we have verified in our numerical simulations, the final result is invariant relatively to the initial orientation of the ion transversal momentum p i, ⊥ ( t = 0) in the x -y plane. This fact is also a consequence of the cylindrical symmetry of both the atomic waveguide and the linear Paul trap. Hence, we can safely assume that along the y -axis the initial energy is zero.

As far as the initial condition for the atom wavefunction is concerned, we use the following ansatz [24]

$$\psi ( { \mathbf r } _ { a }, t = 0 ) & = N \varphi _ { 0 } ( \rho _ { a } ) e ^ { - \frac { ( z _ { a } - z _ { 0 } ) ^ { 2 } } { 2 a _ { z } ^ { 2 } } } e ^ { i k z _ { a } } \\ & = N \varphi _ { 0 } ( \rho _ { a } ) \chi ( z _ { a } - z _ { 0 } ) e ^ { i k z _ { a } } \,. \quad ( 1 4 )$$

Here N is a normalisation constant, ρ a = r a sin θ a with r a ≥ 0 and θ a ∈ [0 , π ), glyph[planckover2pi1] k = √ 2 m E a coll is the initial momentum, z 0 is the position of the initial wave-packet χ , which is far from the ion location, that is, at z a = z 0 the atom and ion do not interact: V ai ( z a = z , r 0 i ) → 0. The longitudinal width of the initial wave-packet (14) is chosen sufficiently broad according to a z glyph[similarequal] 30 -40 R ∗ ( R ∗ = √ 2 µC / 4 glyph[planckover2pi1] is the characteristic length scale of the atom-ion interaction) to satisfy the demand of sufficient monochromaticity of the wave-packet along the z direction D k ( t = 0) = 〈 ψ t ( = 0) ( | k -¯ k ) 2 | ψ t ( = 0) 〉 → 0.

This choice of a z provides insignificant deformation of the envelop χ during scattering due to the small dispersion D k ( ) t of the wave-packet [24]. The above initial condition is interpreted as follows: Initially, the atom and ion are far from each other such that the atom is initially prepared in the ground state of the transverse confinement, that is, in the ground state ϕ 0 of the twodimensional harmonic oscillator as well as in the ground state of an optical dipole trap that is approximated by a harmonic potential of frequency glyph[planckover2pi1] / m a ( a 2 z ). At times t &gt; 0 the atomic longitudinal confinement is suddenly switched off and an initial momentum kick is imparted to the atom (e.g., via a Raman configuration of lasers) such that the atom wave-packet is moving towards the ion trap centre with velocity v 0 = glyph[planckover2pi1] k/m a .

For details on the numerical implementation of the integrators of the atom-ion equations of motion, we refer the interested reader to the appendix A.

## D. Determination of the scattering amplitude

In the course of the collisional process the atom wavepacket splits up into two parts each of them moving in opposite directions z a → ±∞ . Asymptotically we encounter the following behaviour [22, 24, 57]

$$& \psi ( \rho _ { a }, z _ { a }, t \to + \infty ) \underset { z _ { a } \to + \infty } { \longrightarrow } \psi ^ { + } ( \rho _ { a }, z _ { a }, t ) = ( 1 + f ^ { + } ( k ) ) N \varphi _ { 0 } ( \rho _ { a } ) \tilde { \chi } ( z _ { a } - ( z _ { 0 } + v t ) ) e ^ { i k _ { f } z _ { a } }, \\ & \psi ( \rho _ { a }, z _ { a }, t \to + \infty ) \underset { z _ { a } \to - \infty } { \longrightarrow } \psi ^ { - } ( \rho _ { a }, z _ { a }, t ) = f ^ { - } ( k ) N \varphi _ { 0 } ( \rho _ { a } ) \tilde { \chi } ( - z _ { a } - ( z _ { 0 } + v t ) ) e ^ { - i k _ { f } z _ { a } }.$$

Here f ± ( k ) are the atom-ion forward and backward scattering amplitudes in the presence of the external confining potential due to the atomic waveguide and the ion radiofrequency trap. The function ˜( χ z A , t ) describes the atom motion in the longitudinal z -direction, namely the spreading of the initial Gaussian wave-packet.

Our goal is to determine the forward scattering amplitude f + ( k ). Towards that end, we solve first the atom Schr¨dinger equation without the atom-ion interaction, o but in the confining waveguide. The corresponding solution at z a → + ∞ is

$$\psi ^ { ( 0 ) + } ( \rho _ { a }, z _ { a } ) = N \varphi _ { 0 } ( \rho _ { a } ) \, \tilde { \chi } ( z _ { a } - ( z _ { 0 } + v t ) ) \, e ^ { i k z _ { a } } \,, \ \ ( 1 6 ) \quad \text{$\alpha$} \quad \text{$\infty$}$$

whereas at distances z a →-∞ it is identically zero, that is, ψ (0) -( ρ , z a a ) ≡ 0, as a consequence of the fact that there is no scattering centre, i.e. no ion. Because of the sufficient monochromaticity of the wave-packet (see the set of the envelop χ in the initial condition (14)) and the unitarity of the Schr¨dinger equation, it holds o

$$\langle \psi ^ { ( 0 ) + } ( t ) | \psi ( t ) \rangle \underset { t \to + \infty } { \longrightarrow } 1 + f ^ { + } ( k ) \,. \quad \ \ ( 1 7 ) \quad \ \ \ \ g _ { 1 D } ^ { \ }.$$

This relation is used to calculate the amplitude f + ( k ). In order to arrive to the result on the right-hand-side of Eq. (17), we performed the approximation e i ( k -k f ) z a glyph[similarequal] 1 within the region | z a -( z 0 + vt ) | ≤ a z , where the wave packets (15) and (16) do overlap, as a consequence of the smallness of the transmitted momentum during the collision [i.e. 2 π/ k ( -k f ) is much larger than the width of ˜]. χ This approximation has been verified numerically by evaluating the mean atomic energy after the collision ¯ E a ( t → + ∞ ) = lim t → ∞ + 〈 ψ t ( ) | ˆ H ψ t a | ( ) 〉 , which did not exceed the threshold of transverse atomic excitations 3 glyph[planckover2pi1] ω ⊥ .

By exploiting the current conservation law, the transmission, T , and reflection, R , coefficients are defined as

$$\mathfrak { z } \text{ the } \quad T ( k ) = | 1 + f ^ { + } ( k ) | ^ { 2 }, \quad R ( k ) = 1 - | 1 + f ^ { + } ( k ) | ^ { 2 } \,. \ ( 1 8 )$$

We note that the above expression for the reflection coefficient, R k ( ), holds for elastic scattering processes only. For inelastic scattering the current conservation law is in general violated and one has to replace the reflection coefficient by the expression R = | f -( k ) | 2 . This occurs, for instance, when the atom energy is above the energy of the transverse excited state, that is, larger than 3 glyph[planckover2pi1] ω ⊥ . In our study, however, these processes are not treated.

Finally, in the s -wave zero-energy limit f + = f -, the quasi-1D coupling constant is given by [22, 58]

$$g _ { 1 D } = \lim _ { k \to 0 } \frac { \hslash ^ { 2 } k } { m _ { a } } \frac { \Re \{ f ^ { + } ( k ) \} } { \Im \{ f ^ { + } ( k ) \} }. \quad \quad ( 1 9 )$$

The coupling constant g 1 D is the most relevant parameter for analysing confined scattering close to a CIR, where g 1 D → ±∞ [22-24, 27]. However, simultaneously with g 1 D , the position of the CIR is also controlled by checking the condition T → 0.

The above procedure for determining the scattering parameters has been already utilised for ultracold atomic collisions in a waveguide-like atomic trap [24]. The high numerical accuracy of the calculation obtained with this procedure for the transmission and reflection coefficients in wide range of the interactions including the resonant region of the CIR appearance was confirmed by a comparison with simple estimates obtained from direct calculations of the atomic probability to be out of the scattering center after the collision.

We here investigate in detail the scattering dynamics of a ytterbium ion confined in a linear Paul trap and a lithium atom in a waveguide. Atom-ion experiments with such an atomic pair are currently investigated intensively [12, 40, 49]. The atomic trap frequency, ω ⊥ , has been chosen within the range 2 π × (1 -11)kHz, whereas for the Paul trap parameters we have chosen: Ω rf = 2 π × 2 MHz, ω i = 2 π × 63 kHz, a = 0 002 and . q = 0 08 [see also . Eq. (5)]. All numerical simulations have been performed with initial longitudinal atom energy E /k ‖ B = m v / a 2 0 (2 k B ) from the region ∼ 10 nK of low energies.

We have performed simulations for different values of the atom-ion interaction V ai ( r t ( )) by varying it from strong repulsion ( a s glyph[greatermuch] 0) to strong attraction ( a s glyph[lessmuch] 0). Special attention has been focused on the region near the atom-ion CIR. In the zero-energy limit for the atom and under the static approximation for the ion, the ratio a ⊥ /a s [ a ⊥ = √ glyph[planckover2pi1] / m ω ( a ⊥ )] approaches the well known value 1.4603 [22, 25] when R ∗ glyph[lessmuch] a ⊥ . Here R ∗ = √ 2 µC / 4 glyph[planckover2pi1] is the characteristic length scale of the atom-ion interaction (2), where µ denotes the atom-ion reduced mass. The determination of the forward scattering amplitude is performed as follows. For a chosen pair of parameters b and c of the regularised atom-ion potential (3), we calculate the corresponding scattering length a s in free space ( ω ⊥ = 0), which can be easily assessed because of the separation of the center-of-mass and angular part of the atom-ion wave-function. Thereafter, we fix the value of transverse frequency ω ⊥ , i.e. a ⊥ , as well as use the previously determined b and c parameters of the interaction potential (3). Thus, we simulate the time evolution of the atomic wave-packet and the ion trajectory by integrating simultaneously the equations (9) and (10) with the initial conditions (12) and (14). Numerical integration has been performed in the time domain from t = 0 to t = 10 t ⊥ with t ⊥ = 2 π/ω ⊥ defined by the lowest frequency of the problem ω ⊥ glyph[lessmuch] ω , i Ω rf .

An example of such an analysis is illustrated in Fig. 2, where the ion coordinates as a function of time are shown together with the atomic probability density distribution | r a Ψ( r a , t ) | 2 at three different times: Before the collision ( t = 0), at the region of the atom-ion interaction ( t = 5 t ⊥ ), and after the collision ( t = 9 t ⊥ , i.e. t → + ∞ ). The initial conditions and parameters of the atomic trap, that is, z 0 = 70 R ∗ , a z = 30 R ∗ and ω ⊥ /ω ∗ = 0 02, . as well as the atom-ion interaction with a ⊥ /a s = -3 88, . were chosen in such a way that the atom at t = 0 does not interact with the ion [ V ai ( z a = z , r 0 i ) glyph[similarequal] -C / 4 (70 R ∗ ) 4 → 0]. Here ω ∗ = 2 E / ∗ glyph[planckover2pi1] and E ∗ = glyph[planckover2pi1] 2 / [2 µ R ( ∗ ) 2 ]. In this case the ion is performing stable oscillations in the Paul trap with initial conditions (12) and E ⊥ = E ‖ = 0 25 . E ∗ = 4 25 . k B µ K. Here, for the atom-ion pair 6 Li/ 174 Yb + we have: R ∗ glyph[similarequal] 69 77 . nm, that is, z 0 glyph[similarequal] 4 88 . µ m, while ω ∗ = 2 π × 357 16 kHz, . thus ω ⊥ = 2 π × 7 14 kHz and . t ⊥ = 1 ms. As it is shown in Fig. 2, at time t ∼ t ⊥ the ion begins to experience the vicinity of the atom. This is clearly displayed in the y i -component of the ion motion, where from the initial zero value suddenly large oscillations appear. This is due to the fact the ion has been displaced from the trap centre, and therefore it experiences the radiofrequency fields. Similarly, there is an enhancement of the amplitude of the oscillations of the x i -and z i -components of the ion trajectory. The collision occurs approximately within the time window t glyph[similarequal] t ⊥ up to t glyph[similarequal] 5 t ⊥ . Afterwards, when the atom and ion leave the range of the atom-ion interaction, they approach the asymptotic region, that is, the two particle do not interact. The ion coordinates reach a steady-state solution, namely the amplitude of the oscillations is approximately constant, but essentially they exceeded their initial values. This indicates that the ion has 'heated up'. Such a finding has been also identified in the master equation approach [44], and in Monte Carlo simulations for a trapped ion interacting with classical buffer gas [59]. On the other hand, the atomic wave-packet splits up in two parts moving forward and backward. Note also the peak in the atomic density distribution remains near to the origin, that is, at the centre of the Paul trap. This indicates that some part of the initial atomic wave-packet is lost or, in other words, an ionic molecule is formed in the Paul trap due to atomion collision, which is also a consequence of the negative ratio a ⊥ /a s = -3 88, . i.e. attracting atom-ion interaction. Our estimate gives the following value P mol glyph[similarequal] 0 14 . for the probability of creation of molecular ions in such confined atom-ion collisions. Such a phenomenology resembles the situation of resonant molecule formation in atomic confining traps suggested in Ref. [60], where the excess energy is transferred in an excitation of the centerof-mass of the formed molecule. Note that the two-body bound state has been not considered in the study [44], where the Born and Markov approximations have been performed.

Thereafter, we have extracted the scattering parameters f + ( k ), T , and g 1 D as outlined in Sec. II.D. Figures 3 and 4 illustrate the result for two ion energies as well as for three pairs ( b, c ) of the atom-ion potential (3): a ⊥ /a s = 1 544, which corresponds to the case . of the atom-ion CIR obtained in the static ion approximation [25]; a ⊥ /a s = 2 64, i.e. . a rather strong repulsion between atom and ion; a ⊥ /a s = -3 88, that is, consider-. able attraction between the atom and the ion. In Fig. 3 the calculated scattering parameters f + ( k ), T , and g 1 D are presented for the case of the ion being initially ( t = 0) at rest in the centre of the Paul trap with zero initial energy ( E ⊥ = E ‖ = 0). These results demonstrate the efficiency of the computational procedure outlined in Sec. II.D for extracting the scattering parameters for different strengths of the atom-ion interaction (3) including strong repulsion and attraction between atom and ion. We underline that our analysis confirms that if the initial ion energy is zero ( E ⊥ = E ‖ = 0), we find that the position of the CIR coincides rather well with the result

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01151v2-figure-0973dc31-9804-4ce3-99ed-9aad896ba620.jpg)

⊥

FIG. 2: (color online) The calculated evolution in time of the ion trajectory (left column), being initially at the state with E ⊥ = E ‖ = 0 25 . E ∗ = 4 25 . µ K, and the atom probability distribution | r a Ψ( r a , t ) | 2 (right column) for the atom-ion attractive interaction giving ratio a ⊥ /a s = -3 88 and for . ω ⊥ /ω ∗ = 0 02. .

obtained in the static ion approximation [25]. In the plot of the left top panel (i.e. a ⊥ /a s = 1 544) of Fig. 3 it is . shown the behaviour of the scattering parameters, which is characteristic of a CIR, i.e. glyph[Rfractur] [ f + ] → -1, glyph[Ifractur] [ f + ] → 0 and g 1 D → ±∞ [22, 25]. It also shows that away from the CIR position, the effective coupling constant approaches finite values, positive one for the repulsive atom-ion interaction ( a ⊥ /a s = 2 64) and negative one . for the attractive atom-ion potential ( a ⊥ /a s = -3 88). . Since the result at a ⊥ /a s = 1 544 corresponds to . the resonant scattering near the CIR, we observe here the dramatic enhancement of the effective coupling constant g 1 D with respect to the nonresonant g 1 D at a ⊥ /a s = 2 64 . and -3 88. . Strong oscillations in the asymptotic region t → + ∞ observed in g 1 D at a ⊥ /a s = 1 544 with the . time period ∼ 2 π/ (2 ω ⊥ ) corresponds to virtual transitions between the entrance channel and the first closed excited state, which separates by the energy threshold by the amount 2 glyph[planckover2pi1] ω ⊥ . This is perfectly consistent with the physical interpretation of the CIR as a resonance in the first closed transverse channel [58]. Then, we have extended our investigation to the case of the ion oscillating before the collision in the Paul trap with a rather large energy ( E ⊥ = E ‖ = 0 25 . E ∗ = 4 25 . µ K). As it is shown in Fig. 4, if we increase the ion energy the position of the CIR is shifted from the point a ⊥ /a s = 1 544, . which corresponds to the CIR in the case of the ion at rest before the collision. For these parameters, the atomion potential gives a repulsive coupling constant of finite value g 1 D / E R ( ∗ ∗ ) = 2 7 (see low left panel in Fig. 4). . Effective coupling constants g 1 D calculated for the nonresonant cases a ⊥ /a s = 2 64 and = . -3 88 are also shifted . with increasing ion energy with respect to the case of the ion at rest. We note that the oscillation frequency of the coupling constant g 1 D ( ) t in all nonresonant cases presented in Figs. 3 and 4 is considerably smaller than in the resonant case (CIR) and essentially determined by ω i . On the other hand, in the case a ⊥ /a s = -3 88, the . oscillations of g 1 D ( ) t are determined by Ω rf of the micromotion.

We continue by analysing the dependence of the position of the confinement-induced resonance on the initial (mean) ion energy in the Paul trap. To this aim, we first investigate the scenario for the secular time-independent trap

$$U _ { s e c } ( { \mathbf r } _ { i } ) = \frac { m _ { i } } { 2 } \left [ \omega _ { x y } ^ { 2 } ( x _ { i } ^ { 2 } + y _ { i } ^ { 2 } ) + \omega _ { z } ^ { 2 } z _ { i } ^ { 2 } \right ] \quad ( 2 0 ) \quad \text{emor} \\ \text{live}$$

with the frequencies ω z = 2 π × 45 kHz and ω x,y = 2 π × 150 kHz. The positions of the atom-ion CIR were obtained by looking for the positions of the singular points in the coupling constant g 1 D ( E ,E ⊥ ‖ ) (19), namely when the coupling constant diverges. In Fig. 5 we present the calculated dependence of the CIR position a ⊥ /a s on the transversal E ⊥ and longitudinal E ‖ ion energies. The values of a ⊥ /a s given in Fig. 5 at E ⊥ ∼ E ‖ → 0 confirm the result obtained in Ref. [25] in the static ion approximation in the limit R ∗ glyph[lessmuch] a ⊥ for which the atomion CIR position coincides with the well-known result a ⊥ /a s = 1 4603 [22]. . Moreover, it is found that the calculated position of the atom-ion CIR is fixed quite well near the constant value a ⊥ /a s glyph[similarequal] 1 5 in the square domain . E /k ⊥ B , E ‖ /k B ≤ 10 µ K. In other words, in the secular harmonic trap approximation (20), the CIR position is stabilised near the value (1.4603) obtained in the static approximation for the ion (independent of the ion mean energy) if the ion transversal and longitudinal initial energies do not exceed the value of 10 µ K (see the almost flat region in Fig. 5), which is close to the the s-wave threshold energy E ∗ = glyph[planckover2pi1] 2 / [2 µ R ( ∗ ) 2 ] glyph[similarequal] 6 4 . µ K [7].

We proceed further with our analysis of the CIR position by investigating the impact of the time-dependent Paul trap with the confining potential defined in Eq. (5). We have performed calculations for two confining atomic trap frequencies ω ⊥ = 0 02 . ω ∗ and 0 03 . ω ∗ corresponding to ω ⊥ = 2 π × 7 1 kHz and . ω ⊥ = 2 π × 11 kHz, respectively, for the pair 6 Li/ 174 Yb + . The results of the calculation of the CIR positions presented in Figs. 6 and 7 support the approximation of the secular trap in the rectangular domain E ⊥ ≤ 1 µ K, E ‖ ≤ 5 µ K. Outside this domain, the positions a ⊥ /a s calculated for the time-dependent Paul trap (5) start deviate from the values calculated in the secular harmonic approach (20). Herewith, the deviation increases with increasing ion energy. Moreover, out of the mentioned above energy domain we have found a different dependence of the CIR position on the transversal and longitudinal ion mean energies in the Paul trap. In the zero-energy limit for the ion initial energy ( E ,E ⊥ ‖ → 0) we obtain the value for the CIR position coinciding with the result obtained earlier in the static ion approximation [25] as well as the same dependence on the ratio R /a ∗ ⊥ as in the case of the static ion approximation: with decreasing a ⊥ (increasing ω ⊥ ) the position of the CIR shifts towards larger a ⊥ /a s and deviates from the value 1.4603 (see the left top panels of Figs. 6 and 7). However, for larger ion mean energies, outside the domain E ⊥ ≤ 1 µ K, E ‖ ≤ 5 µ K, the CIR position is not anymore captured by the static ion approximation.

The found effect of the ion motion in the radiofrequency field on the CIR position can be qualitatively interpreted as follows. Since the atom moves slowly relative the period of ion oscillations, it feels the averaged effective interaction with the ion. At a rather low ion energy (small amplitude of the ion vibrations), the ion motion only slightly corrects the aforementioned effective atom-ion interaction and slightly changes the relative collisional atom-ion energy, since we observe perfect reflection near the CIR position a ⊥ /a s = 1 46 obtained in . the approximation of a static ion. The CIR position (the point of perfect reflection) starts to deviate considerably from the value calculated in the static ion approximation when the ion energy approaches the s-wave atom-ion collision threshold E ∗ . Therefore, it is natural to assume here a considerable contribution of the p-wave to the scattering amplitude f + . In the case of a comparable contribution of s- and p-waves, total transmission can be observed due to their interference - this is the so-called

1.0

1.0

(t)/[2E * R *

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01151v2-figure-81bc8720-ce55-4354-aa3e-4cf95be99c29.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01151v2-figure-f90cea4a-c7f5-40c3-bfcf-8713ef24f12d.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01151v2-figure-06ca8437-5604-4535-8462-e047365f1856.jpg)

FIG. 3: (color online) The calculated f + ( ), t T t ( ), g 1 D ( ) for the ion being initially at rest (i.e. t with zero energy before the collision with atom) for three different values of the ratio a ⊥ /a s .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01151v2-figure-1301f1f1-6554-46f1-901f-f1d7816901f9.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01151v2-figure-da265cb2-f6f2-4409-960f-49f8e438224b.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01151v2-figure-e5f74136-e4fa-40c7-b4bc-e156d75a490b.jpg)

1.0

(t)/[2E * R *

(t)/[2E * R *

FIG. 4: (color online) The calculated f + ( ), t T t ( ), g 1 D ( ) for an initial ion energy t E ⊥ = E ‖ = 0 25 . E ∗ = 4 25 . k B × µ K ( p ⊥ = p x ) for three different values of the ratio a ⊥ /a s .

FIG. 5: (color online) The calculated position a ⊥ /a s of the atom-ion CIR as a function of the transversal E ⊥ and longitudinal E ‖ kinetic energy of the ion in the secular harmonic trap approximation (20) with frequencies ( ω z = 2 π × 45 kHz and ω x,y = 2 π × 150 kHz).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01151v2-figure-24860791-b244-4029-b83e-975e31ff59db.jpg)

dual CIR predicted by J.I. Kim et. al. [23]. However, the question of dual CIRs in hybrid atom-ion systems demands a separate analysis.

## IV. CONCLUSIONS

We have investigated the conditions for the appearing of atom-ion CIRs in hybrid atom-ion systems. Our analysis has been done for the Li-atom confined in an optical trap situated within a linear Paul trap for a Yb-ion with the realistic parameters and takes into account the motion and micromotion of the ion. Such a choice for the atom-ion pair is motivated by the fact that it is very likely the only pair that permits to attain s-wave collisions in hybrid traps. The shifts of the CIRs due to the ion motion were calculated. We found that the energy of the ion provided by the oscillating radiofrequency fields can affect the resonance position substantially. However, in a broad range of the ion kinetic energies we found that the CIR position is stable even in presence of micromotion. These findings indicate that the intrinsic micromotion of the ion is not detrimental for the occurrence of the resonance and that its position can be controlled by the radiofrequency fields. This provides an additional mean for tuning atom-ion interactions in low spatial dimensions. It also indicates that experimental investigation of atom-ion CIRs in hybrid atom-ion systems is possible. Furthermore, we found a non-negligible probability (about 14%) of forming an atom-ion molecule, thus indicating a pathway for producing such two-body compounds by controlling the confinement of the two atomic species.

The performed study of resonant collisions in the con- fined geometry of hybrid atom-ion traps represents an important advancement in the scattering physics of quantum systems in time-dependent traps. To this end, we have adopted a quantum-semiclassical approach developed in Refs. [45-48]. Our treatment can be applied for further studies of time-dependent problems, e.g. for simulations of two-qubit quantum gates [19, 61] including the effects of micromotion [8]. An extension of the method to hybrid atom-ion systems with comparable atom-ion mass ratios requires further work and it will be carried out in the future. In particular, the functional derivatives (11) will be taken into account. Such studies will be important for understanding more deeply the reliance of the CIR position on the atom-ion mass ratio as well as the impact of the ion micromotion, namely whether a critical mass ratio does exist for the occurrence of CIRs. Indeed, in our study the small atom-ion mass ratio has played a major role to the fact that the micromotion is not too detrimental for the appearance of the CIR. Even more interesting will be to understand whether the s-wave regime can be attained in low spatial dimensions more favorably that in 3D for higher atom-ion mass ratios.

## V. ACKNOWLEDGEMENTS

The work was supported by the Russian Foundation for Basic Research, Grants No. 18-02-00673, by the Cluster of Excellence projects 'The Hamburg Centre for Ultrafast Imaging' of the Deutsche Forschungsgemeinschaft (EXC 1074, Project No. 194651731) and 'CUI: Advanced Imaging of Matter' of the Deutsche Forschungsgemeinschaft (EXC 2056, Project No. 390715994), and by the Polish National Science Center project 2014/14/M/ST2/00015.

## Appendix A: Computational method

In order to integrate the semiclassical atom-ion equations of motion (9) and (10), we applied the splittingup method with a 2D discrete-variable representation (DVR) [62-64]. For an accurate inclusion of the (strong) atom-ion potential (3) in Eq. (9) at the moment of the resonant atom-ion collision, a tailored splitting-up procedure in the 2D-DVR representation has been developed, as we describe below. For the integration of the Hamilton equations of motion, which involve three considerably different scales of frequencies, namely Ω rf , ω i as well as ω ⊥ in the quantum mechanical average 〈 Ψ( r a , t ; r i ) | V ai ( ˆ | r a -r i ( ) ) Ψ( t | | r a , t ; r i ) 〉 , we employed the second order St¨rmer-Verlet method [67]. o

To begin with, the atom wavefunction Ψ( r a , t ) in the 2D-DVR is expanded as [63-66] (for the sake of simplicity, we omit hereafter the parametric dependence on the

a

/

a

a

a

FIG. 6: The dependence of the CIR position a ⊥ /a s on the transversal E ⊥ and longitudinal E ‖ ion energies before the collision for the atomic trap with ω ⊥ = 0 02 . ω ∗ = 2 π × 7 1 kHz. The (blue) arrow in the top left panel (i.e. . E ⊥ = 0) indicates the atomic CIR position 1.4603 at the zero-energy limit [22]. Here the results of the calculations with the time-dependent Paul trap (5) are given as black cycles connected by (black) solid lines. The results obtained in the secular harmonic trap approximation (20) are given as open cycles connected by (red) thin lines.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01151v2-figure-46bf0a77-45d0-4458-87a1-d9b7f4b8dabe.jpg)

a

/

a

FIG. 7: The dependence of the CIR position a ⊥ /a s on the transversal E ⊥ and longitudinal E ‖ ion energies before the collision for the atomic trap with ω ⊥ = 0 03 . ω ∗ = 2 π × 11 kHz. The (blue) arrow in the top left panel (i.e. E ⊥ = 0) indicates the atomic CIR position 1.4603 at the zero-energy limit [22]. Here the results of the calculations with the time-dependent Paul trap (5) are given as black cycles connected by solid lines. The results obtained in the secular harmonic trap approximation (20) are given as open cycles connected by (red) thin lines.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01151v2-figure-0825cc19-58cb-45ed-8634-05829d9f3dfe.jpg)

ion position)

$$\Psi ( r _ { a }, \Omega, t ) = \frac { 1 } { r _ { a } } \sum _ { j = 1 } ^ { N } f _ { j } ( \Omega ) \, \psi _ { j } ( r _ { a }, t ) \quad \ \ ( A 1 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

with the 2D basis defined as

$$f _ { j } ( \Omega ) = \sum _ { \nu = \{ l m \} = 1 } ^ { N } Y _ { \nu } ( \Omega ) ( Y ^ { - 1 } ) _ { \nu j }. \quad \text{ (A2)} \quad \text{ e}.$$

The latter is defined on an angular grid Ω j = ( θ j θ , φ j φ ) of N grid points. The number N is equal to the number of basis functions in the expansion (A1) and the number of terms in the definition (A2). The coefficients ( Y -1 ) νj in Eq. (12) are the elements of the N × N matrix Y -1 inverse to the matrix given by the values Y jν = Y ν (Ω ) j of the polynomials Y ν (Ω) at the grid points Ω . j The construction of the 2D polynomials Y ν (Ω), which slightly deviate from the classical spherical harmonics for large ν ∼ N that are orthogonal on the grid Ω , is described j in detail in Refs. [63-66]. We note that the coefficients ψ j ( r a , t ) in Eq. (A1) define the values of the searching solution Ψ( r a , t ) at the points of the angular grid Ω : j ψ j ( r a , t ) = r a Ψ( r a , Ω j , t ).

With the 2D-DVR, the 3D Schr¨ odinger equation (9) is thus approximated by the system of Schr¨dinger-like o equations (hereafter glyph[planckover2pi1] = 1)

$$i \frac { \partial } { \partial t } \bar { \psi } _ { j } ( r _ { a }, t ) = \sum _ { j ^ { \prime } } ^ { N } [ \hat { H } ^ { ( 0 ) } _ { j } ( r _ { a }, t ) \delta _ { j j ^ { \prime } } + \hat { H } ^ { ( 1 ) } _ { j j ^ { \prime } } ( r _ { a } ) ] \bar { \psi } _ { j ^ { \prime } } ( r _ { a }, t ) \, \quad \text{$\underbrace{\emph{$\tau$ is $i$} \\ \quad \text{$\text{$of the $bs}$} } \\ \text{with $the$Homiliation$concistent$of the $diamond $1$} \quad.$$

with the Hamiltonian consisting of the diagonal

$$\hat { H } _ { j } ^ { ( 0 ) } ( r _ { a } ) = - \frac { 1 } { 2 m _ { a } } \frac { d ^ { 2 } } { d r _ { a } ^ { 2 } } + \frac { 1 } { 2 } \omega _ { \perp } ^ { 2 } r _ { a } ^ { 2 } \sin ^ { 2 } \theta _ { j } + V _ { a i } ( r _ { a }, \Omega _ { j }, \mathbf r _ { i } ( t ) ) \quad \ e q u a t i { \sin \, t h e }$$

(A4)

and offdiagonal

The fact that the 2D-DVR [i.e, the functions (A2)] gives the diagonal representation for ˆ H (0) ( r a , t ) and the Y ν -representation gives the diagonal representations for ˆ H (1) , respectively, has been exploited for constructing an efficient computational algorithm. Actually, for the first and the last steps of Eq. (A6), the exponential operators can be approximated according to

$$& \text{$\mathbb{A}$} ) \quad \exp \left [ - \frac { i } { 2 } \Delta t \hat { H } _ { j } ^ { ( 0 ) } ( r _ { a }, t _ { n } ) \right ] \approx \left [ 1 + \frac { i } { 4 } \Delta t \hat { H } _ { j } ^ { ( 0 ) } ( r _ { a }, t _ { n } ) \right ] ^ { - 1 } \\ \text{$\mathbb{I}_{j}$} ) \quad \times \left [ 1 - \frac { i } { 4 } \Delta t \hat { H } _ { j } ^ { ( 0 ) } ( r _ { a }, t _ { n } ) \right ] + O ( \Delta t ^ { 2 } ) \,, & ( A 7 )$$

Then, we have to solve N independent second order differential equations

$$\left [ 1 + \frac { i } { 4 } \Delta t \hat { H } _ { j } ^ { ( 0 ) } ( r _ { a }, t _ { n } ) \right ] \bar { \psi } \left ( t _ { n } + \frac { 1 } { 2 } \Delta t \right ) = \\ \left [ 1 - \frac { i } { 4 } \Delta t \hat { H } _ { j } ^ { ( 0 ) } ( r _ { a }, t _ { n } ) \right ] \bar { \psi } ( t _ { n } ). & ( A 8 )$$

The intermediate step in Eq. (A6) depending on ˆ H (1) jj ′ ( r a ) is treated in the basis Y ν (Ω ), where the matrix operaj tor ˆ H (1) ν is diagonal. The transformation with help of the simple unitary matrix S jν = λ 1 / 2 j Y jν between 2DDVR (A2) and Y ν -representation provides the efficiency of the computational procedure [63, 66]: computational time increases almost linearly with increasing number N of basis functions.

Simultaneously to the forward in time propagation t n → t n +1 = t n + ∆ t of the atom wave-packet ψ j ( r a , t n ) → ψ j ( r a , t n +1 ) we integrate the Hamilton equations (10), which describe the dynamics of the ion in the Paul trap. To this end, we utilise the second-order St¨rmer-Verlet method [67] o

$$\hat { H } _ { j j ^ { \prime } } ^ { ( 1 ) } ( r _ { a } ) = - \frac { 1 } { 2 m _ { a } r _ { a } ^ { 2 } } \frac { 1 } { \sqrt { \lambda _ { j } \lambda _ { j ^ { \prime } } } } \sum _ { \nu = \{ l m \} = 1 } ^ { N } ( Y ^ { - 1 } 0 _ { j \nu } l ( l + 1 ) ( Y ^ { - 1 } ) _ { \nu j ^ { \prime } } \quad \Phi _ { i } ^ { ( n + } \\ \intertext { a r t e } \text{for $\aleph$ or the while of the Game and rutro} \quad \,.$$

$$\Phi _ { i } ^ { ( n + 1 / 2 ) } = \Phi _ { i } ^ { ( n ) } - \frac { \Delta t } { 2 } \frac { \partial } { \partial \mathbf r _ { i } } H _ { i } ( \mathbf p _ { i } ^ { ( n + 1 / 2 ) }, \mathbf r _ { i } ^ { ( n ) } ) \,,$$

parts. Here λ j are the weights of the Gauss quadratures related with the (angular) grid Ω j [66] and ¯ ( ψ j r a , t ) = √ λ ψ j j ( r a , t ).

To integrate the system of equations (A3) we use the split-operator method that yields the propagation ¯ ( ψ j r a , t n ) → ¯ ( ψ j r a , t n +1 ) with t n → t n +1 = t n + ∆ t according to

$$\bar { \psi } ( t _ { n } + \Delta t ) & \approx \exp \left ( - \frac { i } { 2 } \Delta t \hat { H } ^ { ( 0 ) } \right ) \exp ( - i \Delta t \hat { H } ^ { ( 1 ) } ) \\ & \times \exp \left ( - \frac { i } { 2 } \Delta t \hat { H } ^ { ( 0 ) } \right ) \bar { \psi } ( t _ { n } ) + O ( \Delta t ^ { 3 } ). \ \ ( A 6 ) \quad \text{with} \\ & \mathbf P _ { i } ^ { ( n + 1 }$$

$$\mathbf r _ { i } ^ { ( n + 1 ) } = \mathbf r _ { i } ^ { ( n ) } + \frac { \Delta t } { 2 } \left \{ \frac { \partial } { \partial \mathbf r _ { i } } H _ { i } ( \mathbf p _ { i } ^ { ( n + 1 / 2 ) }, \mathbf r _ { i } ^ { ( n ) } ) \\ + \frac { \partial } { \partial \mathbf r _ { i } } H _ { i } ( \mathbf p _ { i } ^ { ( n + 1 / 2 ) }, \mathbf r _ { i } ^ { ( n + 1 ) } ) \right \} \,,$$

$$\mathbf p _ { i } ^ { ( n + 1 ) } = \mathbf p _ { i } ^ { ( n + 1 / 2 ) } - \frac { \Delta t } { 2 } \frac { \partial } { \partial \mathbf r _ { i } } H _ { i } ( \mathbf p _ { i } ^ { ( n + 1 / 2 ) }, \mathbf r _ { i } ^ { ( n + 1 ) } ). \ ( A 9 )$$

$$( A 6 ) \quad \begin{array} { c c } \text{with} \ p _ { i } ^ { ( n ) } \ = \ p _ { i } \left ( t _ { n } \right ), \ p _ { i } ^ { ( n + 1 / 2 ) } \ = \ p _ { i } \left ( t _ { n } + \frac { \Delta t } { 2 } \right ), \text{ and} \\ \ p _ { i } ^ { ( n + 1 ) } = \ p _ { i } \left ( t _ { n } + \Delta t \right ) \text{ and the same definition for} \ r _ { i } ^ { ( n ) }. \end{array}$$

- [1] R. Cˆ ot´, e V. Kharchenko, and M. D. Lukin, Phys. Rev. Lett. 89 , 093001 (2002).
- [2] J. M. Schurer, A. Negretti, and P. Schmelcher, Phys. Rev. Lett. 119 , 063001 (2017).
- [3] R. Cˆ ot´, Phys. Rev. Lett. e 85 , 5316 (2000).
- [4] R. Mukherjee, Phys. Rev. A 100 , 013403 (2019).
- [5] M. Tomza, K. Jachymski, R. Gerritsma, A. Negretti, T. Calarco, Z. Idziaszek, and P. S. Julienne, Rev. Mod. Phys. 91 , 03500 (2019).
- [6] K.S. Kleinbach, F. Engel, T. Dieterle, R. L¨w, T. Pfau, o and F. Meinert, Phys. Rev. Lett. 120 , 193401 (2018).
- [7] M. Cetina, A.T. Grier, and V. Vuletic, Phys. Rev. Lett. 109 , 253201 (2012).
- [8] L. H. Nguyen, A. Kalev, M. D. Barrett, and B.-G. Englert, Phys. Rev. A 85 , 052718 (2012).
- [9] J. Joger, A. Negretti, and R. Gerritsma, Phys. Rev. A 89 , 063621 (2014).
- [10] H.A. F¨rst and N.V. Ewald, T. Secker, J. Joger, T. Feldu ker, and R Gerritsma, J. Phys. B: At. Mol. Phys. 51 , 195001 (2018).
- [11] M. Tomza, Ch. P. Koch, and R. Moszynski, Phys. Rev. A 91 , 042706 (2015).
- [12] T. Feldker, H. F¨ urst, H. Hirzler, N. V. Ewald, M. Mazzanti, D. Wiater, M. Tomza, and R. Gerritsma, arXiv:1907.10926
- [13] N.V. Ewald, T. Feldker, H. Hirzler, H. F¨ urst, and R. Gerritsma, Phys. Rev. Lett. 122 , 253401 (2019).
- [14] S. Haze, J. Wolf, M. Deiß, L. Wang, G. Raithel, and J.-H. Denschlag, arXiv:1901.11069
- [15] F. Engel, T. Dieterle, T. Schmid, C. Tomschitz, C. Veit, N. Zuber, R. L¨w, T. Pfau, and F. Meinert, Phys. Rev. o Lett. 121 , 193401 (2018).
- [16] T. Secker, R. Gerritsma, A.W. Glaetzle, and A. Negretti, Phys. Rev. A 94 , 013420 (2016).
- [17] T. Secker, N. Ewald, J. Joger, H. F¨rst, T. Feldker, and u R. Gerritsma, Phys. Rev. Lett. 118 , 263201 (2017).
- [18] Z. Idziaszek, T. Calarco, P. S. Julienne, and A. Simoni, Phys. Rev. A 79 , 010702(R) (2009).
- [19] H. Doerk, Z. Idziaszek, and T. Calarco, Phys. Rev. A 81 , 012708 (2010).
- [20] U. Bissbort, D. Cocks, A. Negretti, Z. Idziaszek, T. Calarco, F. Schmidt-Kaler, W. Hofstetter, and R. Gerritsma, Phys. Rev. Lett. 111 , 080501 (2013).
- [21] Casteels, W., J. Tempere, and J. Devreese, J. Low Temp. Phys. 162 , 266 (2011).
- [22] M. Olshanii, Phys. Rev. Lett. 81 , 938 (1998).
- [23] J.I Kim, V.S. Melezhik, and P. Schmelcher, Phys. Rev. Lett. 97 , 193203 (2006).
- [24] V.S. Melezhik, J.I. Kim, and P. Schmelcher, Phys. Rev. A 76 , 053611 (2007).
- [25] V.S. Melezhik and A. Negretti, Phys. Rev. A 94 , 022704 (2016).
- [26] M. Sroczy`ska, T. Wasak, K. Jachymski, T. Calarco, and n Z. Idziaszek, Phys. Rev. A 98 , 012708 (2018).
- [27] S. Shadmehri and V.S. Melezhik, Phys. Rev. A 99 , 032705 (2019).
- [28] E. Haller, M. Gustavsson, M. J. Mark, J. G. Danzl, R. Hart, G. Pupillo, H.-C. N¨gerl, Science a 325 , 1224 (2009).
- [29] J.M. Schurer, R. Gerritsma, P. Schmelcher, and A. Negretti, Phys. Rev. A 93 , 063602 (2016).
- [30] M. Ebgha, S. Saeidian, P. Schmelcher, and A. Negretti,
- arXiv:1902.09594
- [31] A. Michelsen, M. Valiente, N.T. Zinner, and A. Negretti, arXiv:1907.07090
- [32] K. Jachymski, T. Wasak, Z. Idziaszek, P.S. Julienne, A. Negretti, and T. Calarco, Phys. Rev. Lett. 120 , 013401 (2018).
- [33] T. Wasak, K. Jachymski, T. Calarco, and A. Negretti, Phys. Rev. A 97 , 052701 (2018).
- [34] A.S. Dehkharghani, E. Rico, N.T. Zinner, and A. Negretti, Phys. Rev. A 96 , 043611 (2017).
- [35] A. Negretti, R. Gerritsma, Z. Idziaszek, F. SchmidtKaler, and T. Calarco, Phys. Rev. B 90 , 155426 (2014).
- [36] D. Gonz´ alez-Cuadra, P.R. Grzybowski, A. Dauphin, and M. Lewenstein, Phys. Rev. Lett. 121 , 090402 (2018).
- [37] D. Gonzlez-Cuadra, A. Dauphin, P.R. Grzybowski, P. Wjcik, M. Lewenstein, and A. Bermudez, Phys. Rev. B 99 , 045139 (2019).
- [38] G. Pupillo, A. Griessner, A. Micheli, M. Ortner, D.W. Wang, and P. Zoller, Phys. Rev. Lett. 100 , 050402 (2008).
- [39] L. Ratschbacher, C. Sias, L. Carcagni, J. M. Silver, C. Zipkes, and M. K¨hl, Phys. Rev. Lett. o 110 , 160402 (2013).
- [40] H. F¨rst, T. Feldker, N. V. Ewald, J. Joger, M. Tomza, u and R. Gerritsma, Phys. Rev. A 98 , 012713 (2018).
- [41] T. Sikorsky, M. Morita, Z. Meir, A. A. Buchachenko, R. Ben-shlomi, N. Akerman, E. Narevicius, T. V. Tscherbul, and R. Ozeri, Phys. Rev. Lett. 121 , 173402 (2018).
- [42] T. V. Tscherbul, P. Brumer, and A. A. Buchachenko, Phys. Rev. Lett. 117 , 143201 (2016).
- [43] Z. Idziaszek, T. Calarco, and P. Zoller, Phys. Rev. A 76 , 033409 (2007).
- [44] M. Krych and Z. Idziaszek, Phys. Rev. A 91 , 023430 (2015).
- [45] V.S. Melezhik and P. Schmelcher, Phys. Rev. Lett. 84 , 1870 (2000).
- [46] V.S. Melezhik, Hypefine Int. 138 , 351 (2001).
- [47] V.S. Melezhik, J.S. Cohen, and C.Y. Hu, Phys. Rev. A 69 , 032709 (2004).
- [48] V.S. Melezhik and L.A. Sevastianov, Analytical and Computational Methods in Probability Theory, LNCS 10684 , 449 (2017).
- [49] J. Joger, H. F¨rst, N. Ewald, T. Feldker, M. Tomza, and u R. Gerritsma, Phys. Rev. A 96 , 030703(R) (2017).
- [50] S. Haze, R. Saito, M. Fujinaga, and T. Mukaiyama, Phys. Rev. A 91 , 032709 (2015)
- [51] S. Haze, M. Sasakawa, R. Saito, R. Nakai, and T. Mukaiyama, Phys. Rev. Lett. 120 , 043401 (2018).
- [52] R. Saito, S. Haze, M. Sasakawa, R. Nakai, M. Raoult, H. Da Silva, Jr., O. Dulieu, and T. Mukaiyama, Phys. Rev. A 95 , 032709 (2017).
- [53] J. Schmidt, A. Lambrecht, P. Weckesser, M. Debatin, L. Karpa, and T. Schaetz, Phys. Rev. X 8 , 021028 (2018).
- [54] D. Leibfried, R. Blatt, C. Monroe, and D. Wineland, Rev. Mod. Phys. 75 , 281 (2003).
- [55] D.J. Berkeland, J.D. Miller, J.C. Berquist, W.M. Itano, and D.J. Wineland, J. Appl. Phys. 83 , 5025 (1998).
- [56] L.D. Landau and E.M. Lifshitz, in Mechanics (Pergamon, New York, 1976). pp.93-95.
- [57] J.I. Kim, J. Schmiedmayer, and P. Schmelcher, Phys. Rev. A 72 , 042711 (2005).

- [58] T. Bergeman, M.G. Moore, and M. Olshanii, Phys. Rev. Lett. 91 , 163201 (2003).
- [59] R.G. DeVoe, Phys. Rev. Lett. 102 , 063001 (2009).
- [60] V.S. Melezhik and P. Schmelcher, New J. Phys. 11 , 073031 (2009).
- [61] R. Stock, I. H. Deutsch, and E. L. Bolda, Phys. Rev. Lett. 91 , 183201 (2003).
- [62] V.S. Melezhik, AIP Conference Proceedings 1479 (2012) p.1200.
- [63] V.S. Melezhik, A computational method for quantum dynamics of a three-dimensional atom in strong fields, in Atoms and Molecules in Strong External Fields , Eds.
- P. Schmelcher and W. Schweizer (Plenum, New-York and London, 1998) p.89.
- [64] V.S. Melezhik, Phys. Lett. A 230 , 203 (1997).
- [65] V.S. Melezhik and D. Baye, Phys. Rev. C 59 , 3232 (1999).
- [66] V.S. Melezhik, EPJ Web of Conf. 108 , 01008 (2016).
- [67] E. Hairer, Ch. Lubich, and G. Wanner. Geometric Numerical integration. Structure-Preserving Algorithms for Ordinary Differential Equations (Springer-Verlag, Berlin Heidelberg, 2006).Ch.I.
- [68] We have verified this fact by numerical simulations for our parameters as well.
Vladimir S. Melezhik Zbigniew Idziaszek Antonio Negretti
Impact of ion motion on atom-ion confinement-induced resonances in hybrid traps
We investigate confinement-induced resonances in atom-ion quantum mixtures
confined in hybrid traps. Specifically, we consider an ion confined in a
time-dependent radio-frequency Paul trap with linear geometry, while the atom
is constrained to move into a quasi-one-dimensional optical waveguide within
the ion trap. We evaluate the impact of the ion intrinsic micromotion on the
resonance position. Thus, we solve the atom-ion dynamics semiclassically,
namely the atom dynamics is governed by the three-dimensional time-dependent
Schr\"odinger equation, whereas the ion motion is described by the classical
Hamilton equations. We find that the energy of the ion provided by the
oscillating radiofrequency fields can affect the resonance position
substantially. Notwithstanding, the peculiar phenomenology of those resonances
regarding perfect transmission and reflection is still observable. These
findings indicate that the intrinsic micromotion of the ion is not detrimental
for the occurrence of the resonance and that its position can be controlled by
the radiofrequency fields. This provides an additional mean for tuning atom-ion
interactions in low spatial dimensions. The study represents an important
advancement in the scattering physics of compound atomic quantum systems in
time-dependent traps.