arxiv: v2 [cond-mat.quant-gas] 27 Jul 2015

Transcription

1 Fast control of topological vortex formation in BEC by counter-diabatic driving arxiv:57.35v [cond-mat.quant-gas] 7 Jul 5 Shumpei Masuda, Utkan Güngördü, Xi Chen, 3 etsuo Ohmi, 4 and Mikio Nakahara 3,4 James Franck Institute, he University of Chicago, Chicago, Illinois 6637, USA Department of Physics and Astronomy and Nebraska Center for Materials and Nanoscience, University of Nebraska, Lincoln, Nebraska 68588, USA 3 Department of Physics, Shanghai University, 444 Shanghai, People s Republic of China 4 Research Center for Quantum Computing and Department of Physics, Kinki University, Higashi-Osaka, , Japan (Dated: October 3, 8) opological vortex formation has been known as the simplest method for vortex formation in BEC of alkali atoms. his scheme requires inversion of the bias magnetic field along the axis of the condensate, which leads to atom loss when the bias field crosses zero. In this Letter, we propose a scheme with which the atom loss is greatly suppressed by adding counter-diabatic magnetic field. A naive counter-diabatic field violates the Maxwell equations and we need to introduce an approximation to make it physically feasible. he resulting field requires an extra currents, which is experimentally challenging. Finally we solve this problem by applying a gauge transformation so that the counter-diabatic field is generated by controlling the original trap field with the additional control of the bias field. PACS numbers:.3.yy, , Fg, 3.75.Lm Demonstration of quantized vortices in a superfluid is a manifestation of the nonvanishing order parameter. opological structure of vortices often reflects the manifold of the order parameter space. It is, therefore, natural to seek for the method to experimentally demonstrate formation of vortices once BEC of alkali metal atoms has been realized. opological method of vortex formation is the simplest one compared with other proposals, where the spinor structure of the order parameter is fully utilized [ 7]. Although the vortex [6] and monopole [8, 9] in BEC have been observed experimentally, establishment of methods to create such interesting quantum states more accurately in time much shorter than coherence time is important for further experimental study and applications. A condensate in the weak field seeking state (WFSS) is trapped in a quadrupole magnetic field in the xy-plane with a uniform bias field along the z-axis first and, subsequently, the bias field is reversed from B z to B z adiabatically. hen a vortex with the winding number F is formed, were F is the quantum number of the hyperfine spin of the condensate. his theoretical proposal has been demonstrated experimentally by several groups [ 5]. A part of the condensate is lost while vortex formation takes place. his is due to the fact that the gap between the WFSS and the neutral state (NS) and the strong field seeking state (SFSS) is gradually closed near the center of the condensate during the vortex formation and the atoms in the trapping WFSS make transition to the NS and the SFSS through the Majorana flops. It is the purpose of this Letter to show that it is possible to suppress the Majorana flops so that more atoms are left in the trap when a vortex is formed. It is certainly desirable to have more atoms left in the trap since higher atomic density implies (i) easier measurement of physical quantities (ii) fighting against atom loss due to hyperfine spin changing collision and (iii) smaller vortex core size, which is advantageous for forming a vortex in an arbitrary position in the condensate [6]. Motivated by the concept of shortcuts to adiabaticity [7], we use the counter-diabatic formalism (or quantum transitionless algorithm) developed by Demirplak and Rice[8] and formulated by Berry [9, ] and demonstrated experimentally in []. he supplementary counter-diabatic field(cdf) thus obtained, is, however, unphysical in that it violates the Maxwell equation div B = and some approximation must be introduced so that the counterdiabatic control is physically feasible. he approximated CDF turns out to be a quadrupole field, which can be implemented with an ordinary Ioffe bars. However, the CDF is generated by another set of Ioffe bars, which is obtained by rotating the original confining Ioffe bars by π/4. hismeansthattwosetsofioffebarsarerequiredto physically realize this nonadiabatic control, which makes its experimental demonstration rather demanding. o overcome this difficulty, we apply the gauge transformation [] to make the combined magnetic field (the original confining field and the CDF superposed) parallel to the original confining quadrupole field. he modest price one has to pay by this gauge transformation is to modify the time-dependence of the bias field B z (t) from a linear behavior. Experimental demonstration is made much easier with this modification. WefirstregardeachatomwithhyperspinF asanindependent quantum system fixed at r and consider nonadiabatic control thereof. Let B(r,t) = (B (r),b z (t)) be the local magnetic field with the quadrupole field B (r) = B (x, y) and the uniform bias field B z(t) = B z ( t/) for t, B z (t) = B z () for t > where is the inversion time of B z (t). An atom at (r,t)

