Skip to main content

Spatial non-adiabatic passage using geometric phases

Abstract

Quantum technologies based on adiabatic techniques can be highly effective, but often at the cost of being very slow. Here we introduce a set of experimentally realistic, non-adiabatic protocols for spatial state preparation, which yield the same fidelity as their adiabatic counterparts, but on fast timescales. In particular, we consider a charged particle in a system of three tunnel-coupled quantum wells, where the presence of a magnetic field can induce a geometric phase during the tunnelling processes. We show that this leads to the appearance of complex tunnelling amplitudes and allows for the implementation of spatial non-adiabatic passage. We demonstrate the ability of such a system to transport a particle between two different wells and to generate a delocalised superposition between the three traps with high fidelity in short times.

1 Introduction

Adiabatic techniques are widely used for the manipulation of quantum states. They typically yield high fidelities and possess a high degree of robustness. One paradigmatic example is stimulated Raman adiabatic passage (STIRAP) in three-level atomic systems [1–3]. STIRAP-like techniques have been successfully applied to a wide range of problems, and in particular, to the control of the centre-of-mass states of atoms in microtraps. This spatial analogue of STIRAP is called spatial adiabatic passage (SAP) and it relies on coupling different spatial eigenstates via a controllable tunnelling interaction [4]. It has been examined for cold atoms in optical traps [5–12] and for electrons trapped in quantum dots [13, 14]. The ability to control the spatial degrees of freedom of trapped particles is an important goal for using these systems in future quantum technologies such as atomtronics [9, 15, 16] and quantum information processing [17]. SAP has also been suggested for a variety of tasks such as interferometry [11], creating angular momentum [12], and velocity filtering [18]. It is also applicable to the classical optics of coupled waveguides [19, 20].

However, the high fidelity and robustness of adiabatic techniques comes at the expense of requiring long operation times. This is problematic as the system will therefore also have a long time to interact with an environment leading to losses or decoherence. To avoid this problem, we will show how one can speed-up processes that control the centre-of-mass state of quantum particles and introduce a new class of techniques which we refer to as spatial non-adiabatic passage. The underlying foundation for these are shortcuts to adiabaticity (STA) techniques, which have been developed to achieve high fidelities in much shorter total times, for a review see [21, 22]. Moreover, shortcuts are known to provide the freedom to optimise against undesirable effects such as noise, systematic errors or transitions to unwanted levels [22–31].

Implementing the STA techniques for spatial control requires complex tunnelling amplitudes. However, tunnelling frequencies are typically real. To solve this, we show that the application of a magnetic field to a triple well system containing a single charged particle (which could correspond to a quantum dot system [32–37]) can achieve complex tunnelling frequencies through the addition of a geometric phase. This then allows one to implement a counter-diabatic driving term [21, 22, 38–40] or, more generally, to design dynamics using Lewis-Riesenfeld invariants [41].

The paper is structured as follows. In the next section, we present the model we examine, namely a charged particle in a triple well ring system with a magnetic field in the centre. In Section 3, we introduce the spatial adiabatic passage technique in a three-level system and show that making one of the couplings imaginary allows the implementation of transitionless quantum driving. We then show, in Section 3.3, how to create inverse-engineering protocols in this system using Lewis-Riesenfeld invariants. Results for two such protocols, namely transport and generation of a three-trap superposition, are given in Section 4. Section 5 presents a more realistic one-dimensional continuum model for the system, where the same schemes are implemented. Finally, in Section 6, we review and summarise the results.

2 System model

We consider a charged particle trapped in a system of three localised potentials, between which the tunnel coupling can be changed in a time-dependent manner. In order to have coupling between all traps, they are assumed to be arranged along a ring and a magnetic field exists perpendicular to the plane containing the traps, see Figure 1. The particle will initially be located in one of the traps and we will show how to design spatial non-adiabatic passage protocols where a specific final state can be reached within a finite time and with high fidelity. Such a model could, for example, correspond to an electron trapped in an arrangement of quantum dots, where gate electrodes can be used to change the tunnelling between different traps [42]. Another option would be to use ion trapping systems [43], where ring configurations have been recently demonstrated [44–46]. In these systems, tunnelling of an ion has already been observed (and controlled by manipulating the radial confinement), as well as the Aharonov-Bohm phase [47] acquired due to the presence of an external magnetic field [44].

Figure 1
figure 1

Diagram of the system consisting of three coupled quantum wells and a localised magnetic field in the centre. The basis states and the couplings strengths used in the three-level approximation are indicated. The coordinate system for the continuous model in Section 5 is also shown. The distance between two traps along the ring is defined as l, so that the total circumference of the ring is 3l.

Let us start by considering the single-particle Schrödinger equation

$$ i\hbar\frac{\partial\psi}{\partial t} = \frac{1}{2m} (- i \hbar \nabla- q \vec{A} )^{2} \psi+ V \psi, $$
(1)

where m and q are the mass and charge of the particle, respectively, and V corresponds to the potential describing the trapping geometry. We assume that the vector potential is originating from an idealised point-like and infinitely long solenoid at the origin (creating a magnetic flux \(\Phi_{B}\)) and it is therefore given by \(\vec {A} = \frac{\Phi_{B}}{2 \pi r} \hat{e}_{\varphi}\) (for \(\vec{r} \neq0\)). Here r, φ, z are cylindrical coordinates and \(\hat{e}_{\varphi}\) is a unit vector in the φ direction.

At low energies such a system can be approximated by a three-level (3L) model, where each basis state, \(|j\rangle\), corresponds to the localised ground state in one of the trapping potentials (see Figure 1). These states are isolated when a high barrier between them exists, but when the barrier is lowered the tunnelling amplitude \(\Omega_{jk}\) between states \(|j\rangle\) and \(|k\rangle\) becomes significant.

The presence of the magnetic field leads to the particle acquiring an Aharonov-Bohm phase [47] whenever it moves (tunnels) between two different positions (traps). This phase is given by \(\phi _{j,k} = \frac{q}{\hbar}\int_{\vec{r}_{j}}^{\vec{r}_{k}}\vec{A}(\vec {r})\cdot d\vec{r}\), where \(\vec{r}_{j}\) is the position of the jth trap, and for consistency, we always chose the direction of the path of the integration to be anti-clockwise around the pole of the vector potential (at \(\vec{r} = 0\)). The effects of this phase on the tunnelling amplitudes is given through the Peierls phase factors [48–50], \(\exp (i \phi_{j,k} )\), and the Hamiltonian for the 3L system can be written as

