Abstract
We study the thermodynamic properties of a spin-1 Bose gas across the Bose–Einstein condensation transition. We present the theoretical description of the thermodynamics of a trapped ideal spin-1 Bose gas and we describe the phases that can be obtained in this system as a function of the temperature and of the populations in the different spin components. We propose a simple way to realize a 'synthetic magnetization' that can be used to probe the entire phase diagram while keeping the real magnetization of the system fixed. We experimentally demonstrate the use of such method to explore different phases in a sample with zero total magnetization. Our work opens up new perspectives to study isothermal quenching dynamics through different magnetic phases in spinor condensates.
Export citation and abstract BibTeX RIS
Original content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.
1. Introduction
Spinor Bose gases and spinor Bose–Einstein condensates (BECs) are characterized by the fact that their constituent particles have an internal degree of freedom: their spin. For example, in alkali atoms if the total spin of the atoms is F and mF denotes its magnetic quantum number, the different Zeeman states of one hyperfine manifold coexist in such systems. Concerning spinor BECs, the combination of magnetic ordering and superfluidity makes them interesting systems to study phenomena involving spontaneous symmetry breaking [1, 2], spin superfluidity [3], vortex dynamics [4], or collective magnetic excitations [5]. Of particular interest are the understanding of spin dynamics and the characterization of the ground states properties of these systems, which are determined by collisional processes [6]. Collisions between the different internal states of the atoms allow spin-changing collisions that have been studied in detail in [7–10]. These collisions can notably be employed to generate spin squeezing [11–13], that can be used to overcome the quantum shot noise limit [14–18]. Spinor dynamics was also studied in two dimensional systems [15], and in the presence of periodic potentials and across the superfluid to Mott insulator transition [19].
Here, we focus on spin-1 bose gases, and in particular on alkali atoms in the hyperfine F = 1 state, where the three magnetic Zeeman substates mF = 1, 0, −1 coexist. Spin-1 gases can display ferromagnetic or antiferromagnetic character depending on the sign of the spin-dependent contact interaction term , where aF are the s-wave scattering lengths for the two allowed spin collisional channels F = 0 and F = 2, and m is the atomic mass [20]. The rotational symmetry of s-wave collisional processes (provided that dipolar interactions are negligible) implies that the total magnetization of the system, defined as with N±1 the populations in the mF = ±1 Zeeman substates, is a conserved quantity. Extensive work has been done to study the phases and mean-field ground states of both ferromagnetic [20–24] and antiferromagnetic [25–27] spin-1 condensates.
While the ground state properties of spinor Bose gases have attracted substantial interest, their finite temperature behavior has not been investigated thoroughly. The additional internal degree of freedom makes these systems richer than single component Bose gases, and a large number of different thermodynamic phases can be observed. To date, only the thermodynamics of antiferromagnetic systems has been studied, where step-wise condensation of the spin components was observed as a function of the initial magnetization and the external magnetic field [28].
In this work, we study the thermodynamics of a spin-1 Bose gas using the relevant non-interacting theory for an ideal trapped gas, extending the work presented in [29]. We classify the different magnetic phases of this system, and demonstrate that it is possible to induce a 'synthetic magnetization' by exploiting a spin-dependent trapping potential. We use such method to experimentally explore the phase diagram of a symmetric polar (SP) Bose gas of 87Rb, which is characterized by zero total magnetization and equal population of the three spin states . We show that our method can be used to realize highly magnetized condensates while keeping the total magnetization of the system to zero.
Due to their rich phase diagram, spinor Bose gases have recently been promoted as an optimal system to study non-equilibrium dynamics. For example thermal quenches were used to cross over the BEC transition [30, 31], and microwave dressing allowed to operate selectively on the Zeemen energy levels [32]. Our work opens new possibilities in performing isothermal quenches across different phases using the synthetic magnetization, a method that could be exploited in future experiments to study out-of-equilibrium physics in spinor systems.
This paper is organized as follows: in section 2 we present the non-interacting model for a spin-1 Bose gas, highlighting the presence of up to three critical temperatures and classifying the corresponding magnetic phases. In section 3, we describe our method to experimentally generate the synthetic magnetization, we present the details of our experiment and the results on the experimental exploration of the phase diagram of the 87Rb SP Bose gas. Finally, section 4 is devoted to the conclusions.
2. Theory of ideal spin-1 Bose gases
In this section, we present the theory describing the condensation dynamics of an F = 1 ideal Bose gas within the grand canonical ensemble formalism. We extend the theory presented in [29] to the more general situation when not only the total atom number and the magnetization is fixed but also the number of atoms in the mF = 0 Zeeman component. We will give analytical expressions for the critical temperatures and we will classify the different magnetic phases that can be realized with this kind of systems.
Let us consider an ideal, trapped, dilute spin-1 Bose gas in the presence of a magnetic field. In the case of alkali atoms, the effect of a non-zero magnetic field B along the direction, which sets the quantization axis, can be expressed analytically through the Breit–Rabi formula [23]. The contribution to the total energy of the system can be decomposed into linear and quadratic parts , with , and /G2 for 87Rb atoms, where gJ, gI are the gyromagnetic ratios of the electron and nucleus, Ehfs is the hyperfine energy splitting for zero magnetic field, μB is the Bohr magneton and we skipped constant terms. The linear contribution is irrelevant as it is proportional to the magnetization M which is a constant of motion. The quadratic part is of the main importance in the lowest order approximation, even for a realistic system composed of interacting atoms. The presence of the spin-mixing collisional processes makes the linear part of this effect irrelevant for the dynamics. In other words, the chemical equilibrium required by the spin-changing collisional processes implies that the effective chemical potentials of the individual species in the condensate are constrained by the relation μ1 + μ−1 = 2μ0. The consequence of this is that condensing at a fixed magnetization has the same effect as condensing under the effect of an effective external magnetic field. Therefore, the applied magnetic field can be viewed as an effective magnetization of the sample.
Given these preliminary considerations, the Hamiltonian of a trapped spin-1 ideal Bose gas can be written as
where , and . In the above Hamiltonian, the chemical potential μ, the linear Zeeman shift η, and γ are Lagrange multipliers that enforce all the constrains in our system i.e the conservation of the total number of atoms , the magnetization , and the population of the N0 state . The Zeeman energy is included in the Lagrange multipliers. In other words, the Lagrange multipliers are shifted by the non-zero magnetic field, i.e. and . The energy spectra for the three spin components are therefore
Assuming equal trapping energies for all three Zeeman components, it is clear that a state with the lowest energy will be determined by the interplay between the external magnetic field parametrized by γ and the effect of fixed magnetization parametrized by η. We will distinguish two limits for non-negative magnetization: (a) γ > η, when E0 is the lowest energy state; (b) γ < η, when the lowest energy state is E1. Using the above Hamiltonian, we can write the grand canonical partition function of the system as
and with kB being the Boltzmann constant and T the temperature. The particular fugacities are . According to the rules of grand canonical formalism, we can derive the ensemble average population
in which the effective fugacities are . In the thermodynamic limit, following the standard derivation we can write the condensate fraction in the mF Zeeman component as
while the number of thermal atoms in each component can be expressed
where and g3(x) is the Bose function.
2.1. In the large magnetic field limit when γ > η
First, we will focus on the limit of γ > η, i.e. when the energy associated to the quadratic Zeeman effect dominates over the mean-field energy associated to the fixed magnetization. In this regime it is the mF = 0 component that condenses first, leading mathematically to μ → −γ, z0→1 and , . The number of thermal atoms in each Zeeman state reads
where ζ(3) is the Euler–Riemann zeta function and the number of condensed atoms are N±1,0 → 0 and . Following the arguments in [29], the first critical temperature T0,c for the mF = 0 state is defined by which gives
because the number of atoms in the mF = 0 component is fixed.
The second phase transition occurs when η → γ. At this point it is the mF = 1 component which condenses, leading to the following relations for the number of thermal atoms:
and the number of condensed atoms N0,0 ≫ 1, N−1,0 ≫ 1 and . The second critical temperature T1,c is defined as
The third phase transition occurs when γ → 0. Assuming that γ = q + γz as above, the transition is possible when γz → −q and then one can define the third critical temperature T−1,c for the mF = −1 as
Notice that when the value of magnetization is zero then T1,c = T−1,c and both mF = 1 and mF = −1 components condense at the same temperature. Moreover, in the symmetric case for N0 = N/3 the three components condense simultaneously. In addition, the second and third phase transition can be defined also when η → −γ and then the role of mF = 1 and mF = −1 components exchange. This case corresponds to negative values of magnetization.
2.2. In the low magnetic field limit when
In the case when the effect of fixed magnetization dominates it is the mF = 1 component that condenses first, leading to μ → −η and the mathematical relations among fugacities z1 → 1, , , the number of thermal atoms
the condensate fractions and the sharp grow of the N1,0(T) value above the critical point. The first critical temperature T1,c for the mF = 1 state can be defined as
The magnetization of condensed atoms is zero at T = T1,c, but it starts grow up above T1,c. The second phase transition occurs when η → γ and N0,0 ≫ 1. The second critical temperature in this situation, T0,c can be defined from the constraint of the number of atoms in the mF = 0 component
At T0,c the magnetization of condensed atoms is already non-zero but still some thermal atoms contribute in order to take into account its fixed value
The third phase transition takes place when γ → 0 as we consider the limit q → 0. Now, the mF = −1 starts to condense. The third critical temperature T−1,c can be defined as
What is more interesting, one can show that at T−1,c the value of magnetization is determined by condensed atoms only as the contribution of thermal atoms compensate each other :
An evidence of the third transition is the relative magnetization of condensed atoms equal to one, which is a characteristic feature of the low magnetic field region. The behavior of the three normalized critical temperatures , with as a function of N0/N and M/N is reported in figure 1.
2.3. Magnetic BEC phases
From the theoretical model just described, a F = 1 spinor Bose gas with fixed magnetization features one, two, or three critical temperatures depending on the balance between the Zeeman populations of the sample. The presence of three different critical temperatures gives rise to a number of different phases of the condensed part of the spinor gas. In particular, for we have that only the mF = +1 component is condensed, therefore we label this phase as magnetized BEC (I). When instead both the mF = 1 and mF = −1 components are condensed, and we label this phase as axisymmetric BEC (II). The case in which corresponds to a situation when the mF = 0 and mF = 1 components are condensed,therefore this phase corresponds to transverse magnetized BEC (III). When only the mF = 0 state is condensed and we have a polar BEC (IV), while when all the spin components are condensed and therefore this phase corresponds to a spinor BEC (V). In figure 2 we report the different phases as a function of the total magnetization of the system for the three different cases of γ > η, γ = η and γ < η. In case of zero magnetization (η = 0) we have that . The SP state (η = γ = 0, indicated as a dot in figure 2) is the only one for which the three critical temperatures coincide and there is a direct transition from normal gas to spinor BEC.
Download figure:
Standard image High-resolution image3. Exploring the phase diagram
Experimentally, it is possible to access the different thermodynamic phases of a spin-1 Bose gas by adjusting the populations in the three states and changing the temperature of the sample, as it was done in [28]. A different method, that we employ in this work, is to use a spin-selective trapping potential that induces a 'synthetic magnetization'. This allows us to control the thermodynamic properties of the system and explore the phase diagram without the need of changing the populations in the three spin components. This method opens up the possibility of performing isothermal quenches across the different thermodynamic phases, with the additional benefit that the process can be reversible.
Details about our experimental sequence and methods can be found in [33]. In brief, we load 87Rb atoms from a 3D MOT into a bichromatic crossed dipole trap made by two lasers with wavelengths λ1 = 1064 nm and λ2 = 1550 nm. At the beginning of the evaporation, the atoms are evenly distributed among the three Zeeman states of the F = 1 ground state. This means that the total magnetization M is zero (within our experimental error bars). Since no coherences are present in the system, this state is the SP = (N/3, N/3, N/3) state. We then start the evaporative cooling process, that equally removes atoms from the three Zeeman states, thus preserving the magnetization and the symmetry of the state at every temperature. In other words, in our experiment we preserve the SP state at every temperature. We stop the evaporation at different times, corresponding to different final temperatures in the range 250−0 nK (with 0 nK we indicate a BEC with no measurable thermal component). The trapping frequencies at the end of the evaporation, i.e. when we have a pure BEC, are 2π × (284; 284; 60) Hz. For each temperature, we let the system thermalize and equilibrate for 5 s, longer than reported in previous experimental works [24, 28]. This time is needed to ensure that any spin dynamics in the system has evolved towards its equilibrium state. To detect each component separately, we switch off the trapping potential and we let the atomic cloud fall during a time of flight of 30 ms. During this time, we apply a magnetic field gradient that spatially separates the three Zeeman substates. This allows to image the three clouds independently and to fit each of them using independent routines, extracting the temperature, the number of atoms and the condensate fraction.
3.1. Synthetic magnetization
As it can be observed in figure 2, the SP state would normally feature a single critical temperature and the sample would undergone a one-step transition from normal gas to spinor condensate (V). However, for trapped samples, it is possible to induce a 'synthetic magnetization' also for the SP state by selectively acting on the external trapping potential of the three spin components.
Let us consider the Hamiltonian for an SP state (γ = η = 0) with spin-selective trapping potentials:
A way to achieve such configuration is to use a dipole trap with elliptical polarization, indeed for alkali atoms, in case of large detunings and negligible saturation, the dipole potential is given by [34]
with c the speed of light, Γ the atom decay rate, I(r) the intensity profile of the laser, gF the Landé factors, and the frequencies of the D1 and D2 atomic transition, Δ1,F and Δ2, F the corresponding detunings of the laser light and P the laser polarization (P = 0 for linear polarizations and P = ±1 for σ± polarizations). Clearly, in case of σ polarized light the same laser beam produces a different potential for the three spin states, while a linearly polarized or unpolarized light produces a potential that is not state dependent. In case the polarization of the light is either σ+ or σ−, the corresponding energy spectra are
where it is apparent that the difference in the three trapping potentials acts as if there was a synthetic magnetization in the system. Indeed for large detunings the shift of the trapping potential for the mF = 1 state with respect to the potential for the mF = 0 state is with very good approximation opposite to the shift for the mF = −1 state. Therefore, by controlling the polarization of the light from circular to linear—or not defined—we can control the synthetic magnetization and access different regions of the phase diagram. In our experiment, we achieve this by rotating the quantization axis of the system.
In the absence of any compensation field, in our setup there is a small horizontal magnetic field G. Our trapping lasers propagate also in the horizontal plane and the one at 1064 nm have an excess of σ+ polarization of ≃15%. At the beginning of the evaporation, we adiabatically ramp the current in a pair of vertical coils arranged in Helmholtz configuration, that sets the magnitude of the vertical magnetic field Bz. Therefore, by increasing the current in the Helmholtz coils we rotate the direction of the magnetic field. Accordingly, the quantization axis of the system rotates from horizontal to vertical (see the inset in figure 3), and the polarization of the light changes from elliptical to undefined. In other words, P in equation (19) goes from a finite value to zero and therefore U goes from being spin-dependent to be spin-independent. In figure 3, we plot the difference as a function of the applied magnetic field Bz, calculated using equation (19) with our experimental parameters. The difference in trap depth goes from a maximum of ≃25 nK when the magnetic field is horizontal, to zero when the quantization axis is almost vertical. In our experiment, the dipole potential becomes state independent when Bz ≥ 3 G.
Download figure:
Standard image High-resolution image3.2. Experimental results
In figure 4 we show the measured condensed fractions of the three spin components as a function of the temperature of the system and the applied magnetic field Bz. At 210 nK, which is the highest temperature shown, we have ≃5 × 104 atoms in each spin component. As we proceed with the forced evaporation towards lower temperatures, the number of atoms progressively decreases. At 50 nK we have ≃1.5 − 2 × 104 atoms in each spin component. We observe that for low magnetic fields the three components condense at different temperatures, while for higher values of Bz the three critical temperatures coincide. At ≃3 G (white arrow) we observe an anomalously low number of atoms in the mF = 1 component. The origin of such feature is not completely clear, however our trapping laser at 1064 nm is intrinsically modulated in amplitude at some specific frequencies. The spectrum contains a peak at ≃20 MHz that at this value of the magnetic field could induce two-photon transitions to an s-wave bound state located 24.37 MHz from the atomic threshold [35]. The process is similar to the one studied in [36] and its detailed study will be the subject of future works.
Download figure:
Standard image High-resolution imageFor low values of Bz, where we observe three different critical temperatures, we have the highest difference between the three trapping potentials and therefore the highest value of the synthetic magnetization. As the temperature is decreased, the first atoms to condense are those in the mF = 1 state, therefore realizing the magnetized BEC phase (I). In this phase the magnetization of the condensate fraction is indeed . Further decreasing the temperature, the atoms in the mF = 0 condense, realizing a transverse magnetized BEC (III), where 1 > Mc > 0. At lower temperatures, also the atoms in the mF = −1 condense and we enter the spinor condensate phase (V). The difference between T1,c and T0,c is larger than the difference between T0,c and T−1,c even if the difference in the trapping potentials is the same. This is due to the mean-field potential exerted by the atoms already condensed that 'flattens' the potential for the non-condensed atoms, an effect not included in the non-interacting theory.
As we increase the field, we rotate the quantization axis of the system and we therefore reduce the synthetic magnetization. We observe that the critical temperature T1,c progressively decreases until T1,c = T0,c = T−1,c at ≃3 G. This corresponds to the situation in which the quantization axis of the system is almost completely vertical, and the trapping potential becomes effectively spin-independent. In this regime we observe the direct transition from the normal gas to spinor BEC.
It is important to remark that for every value of Bz the total real magnetization of our spinor gas is always zero, as shown in figure 5(a). Within our error bars, the system as a whole, i.e. accounting for the condensed and non-condensed parts of the system, remains unmagnetized and the SP state is conserved. As discussed, the use of the spin-selective potential allows us to generate a synthetic magnetization that manifests itself in the onset of different magnetic phases for the condensed part of the sample, as shown in figure 5(b). By changing the direction of the external magnetic field from horizontal to vertical and by controlling the temperature of the sample, we are able to explore the whole phase diagram reported in the central panel of figure 2. Increasing Bz corresponds indeed to decreasing the (synthetic) magnetization and therefore to moving from right to left in figure 2(b).
Download figure:
Standard image High-resolution image4. Conclusions and outlook
In this work, we presented the non-interacting thermodynamic theory of a trapped spin-1 Bose gas. We classified the different magnetic phases of the condensed part of the system, and we derived analytic expressions for the critical temperatures of the three Zeeman substates . We proposed a method to induce a synthetic magnetization in these systems using a spin-dependent trapping potential, and we presented a simple way to control such potential combining static magnetic fields and optical traps made with elliptically polarized light. In addition, we demonstrated experimentally that controlling the synthetic magnetization we were able to explore the phase diagram of a spin-1 Bose gas using an atomic sample with total zero magnetization. The extension of the technique presented in this work opens new exciting possibilities to study out-of-equilibrium physics in ferromagnetic spinor gases after a sudden (isothermal) quench of the synthetic magnetization, a situation that remains unexplored in the burgeoning field of quenched spinor BECs, that is attracting increasing interest both theoretically [37, 38] and experimentally [39].
Acknowledgments
We are very grateful to Jeremy Hutson for performing the bound state calculations for our system. We thank J Sánchez Claros for his help during the measurement campaign and the people of the Cold Atoms group at the University of Birmingham for useful discussions and suggestions.