2 interacts with the magnetic field through the Hamiltonian H B (r,t) = γb(r,t) F, () where γ = µ B g F, µ B is the Bohr magneton, g F is the g-factor for the hyperfine spin F and F = (F x,f y,f z ) is (F + )-dimensional irreducbile representation of the Lie algebra su(). We take F = in the following for definiteness. For a given B(r, t), the Hamiltonian () has three normalized eigenstates WFSS = (B B z )e iφ B e iφ, B B +B z NS = B e iφ Bz e iφ, () B SFSS = B B (B +B z )e iφ B e iφ, B B z where B = B (r)+b z (t) and φ is the azimuthal angle. he corresponding eigenvalues are γ B,, and γ B, respectively, where γ = µ B /. he WFSS with the highestenergy γ B isthe onlystate thatcan bemagnetically trapped. We start with a pure WFSS at t = and want to keep as many atoms in the WFSS as possible at t =, when the vortex formation is complete, in spite of vanishing gaps among the three states at the origin (r = ) at t = /. his is realized by adding a CDF on the atom. Here we briefly summarize the counter-diabatic approach to nonadiabatic quantum control and apply it to our hyperfine spin control. Let H (t) be a time-dependent Hamiltonian and n(t) be an instantaneous eigenvector with the eigenvalue E n such that H (t) n(t) = E n (t) n(t). In the adiabatic approximation, the solution to the time dependent Schrödinger equation is ψ n (t) = e iγn(t) n(t), where γ n (t) = t dt E n (t ) i t dt n(t ) t n(t ). (3) Let U(t) = n ψ n(t) n() be the time-evolution operatorsuchthatu(t) : n() ψ n (t). heoperatoru(t) defines a Hamiltonian H(t) = i ( t U(t))U (t) for an arbitrary time evolution, not necessarily adiabatic, namely ψ n (t) is the exact solution of i t ψ n (t) = H(t) ψ n (t), (4) where ψ n () = n() is satisfied by definition. If we write H(t) = H (t) + H CD (t), the counter-diabatic Hamiltonian H CD (t) is written as H CD (t) = i n t n(t) n(t), (5) when the parallel condition, n(t) t n(t) =, is satisfied. Let us apply the counter-diabatic scheme to our hyperfine spin system by taking H and n as H (r,t) = H B (r,t) and the instantaneous eigenstates in Eqs. (). Note that the coordinate r here is just a parameter specifying the position of the atom in the condensate, and the geometric phase in Eq. (3) vanishes. By substituting Eq. () into (5), the counter-diabatic Hamiltonian is obtained as where H CD = γb CD F, (6) B CD (r,t) = γ B z ()B B (y,x,) (7) (r,t) is the CDF. his magnetic field looks as if it could be generated by four Ioffe bars that are obtained by rotating the original Ioffe bars by an angle π/4 around the z-axis. An important remark is in order here. When the counter-diabatic scheme is applied to a quantum system, it produces a counter-diabatic potential (the magnetic field in our case) that prevents the quantum system from escaping from an adiabatic time evolution. here is no guarantee, however, that the counter-diabatic potential will be physically feasible. In fact, B CD in Eq. (7) does not satisfy divb CD =. his is not unexpected since we did not consider the Maxwell equation while we obtained B CD. he way out of this problem is to fix the coordinate r to r in B (r,t) in the denominator of B(r,t) so that B CD (r,t) (y,x,) and divb CD vanishes. In our calculation below, we show that the CDF increases the number of atoms forming the vortex in a wide range of the inversion time by using several values of r. he performance may be further improved if we optimize the number of atoms left in the trap by varying r, which will be worked out in our future publication. We call the combined magnetic field as B = B +B CD. Now we solve the Gross-Pitaevskii equation (GPE) with the designed magnetic field B; i t Ψ m (r,t) = {h mn +g n δ mn Ψ p p +g s (Ψ l (F α ) lp Ψ p )(F α ) mn }Ψ n, where l,m,n,p {,,}, α {x,y,z}, h mn = α l,p ( M µ ) δ mn +B mn, B = γ B F, (8) with the mass of the atom M, and the chemical potential µ (the eigenvalue of the GPE at t = ). o begin with, we need to find the initial condition to solve the GPE (8). Assuming that the initial state is in the WFSS