$$ H = -\frac{\hbar}{2} \begin{pmatrix} 0 & \Omega_{12}e^{i\phi_{1,2}} & \Omega_{31}e^{-i\phi_{3,1}} \\ \Omega_{12}e^{-i\phi_{1,2}} & 0 & \Omega_{23}e^{i\phi_{2,3}} \\ \Omega_{31}e^{i\phi_{3,1}} & \Omega_{23}e^{-i\phi_{2,3}} & 0 \end{pmatrix} . $$
(2)

Here the \(\Omega_{jk}\) are the coupling coefficients in the absence of any vector potential. The total phase around a closed path containing the three traps is then given by

$$ \Phi\equiv\phi_{1,2}+\phi_{2,3}+\phi_{3,1} = \frac{q}{\hbar} \oint\vec {A}(\vec{r})\cdot d\vec{l} = \frac{q}{\hbar} \Phi_{B}, $$
(3)

and is non-zero due to the pole of the vector potential A⃗ at the origin.

To simplify the Hamiltonian (2) one can use the following unitary transformation, which only employs local phases,

$$ U= \begin{pmatrix} 1 & 0 & 0\\ 0 & e^{-i\phi_{1,2}} & 0\\ 0 & 0 & e^{-i (\phi_{1,2}+\phi_{2,3} )} \end{pmatrix} , $$
(4)

and transforms the Hamiltonian as

$$ H \rightarrow U^{\dagger} H U= - \frac{\hbar}{2} \begin{pmatrix} 0 & \Omega_{12} & \Omega_{31}e^{- i \Phi}\\ \Omega_{12} & 0 & \Omega_{23}\\ \Omega_{31}e^{i \Phi} & \Omega_{23} & 0 \end{pmatrix} , $$
(5)

so that two of the tunnelling amplitudes become real-valued.

A case of particular interest is when \(\Phi=\pi/2\), i.e., when the magnetic flux is \(\Phi_{B} = \pi\hbar/ 2 q\). In this case the Hamiltonian becomes

$$ H = -\frac{\hbar}{2} (\Omega_{12} K_{1}+ \Omega_{23} K_{2}+\Omega_{31} K_{3} ) , $$
(6)

where each \(K_{j}\) is a spin 1 angular momentum operator defined as

$$ K_{1}= \begin{pmatrix} 0 & 1 & 0\\ 1 & 0 & 0\\ 0 & 0 & 0 \end{pmatrix} ,\qquad K_{2}= \begin{pmatrix} 0 & 0 & 0\\ 0 & 0 & 1\\ 0 & 1 & 0 \end{pmatrix} ,\qquad K_{3}= \begin{pmatrix} 0 & 0 & -i\\ 0 & 0 & 0\\ i & 0 & 0 \end{pmatrix} , $$
(7)

satisfying \([K_{j}, K_{k}] = i \epsilon_{jkl} K_{l}\) and \(\epsilon_{jkl}\) is the Levi-Civita symbol [51]. This means that the tunnel coupling between \(|3\rangle\) and \(|1\rangle\) becomes purely imaginary. We will show in the next section that this allows for the implementation of spatial non-adiabatic passage processes by either applying a transitionless quantum driving protocol or by using Lewis-Riesenfeld invariants.

3 Processes in the three-level approximation

3.1 Adiabatic methods

A series of spatial adiabatic passage (SAP) techniques have been developed in recent years, which allows one to manipulate and control the external degrees of freedom of quantum particles in localised potentials with high fidelity [4]. The standard SAP protocol for the transport of a single particle in a triple well system [5, 13] is the spatial analogue of the quantum-optical STIRAP technique [1–3]. It involves three linearly arranged, degenerate trapping states, \(|j\rangle\) with \(j = 1, 2\mbox{ and }3\), that can be coupled through tunnelling by either changing the distance between the traps or lowering the potential barrier between them. The system in the 3L approximation is described by the Hamiltonian

$$ H_{0} = -\frac{\hbar}{2} (\Omega_{12} K_{1}+\Omega_{23} K_{2} ), $$
(8)

which has a zero-energy eigenstate of the form

$$ |\lambda_{0}\rangle = \cos\theta|1\rangle - \sin \theta|3\rangle\quad \text{with } \tan\theta= \Omega_{12}/ \Omega_{23} . $$
(9)

This state is often called the dark state and SAP consists of adiabatically following \(|\lambda _{0}\rangle\) from \(|1\rangle\) (at \(t=0\)) to \(-|3\rangle\) (at a final time \(t=T\)), effectively transporting the particle between the outer traps one and three. This corresponds to changing θ from 0 (\(\Omega_{23} \gg \Omega_{12}\)) to \(\pi/2\) (\(\Omega_{23} \ll\Omega_{12}\)). Hence in the case of ideal adiabatic following, trap two (located in the middle) is never populated.

3.2 Transitionless quantum driving

The main drawback of SAP is that it requires the process to be carried out adiabatically and therefore slowly compared to the energy gap [4]. If this requirement is not met, unwanted excitations will lead to imperfect transport. One way to specifically cancel possible diabatic transitions in STIRAP was discussed in [52] and a general approach for recovering adiabatic dynamics in a non-adiabatic regime is to use shortcuts to adiabaticity, such as transitionless quantum driving [38–40]. This technique consists of adding a counter-diabatic term to the original Hamiltonian, whose particular form is given as

$$ H_{\mathrm{CD}} = i \hbar\sum_{n} \bigl(| \partial_{t} \lambda _{n}\rangle \langle\lambda _{n}| - \langle\lambda_{n}|\partial_{t} \lambda_{n}\rangle |\lambda _{n}\rangle \langle\lambda _{n}| \bigr), $$
(10)

where the \(|\lambda_{n}\rangle\) are the eigenstates of \(H_{0}\). For the reference Hamiltonian in Eq. (8) this gives [40]

$$ H_{\mathrm{CD}} =- \frac{\hbar\Omega_{31}(t)}{2} K_{3}, \quad \text{with } \Omega_{31}(t) = 2 \dot{\theta}(t) =2 \biggl( \frac{\Omega_{23} \dot {\Omega}_{12} - \Omega_{12} \dot{\Omega}_{23}}{\Omega_{12}^{2} + \Omega _{23}^{2}} \biggr). $$
(11)

