APS/1-QED Population and phase dynamics of F = 1 spinor condensates in an external magnetic field arxiv:cond-mat/4869v1 [cond-mat.stat-mech] Aug 4 D. R. Romano and E. J. V. de Passos Instituto de Física, Universidade de São Paulo Caixa Postal 6618, CEP 515-97, São Paulo, São Paulo, Brasil Dated: March 19, 18) We show that the classical dynamics underlying the mean-field description of homogeneous mixtures of spinor F = 1 Bose-Einstein condensates in an external magnetic field is integrable as a conseuence of number conservation and axial symmetry in spin space. The population dynamics depends only on the uadratic term of the Zeeman energy and on the strength of the spin-dependent term of the atom-atom interaction. We determine the euilibrium populations as function of the ratio of these two uantities and the miscibility of the hyperfine components in the ground state spinors are thoroughly discussed. Outside the euilibrium, the populations are always a periodic function of time where the periodic motion can be a libration or a rotation. Our studies also indicates the absence of metastability. PACS numbers: 5..Jp,.75.Hh I. INTRODUCTION The recently realized trapping of sodium atoms by purely optical means[1] opens up the possibility of studying spinor condensates in which the spin degrees of freedom are not frozen []. Many authors have investigated, in the framework of Bogoliubov theory, the ground state configurations and the low-lying collective excitations of homogeneous mixtures of F = 1 spinor condensates in the absence [] and in the presence of an external magnetic field [, 4]. These studies predicted a variety of new phenomena such as the existence of spin domains in the ground state [, 5] and the propagation of spin waves [, 4]. The mean-field dynamics inherent in these works is known to be euivalent to a classical dynamics whose degrees of freedom are the phase and population of the three hyperfine levels [, 4]. In this paper we show that this classical dynamics is integrable as a conseuence of number conservation and axial symmetry in spin space. By a proper canonical transformation it reduces to a dynamics involving only one degree of freedom. The determination of the euilibrium points reveals a rich structure in phase-space. Contour plots of the constant energy surfaces show that the population is a periodic function of time, where the periodic motion can be a libration or a rotation. Our studies differs from references [,, 4] by taking explicitly into account the constraint of the axial symmetry in spin space. This allow us to make a detailed discussion of the properties of the euilibrium configurations and of the population and phase dynamics of the spinor condensate, which complements previous works [,, 4]. II. THE EQUATIONS OF MOTION The hamiltonian of our system of F = 1 homogeneous mixture of bosonic atoms in the presence of an external magnetic field is eual to [,, 4] Ĥ = ǫ k a α a k α k p α Ŝz α a α a k α k α, k α, k + α Ŝ z α a α a k α k with: + c V + c V α, k a α k a 1+ β a k β k a α k1 1) k1, k,α,β α ˆ S α β ˆ S β a α k a 1+ β a k β k a α k1 k1, k,α,β α,β ǫ k = k m In ), a α creates an atom in the hyperfine level α, k α = 1,, 1, with momentum k, p and are the intensities of the linear and uadratic terms of the Zeeman energy [] and c and c are, respectively, the strengths of the spin-independent and spin-dependent terms of the atom-atom interaction [,, 4]. The hamiltonian ) is number conserving and axially symmetric in spin space, ] [Ĥ, ˆN] i =, [Ĥ,e φŝz =
In the mean-field theory we suppose that the condensate is a coherent combination of atoms in the p = state, z = e 1 α zα e α zαa α ) where is the vacuum. The complex numbers z α are the condensate wave-functions for the atoms in the hyperfine level α. To find the time evolution of z α, we use the time-dependent variational principle [, 6] δs = δ which reduces to α i z z z Ĥ z dt = ) δs = δ i zα z α H z, z ) dt = 4) where the hamiltonian H z, z ) is given by H z, z ) = z Ĥ z = p z 1 z 1 ) + z 1 + z 1 ) + c z α V α + c z1 z 1 ) V + z z 1 + z 1 ) + z 1z 1z +z 1 z 1 z ) ) 5) Imposing that the action is stationary with respect to variations of z α, we get Hamilton euations of motion in complex coordinates i ż α = H zα, i żα = H 6) z α To take advantage of number conservation and axial symmetry in spin space we perform two canonical transformations. The first one introduces, as canonical variables in phase-space, the population and phase of each hyperfine component by the transformation such that the action 4) reduces to S = t t 1 z α = N α e iθα 7) α ) N α θ α H N, θ) dt 8) The second canonical transformation is given by θ 1 = θ 1 +θ +θ 1 ) N 1 = N 1 +N +N 1 θ = θ θ 1 +θ 1 ) N = N 1 N 1 +N 1 )9) θ = θ 1 θ 1 N = N 1 N 1 ) Two of these variables are, respectively, the mean number of atoms N 1 = z ˆN z 1) andone-halfthemeanvalueofthecomponentofthetotal hyperfine spin in the direction of the magnetic field N = z Ŝz z 11) In terms of these new canonical variables the action becomes eual to S = t t 1 N α θα H N 1, N, N, θ ) dt 1) α where the Hamiltonian is ciclic in the coordinates θ 1 and θ, H N 1, N, N, θ ) = p N + N 1 N + c V N 1 + c V N + 1 N 1 + N ) N 1 N + N 1 N 4 N 1 N 1 + N cos θ 1) As stated before, the dynamics follows from the condition that the action 1) is stationary, which leads to Hamilton euations of motion in these new canonical variables, θ 1 = H N 1 N1 = θ = H N N = H θ 14) θ = H N N = The propertythat the hamiltonian is ciclic in θ 1 and θ has two important conseuences. One is that the meanvalue of the total number of particles and of the component ofthe total spin ofthe condensate in the directionof the magnetic field are constants of the motion, N1 = N,
N = z S z z. Theotheristhatthe dynamicsinvolves only one degree of freedom θ = H N N = H θ 15) To simplify the euations of motion 15) we define new variables which are the old variables divided by the number of particles, n α = N α N, to write the euations of motion as 1 n = c ρ n 4 n + n sin θ θ = +c ρ n 1 ) + ) n n 1 +4 n cos θ 16) n ) 4 n where ρ = N V is the density of the condensate. The euations 16) are the analog of the mean-field classical euations of motion that describes atoms in two states coupled by a Josephson-type term [7]. From euations 16), the following general properties of the euations of motion emerges: 1) The population dynamics is independent of the strength of the linear term of the Zeeman energy and of the spin-independent component of the atom-atom interaction. Indeed, only the phase θ depends on p and only the phase θ 1 depends on c. ) By a proper choice of time scale, τ = c, the population dynamics, in the limit =, is independent of the magnitude of c, depending only on its sign. In general, that is when, the population dynamics depends on the ratio c and on the sign of c ρ. III. PROPERTIES OF THE EQUATIONS OF MOTION A. Euilibrium configurations From the euations of motion, euations 16), we see that the euations which determine the euilibrium configurations are n ) 1 4 n + n sin θ = c + n 1 ) + ρ + ) n n 1 +4 n cos θ = n ) 4 n 17a) 17b) In euations 16) and 17a,17b) the constant of the motion n is defined in the interval 1 < n < 1 and the dynamic variable n in the interval 1 < n < n. Inourdiscussionofthesolutionsoftheeuilibrium euations we consider separately the cases n and n =. a) n We have solutions which depends on the phase θ and solutions which are independent of θ a1) Solution which depends on θ with cos θ = 1. In this case the euilibrium value of n is given by euation 17b) with cos θ = 1. This euation have one solution in the interval < c < 1 + ρ 1 n ). When c the euilibrium value of n ρ approaches theupperboundary, n = n, forwhichthefraction ofatomsoccupyingthehyperfinelevelsaren 1 = n + n, n = 1 n, n 1 = n n. On the other hand when c = 1 + ρ 1 n ), n is at the lower boundary,, in which case the fraction of atoms occupying n = 1 the hyperfine levels are n 1 = 1 1+ n ), n =, n 1 = 1 1 n ). When we neglect the uadratic term of the Zeeman energy, that is =, the euilibrium value of n is, n = 1 1 n ) ) and the occupation fractions are n 1 = 1 4 1+ n ), n = 1 n), n 1 = 1 4 1 n ). a) Solution which depends on θ with cos θ = 1 The euilibrium value of n is now given by euation 17b) with cos θ = 1. This euation has one solution in the interval c > 1 ρ 1 n ). When c = ρ 1 1 n ) n is at the lower boundary n = 1. On the other hand when c ρ it approaches the upper boundary n = n. a) Solution which does not depend on the phase θ In this case n is at the lower boundary n = 1, independently of the value of c ρ. These properties are illustrated in FIG.1 where we plot the euilibrium values of n as a function of c for n ρ = 1 4. b) n = In this case the euilibrium euations, 17a,17b), reduces to
4 TABLE I: Euilibrium configurations for different parameter domains. n, antiferromagnetic case, c ρ >. < 1 c 1 n ) a1) maximum; a) minimum 1 1 n ) < < 1+ c 1 n ) a1) maximum; a) minimum; a) undefined a > 1+ c 1 n ) a) minimum; a) maximum a The classification undefined means that when we leave the line defining the corresponding boundary the energy increases or decreases depending of the value of the phase θ. TABLE II: Euilibrium configurations for different parameter domains. n, ferromagnetic case, c ρ <. > c 1 ) 1 n ) a1) minimum; a) maximum 1+ ) 1 n ) < < c 1 ) 1 n ) a1) minimum; a) maximum; a) undefined < c 1+ ) 1 n ) a) maximum; a) minimum n n = 1 1 c ρ..1 b) Solution which depends on the phase θ with cos θ = 1 replacements - -1 1 -.1 c ρ In this case there is a solution only at = and it is the line cos θ = 1, 1 < n <. When n =, there are two solutions which does not depend on the phases. -. -. b) One is the lower boundary n = 1 in which case the occupation fractions are n 1 = n 1 = 1, n =. FIG. 1: Euilibrium values of n as function of for c ρ n = 1. Solid curve is a1), dashed is a), and the lower 4 boundary n = 1 is a). The horizontal straight lines are the boundaries of n. 1 n + n sin θ = c ρ + n 1 ) 1+cos θ = ) Again, we have two classes of solutions, dependent and independent of the phase θ. b1) Solution which depends on the phase θ with cos θ = 1 Theeuilibriumvalueof n is n = 1 1 c ρ inthe interval < c ρ <. In this case the occupation fraction of the hyperfine levels are n 1 = n 1 = 1 4 1+ c ρ, b4) The other is the upper boundary n = in which case the occupation fractions are n 1 = n 1 =, n = 1. The solutions b) and b4) exist for any value of c. ρ Thepropertiesofthe n = euilibriumconfigurations are illustrated in FIG., where we plot the euilibrium values of n as a function of c. ρ We would like to point out that it can be shown that the phase dependent solutions of the euilibrium euations are the roots of the third order euation, 4 c ρ f + ) f n ) 1+ ) f c ρ c ρ c ρ + n ) 1+ c ρ ) + n ) ) = f n 18) with n = f, wheref isthefractionofatomsoutside the m f = hyperfine level, f = n 1 +n 1. We can find analytical expressions for the roots of euation 18) but they are not particularly illuminating and we will not write them here.
5 TABLE III: Euilibrium configurations for different parameter domains. n =, antiferromagnetic limit, c ρ >. c < b) minimum; b4) maximum < c < b1) maximum; b) minimum; b4) undefined = b1) maximum; b),b),b4) degenerate minimum < c c < b1) maximum; b) undefined; b4) minimum > b) maximum; b4) minimum TABLE IV: Euilibrium configurations for different parameter domains. n =, ferromagnetic limit, c ρ <. c > b) maximum; b4) minimum < c < b1) minimum; b) maximum; b4) undefined = b1) minimum; b),b),b4) degenerate maximum < c c < b1) minimum; b) undefined;b4) maximum < b) minimum; b4) maximum replacements.6.4. - -1 1 -. n c ρ FIG.: Euilibriumvaluesof n as functionof c for n =. ρ The dashed straight line is b1). The vertical straight line = is b). The horizontal straight lines are the boundaries of n. A summary of our discussion is displayed in the tables I,II,III,IV), where we consider separately the antiferromagnetic and the ferromagnetic limits. From the tables we can easily determine what are the ground state configurations for the different parameter domains and in the two limits considered above. Ground state configurations have been investigated in reference [], in the absence of the magnetic field p = = ). They observe that the degenerate) ground state in the antiferromagnetic limit are the polar states and in the ferromagnetic limit, the ferromagnetic states. The polar states of [] have n = and are the states b),b), and b4) at = see Table III). When n =, the ground state in the ferromagnetic limit is the state b1) which is eual to the ferromagnetic state of [] with n =. When n, the states a1) at = are eual to the ferromagnetic states of [] and are the ground states in the ferromagnetic limit see Table II). On the other hand, the ground state in the antiferromagnetic limit is the state a), which is not a polar state. Reference [4] investigate the ground state configurations in the presence of the magnetic field and neglecting the uadratic term of the Zeeman energy. As shown in this paper, in this case the ground state configurations coincide with the ground state configurations in the absence of the magnetic field, a fact overlooked in [4]. Besides a state with maximum value of n, n = 1, [4] identify only the spinor a), which is the ground state just in the antiferromagnetic limit see Table I). The general case when we consider both the linear and uadratic terms of the Zeeman energy have been studied in reference []. Our approach differs from [], in the sense that we take explicitly into account the constraints of axial symmetry which simplifies considerably the discussion. B. Miscibility One uestion that we can adress is the miscibility of the hyperfine components in the ground state spinors []. Consider first the case n and the antiferromagnetic limit. For < 1 1 n ), the ground state
6 Sfrag replacements Energy.7.6.5 θ 1 -. -. -.1.1 n PSfrag replacements -. Energy -.4 θ 1 -...6.4 n.6.1.4 n -.1 n. -. Sfrag replacements Energy -..5 1 1.5.5 θ PSfrag replacements Energy -..5 1 1.5.5 θ FIG. : Top: energy surface as function of θ and n. Bottom: contour plot of the energy surface in the θ n plane. In these plots n = 1 and 4 c =.4, in the antiferromagnetic limit. The energy is in units of c. Darker colors means lower energy. FIG. 4: Top: energy surface as function of θ and n. Bottom: contour plot of the energy surface in the θ n plane. In this plots n = and c =, in the ferromagnetic limit. The energy is in units of c. spinor is a), where only the m f = ±1 components are miscible, the population of the m f = component being null. For c > 1 1 n ), the ground state spinor is a), in which case the three components are generally miscible. Actually the miscibility of the hyperfine components depends on the ratio c. When c approaches its lowest value, n is near the lower boundary, n = 1, therefore the m f = component is practically immiscible with the m f = ±1 components. When c increases, a) approaches the upper boundary n = n, in which case the m f = component mixes predominantly with the m f = +1 1) component depending on the sign of n, n > < ), the population of the third component being negligible. The next case is n and the ferromagnetic limit. For c < 1+ ) 1 n ), the ground state spinor is a), where only the m f = ±1 components are miscible, the population of the m f = com- c > ponent being null. On the other hand, when 1+ ) 1 n ) the ground state spinor is a1). As in the antiferromagnetic limit, in this case the three hyperfine components are generally miscible, the degree of c miscibility depending on the ratio. When this ratio approches its lowest value, n is near the lower boundary, n = 1, and the m f = component is practically immiscible with the m f = ±1 component. On the other hand, when c increasesapproachig+, n is near the upper boundary, n = n, in which case the m f = component mixes predominantly with the m f = +1 1) component depending on the sign of n, n > < ), the population of the third component being negligible. The only case left is when n =. Consider first the antiferromagnetic limit. For c < the ground state spinor is b), for which n is at the lower boundary, n = 1. In this spinor we have eual population of the m f = ±1 components, the population of the m f = component being null. On the other hand, for c >,
7 replacements θ.15 -. n t PSfrag replacements 5 1 15 5 n θ 1 8 4 5 1 15 5 t contour curves H N =constant in the phase space plane θ n. Examples are shown in FIG. and FIG.4. We see that n t) is always a periodic function of time. The motion can be a libration when θ is a limited function of time and a rotation in which case θ always increases decreases) as a function of time. Examples of these behaviors are displayed in FIG.5. a) b) V. SUMMARY AND CONCLUSIONS replacements θ.5 -.1 n PSfrag replacementst n 5 1 15 5 c) θ 1.7 1.5 5 1 15 5 FIG. 5: a) n,b) θ as function of time, in the case of rotation; c) n, d) θ as funtion of time, in the case of libration. Notice that n is always periodic. The time is in units of c and the angles, in radians. The parameters are eual to the ones in FIG.. The initial conditions are: n =.159 and θ = π for rotation, n =.1 and θ = π for libration. the ground state spinor is b4) for which n is at the upper boundary, n =. In this spinor all atoms are in the m f = state. The ferromagnetic limit is richer than the previous one. For c < the ground state spinor is b). However, for < c <, the ground state spinor is b1), where the three components are generally miscible. Actually, it changes from a spinor where practically only the m f = ±1 are miscible, the population of the m f = component being negligible, near the lowest value of c to one where practically only the component m f = is populated near the highest value of c. For ρ c > the ground state spinor is b4). IV. DYNAMICS The ualitative features of the population dynamics can be easily visualized if we make portraits of the d) t In summary we have studied in this paper the classical dynamics that underlies the mean-field description of an homogeneous mixture of spinor F = 1 condensate in an external magnetic field. As a conseuence of number conservation and axial symmetry in spin space this dynamics is integrable. The euations of motion show that the population dynamics depends only on the uadratic term of the Zeeman energy and on the strength of the spin-dependent component of the atom-atom interaction. For a fixed mean-value of the component of the condensate spin in the direction of the magnetic field we determine the euilibrium configurations as a function of the ratio c ρ. We also make a detailed discussion of the miscibility of the three hyperfine components in the ground state spinor as a function of c ρ. Our studies revealtheabsenceofmetastabilityinthesensethatthereis no two local minima configurations in the parameter domain. We have shown that outside the euilibrium, the populations are always a periodic function of time, where the periodic motion can be a libration or a rotation. In the first case the phase θ is always limited whereas in the second case it always increase decrease), as function of time. Finally we would like to remark that the restriction to a homogeneous mixture is not only of academic interest [, 4]. Besides being a guide to what happens in the case of trapped condensates, classical Hamiltonians of the type considered in this paper emerges in a mean-field description of the uantum single-mode approximation for spinor condensates in a trap [8]. DRR would like to acknowledge financial support from CNP, and EJVP partial support from CNP. [1] D. M. Stamper-Kurn, M. R. Andrews, A. P. Chikkatur, S. Inouye, H. -J. Miesner, J. Stenger, and W. Ketterle, Phys. Rev. Lett. 8, 7 1998) [] J. Stenger, S. Inouye, D. M. Stamper-Kurn, H. -J. Miesner, A. P. Chikkatur, and W. Ketterle, Nature 96, 45 1998) [] T. -L. Ho, Phys. Rev. Lett. 81, 74 1998) [4] T. Ohmi, K. Machida, J. Phys. Soc. Jpn. 67, 18 1998) [5] T. Isoshima, K. Machida, and T. Ohmi, Phys. Rev. A 6, 4857 1999) [6] F. Dalfovo, S. Giorgini, L. P. Pitaevskii, S. Stringari, Rev. Mod. Phys., 71, 46 1996)
8 [7] A. Smerzi, S. Fantoni, S. Giovanezzi, and S. R. Shenoy, Phys. Rev. Lett. 79, 495 1997). S. Ragavhan, A. Smerzi, S. Fantoni, ands. R.Shenoy, Phys. Rev. A 59, 6 1999) [8] C. K. Law, H. Pu, and N. P. Bigelow, Phys. Rev. Lett., 81, 557 1998)