3 B CD B 3 Ψ() = f(r) WFSS(), the condensate wave function is obtained by solving the stationary GPE ( M[ d r d ) β r dr dr ] 4r (7 8cosβ +cosβ) f(r) +γ B(r,)f(r)+gf 3 (r) = µf(r), (9) where g = g n + g s, B(r,t) = {[B (r) + B CD (r,t)] + Bz(t)} /, and the WFSS is parametrized as WFSS = cosβcosφ isinφ cosβsinφ icosφ, sinβ with [ ] β = tan B (r)+b CD (r,t). B z (t) Here we note that f(r) and B (r) + B CD (r,) are in fact the functions of r only. It turns out that B CD is negligibly small at t = and it can be ignored safely in solving Eq. (9). he time-dependent GPE is solved numerically and we summarize the results below. hroughout our calculation, we have taken the parameters of 3 Na atoms [3], g n =.378a 3 HO ω, g s = and B z () = G, B = 3 G/cm for which the harmonic oscillator length is found to be a HO 9.4 µm and ω erg. Lengths and energies are scaled by these parameters. In addition, we find ω L (r =,t = ) rad/s. he time scale τ = π/ω L (,).43 µs is a reasonable measure of adiabaticity. he chemical potential is found to be µ ω L (,) 3.66 in dimensionless units. Figure shows the time dependence of the ratio B CD / B for log (/τ) =.4 and r /a HO =, which is independent of r and depends only on t. he inset is the profile of the condensate at t =, which is obtained by solving Eq. (9) numerically. he units of time, length and f are τ, a HO and a 3/ HO, respectively. Now the GP equation (8) is solved to find interesting quantities. Let us show first how much atoms are left in the trap after the vortex formation takes place. o mimic the trap loss, we multiply the order parameter by h(r) =.5[ tanh((r r )/λ)] with r = 3a HO and λ = a HO at each time step of the numerical simulation [3]. In Fig. we show the snapshots of the amplitude of the order parameter n Ψ n for t with the CDF with the same parameter set as Fig.. he condensate in the trapping region at t = is mostly in WFSS, while SFSS and NS are gradually removed from the trap in the domain r > r due to h(r). he zero of the order parameter at r =, t = is a manifestation of the vortex of the WFSS. Fig. 3 shows the fraction of atoms N(t)/N() left in the trap at t. Since we replacedan unphysical B CD with a physicalb CD satisfying div B CD =, the result depends on the parameter r 4.5 f r 5 FIG.. (Color online) ime dependence of the ratio B CD / B, which is independent of r and depends only on t. Inset shows the bound state wave function obtained by solving Eq. (9) for 3 Na. r 7 FIG.. (Color online) Snapshots of the amplitude of the order parameter for t with the same parameter set as that in Fig.. ime intervals between snapshots are equal with each other. he dashed and thick solid curves correspond to t = and t =, respectively. he unit of n Ψn is a 3 HO..7 N(t) N() 4 log (/τ) r a HO.5.5 FIG. 3. (Color online) Fraction N(t)/N() of atoms left in the trap long time after a vortex is formed. Since t, the condensate is in the pure WFSS. shows the fraction without B CD while other symbols show the fractions with B CD for different r.