We will see in Section 4.1 how this exact same scheme can also be obtained using Lewis-Riesenfeld invariants.

Shortcuts to adiabaticity have been studied in the context of STIRAP [40, 53], i.e., population transfer between internal levels. Its spatial analogue is more challenging as it requires that the additional tunnelling coupling between sites one and three is imaginary (see the definition of \(K_{3}\) in Eq. (7)). However, the system we have presented here is ideal for this, as the system Hamiltonian Eq. (6) is already equal to the total Hamiltonian \(H_{0} + H_{\mathrm{CD}}\). Other methods to implement the imaginary coupling could be, for example, the use of artificial magnetic fields [54] or angular momentum states [55].

A heuristic but not rigorous explanation of why the coupling needs to be imaginary can be obtained by examining the two ‘paths’ the particle can take to move from trap one to trap three. The first is via SAP and leads to \(|1\rangle \to- |3\rangle\) whereas the second is via the direct coupling the shortcut introduces, which leads to \(|1\rangle \to i e^{i \Phi} |3\rangle\). One can then immediately see that for constructive interference of these two terms the phase needs to have the value \(\Phi= \pi/2\), which corresponds to the required imaginary coupling between states \(|1\rangle\) and \(|3\rangle\). It is also interesting to note that the coupling between traps one and three in the shortcut has the form of a π-pulse

$$ \int_{0}^{T} \Omega_{31}(t) \, dt = 2 \int_{0}^{T} \dot{\theta}(t) \, dt = 2 \bigl[ \theta(T) - \theta(0) \bigr] = \pi. $$
(12)

3.3 Invariant-based inverse engineering

Another method of designing shortcuts to adiabaticity is by means of inverse-engineering using Lewis-Riesenfeld (LR) invariants [41, 56]. In this section we will briefly review these methods and then apply them to our particular system to both transport the particle and create a superposition between the three wells.

A LR invariant for a Hamiltonian \(H(t)\) is a Hermitian operator \(I(t)\) satisfying [41]

$$ \frac{\partial I}{\partial t}+\frac{i}{\hbar} [H,I ]=0. $$
(13)

Since \(I(t)\) is a constant of motion it can be shown that it has time-independent eigenvalues. It can be further shown that a particular solution of the Schrödinger equation,

$$ i\hbar\partial_{t} \bigl\vert \psi(t) \bigr\rangle = H(t) \bigl\vert \psi(t) \bigr\rangle , $$
(14)

can be written as

$$ \bigl\vert \psi_{k}(t) \bigr\rangle =e^{i \alpha_{k}(t)} \bigl\vert \phi_{k}(t) \bigr\rangle , $$
(15)

where the \(|\phi_{k}(t)\rangle\) are the instantaneous eigenstates of \(H(t)\) and

$$ \alpha_{k}(t)=\frac{1}{\hbar} \int_{0}^{t} \bigl\langle \phi_{k}(s) \bigr\vert \bigl[i\hbar \partial_{s}-H(s) \bigr] \bigl\vert \phi_{k}(s) \bigr\rangle \, ds $$
(16)

are the LR phases. Hence a general solution to the Schrödinger equation can be written as

$$ \bigl\vert \psi(t) \bigr\rangle =\sum_{k} c_{k} \bigl\vert \psi _{k}(t) \bigr\rangle , $$
(17)

where the \(c_{k}\) are independent of time.

The idea behind inverse engineering using LR invariants is not to follow an instantaneous eigenstate of the \(H(t)\) as one would in the adiabatic case, but rather follow an eigenstate of \(I(t)\) (up to the LR phase). To guarantee that the eigenstates coincide at the beginning and the end of the process, it is necessary that the invariant and the Hamiltonian commute at these times, i.e.,

$$ \bigl[I(0),H(0) \bigr]= \bigl[I(T),H(T) \bigr]=0. $$
(18)

One is then free to choose how the state evolves in the intermediate time and once this is fixed, Eq. (13) determines how the Hamiltonian should vary with time to achieve those dynamics.

A LR invariant for a three-level system described by Eq. (6) can be written as

$$ I = -\sin\beta\sin\alpha K_{1}-\sin\beta\cos\alpha K_{2}+ \cos\beta K_{3} , $$
(19)

where α and β are time dependent functions which must fulfil the following relations (imposed by Eq. (13))

$$\begin{aligned}& \dot{\alpha} = \frac{\Omega_{12} \sin\alpha+ \Omega_{23} \cos \alpha }{2 \tan\beta} + \frac{\Omega_{31}}{2}, \end{aligned}$$
(20)
$$\begin{aligned}& \dot{\beta} = \frac{1}{2} (\Omega_{23} \sin \alpha- \Omega_{12} \cos \alpha). \end{aligned}$$
(21)

The eigenstates of this invariant are

$$\begin{aligned}& \bigl\vert \phi_{0}(t) \bigr\rangle = \left ( \begin{matrix} -\sin\beta\cos\alpha\\ -i\cos\beta\\ \sin\beta\sin\alpha \end{matrix} \right ), \end{aligned}$$
(22)
$$\begin{aligned}& \bigl\vert \phi_{\pm}(t) \bigr\rangle = \frac{1}{\sqrt{2}}\left ( \begin{matrix} \cos\beta\cos\alpha\pm i\sin\alpha\\ -i\sin\beta\\ -\cos\beta\sin\alpha\pm i\cos\alpha \end{matrix} \right ), \end{aligned}$$
(23)

with respective eigenvalues \(\mu_{0}=0\) and \(\mu_{\pm}=\pm1\). One solution of the time-dependent Schrödinger equation is then given by \(|\Psi(t)\rangle = |\phi_{0} (t)\rangle\) as the corresponding LR phase is zero in this case. Note that this invariant is a generalisation of the invariant considered in [57] where a third coupling \(\Omega _{31}\) was not taken into account.

After fixing the boundary conditions using Eq. (18), one is free to choose the functions \(\alpha(t)\) and \(\beta(t)\). Moreover, in this case, one is also free to directly choose the function \(\Omega_{31}\). By inverting Eqs. (20) and (21), the other coupling coefficients are then given by

