Two-flavor mixture of a few fermions of different mass in a one-dimensional harmonic trap

A system of two species of fermions of different mass confined in a one-dimensional harmonic trap is studied with an exact diagonalization approach. It is shown independently on the number of particles that a mass difference between fermionic species induces a separation in the lighter flavor system. The mechanism of emerging of separated phases is explained phenomenologically and confirmed with the help of a direct inspection of the ground-state of the system. Finally, it is shown that the separation driven by a mass difference, in contrast to the separation induced by a difference of populations, is robust to the interactions with thermal environment.


Introduction
Development of cooling and trapping techniques of a few quantum particles (bosons as well as fermions) have opened completely a new area of modern quantum engineering [1][2][3][4][5]. In todays experiments, it is possible to rigorously control not only the strength of mutual interactions between atoms but also the precise number of quantum particles of different flavours confined in a trap of almost arbitrary shape [5][6][7]. These possibilities are strongly related to the fundamental theoretical questions about the properties of a mesoscopic number of particles which up to now had not been answered in a satisfactory way. One of these questions is about the transition from two-, three-body physics [8][9][10] that are quite well understood with standard approaches, to macroscopic systems, that are well described with the help of sophisticated methods of condensed matter physics. A first attempt to give an experimental answer to this question was given in [11] where the emergence of the Fermi sea with increasing number of fermions was studied in a few-body system. In the attractive regime the effective appearance of two-body pairing as a function of the system's size was observed [12] and a theoretical description of this phenomenon was given [13][14][15][16]. Recently, many different results in the area of ultra-cold few-body systems were reported both in experimental [17][18][19][20][21][22][23][24] and theoretical  works.
In a recent paper [48], a theoretical analysis of phase separation in a one-dimensional imbalanced mixture of two kinds of fermions was made. It was shown that strong interactions with a single impurity dipped in a mesoscopic Fermi sea lead to a separation of density in the bulk flavour. This separation is directly induced by the imbalance in the number of fermions of the two types. In particular it can not be observed in a balanced system. However, it is known that in balanced systems consisting of flavours of different atomic mass, separation is also possible. For example, in the recent paper [49], density separation driven by a mass difference was shown in a system confined in a uniform trap. Since in experiments with ultra-cold fermionic mixtures it is possible to confine atoms with different masses [6,7,21,23,32], therefore it is interesting to study the properties of such a system in detail. In this paper, we study the separation mechanism in a mixture of two-flavour fermionic atoms of different mass confined in a one-dimensional harmonic trap. The whole analysis is performed with an exact diagonalisation of the many-body Hamiltonian.
At this point it is worth noting that Hubbard-like models with imbalanced mass of particles were also intensively studied [50,51].

The Model
In this paper we study a two-flavour system of a few ultra-cold fermions confined in a one-dimensional harmonic trap of frequency ω and interacting via a short-range deltalike potential. The fermions of different flavours are fundamentally distinguishable and in principle may have different masses m ↑ = m ↓ . If one expresses all quantities in the natural units of the harmonic oscillator related to one of the flavours (in the following we choose σ =↓) then the Hamiltonian of the system in the second-quantized form reads: whereΨ σ (x) annihilates a fermion with spin σ at a point x and the dimensionless interaction constant g = (m ↓ / 3 ω) 1/2 g 1D is related to the effective one-dimensional interaction g 1D between fermions of opposite spins. Experimentally this effective onedimensional model is obtained by applying very strong harmonic confinement in two remaining perpendicular directions [17]. Then it can be assumed that the dynamics in these directions is frozen and particles occupy only the ground state of the trap. A direct relation between effective interaction g 1D and the s-wave scattering length a s which can be tuned with Feshbach resonance was discussed with all the details in [28]. The single-particle Hamiltonians H σ take into account the mass ratio α = m ↑ /m ↓ and have the following forms It is worth noticing that in principle both flavours of fermions can be confined in traps with different frequencies.
Here we assume that the experimental parameters are tuned such that the frequencies are equal. At the level of the single-particle Hamiltonian the only difference between flavours is related to the mass difference. This assumption has some consequences on the properties of the system. For example, the single-particle excitation energy is flavour-independent and it is equal to ω.
Since particles of different flavours are distinguishable, the field operators for opposite flavours fulfil standard commutation relations Ψ ↑ (x),Ψ † ↓ (x ) = Ψ ↑ (x),Ψ ↓ (x ) = 0. However, the field operators for the same flavour obey the canonical fermionic anti-commutation relations Ψ σ (x),Ψ † σ (x ) = δ(x − x ) and Ψ σ (x),Ψ σ (x ) = 0. Obviously, the Hamiltonian (1) commutes separately with the number of fermions of a given flavourN σ = dxΨ † σ (x)Ψ σ (x). Therefore, further analysis will always be performed under the assumption that the number of fermions in a given flavour is known and well defined experimentally. This assumption is reasonable since in ongoing experiments with few fermion systems the number of particles is controlled to a very high accuracy [5,11].