4 4 N W t ( ) 8 ratio r a HO.5.5 log (/τ) 4 FIG. 4. (Color online) Number of atoms N(t) for t in the control with the CDF normalized by N(t) in the control without the CDF. he ratio is larger for small. ( ) N works reasonablly well. Now an important observation is in order. In our original proposal, we need to prepare two sets of Ioffe bars to produce the original confining quadrupole field and the counter-diabatic magnetic field. Clearly this is demanding for experimentalists. All topological vortex formation experiments so far were conducted with a single set of Ioffe bars. Even if one could build a trap with two sets of Ioffe bars, aligning their centers exactly at the same place is practically impossible. o circumvent this problem, we introduce time-dependent gauge transformation so that the combined field B +B CD is rotated and the resulting field is in parallel to the original B. Let B α B = tan B () be the anglebetween B and B. If B is rotatedby α B, it becomes parallel to B and hence it can be generated by a single set of Ioffe bars. his rotation is implemented by the unitary transformation U(α B ) = e iαbfz. () Now the Zeeman term of the Hamiltonian is transformed as U(α B )γ B FU (α B ) = γ B F, () FIG. 5. (Color online) Instantaneous fraction N W(t)/N() of atoms in WFSS during vortex formation for log (/τ) =.4 and r = a HO. he dashed and solid curves correspond to the controls without and with the CDF, respectively. in the denominator of Eq. (7). Observe the prominent improvement in the ratio N(t)/N() for small. his is more clearly seen by plotting N(t) in the control with the CDF normalized by N(t) in the control without the CDF. Figure 4 shows the ratios with the same r as Fig. 3. he ratio is almost for = τ and r =.5a HO. Next, we show the CDF really suppresses the transitions WFSS NS and WFSS SFSS by looking at the projected atom numbers. For this purpose, we define the projection operators Π W = WFSS WFSS,Π N = NS NS and Π S = SFSS SFSS and evaluate the numbers of atoms in these states as N W = Ψ Π W Ψ,N N = Ψ Π N Ψ,N S = Ψ Π S Ψ. Figure 5 shows the time dependence of N W in the controls without CDF (dashed curve) and with CDF (solid curve). he figure clearly shows that the nonadiabatic transitions at t / are greatly suppressed and transitions take place only at around t = /. his means that our approximation introduced in the denominator of (7) is not justifiable at t /, where the original B CD diverges at r =. Otherwise, our simple approximation where B = B (x, y,). Of course, this is not the whole story and there is a price we need to pay. he Hamiltonian in the rotating frame acquires a gauge term iu t U = α B F z. (3) his term, proportional to F z, works as a magnetic field in the z-direction and the bias field B z (t) = B z ()( t/) is replaced by ( B z (t) = B z () t ) α B. (4) Note that α B depends on time but not on space coordinates. Numerical caluculation shows that the number of atoms in WFSS in the control with B and B z is exactly the same as that for two sets of Ioffe bars for t. Figure 6 shows B z (t) and B (t) / B for 3 Na with the parameters used in Fig. 5. In summary, we have proposed a method to suppress nonadiabatic transitions while topological vortex formation takes place in BEC of alkali atoms. he counterdiabatic field is generated by a set of Ioffe bars, which is obtained by rotating the original Ioffe bars producing the confining quadrupole field by π/4. Our numerical calculation demonstrates that nondaiabatic transtions are suppressed for any inverstion time and, in particular, it is most impressive for a small. We can further improve this scheme by applying a gauge transformation to a rotating frame so that the combined field B is parallel to