$$\begin{aligned}& \Omega_{12} = 2 \dot{\alpha}\sin\alpha\tan\beta- 2 \dot{\beta}\cos \alpha- \Omega_{31} \sin\alpha\tan\beta, \end{aligned}$$
(24)
$$\begin{aligned}& \Omega_{23} = 2 \dot{\alpha}\cos\alpha\tan\beta+ 2 \dot{\beta}\sin \alpha- \Omega_{31} \cos\alpha\tan\beta. \end{aligned}$$
(25)

4 Examples of spatial non-adiabatic passage schemes

In the following we will discuss two examples of spatial non-adiabatic passage derived from LR invariant based inverse engineering in the 3L approximation. The first one is the transport between two different traps, which is shown to be equivalent to the transitionless quantum driving method from Section 3 in some cases. The second scheme will create an equal superposition of the particle in all three traps.

4.1 Transport

The first example of control we examine is the population transfer determined by

$$ \bigl\vert \Psi(0) \bigr\rangle =\vert 1 \rangle \rightarrow \vert \Psi_{\mathrm{target}} \rangle = \bigl\vert \Psi(T) \bigr\rangle =- \vert 3 \rangle, $$
(26)

which was considered in the optical regime in [40]. This can be achieved by choosing auxiliary functions that fulfil the boundary conditions

$$ \beta(0)= \beta(T)= - \frac{\pi}{2},\qquad \alpha(0)=0, \quad \text{and} \quad \alpha(T)=\frac{\pi}{2}. $$
(27)

The experimentally required tunnelling frequencies are then explicitly given by Eqs. (24) and (25).

For the special choice of \(\beta(t) = -\pi/2\), one can show that \(\langle2|\Psi(t)\rangle = 0\) for all times, i.e. trap two is never occupied during the process. This choice then results in

$$ \tan\alpha= \frac{\Omega_{12}}{\Omega_{23}} \quad \text{and}\quad \Omega_{31} = 2\dot{\alpha}. $$
(28)

By identifying α with θ (see Eq. (9)) one can immediately see that this is the same pulse as in the STA scheme derived in Section 3.2.

The transport scheme can be implemented by the choosing the counterintuitive SAP pulses \(\Omega_{12}\) and \(\Omega_{23}\) to have a Gaussian profile [4]

$$\begin{aligned}& \Omega_{12}(t) = \Omega_{0} \exp \bigl[-100 (t/T - 1/2 )^{2} \bigr], \end{aligned}$$
(29)
$$\begin{aligned}& \Omega_{23}(t) = \Omega_{0} \exp \bigl[-100 (t/T - 1/3 )^{2} \bigr], \end{aligned}$$
(30)

and then calculating \(\Omega_{31}\) from Eq. (28). The resulting pulses and associated dynamical populations are shown in Figure 2. As expected the system follows exactly the dark state, transferring the population between states \(|1\rangle\) and \(|3\rangle\) without populating state \(|2\rangle\).

Figure 2
figure 2

Spatial non-adiabatic passage transport in the 3L approximation. \(T/\tau =100\) for \(\Omega _{0} \tau = 0.25\). (a) Modulus of the tunnelling amplitudes. (b) Evolution of the populations \(P_{i}=|\langle i|\Psi (t)\rangle|^{2}\). The time unit τ is defined as \(\tau = m l^{2}/\hbar \).

The fidelity of the transport process as a function of the total time and the phase Φ generated by the magnetic field is shown in Figure 3(a). Transport can be seen to occur with perfect fidelity for any value of the total time if the phase takes the appropriate value \(\Phi= \pi/2\). It can also be seen that the shortcut is successful for any value of the phase in the limit of very short or very long times. The latter one is not surprising, as \(\Omega_{31}\) can be neglected in the adiabatic limit, and hence its phase becomes irrelevant. A similar effect occurs for short total times, where the roles are reversed. In this limit \(\Omega_{31}\) is the largest of all three couplings, and hence the phase relation between it and the other couplings becomes inconsequential. As \(\Omega_{31}\) is a π pulse, perfect population transfer in this regime can be achieved regardless of the phase.

Figure 3
figure 3

Transport process \(\pmb{|1\rangle \to -|3\rangle}\) in the 3L approximation. (a) Fidelity as a function of the total time and the total magnetic phase traversing the system. The green contour line is defined by \(P_{3}=99\%\). (b) Probabilities of population in each of the traps for \(T/\tau =48\) (indicated by a dashed white line in (a)) as a function of the total magnetic phase traversing the system. The dashed black line indicates the optimal value of the phase \(\Phi =\pi /2\).

However, in order to maintain this pulse area, a strong coupling is required for very short processes, as the strength of \(\Omega_{31}\) is inversely proportional to T. This sets a bound on how fast this scheme can be implemented, as any physical implementation will have a maximum tunnelling amplitude. Setting the maximum value of \(\Omega_{31}\) to \(0.25/\tau\), the minimum process times T to achieve fidelities above 99% are approximately 880Ï„ for SAP and 100Ï„ for the shortcut scheme. These times are similar to the ones achievable in a spin-dependent transport scheme recently presented by Masuda et al. [58], however the setup in their work requires four traps and a constant and an AC magnetic field.

It is worth noting that this system also allows for the possibility of measuring the magnetic flux \(\Phi_{B}\), as the amount of transferred population oscillates as a function of the total phase Φ, which is directly related to the magnetic flux as \(\Phi=\frac{q}{\hbar}\Phi_{B}\). As an example we show the occupation probabilities for \(T/\tau= 48\) in each trap at the end of the process as a function of the phase in Figure 3(b). One can see that the populations strongly depend on the phase and over a large range of values one can therefore determine the magnetic flux. The exact relationship between the probabilities and the magnetic flux differs for different total times T.

4.2 Creation of a three-trap superposition

The second scheme we discuss highlights the generality of the LR invariant based method. In this scheme we create an equal superposition state between the particles being in all three traps, which means that the initial and target states are

$$ \bigl\vert \Psi(0) \bigr\rangle = \vert 1 \rangle \rightarrow \vert \Psi_{\mathrm{target}} \rangle = \bigl\vert \Psi (T) \bigr\rangle = \frac{1}{\sqrt{3}} \bigl( \vert 1 \rangle- i \vert 2 \rangle - \vert 3 \rangle\bigr). $$
(31)

This can be realised by imposing the boundary conditions