Method
The properties of the system described by the Hamiltonian (1) are studied with the exact diagonalisation approach what was previously adopted for equal mass fermionic mixtures [14,15,52], in both the repulsive and attractive case. First, we decompose the field operatorsΨ σ (x) onto the basis of the eigenfunctions of the appropriate singleparticle Hamiltonian, i.e.,Ψ σ (x) = i φ iσ (x)â iσ . Note, that due to the mass difference, the basis functions of different flavours have different spatial spread and they are related by the scaling relation φ i↑ (x) = α 1/4 φ i↓ ( √ α x). In this language, the Hamiltonian (1) has the following form wheren iσ =â † iσâ iσ . The single-particle energy E i = (i + 1/2) and the two-body interaction term U ijkl can be calculated directly Since the wave functions φ iσ (x) form a complete set in the space of single-particle states, the Hamiltonian (3) is fully equivalent to the original Hamiltonian (1). In practice, we have to cut the decomposition of the field operators at some sufficiently large level index N max . For appropriately chosen cut-off, the results do not depend on its value. For example, for the parameters N ↑ = 3, N ↓ = 1, α = 1 and g = 10 it is sufficient to take N max = 18 single-particle eigenstates in the decomposition of the field operators. Then, the dimension of the many-body Hilbert space is dim(H) = 14688. After obtaining the matrix form of the Hamiltonian (1) in the Fock basis for a given number of fermions and interaction strength g, and after subsequent diagonalisation, we  (1) for N ↑ = 3 and N ↓ = 1 and strong interaction g = 10 as a function of the mass ratio α. It is seen that for increasing α, the energy of the two many-body eigenstates belonging to the ground-manifold rapidly split and the ground-manifold has only a two-fold degeneracy.
find the M lowest many-body eigenstates |G i and corresponding eigenenergies E i of the Hamiltonian.