5 B z B z () FIG. 6. (Color online) ime dependenceof B z(t)/b z()(solid curve) and B z(t) (dashed line) for the same parameters used in Fig. 5. he vertical and horizontal axes of the inset are B (r,t) / B (r) and time with the range /4 t 3/4, respectively. B. his requires modulation of B z (t) from linear timedependence. We believe our proposal is experimentally feasible by simple modifications of the existing setup. MN would like to thank Yuki Kawaguchi and akeshi Kuwamoto for useful discussions. he work of SM is partly supported by JSPS postdoc fellowship. UG is supported in part by the NSF under Grant No. Phy- 456 and nsf-epscor 494. XC s work was partially supported by the NSFC (47493 and 6768), the Shuguang and Pujiang Program (4SU35 and 3PJ43), the Specialized Research Fund for the Doctoral Program (3383), the Program for Eastern Scholar. MN s work is supported in part by a opological Quantum Phenomena Grant-in Aid for Scientific Research on Innovative Areas (No. 33) from the Ministry of Education, Culture, Sports, Science and echnology (MEX) of Japan. MN, O and UG are also grateful to JSPS for partial support from Grants-in- Aid for Scientific Research (Grant No. 644). [] M. Nakahara,. Isohima, K. Machid, S.-I. Ogawa, and. Ohmi, Physica B 84, 7 (). []. Isoshima, M. Nakahara,. Ohmi, and K. Machida, Phys. Rev. A 6, 636 (). [3] S.-I. Ogawa, M. Möttönen, M. Nakahara,. Ohmi, and H. Shimada, Phys. Rev. A 66, 367 (). [4] M. Möttönen, N. Matsumoto, M. Nakahara, and. Ohmi, J. Phys.: Condens. Matter 4, 348 (). [5] M. Möttönen, V. Pietilä, and S. M. M. Virtanen, Phys. Rev. Lett. 99, 546 (7). [6] For a review, see V. Pietila, M. Möttönen and M. Nakahara, in Electromagnetic, Magnetostatic, and Exchange- Interaction Vortices in Confined Magnetic Structures, ed. by E. O. Kamenetskii, ransworld Research Network, 97 (8). [7] P. Kuopanportti, B. P. Anderson, and M. Möttönen, Phys. Rev. A 87, 3363 (3). [8] M. W. Ray, E. Ruokokoski, S. Kandel, M. Möttönen, and D. S. Hall, Nature 55, 657 (4). [9] M. W. Ray, E. Ruokokoski, K. iurev, M. Möttönen, and D. S. Hall, Science 348, 544 (5). [] A. E. Leanhardt, A. Görlitz, A. P. Chikkatur, D. Kielpinski, Y.Shin,D.E.Pritchard, andw.ketterle.phys.rev. Lett. 89, 943 (). [] A. E. Leanhardt, Y. Shin, D. Kielpinski, D. E. Pritchard, and W. Ketterle, Phys. Rev. Lett. 9, 443 (3). [] M. Kumakura,. Hirotani, M. Okano, Y. akahashi, and. Yabuzaki, Phys. Rev. A 73, 6365 (6). [3]. Isoshima, M. Okano, H. Yasuda, K. Kasa, J. A. M. Huhtamäki, M. Kumakura, and Y. akahashi, Phys. Revl. Lett. 99, 43 (7). [4]. Kuwamoto, H.Usuda, S.ojo, and.hirano, J.Phys. Soc. Jpn. 79, 344 (). [5] H. Shibayama, Y. Yasaku, and. Kuwamoto, J. Phys. B: At. Mol. Opt. Phys. 44, 753 (). [6] Y. Kawaguchi, private communication. [7] E. orrontegui, S. Ibáñez, S. Martnez-Garaot, M. Modugno, A. del Campo, D. Guéry-Odelin, A. Ruschhaupt, X. Chen, and J. G. Muga, Adv. Atom. Mol. Opt. Phys. 6, 7 (3). [8] M. Demirplak and S. A. Rice, J. Phys. Chem. A 7, 9937 (3). [9] M Berry, J. Phys. A: Math. Gen. 4, (9). [] X. Chen, I. Lizuain, A. Ruschhaupt, D. Guéry-Odelin, and J. G. Muga, Phys. Rev. Lett. 5, 33 (). [] M. G. Bason, M. Viteau, N. Malossi, P. Huillery, E. Arimondo, D. Ciampini, R. Fazio, V. Giovannetti, R. Mannella and O. Morsch, Nat. Phys. 8, 47 (). [] S. Ibáñez, X. Chen, E. orrontegui, J. G. Muga, and A. Ruschhaupt, Phys. Rev. Lett. 9, 43 ().