$$\begin{aligned}& \beta(0)=-\frac{\pi}{2},\qquad \beta(T)=-\arctan\sqrt{2}, \end{aligned}$$
(32)
$$\begin{aligned}& \alpha(0)=0,\qquad \alpha(T)=\frac{\pi}{4}, \end{aligned}$$
(33)

on the auxiliary functions. A simple ansatz which fulfils these boundary conditions is a fourth order polynomial for \(\beta(t)\) and third order polynomials for \(\alpha (t)\) and \(\Omega_{31}(t)\). The pulses are then obtained from Eqs. (24) and (25) and their form is shown in Figure 4(a). From Figure 4(b) it can be seen that this choice creates the target state at the final time with perfect fidelity.

Figure 4
figure 4

Spatial non-adiabatic superposition scheme \(\pmb{\vert 1 \rangle \to\frac{1}{\sqrt{3}} ( \vert 1 \rangle - i \vert 2 \rangle - \vert 3 \rangle)}\) in the 3L approximation. \(T/\tau =400\). Sub-figures are the same as in Figure 2 and the fidelity shown in (b) is defined as \(F=\vert \langle \Psi _{\mathrm{target}}|\Psi (t)\rangle \vert ^{2}\).

5 Spatial non-adiabatic passage in the continuum model

While the 3L approximation discussed above gives a clear picture of the physics of the system, it does not include effects such as excitations to higher energy states that can occur during the process. We will therefore in the following test the approximation by numerically integrating the full Schrödinger equation in real space. For this, we will consider traps that are narrow enough to limit the system dynamics to an effectively one-dimensional setting along the azimuthal coordinate, \(x = \varphi R\), i.e., around a circle of radius R, see Figure 1. Moreover, we will assume that the magnetic field is characterised by a vector potential in the azimuthal direction, \(\vec{A} = A \hat {e}_{\varphi}\).

We are therefore dealing with a one-dimensional system of length \(2 \pi R\) with periodic boundary conditions, whose dynamics are described by the following Schrödinger equation

$$ i \hbar\frac{\partial\psi}{\partial t} = \frac{1}{2m} \biggl(- i \hbar \frac{\partial}{\partial x} - q A \biggr)^{2} \psi+ V(x) \psi. $$
(34)

We assume a constant vector potential throughout the dynamical part of the protocols, as any time-varying vector potential would produce an unwanted force due to the electric field \(\vec{E} = - \partial_{t} \vec{A}\).

In order to be able to apply a well-defined phase we model the trapping sites as highly localised point-like potentials of depth \(\epsilon_{j}\) at the positions \(x_{j} = j l - l/2\) (see Figure 5). They are separated by square barriers of heights \(V_{jk}(t)\) (and length l), giving a total potential

$$ V(x,t) = - \sum_{j=1}^{3} \epsilon_{j}(t) \delta(x - x_{j})+ \textstyle\begin{cases} V_{31}(t) & \text{if } 0 < x < x_{1}, \\ V_{12}(t) & \text{if } x_{1} < x < x_{2}, \\ V_{23}(t) & \text{if } x_{2} < x < x_{3}, \\ V_{31}(t) & \text{if } x_{3} < x < 3l . \end{cases} $$
(35)

Since point-like potentials are difficult to implement numerically, in the simulations below they are implemented as narrow Gaussians. It is important to note that this model is not designed to give realistic estimates for the fidelities or exactly reproduce the dynamics of the 3L approximation. It is a toy model to validate the basic underlying processes and show that our schemes also make sense in the continuum.

Figure 5
figure 5

Schematic of the potential used in the numerical simulations (black line) with the localised states in each trap (coloured areas). The Gaussian shape of the traps is exaggerated here for clarity.

As mentioned above, the tunnelling amplitudes \(\Omega_{jk}(t)\) in the 3L approximation are related to the barrier heights \(V_{jk}(t)\) of the continuum model, see the Appendix. However, changing the barrier heights in order to achieve tunnelling will also affect the energies of the localised states in the neighbouring traps. Therefore, in order to reproduce the resonance of the 3L approximation (where the diagonal elements of the Hamiltonian are always zero) in the continuum model, the depths of the delta potentials \(\epsilon_{j}\) have to be adjusted as the barriers heights change, see Figure 5. Finally, to map the barrier heights \(V_{jk}\) and trap depths \(\epsilon _{j}\) parameters of the continuum model to the tunnelling amplitudes \(\Omega_{jk}\) of the 3L approximation, we numerically calculate the overlaps of neighbouring delta-trap eigenstates.

Results for transport of a particle using the shortcut scheme described in Section 4.1 are shown in Figure 6 and the barrier heights and trap depths used to match the pulses given in Figure 2 are shown in Figure 6(a), (b). The probability density during the process can be seen in Figure 6(c) and the populations in each trap are given in Figure 6(d). While the process is not perfect, one can see that the particle is transported to the final trap with a fidelity of 87%. The effect of the magnetic field can be seen in Figure 6(e), (f), where we show results for the same process but with an inverted magnetic field (using a total phase of \(\Phi= -\pi/2\)). In this case the interference between the adiabatic and shortcut paths is destructive, and almost no population ends up in the final trap.

Figure 6
figure 6

Spatial non-adiabatic transport process in the continuum model. \(T/\tau =100\). (a), (b) Barrier heights and trap depths obtained by mapping the couplings in Figure 2(a). (c) Evolution of the particle density \(|\psi (x,t)|^{2}\). (d) Corresponding populations \(P_{i}=\vert \langle i|\Psi (t)\rangle \vert ^{2}\) in each trap and of the target state. (e), (f) are the same as (c), (d) but with the magnetic flux flowing in the opposite direction. The width of the Gaussian traps is 10−4 l.

The results for the creation of the superposition state discussed in Section 4.2 are shown in Figure 7. The observed dynamics are very similar to the one in the 3L approximation and the process reaches a final fidelity of the target state of 91%.

Figure 7
figure 7

Same as Figure  6 (a)-(d) but for the spatial non-adiabatic superposition scheme given in Eq. ( 31 ) in the continuum model. \(T/\tau =400\). \(F=\vert \langle \psi _{\mathrm{target}}|\psi (t)\rangle \vert ^{2}\) is the fidelity of the process.

Since the continuum model has many more degrees of freedom than the 3L model, it is not surprising that the fidelities obtained are lower. Nevertheless, the basic functioning of our spatial non-adiabatic techniques is clearly established from the calculations shown above. Optimising the fidelity in the continuum is an interesting task which, however, goes beyond the scope of the current work.