Spectrum of the Hamiltonian
The first significant consequence of the mass difference between flavours is visible in the spectrum of the Hamiltonian. In the standard situation, when α = 1, the Hamiltonian of the system is symmetric under exchange of the flavours. As a matter of fact, this symmetry, together with the Pauli exclusion principle, leads directly to the appearance of a degeneracy in the manifolds of many-body eigenstates of the Hamiltonian in the limit of infinite repulsion [52,53]. Each degenerate manifold is spanned by states with the same spatial distribution of particles but with different characters of the irreducible representation of the exchange-symmetry group mentioned above. For large but finite repulsions the spectrum of the Hamiltonian is not degenerate but it factorises into distant manifolds of quasi-degenerate eigenstates. The mechanism that lifts the degeneracy and its experimental consequences have been broadly discussed recently [39,52]. It is worth noticing that the case of N ↑ = N ↓ = 1 is solvable analytically for any value of the mass ratio α. The two-body problem is brought to a single-particle one by introducing centre of mass and relative distance coordinates. In these coordinates, the Hamiltonian of the system is decomposed into two single-particle problems for which analytical expressions for eigenenergies and eigenfunctions are known [54]. Many properties of two particles confined in harmonic trap were recently studied in fermionic as well as in bosonic systems [17,[55][56][57][58][59].
Whenever α = 1, the exchange-symmetry is violated since the interaction part of the many-body Hamiltonian, in contrast to the single-particle one, is not symmetric under exchange of particles. This directly leads to different spectral properties of the Hamiltonian in the limit of large repulsion. To illustrate this mechanism, in Fig. 1a we compare spectra of the Hamiltonian (1) with α = 1 and α = 40/6 (corresponding to the mass ratio in the mixture of lithium and potassium atoms [6,7,21,23,32]) calculated for a fixed number of fermions N ↑ = 3 and N ↓ = 1 as a function of interaction g. As it can be seen, for large interactions the degeneracy manifolds do split. Our numerical results suggest that this mechanism is present for different N ↑ and N ↓ . Unfortunately, due to the numerical complexity, our method enables us to study only those configurations for which N ↑ + N ↓ ≤ 10. However, for all of them the same mechanism is present if the mass ratio α is large enough. Therefore, we believe that the mechanism is very general and it is present for any number of particles. It is quite natural that for different pairs of N σ , the degeneracy may be lifted into sub-manifolds with different numbers of eigenstates. To give a better explanation of how this degeneracy is lifted we show in Fig. 1b the spectrum of the Hamiltonian (1) for g = 10 as a function of the mass ratio α. For α = 1 the ground-manifold is spanned by four many-body eigenstates. However, for a larger α the eigenenergy of the two states rapidly grows. Due to changes of the shapes of the corresponding single-particle wavefunctions φ iσ , the contact interaction between flavours becomes significant. It is seen that, the same mechanism of quasi-degeneracy splitting repeats also for higher manifolds.

Properties of the ground-state
The mass difference between fermionic species significantly changes not only the properties of the spectrum of the Hamiltonian but also some properties of many-body states. First, let us discuss changes of the many-body ground-state |G 0 . For α = 1 the properties of the ground-state were studied recently in [48]. It was shown that repulsive interactions between species lead to phase separation whenever the number of particles is imbalanced, i.e. N ↓ = N ↑ . In particular, the authors studied the interaction of a single impurity dipped in the Fermi sea of the remaining species (N ↑ > 1 and N ↓ = 1). They found that for strong interaction, the bulk flavour is pushed out of the centre of the trap and two spin domains precursors appear. The remaining impurity is located in the middle with an almost unchanged density profile. This observation can be visualized by plotting the single-particle densities defined as follows Note, that the density ρ σ is normalized to N σ . In the upper panel of Fig. 2 we show density profiles for N ↑ = 3 and N ↓ = 1 obtained with the exact diagonalisation approach for different g. In this way we reconstruct the corresponding result of [48] by using an exact diagonalisation method which is accurate near g = 0 for weak interaction. The situation is qualitatively different when a difference of mass between flavours is introduced. As an illustration, we repeat all the calculations for α = 40/6 (lower Single-particle densities ρ ↑ (x) (thin blue line, heavier flavour) and ρ ↓ (x) (thick red line, lighter flavour) calculated in the ground-state of the system for different interactions g and mass ratios α (N ↑ = 3 and N ↓ = 1). For weak interactions, the single-particle densities reproduce the density profile for ideal gasses in both species. When interaction is strong enough, a phase separation takes place in the system. For equal masses, the separation occurs in the species with the larger number of atoms (as predicted in [48]). On the contrary, for α = 40/6 the lighter species undergoes separation. We have checked that this behaviour is independent on the number of fermions in individual flavours (see main text and Fig. 3). panel in Fig. 2). Obviously, in the limit of vanishing interaction the only difference comes from the fact that the heavier impurity is better localised in the middle of the trap. However, for strong enough interactions the separation is not present in the bulk subsystem, but a significant separation of the lighter impurity appears. Note that in this case, the separation of the lighter flavour means that the impurity is in the "Schrödinger cat"-like state, i.e. the probability of finding the particle is non-zero in the two, spatially separated regions of space. It was reported recently that this sort of split state can also be seen for impurities in a bosonic bath [60]. Also a quite similar effect of the phase separation driven by the mass difference was studied recently for spin-1/2 fermions in a uniform trap [49].
The mechanism of separation induced by the mass imbalance is fundamentally different form separation driven by interactions. For example, for equal masses α = 1, the separation is present in the ground-state only for different number of particles (N ↑ = N ↓ ) and it appears always in the flavour with larger N σ . It can be viewed as the consequence of the Pauli exclusion principle and the flavour-exchange symmetry described above. Interactions favour such configurations in which different fermions do not overlap, and as a consequence particles are promoted to higher single-particle levels respecting the Pauli principle for fermions with the same spin. At the same time, the ground-state of the system has to remain in the symmetric configuration under exchange of two particles with opposite spins (this is related to the antiferromagnetic structure of the ground state, see [47]). These two requirements force the system to excite favours Position (osc. u.) with larger N σ to higher states than the particles of the minority. In consequence, the probability of finding a minority fermion in the centre of the trap decreases slower than for majority fermions. This is the reason why this kind of separation requires an imbalanced particle number. It is worth noticing that this argumentation is valid only for the ground-state of the system which in the limit of infinite interactions becomes degenerate with states of other exchange-symmetries. Those other states may not manifest separation in any of the flavours [47]. For α = 1, the separation has a different origin. In this case, the singleparticle states φ iσ are not symmetric under flavour-exchange and as a consequence, the interaction energies between the two particles occupying different levels of the harmonic oscillator are not symmetric either. This fact leads directly to the separation in the flavour with smaller mass. To explain this mechanism, let us concentrate on the case N ↑ = 3 and N ↓ = 1. In the limit of vanishing interactions, the ground-state of the system has the form |vac . When interactions are present then particles of different spins try to avoid each other to minimize the interaction energy and as a consequence they change their Fock-space configuration. The interaction energy can be minimized by exciting particles to higher single-particle states. In the model studied here, the energy cost of such an excitation does not depend on the flavour. However, due to the violation of the flavour-exchange symmetry in the interaction sector, the energy gain is larger when lighter particles (i.e. particles with a larger spread of the wavefunction) are promoted. In consequence, the Fock states of the form k↓ |vac with k > 0 start to contribute and dominate the exact many-body ground-state of the system. Direct inspection of the ground-state obtained from the exact diagonalisation method proves indeed that this simple heuristic argumentation is correct. For example, for g = 10 and α = 40/6 the projection of the true ground state |G 0 onto the non-interacting ground-state is | 210; 0|G 0 | 2 < 10 −4 . At the same time, the cumulative probability of finding all three heavy particles in the lowest states, i.e. k | 210;k|G 0 | 2 ≈ 0.9. This means that the energy of the groundstate is minimised mainly due to the mobility of the lighter particle, i.e. the particle with the larger spreading of single-particle states.
We have checked that the mechanism described above is very general and does not depend on the number of fermions of different flavours. In Fig. 3 we present the singleparticle densities for different number of fermions N ↑ and N ↓ obtained in the strong interaction regime g = 10. In all plots, the mass ratio is α = 40/6. It is seen that the separation induced by a mass difference always occurs in the lighter component (red thick line). Our results performed for all configurations with N ↑ + N ↓ ≤ 10 confirm that the lighter component is always divided into two domains regardless of the number of fermions.
It is worth noticing that for N ↑ = 1 and N ↓ = 2 our results obtained in the strong interaction limit are consistent with the density profiles obtained in [39] where a variational approach was used for this three-body problem.

Non-zero temperatures
The method used allows us to study the properties of the system also at non-zero temperatures, i.e. in a mixed state of the formρ T = Z −1 i p i |G i G i |, where p i = exp(−E i /k B T ) and Z = i p i is the partition function. A discussion of whether the one-dimensional system described by the Hamiltonian (1) reaches the thermal stateρ T when interacting with the thermal bath is beyond the scope of this article. However, our aim is to discuss the influence of the interaction with an environment on the phenomenon of phase separation. From this point of view, the temperature T can be viewed as an effective parameter which controls strength of the interaction with an external environment.
After defining the mixed stateρ T , the whole analysis can be performed in analogy to the zero temperature case provided that the single-particle densities ρ σ (x) are calculated  4. Single-particle densities ρ ↑ (x) (thin blue line, heavier flavour) and ρ ↓ (x) (thick red line, lighter flavour) calculated in the thermal stateρ T with N ↑ = 3 and N ↓ = 1 for g = 10 as a function of temperature T . For α = 1 the phase separation induced by the imbalanced number of fermions is instantly destroyed at non-zero temperature. In contrast, for α = 40/6, when the phase separation is driven by a mass difference between fermions, the single-particle density of lighter flavour is almost unchanged and it is very robust to the mixing with higher many-body states.
In practice, we cut the number of exact many-body eigenstates at the number M (defined in Sec. 3) which is temperature dependent. For a given temperature, this number is determined from the condition that the final single-particle densities do not change significantly when M is increased. It is natural that in the limit T → 0 we reproduce the previous results obtained for ρ σ (x) defined by eq. (5). In Fig. 4 we visualize the influence of non-zero temperatures on the phase separation for a number of particles equal to N ↑ = 3, N ↓ = 1 and strong interactions g = 10. As seen in the upper panel for α = 1, the phase separation is destroyed almost immediately when the higher eigenstates of the Hamiltonian contribute. Density profiles are becoming smoother and spatial density fluctuations vanish. In contrast, for α = 40/6, the separation induced by a large mass imbalance is almost insensitive to the temperature, i.e. the lighter particle remains in the Schrödinger cat-like state. The only visible consequence of higher temperature is a smoothing of the bulk flavour density. This result confirms that the mechanism of the separation induced by a mass imbalance between flavours has a completely different origin than the separation induced by interactions for a system with flavour-exchange symmetry.

Conclusions
We have studied the properties of a two-component mixture of fermions in a onedimensional harmonic trap for different mass ratios between atoms. We show that the mass difference induces a phase separation in the lighter flavour independently of the distribution of fermions between the flavours. The mechanism observed is general and is forced by the interaction energy cost when particles are excited to higher states of the harmonic oscillator. Finally, we show that in contrast to the standard separation induced by an imbalanced number of fermions, the separation driven by the mass difference is robust to non-zero temperatures. We believe that our predictions can be confirmed in upcoming experiments with ultra-cold mixtures of fermions.