Dynamics of a Bose-Einstein Condensate on Changing Speeds of an Atomchip Trap Potential
Dynamics of a Bose-Einstein Condensate on Changing Speeds of an Atomchip Trap Potential
Journal of the Optical Society of Korea. 2014. Dec, 18(6): 633-6388
Copyright © 2014, Optical Society of Korea
  • Received : July 07, 2014
  • Accepted : October 10, 2014
  • Published : December 25, 2014
Export by style
Cited by
About the Authors
Seung Jin Kim
Department of Physics Education, Korea National University of Education, Cheongju 363-791, Korea
Jae June Noh
Department of Physics Education, Korea National University of Education, Cheongju 363-791, Korea
Min Seok Kim
Department of Physics Education, Korea National University of Education, Cheongju 363-791, Korea
Jin Seung Lee
Agency for Defense Development, Daejon 305-600, Korea
Hoon Yu
Department of Physics and JILA, University of Colorado, Boulder, Colorado, 80309-0440, USA
Jung Bog Kim
Department of Physics Education, Korea National University of Education, Cheongju 363-791, Korea
We report experimental behaviors of condensed 87 Rb atoms responding to changes in the trap potential of the atomchip. The two-types of adiabatic and non-adiabatic overall changes were implemented by changing the ramp-down speed of the chip-wire current, which can dominantly modify the one-axis magnetic field gradient. Under the adiabatic process, a pure condensate stayed in the initial spin state and collectively oscillated with both monopole and dipole modes, while an atomic cloud above the critical temperature exhibited sound waves in a dense ultracold gas. On the other hand, Bose-Einstein condensate atoms with non-adiabatic perturbation were split into spatially different positions by spin states through spin-flip. We investigated the split ratio among spin states depending on final evaporation frequency. Potential changes, of course, cause collective oscillations regardless of the changing process.
Since the first realization of Bose-Einstein condensation in dilute gases confined in a magnetic trap [1 , 2] , intensive studies of the dynamical behavior of mesoscopic quantum systems have been performed [3 - 14] , and it was expected to explain, for instance, the superfluidity of liquid helium [15] . The study of collective excitations (normal modes) is one of the standard methods used to characterize the behavior of highly-degenerate quantum systems. Since the first measurements of the lowest normal mode frequencies [3 , 4] verified that the Gross-Pitaevskii equation (GPE) correctly predicted the dynamic behavior of condensate under excitation, the spectroscopy of the excited condensate has been extended to a variety of studies, e.g., frequency shifts [8 , 16] , mode coupling [8 , 16 , 17] and damping [6 , 18] .
In the hydrodynamic limit, which is applied to condensates with a large number of atoms, Stringari found the dispersion law for the elementary excitations of a dilute atomic Bose gas confined in both a harmonic trap and repulsive interacting forces [9] . For a cylindrically-symmetric potential,
PPT Slide
Lager Image
where M is the atomic mass and ωi is the angular frequency of oscillation along the i -axis, the m = 0 modes’ eigenfrequencies are given by [9] :
PPT Slide
Lager Image
Here - and + signs are to represent monopole and quadrupole mode, respectively. This prediction has been successfully testified by experimental results [3 , 4] .
Since low-lying collective modes have an energy on the order of ħωi for a harmonically confined Bose-Einstein Condensate (BEC), where ħ is the Planck constant divided by 2 π , in previous studies collective excitation in a BEC was performed by applying a series of sinusoidal modulations near the resonance frequency on the confining potential [3 - 6 , 8] . That is, the energy to excite condensate was provided from external sinusoidal perturbation.
From these previous results, we examined a possibility that a decrease of trap frequency also excites atoms with a contribution of the residual energy since modification of the trap geometry is the basic essential for the dynamical control of the ultracold atoms. The excitation behavior according to the adiabaticity of the deformation is also a topic of our question. For this, the atomchip is a good experimental apparatus because it can easily modify trap geometry at a micrometer scale within a wide range of changing times [19] . In this paper, we report on the observation of the dynamical behaviors of ultracold 87 Rb atoms on an atomchip under trap deformation. From the adiabatic deformation, we observed low-lying collective excitations and sound waves with pure condensate and a thermal cloud before the onset of a condensate, respectively. The frequencies were compared with theoretical predictions from the GPE in the hydrodynamic limit and to the sound wave. We were able to separate spatially the condensate into groups by different spin states by Majorana spin-flip through the non-adiabatic deformation. In addition, we investigated the population ratio of the different spin components by varying the final frequency in the rf evaporation step for the BEC.
The detailed apparatus and procedures we use for creating the BEC are described elsewhere [20] . In summary, our experimental setup consists of a compact double magnetooptical trap (MOT) apparatus with an anodically bonded atomchip which produces magnetically trapped 87 Rb Bose-Einstein condensates in the F = 2, m F = 2 state. An atomic flux of high density is generated from a 2D + MOT in the isolated 2D + MOT chamber [21] . The atoms pass through a 750 μm hole and are trapped in a 3D-MOT. The atoms in the 3D-MOT are further cooled below the Doppler limit through a compressed-MOT and polarization gradient cooling. By using an external Z-wire trap, the atoms are transported from the initial 3D-MOT position to near the atomchip surface after optical pumping. About 9 × 10 7 atoms can be trapped in an atomchip magnetic trap by the currents through the atomchip wires. After two-step rf forced evaporation, a BEC can be manipulated. Resonant absorption images were taken with the 100 μs illumination time, and the optical resolution was 3.5 μm.
The 26 × 26 mm 2 sized atomchip is formed by depositing 10 μm-thick copper conductor onto a 450 μm-thick silicon substrate ( Fig. 1(a) ). The central Z-shaped wire (called the Z-wire, 100 μm-width) is used with external bias fields to create an Ioffe-Pritchard trap. Higher trap frequency along the x-axis can be obtained with a dimple trap that is created by running a current through the wire that is perpendicular to the Z-wire (called dimple-wire, 100 μm-width). The two wires beside the Z-wire are designed for an rf-induced potential experiment and are not used for this study. To investigate the motion of the condensate under a trap deformation, the magnetic field of the dimple trap is modified by a ramp-down of the dimple current I d . Since the dimple current makes an x-directional field at the center of the trap, the x-directional trap frequency can be changed dominantly with the current ( Fig. 1(b) ). Fig. 1(c) is the measured anisotropy ( λ = ωx / ω ) of the deformed trap as a function of the final Id in the range of our performance. More reduction of the current makes the trap more anisotropic [22] .
PPT Slide
Lager Image
(a) Schematic of the atomchip used for magnetic trapping potentials. The central Z-shaped wire (B-G) carries a current and is used together with an external bias field along the y-axis to produce an Ioffe-Pritchard trap. The magnetic field gradient along the x-axis can be controlled with the current through the central straight wire (D-E), the dimple wire. (The two Z-shaped wires beside the Z-wire are not used for this study.)
For this study, initially, a BEC was prepared after rf forced evaporative cooling within 2.3 s which is sufficient to damp the atoms’ motions induced by the atomic cloud transfer. The prepared BEC typically contains 5.5 × 10 4 atoms without recognizable thermal fraction and is positioned at 180 μm below the atomchip surface with the Z-wire current I z = 2.98 A, the x-axis bias field B x = 6.86 G, the y-axis bias field B y = 33.10 G and the I d = 0.98 A.
Trap deformation for investigating dynamical behavior was performed in two ways, adiabatic and non-adiabatic changes, by changing the ramp-down duration of dimple current, which is determined by the speed changing trap potential.
- 3.1. Adiabatic Changing – Collective Excitations
For the adiabatic deformation, the dimple current I d was ramped down from 980 mA to its final value within 10 ms, corresponding to a reduction of the trap frequency. The 10 ms ramping time was sufficient to guarantee adiabaticity for our experiments because the Larmor frequency at the center of the dimple trap can be estimated as > 2 MHz by the trap bottom (~ 700 kHz corresponding to ~ 1 G). After the deformation, the deformed trap was maintained for a variable time, ranging from 0 ms to 60 ms, to strobe the free dynamical behavior in the deformed trap. The atoms were observed along the y-direction by absorption imaging after a sudden switch off of the magnetic trap and 10 ms time-of-flight. There was no trap loss during the trap maintenance time over which the delay was varied.
The condensate was excited by the deformation of the trap and two modes of collective excitation were identified from the absorption images: dipole and monopole m = 0 modes. The x and z-directional center-of-mass oscillations (dipole modes) were excited because change of the dimple current displaced the minimum point of the trap slightly due to asymmetry of the trapping Z-wire. From this center-of-mass oscillation, accurate values of trap frequencies can be obtained because dipole mode excitations are not affected by the interatomic forces, nor by the temperature or the statistics [23] . Figure 2(a) shows the shape oscillation of a condensate under the deformation with final I d = 47 mA. The condensate exhibited sinusoidal width oscillation with large amplitude along the x-axis. In addition, the z-directional width oscillated with in-phase and relatively small amplitude ( Fig. 2(b) ). Such width-oscillation behaviors represent monopole m = 0 mode (breathing mode) [9 , 10] . Note that the widths were convoluted since the condensate was observed after ballistic expansion. The striations parallel to the z-direction is the optical diffraction fringe formed by the narrowed condensate. The widths of the clouds were extracted from a Thomas-Fermi fit of the absorption profiles along the x and z-axis at the center. Beyond the optical diffraction limit, which is when the striations occurred, the widths were calibrated from the diffraction spacing. For comparison with the theoretical expectation, the x-directional oscillating widths were fitted with an exponentially decaying sine function, σ ( t ) = σ (0) + A exp(– t / τ )sin(2 π ( t φ ) / T ), and the frequency of the breathing mode and the decay rate were determined as 127.9±1.4 Hz and 28.2±5.1 s -1 , respectively. The expectation value of the breathing mode followed by Eq. (2) was ω /2 π = 126.4 Hz because the condensates were in the hydrodynamic limit Na / aHO =150≫1 and the trap frequencies of the deformed trap were determined as ωx /2 π = 80.2 ± 1.8 Hz and ω /2 π = 367.3 ± 0.5 Hz from the dipole mode oscillation. The experimentally determined frequency is consistent with the theoretically expected value within an error bound. The same deformation experiment on a thermal cloud at T / Tc ≈ 1.5 was performed. Fig. 2(c) shows the width oscillation and the frequency was 152.5±7 Hz. This frequency is also consistent with 2( ωx /2 π ) = 160.4 ± 3.5 Hz, which is expected for the non-interacting Bose gas [9] . These width oscillations in a thermal cloud can be understood by sound waves. In addition, rapid diffusion of the cloud is identified from the fit line in Fig. 2(c) , which indicates the exponentially decaying sine function for accounting for diffusion due to heat,
PPT Slide
Lager Image
. We consider the fast diffusion of the cloud is due to the finite depth of the deformed trap and collision induced thermalization of the excitation. From those results, which were consistent with the theoretical expectation, we were able to conclude that the adiabatic deformation of the trap, not by sinusoidal modulation but by ramping, induces collective excitations on the BEC.
PPT Slide
Lager Image
The width oscillations induced by the adiabatic deformation with final Id = 47.0 mA. (a) Shape oscillation of a pure BEC. After the deformation, the condensate was allowed to freely oscillate in the deformed trap for a variable time. The absorption images were taken along the y-axis after a 10 ms time-of-flight. Each absorption image has a 360 × 360 μm2 size and the trap maintenance time below. (b) The widths oscillation of pure BEC. The x and z-directional widths oscillate in the same phase with large and small amplitudes, respectively: monopole m = 0 mode. From fitting an exponentially decaying sine function, the frequency and the damping rate were determined as 127.9±1.4 Hz and 28.2±5.1 s-1, respectively. (c) The x-directional width oscillation of the thermal cloud at T/Tc≈1.5. The determined frequency is 152.5±7 Hz and the width shows rapid diffusion.
After the confirmation of the collective monopole m = 0 mode on the condensate, the experiment was done with final I d range of 0 to 94 mA. We were not able to include the data with the final I d over 94 mA because the amplitudes of the dipole mode excitations were too small to determine the trap frequencies. The frequencies of the monopole m = 0 mode for various deformed trap geometries were normalized by the corresponding x-directional trap frequency and plotted as a function of the trap anisotropy parameter 1/ λ as shown in Fig. 3(a) . The solid line is the theoretical prediction of the mode frequencies in the hydrodynamic limit as given in Eq. (2). Despite the discrepancy, the measured frequencies demonstrate a similar tendency to the anisotropy and are well below 2( ωx /2 π ) which is expected for noninteracting gases. The mode frequencies higher than the prediction for the highly anisotropic trap can be understood as the effect of anharmonic shift due to the large response amplitudes [5] . Fig. 3(b) presents the decay rate γ as a function of the breathing mode frequency. Note that reduction of the final dimple current causes not only a decrease in the mode frequency but also an increase in the response amplitudes. However, we conclude that the decay rate is dominantly affected by the mode frequency because a damping rate of the m = 0 condensate mode tends to increase with the response amplitude [5] . The damping rate dependent on mode frequency can be understood with the fact that the atom-atom collision rate is proportional to the velocity of the atoms.
PPT Slide
Lager Image
(a) Anisotropy-dependent monopole m = 0 mode excitation spectrum. The data were obtained from the adiabatic deformation experiments with the final Id range of 0 to 94 mA. Solid line indicates theoretical prediction from Eq. (2). (b) The damping rates versus the frequency of monopole m = 0 mode excitation. The damping rate quickly decreases with the mode frequency. The solid line is a guide to the eye.
- 3.2. Nonadiabatic Changing – Majorana Transitions
To investigate the dynamics under a nonadiabatic deformation, we ramped down the dimple current to 47 mA within 20 μs, ~50 T/s, which is sufficiently fast to induce Majorana spin-flips. A variable time delay was introduced between the deformation and the observation of the cloud. The absorption imaging was performed after a 2 ms time-of-flight. Figure 4(a) is the absorption image taken with a 2.8 ms time delay and shows the three spatially separated components of ultracold atoms with different spin states: m F = 2, m F = 1 and m F = 0. During the deformation, the BEC prepared initially in the m F = 2 state underwent a fast change in the magnetic field and was spin-flipped by Majorana transition [24 , 25] . After the transition, the spin-flipped m F = 0 component freely fell by gravity while the m F = 2 and the m F = 1 state components oscillated in the deformed trap with different spring constants ( Fig. 4(b) ). The oscillation frequencies of the two components have a ratio of
PPT Slide
Lager Image
which is highly consistent with the expectation
PPT Slide
Lager Image
. Note that there were no detectable m F = -2 and m F = -1 state atoms because they experienced outward forces from the deformed trap and rapidly dissipated. In addition, we investigated the population ratio of the m F = 1 atoms to the sum of the m F = 1 and the m F = 2 atoms as a function of the final frequency used in the rf evaporation ( Fig. 4(c) ). Each data point represents 30 measurements. The ratio shows an obvious increasing tendency as the final frequency decreases. For now, we do not have a clear explanation for this. However, from the fact that the population of the m F = 2 state sharply decreases with the frequency while the population of the m F = 1 state remains nearly constant, an inhomogeneous trapping field can be the most plausible reason for the tendency. Some other possible reasons need to be studied, such as temperature dependency of the transition.
PPT Slide
Lager Image
Separation of different-spin ultracold atoms under the nonadiabatic deformation. (a) The 2 ms time-of-flight absorption image with a trap maintaining time of 2.8 ms after the deformation. Different spin components can be found at different positions since the mF = 0 atoms freely fall by gravity while the mF = 1 and the mF = 2 components oscillate in the deformed trap with different frequencies. (b) The evolution of the dipole mode oscillation with time for the mF = 2 and the mF = 1 atoms. The mode frequency ratio of the mF = 2 to the mF = 1 is 1.43 . (c) Population ratio of the spin-flipped mF = 1 atoms to the sum of the mF = 1 and the mF = 2 atoms as a function of the final frequency used in the rf evaporation. The vertical dashed line marks the observed transition temperature.
In summary, we studied the dynamical behaviors of ultracold atoms on an atomchip under a trap deformation. With an adiabatic deformation, a BEC shows collective excitations, m = 0 monopole and dipole mode oscillations in the deformed trap. The width-oscillation frequency difference between the pure condensate and the non-condensed atomic cloud was clear evidence for the fact that the trap deformation was able to induce collective excitation on the BEC. The measured breathing mode frequencies in the range of the measurements were consistent with the prediction of the mean-field theory in the hydrodynamic limit. Under a nonadiabatic trap deformation that induced a Majorana transition, ultracold atoms were able to be separated spatially due to their trap frequency difference in the deformed trap.
Although collective excitation is not a new topic, the collective behaviors of the BEC on an atomchip is an interesting subject, for example, nonlinear dynamics [26 , 27] , because the trap geometry can be changed with a high aspect ratio, wide range of time, and high anharmonicity, etc. We believe this data, along with further theoretical study, will help in dynamical control of ultracold atoms on an atomchip.
We are grateful to Prof. Dana. Z. Anderson in JILA for manufacturing the atomchip and the vacuum chamber. This research was supported by a grant to the Atomic Interferometer Research Laboratory for the National Defense funded by DAPA/ADD and by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (2012R1A1A2041856).
Anderson M. H. , Ensher J. R. , Matthews M. R. , Wieman C. E. , Cornell E. A. 1995 “Observation of Bose-Einstein condensation in a dilute atomic vapor,” Science 269 198 - 201
Davis K. B. , Mewes M. O. , Andrews M. R. , van Druten N. J. , Durfee D. D. , Kum D. M. , Ketterle W. 1995 “Bose-Einstein condensation in a gas of sodium atoms,” Phys. Rev. Lett. 75 3969 - 3973
Jin D. S. , Ensher J. R. , Matthews M. R. , Wieman C. E. , Cornell E. A. 1996 “Collective excitations of a Bose-Einstein condensate in a dilute gas,” Phys. Rev. Lett. 77 420 - 423
Mewes M.-O. , Andrews M. R. , van Druten N. N. , Kurn D. M. , Durfee D. S. , Townsend C. G. , Ketterle W. 1996 “Collective excitations of a Bose-Einstein condensate in a magnetic trap,” Phys. Rev. Lett. 77 988 - 991
Jin D. S. , Matthews M. R. , Ensher J. R. , Wiemann C. E. , Cornell E. A. 1997 “Temperature-dependent damping and frequency shifts in collective excitations of a dilute Bose-Einstein condensate,” Phys. Rev. Lett. 78 764 - 767
Stamper-Kurn D. M. , Miesner H.-J. , Inouye S. , Andrews M. R. , Ketterle W. 1998 “Collisionless and hydrodynamic excitations of a Bose-Einstein condensate,” Phys. Rev. Lett. 81 500 - 503
Maragò O. M. , Hopkins S. A. , Arlt J. , Hodby E. , Hechenblaikner G. , Foot C. J. 2000 “Observation of the scissors mode and evidence for superfluidity of a trapped Bose-Einstein condensed gas,” Phys. Rev. Lett. 84 2056 - 2059
Hechenblaikner G. , Maragò O. M. , Hodby E. , Arlt J. , Hopkins S. , Foot C. J. 2000 “Observation of harmonic generation and nonlinear coupling in the collective dynamics of a Bose-Einstein condensate,” Phys. Rev. Lett. 85 692 - 695
Stringari S. 1996 “Collective excitations of a trapped Bose-condensed gas,” Phys. Rev. Lett. 77 2360 - 2363
Edwards M. , Ruprecht P. A. , Burnett K. , Dodd R. J. , Clark C. W. 1996 “Collective excitations of atomic Bose-Einstein condensates,” Phys. Rev. Lett. 77 1671 - 1674
Singh K. G. , Rokhsar D. S. 1996 “Collective excitations of a confined Bose condensate,” Phys. Rev. Lett. 77 1667 - 1670
Pérez-García V. M. , Michinel H. , Cirac J. I. , Lewenstein M. , Zoller P. 1996 “Low energy excitations of a Bose-Einstein condensate: A time-dependent variational analysis,” Phys. Rev. Lett. 77 5320 - 5323
Pérez-García V. M. , Michinel H. , Cirac J. I. , Lewenstein M. , Zoller P. 1997 “Dynamics of Bose-Einstein condensates: Variational solutions of the Gross-Pitaevskii equations,” Phys. Rev. A 56 1424 - 1432
You L. , Hoston W. , Lewenstein M. 1997 “Low-energy excitations of trapped Bose condensates,” Phys. Rev. A 55 R1581 - R158
Bogoliubov N. N. 1947 “On the theory of superfluidity,” J. Phys. (USSR) 11 23 - 32
Dalfovo F. , Minniti C. , Pitaevskii L. P. 1997 “Frequency shift and mode coupling in the nonlinear dynamics of a Bose-condensed gas,” Phys. Rev. A 56 4855 - 4863
Hechenblaikner G. , Morgan S. A. , Hodby E. , Maragò O. M. , Foot C. J. 2002 “Calculation of mode coupling for quadrupole excitations in a Bose-Einstein condensate,” Phys. Rev. A 65 033612 -
Pitaevskii L. , Stringari S. 1997 “Landau damping in dilute Bose gases,” Phys. Lett. A 235 398 - 402
Fortágh J. , Zimmermann C. 2007 “Magnetic microtraps for ultracold atoms,” Rev. Mod. Phys. 79 235 - 289
Yu H. , Kim S. J. , Moon Y. L. , Lee J. S. , Kim J. B. 2013 “Compact experimental apparatus for producing high-repetition-rate87Rb Bose-Einstein condensation on atom chip,” J. Korean Phys. Soc. 63 900 - 906
Dieckmann K. , Spreeuw R. J. C. , Weidemuller M. , Walraven J. T. M. 1998 “Two-dimensional magneto-optical trap as a source of slow atoms,” Phys. Rev. A 58 3891 - 3895
For the clarity, change of the dimple current induce trap rotation in the x-y plane which can be calculated as θxy ≈ arctan(Id/Iz) . However, there was no noticeable rotational-behavior of the condensate because Id≪Iz and the directions of the imaging are parallel to the rotation plane
Cohen-Tannoudji C. , Guéry-Odelin D. 2011 Advances in Atomic Physics World Scientific Singapore
Majorana E. 1932 “Majorana flops,” Nuovo Cimento 9 43 -
Rabi I. I. 1937 “Space quantization in a gyrating magnetic field,” Phys. Rev. 51 652 - 654
Li G.-Q. , Fu L.-B. , Xue J.-K. , Chen X.-Z. , Liu J. 2006 “Collective excitations of a Bose-Einstein condensate in an anharmonic trap,” Phys. Rev. A 74 055601 -
Ott H. , Fortágh J. , Kraft S. , Günther A. , Komma D. , Zimmermann C. 2003 “Nonlinear dynamics of a Bose-Einstein condensate in a magnetic waveguide,” Phys. Rev. Lett. 91 040402 -