6 Conclusions and outlook

We have shown how complex tunnel frequencies in single-particle systems allow one to develop spatial non-adiabatic passage techniques that can lead to fast and robust processes for quantum technologies. In particular, we have discussed the case of a single, charged particle in a microtrap environment. The complex tunnelling couplings are obtained from the addition of a constant magnetic field, and have allowed us to generalise adiabatic state preparation protocols beyond the usual spatial adiabatic passage techniques [4]. This demonstrates that non-adiabatic techniques can be as efficient as their adiabatic counterparts, without requiring the long operation times.

In particular, we have discussed the implementation of the counter-diabatic term for spatial adiabatic passage transport via a direct coupling of all the traps. This was, in a second step, generalised to a flexible and robust method for preparing any state of the single-particle system by using Lewis-Riesenfeld invariants. As an example, we have shown that an equal spatial superposition state between the three wells can be created on a short time scale. Finally, we have presented numerical evidence that spatial non-adiabatic processes work also in a one-dimensional toy model by introducing a mapping between the discrete three-level approximation and a continuum model.

While in this work we have focused on a three-trap system, an interesting extension would be to investigate similar schemes in larger systems, or in different physical settings (for example, superconducting qubits [59]). Often, if the transitionless quantum driving technique is directly applied to complex quantum systems, the additional counter-adiabatic terms become very complicated, hard to implement or even unphysical. Nevertheless, the steps outlined in our work (using a few-level approximation, applying the shortcut technique, and then mapping everything back to a continuous model) can in principle be applied to any trap configuration. These steps might lead to schemes which are much easier to implement experimentally than the direct application of the transitionless quantum driving. However, each of these generalised configuration would need to be studied on an individual basis.

It would also be very interesting to see the effect of interactions in this system. For very strong interactions such that double occupancy of a site is suppressed and a single empty site is present, one might expect to observe similar dynamics but for the empty site [9]. In this case, spatial non-adiabatic ideas can be straightforwardly transferred. For intermediate interaction strengths (but stronger than the tunnelling couplings), repulsively-bound pair processes have been shown to dominate the dynamics and single-particle-like dynamics can be recovered for the pair [10, 60, 61]. In this case the presented techniques might be extended for a particle pair.

Finally, it is also worth noting that these complex tunnelling couplings we introduce can be used to implement techniques based on composite pulses [62].

References

  1. Bergmann K, Theuer H, Shore BW. Rev Mod Phys. 1998;70:1003.

    Article  ADS  Google Scholar 

  2. Bergmann K, Vitanov NV, Shore BW. J Chem Phys. 2015;142:170901.

    Article  ADS  Google Scholar 

  3. Vitanov NV, Rangelov AA, Shore BW, Bergmann K. Rev Mod Phys. 2017;89:015006.

    Article  Google Scholar 

  4. Menchon-Enrich R, Benseny A, Ahufinger V, Greentree AD, Busch T, Mompart J. Rep Prog Phys. 2016;79:074401.

    Article  ADS  Google Scholar 

  5. Eckert K, Lewenstein M, Corbalán R, Birkl G, Ertmer W, Mompart J. Phys Rev A. 2004;70:023606.

    Article  ADS  Google Scholar 

  6. Eckert K, Mompart J, Corbalán R, Lewenstein M, Birkl G. Opt Commun. 2006;264:264.

    Article  ADS  Google Scholar 

  7. McEndoo S, Croke S, Brophy J, Busch T. Phys Rev A. 2010;81:043640.

    Article  ADS  Google Scholar 

  8. Gajdacz M, Opatrný T, Das KK. Phys Rev A. 2011;83:033623.

    Article  ADS  Google Scholar 

  9. Benseny A, Fernández-Vidal S, Bagudà J, Corbalán R, Picón A, Roso L, Birkl G, Mompart J. Phys Rev A. 2010;82:013604.

    Article  ADS  Google Scholar 

  10. Benseny A, Gillet J, Busch T. Phys Rev A. 2016;93:033629.

    Article  ADS  Google Scholar 

  11. Menchon-Enrich R, McEndoo S, Busch T, Ahufinger V, Mompart J. Phys Rev A. 2014;89:053611.

    Article  ADS  Google Scholar 

  12. Menchon-Enrich R, McEndoo S, Mompart J, Ahufinger V, Busch T. Phys Rev A. 2014;89:013626.

    Article  ADS  Google Scholar 

  13. Greentree AD, Cole JH, Hamilton AR, Hollenberg LCL. Phys Rev B. 2004;70:235317.

    Article  ADS  Google Scholar 

  14. Fountoulakis A, Paspalakisa E. J Appl Phys. 2013;113:174301.

    Article  ADS  Google Scholar 

  15. Seaman BT, Krämer M, Anderson DZ, Holland MJ. Phys Rev A. 2007;75:023615.

    Article  ADS  Google Scholar 

  16. Pepino RA, Cooper J, Anderson DZ, Holland MJ. Phys Rev Lett. 2009;103:140405.

    Article  ADS  Google Scholar 

  17. Jaksch D, Briegel H-J, Cirac JI, Gardiner CW, Zoller P. Phys Rev Lett. 1999;82:1975.

    Article  ADS  Google Scholar 

  18. Loiko Y, Ahufinger V, Menchon-Enrich R, Birkl G, Mompart J. Eur Phys J D. 2014;68:147.

    Article  ADS  Google Scholar 

  19. Longhi S. Phys Rev E. 2006;73:026607.

    Article  ADS  Google Scholar 

  20. Longhi S, Della Valle G, Ornigotti M, Laporta P. Phys Rev B. 2007;76:201101.

    Article  ADS  Google Scholar 

  21. Torrontegui E, Ibáñez S, Martínez-Garaot S, Modugno M, del Campo A, Guéry-Odelin D, Ruschhaupt A, Chen X, Muga JG. Adv At Mol Opt Phys. 2013;62:117.

    Article  ADS  Google Scholar 

  22. Ruschhaupt A, Muga JG. J Mod Opt. 2013;61:828.

    Article  ADS  Google Scholar 

  23. Ruschhaupt A, Chen X, Alonso D, Muga JG. New J Phys. 2012;14:093040.

    Article  Google Scholar 

  24. Daems D, Ruschhaupt A, Sugny D, Guérin S. Phys Rev Lett. 2013;111:050404.

    Article  ADS  Google Scholar 

  25. Lu XJ, Chen X, Ruschhaupt A, Alonso D, Guérin S, Muga JG. Phys Rev A. 2013;88:033406.

    Article  ADS  Google Scholar 

  26. Kiely A, Ruschhaupt A. J Phys B, At Mol Opt Phys. 2014;47:115501.

    Article  ADS  Google Scholar 

  27. Guéry-Odelin D, Muga JG. Phys Rev A. 2014;90:063425.

    Article  ADS  Google Scholar 

  28. Lu XJ, Muga JG, Poschinger UG, Schmidt-Kaler F, Ruschhaupt A. Phys Rev A. 2014;89:063414.

    Article  ADS  Google Scholar 

  29. Zhang Q, Chen X, Guéry-Odelin D. Phys Rev A. 2015;92:043410.

    Article  ADS  Google Scholar 

  30. Kiely A, Benseny A, Busch T, Ruschhaupt A. J Phys B, At Mol Opt Phys. 2016;49:215003.

    Article  ADS  Google Scholar 

  31. Zhang Q, Muga JG, Guéry-Odelin D, Chen X. J Phys B, At Mol Opt Phys. 2016;49:125503.

    Article  ADS  Google Scholar 

  32. Hsieh C-Y, Shim Y-P, Korkusinski M, Hawrylak P. Rep Prog Phys. 2012;75:114501.

    Article  Google Scholar 

  33. Domínguez F, Platero G, Kohler S. Chem Phys. 2010;375:284.

    Article  ADS  Google Scholar 

  34. Huneke J, Platero G, Kohler S. Phys Rev Lett. 2013;110:036802.

    Article  ADS  Google Scholar 

  35. Jong LM, Greentree AD. Phys Rev B. 2010;81:035311.

    Article  ADS  Google Scholar 

  36. Mousolou VA. Europhys Lett. 2017;117:10006.

    Article  Google Scholar 

  37. Zeng Q-B, Chen S, Lü R. arXiv:1608.00065 [quant-ph].

  38. Demirplak M, Rice SA. J Phys Chem A. 2003;107:9937.

    Article  Google Scholar 

  39. Berry MV. J Phys A. 2009;42:365303.

    Article  MathSciNet  Google Scholar 

  40. Chen X, Lizuain I, Ruschhaupt A, Guéry-Odelin D, Muga JG. Phys Rev Lett. 2010;105:123003.

    Article  ADS  Google Scholar 

  41. Lewis HR, Riesenfeld WB. J Math Phys. 1969;10:1458.

    Article  ADS  Google Scholar 

  42. Braakman FR, Barthelemy P, Reichl C, Wegscheider W, Vandersypen LMK. Nat Nanotechnol. 2013;8:432.

    Article  ADS  Google Scholar 

  43. Seidelin S, Chiaverini J, Reichle R, Bollinger JJ, Leibfried D, Britton J, Wesenberg JH, Blakestad RB, Epstein RJ, Hume DB, Itano WM, Jost JD, Langer C, Ozeri R, Shiga N, Wineland DJ. Phys Rev Lett. 2006;96:253003.

    Article  ADS  Google Scholar 

  44. Noguchi A, Shikano Y, Toyoda K, Urabe S. Nat Commun. 2014;5:3868.

    Article  ADS  Google Scholar 

  45. Tabakov B, Benito F, Blain M, Clark CR, Clark S, Haltli RA, Maunz P, Sterk JD, Tigges C, Stick D. Phys Rev Appl. 2015;4:031001.

    Article  ADS  Google Scholar 

  46. Yoshimura B, Stork M, Dadic D, Campbell WC, Freericks JK. EPJ Quantum Technol. 2015;2:2.

    Article  Google Scholar 

  47. Aharonov Y, Bohm D. Phys Rev. 1959;115:485.

    Article  ADS  MathSciNet  Google Scholar 

  48. Graf M, Vogl P. Phys Rev B. 1995;51:4940.

    Article  ADS  Google Scholar 

  49. Ismail-Beigi S, Chang EK, Louie SG. Phys Rev Lett. 2001;87:087402.

    Article  ADS  Google Scholar 

  50. Cehovin A, Canali CM, MacDonald AH. Phys Rev B. 2004;69:045411.

    Article  ADS  Google Scholar 

  51. Carroll CE, Hioe FT. J Opt Soc Am B. 1988;5:1335.

    Article  ADS  Google Scholar 

  52. Unanyan RG, Yatsenko LP, Bergmann K, Shore BW. Opt Commun. 1997;139:48.

    Article  ADS  Google Scholar 

  53. Du YX, Liang ZT, Li YC, Yue XX, Lv QX, Huang W, Chen X, Yan H, Zhu SL. Nat Commun. 2016;7:12479.

    Article  ADS  Google Scholar 

  54. Dalibard J, Gerbier F, Juzeliūnas G, Öhberg P. Rev Mod Phys. 2011;83:1523.

    Article  ADS  Google Scholar 

  55. Polo J, Mompart J, Ahufinger V. Phys Rev A. 2016;93:033613.

    Article  ADS  Google Scholar 

  56. Chen X, Ruschhaupt A, Schmidt S, del Campo A, Guéry-Odelin D, Muga JG. Phys Rev Lett. 2010;104:063002.

    Article  ADS  Google Scholar 

  57. Chen X, Muga JG. Phys Rev A. 2012;86:033405.

    Article  ADS  Google Scholar 

  58. Masuda S, Tan KY, Nakahara M. arXiv:1612.08389 [cond-mat.mes-hall].

  59. Roushan P, Neill C, Megrant A, Chen Y, Babbush R, Barends R, Campbell B, Chen Z, Chiaro B, Dunsworth A, Fowler A, Jeffrey E, Kelly J, Lucero E, Mutus J, O’Malley PJJ, Neeley M, Quintana C, Sank D, Vainsencher A, Wenner J, White T, Kapit E, Neven H, Martinis J. Nat Phys 2017;13:146.

    Article  Google Scholar 

  60. Bello M, Creffield CE, Platero G. Sci Rep. 2016;6:22562.

    Article  ADS  Google Scholar 

  61. Bello M, Creffield CE, Platero G. Phys Rev B. 2017;95:094303.

    Article  Google Scholar 

  62. Torosov BT, Vitanov NV. Phys Rev A. 2011;83:053420.

    Article  ADS  Google Scholar 

Download references

Acknowledgements

This work has received financial support from Science Foundation Ireland under the International Strategic Cooperation Award Grant No. SFI/13/ISCA/2845 and the Okinawa Institute of Science and Technology Graduate University. We are grateful to David Rea for useful discussion and commenting on the manuscript.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Albert Benseny.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

AB and AK performed the calculations. All authors participated in the design and discussion of the project, and contributed equally to the writing of this paper. All authors read and approved the final manuscript.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Appendix:  Mapping between three-level approximation and continuum model

Appendix:  Mapping between three-level approximation and continuum model

In this appendix we give more details on how to connect the parameters of the three-level (3L) approximation and the continuum model. For clarity, we set \(\hbar=m=1\) in the following. Let us first recall the eigenfunctions of a single asymmetric delta potential given by

$$ V(x)=-\epsilon\delta(x)+ \textstyle\begin{cases} V_{\mathrm{L}},& x< 0, \\ V_{\mathrm{R}},& x \geq0. \end{cases} $$
(36)

This potential has only one bound state (as long as \(2 \epsilon^{2} > |V_{\mathrm{L}}- V_{\mathrm{R}}|\)) which is of the form

$$ \psi(x)= \textstyle\begin{cases} \phi_{\mathrm{L}}(x), & x< 0, \\ \phi_{\mathrm{R}}(x), & x \geq0, \end{cases} $$
(37)

where

$$\begin{aligned}& \phi_{\mathrm{L}}(x) = \exp \biggl[\frac{ (2 \epsilon ^{2}+V_{\mathrm{L}}-V_{\mathrm{R}} )x}{2 \epsilon} \biggr] \frac{\sqrt{4 \epsilon^{4}- (V_{\mathrm{L}}-V_{\mathrm{R}} )^{2}}}{2 \epsilon^{3/2}}, \end{aligned}$$
(38)
$$\begin{aligned}& \phi_{\mathrm{R}}(x) = \exp \biggl[-\frac{ (2 \epsilon ^{2}+V_{\mathrm{R}}-V_{\mathrm{L}} )x}{2 \epsilon} \biggr] \frac{\sqrt{4 \epsilon^{4}- (V_{\mathrm{L}}-V_{\mathrm{R}} )^{2}}}{2 \epsilon^{3/2}}, \end{aligned}$$
(39)

with an energy

$$ E=-\frac{4 \epsilon^{4}+ (V_{\mathrm{R}}-V_{\mathrm{L}} )^{2}-4 \epsilon ^{2} (V_{\mathrm{L}} +V_{\mathrm{R}} )}{8 \epsilon^{2}}. $$
(40)

In our work we use these eigenstate as the localised basis states in each of the three delta trap potentials. For example, the basis state \(\psi_{1}(x)\) for the first trap can be constructed from the substitutions \(\epsilon\rightarrow\epsilon_{1}\), \(V_{\mathrm{L}}\rightarrow V_{13}\), \(V_{\mathrm{R}}\rightarrow V_{12}\) and \(x \rightarrow x-l/2\) (see Figure 5) and the states \(\psi_{2}(x)\) and \(\psi _{3}(x)\) for the other two wells can be obtained in a similar manner. While choosing a basis for the system this way does not necessarily lead to an orthogonal basis set, we have checked numerically that the states are approximately orthogonal at all times during our simulations. This allows us to approximate the Hamiltonian associated with Eq. (34) as Eq. (6).

The couplings constants between each pair of neighbouring two traps can be determined by calculating the overlap between the two respective trap states in the barrier region between them, i.e.,

$$\begin{aligned}& \Omega_{12} \approx-2 \int_{l/2}^{3l/2} \psi_{1}(x) \biggl[- \frac {1}{2}\, \partial_{x}^{2} \psi_{2}(x)+V_{12} \psi_{2}(x) \biggr] \,dx, \end{aligned}$$
(41)
$$\begin{aligned}& \Omega_{23} \approx-2 \int_{3l/2}^{5l/2} \psi_{2}(x) \biggl[- \frac {1}{2}\, \partial_{x}^{2} \psi_{3}(x)+V_{23} \psi_{3}(x) \biggr] \,dx, \end{aligned}$$
(42)
$$\begin{aligned}& \Omega_{31} \approx-2 \int_{5l/2}^{7l/2} \psi_{1}(x-3l) \biggl[- \frac {1}{2}\, \partial_{x}^{2} \psi_{3}(x)+V_{31} \psi_{3}(x) \biggr] \,dx. \end{aligned}$$
(43)

Similarly, the on-site energies (or diagonal elements) are approximated by considering only the regions for which the basis states are significant, e.g.,

$$\begin{aligned} E_{2} \approx& \int_{l/2}^{3l/2} \psi_{2}(x) \biggl[- \frac {1}{2}\, \partial _{x}^{2} \psi_{2}(x)+V_{12} \psi_{2}(x) \biggr] \, dx \\ &{}+ \int_{3l/2}^{5l/2} \psi_{2}(x) \biggl[- \frac{1}{2}\, \partial_{x}^{2} \psi _{2}(x)+V_{23} \psi_{2}(x) \biggr] \, dx \end{aligned}$$
(44)

and correspondingly for the energies \(E_{1}\) and \(E_{3}\). We can then tabulate sets of values for the barrier heights, \(\{ V_{12},V_{23},V_{13}\}\), trap depths, \(\{ \epsilon_{1},\epsilon _{2},\epsilon _{3}\}\), and tunnelling amplitudes, \(\{ \Omega_{12},\Omega_{23},\Omega _{13}\}\), such that the energies all match a fixed reference value, i.e., \(E_{1}=E_{2}=E_{3}=E_{0}\) where \(E_{0}\) is fixed to some constant value. Since for a given protocol the required tunnelling amplitudes are known, we can finally numerically invert the table in order to determine how the barrier heights and trap depths have to vary in time.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Benseny, A., Kiely, A., Zhang, Y. et al. Spatial non-adiabatic passage using geometric phases. EPJ Quantum Technol. 4, 3 (2017). https://doi.org/10.1140/epjqt/s40507-017-0056-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1140/epjqt/s40507-017-0056-x

Keywords