# Exact density profiles and symmetry classification for strongly interacting multi-component Fermi gases in tight waveguides

Jean Decamp, Pacome Armagnat, Bess Fang, Mathias Albert, Anna Minguzzi, Patrizia Vignolo

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Jean Decamp, Pacome Armagnat, Bess Fang, Mathias Albert, Anna Minguzzi, et al.. Exact density profiles and symmetry classification for strongly interacting multi-component Fermi gases in tight waveguides. New Journal of Physics, 2016, 18, pp.055011. 10.1088/1367-2630/18/5/055011. hal-01339365

HAL Id: hal-01339365

https://hal.sorbonne-universite.fr/hal-01339365

Submitted on 29 Jun 2016

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Exact density profiles and symmetry classification for strongly interacting multi-component Fermi gases in tight waveguides

This content has been downloaded from IOPscience. Please scroll down to see the full text.

2016 New J. Phys. 18 055011

(http://iopscience.iop.org/1367-2630/18/5/055011)

View the table of contents for this issue, or go to the journal homepage for more

Download details:

IP Address: 134.157.80.213

This content was downloaded on 29/06/2016 at 15:14

Please note that *terms and conditions* apply.
<figure><img src="image_3.png" /></figure>

CrossMark

OPEN ACCESS

RECEIVED

29 February 2016

REVISED

5 April 2016

ACCEPTED FOR PUBLICATION

3 May 2016

PUBLISHED

18 May 2016

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.

<figure><img src="image_4.png" /></figure>

PAPER

Exact density profiles and symmetry classification for strongly
interacting multi-component Fermi gases in tight waveguides

Jean Decamp¹, Pacome Armagnat²,³, Bess Fang³, Mathias Albert¹, Anna Minguzzi⁴,⁵ and Patrizia Vignolo¹,⁷

¹ Université de Nice-Sophia Antipolis, Institut Non Linéaire de Nice, CNRS, 1361 route des Lucioles, F-06560 Valbonne, France

² Université Grenoble-Alpes, INAC-PHELIQS, F-38000 Grenoble, France

³ CEA, INAC-PHELIQS, F-38000 Grenoble, France

⁴ LNE-SYRTE, Observatoire de Paris, PSL Research University, CNRS, Sorbonne Universités, UPMC Univ. Paris 06, 61 avenue de l'Observatoire, F-75014 Paris, France

⁵ Université Grenoble-Alpes, LPMMC, BP166, F-38042 Grenoble, France

⁶ CNRS, LPMMC, BP166, F-38042 Grenoble, France

⁷ Author to whom any correspondence should be addressed.

E-mail: patrizia.vignolo@inln.cnrs.fr

Keywords: one-dimensional systems, strongly interacting Fermi gases, Lieb–Mattis theorem

Abstract

We consider a mixture of one-dimensional strongly interacting Fermi gases with up to six
components, subjected to a longitudinal harmonic confinement. In the limit of infinitely strong
repulsions we provide an exact solution which generalizes the one for the two-component mixture.
We show that an imbalanced mixture under harmonic confinement displays partial spatial separation
among the components, with a structure which depends on the relative population of the various
components. Furthermore, we provide a symmetry characterization of the ground and excited states
of the mixture introducing and evaluating a suitable operator, namely the conjugacy class sum. We
show that, even under external confinement, the gas has a definite symmetry which corresponds to the
most symmetric one compatible with the imbalance among the components. This generalizes the
predictions of the Lieb–Mattis theorem for a Fermionic mixture with more than two components.

Introduction

Ultracold atomic gases made with rare-Earth elements cooled to quantum degeneracy and subjected to two-
dimensional optical lattices provide a beautiful realization of the model of one-dimensional multicomponent
Fermi gases with strong and equal intercomponent repulsion among the species [1].

In the absence of external harmonic confinement in the longitudinal direction, the system of multi-
component Fermions with intercomponent delta interactions is a generalization of the Yang–Gaudin
Hamiltonian [2–4] and can be solved by nested Bethe Ansatz [5]. The homogeneous system may also be
described at low energy as a multi-component Luttinger liquid. This model has been extensively studied in the
context of electronic multichannel systems [6, 7] and exotic condensed matter materials [7–11].

The experiments with ultracold atoms are characterized by the presence of an external longitudinal
confinement, which can be well approximated as harmonic. In this case, the Bethe Ansatz solutions do not apply,
however, a special integrable case is provided by the limit of infinitely strong repulsions among the species. In
this regime, corresponding to the Tonks-Girardeau regime for ultracold bosonic atoms, Fermions belonging to
different components further Fermionize, as has been experimentally shown in [12]: the wavefunction vanishes
at contact and they can be mapped onto a noninteracting Fermi gas in the same external confinement with
particle number corresponding to the total number of Fermions in the mixture, following the same idea as the
original Girardeau solution for bosons [13]. An additional difficulty in the multicomponent case is that the
manifold displays a large degeneracy [14]. This is associated to the arbitrariness in fixing the relative phase once
two Fermions belonging to different components get in contact and then exchange their mutual position. The
degeneracy is however broken at finite interactions, where one expects a unique ground state [15, 16]. The
ground state branch can be obtained by performing a strong-coupling 1/g expansion, g being the interaction
strength among the Fermions. This provides a unique way to identify in the degenerate manifold what would be the ground and excited states at finite interactions [16, 17].

In this work, we explore this type of solution for multicomponent Fermi mixtures with a number of components that ranges from two to six as in the ¹⁷³Yb experiment of [1]. We study in particular the density profiles of the imbalanced mixtures, thus generalizing the works of [18, 19]. As main result, we find a complex, inhomogeneous spatial structure. We then explore the symmetry property of both the ground and excited state branches of the solution, introducing and evaluating suitable operators for the mixture, i.e. the transposition class-sum and the three-cycle class sum operators [20, 21]. This allows us to test (and generalize) the Lieb–Mattis theorem [22] for the trapped multicomponent mixtures, showing that the ground state carries the most symmetric configuration allowed by the imbalance among the components.

## 1. Model

We consider a system of N Fermions of equal mass m, divided in r species with population $N_1, N_2, ..., N_r$. We assume that all components are subjected to the same harmonic potential $V(x) = m\omega^2x^2/2$, as is the case of Fermions in optical traps. The Fermions belonging to different species interact with each other via the contact potential $v(x - x') = g\delta(x - x')$, where g is the interaction strength, and $\delta(x)$ is the Dirac delta function. The total Hamiltonian reads

$$ \mathcal{H} = \sum_{i=1}^{N} \left[ -\frac{\hbar^2}{2m} \frac{\partial^2}{\partial x_i^2} + \frac{1}{2} m \omega^2 x_i^2 \right] + g \sum_{i<j} \delta(x_i - x_j), \quad (1) $$

where $x_1, ..., x_N$ are the coordinates of the Fermions. The effect of contact interactions can be replaced by a cusp condition on the many-body wavefunction:

$$ \frac{\hbar^2}{2mg} \left[ \left( \frac{\partial \Psi}{\partial x_i} - \frac{\partial \Psi}{\partial x_j} \right)_{x_i - x_j \to 0^+} - \left( \frac{\partial \Psi}{\partial x_i} - \frac{\partial \Psi}{\partial x_j} \right)_{x_i - x_j \to 0^-} \right] = \Psi(x_i = x_j). \quad (2) $$

In the following, we will focus on the impenetrable limit $g \to \infty$. In this case the cusp condition imposes that the many-body wavefunction vanishes when two particles meet, i.e. $\Psi(x_j = x_\ell) = 0$ for each pair $\{j, \ell\}$. This condition is exactly satisfied by the fully antisymmetric solution $\Psi_A(x_1, ..., x_N)$ of $N = N_1 + N_2 + ... + N_r$ noninteracting Fermions in the same confining potential, corresponding to a Slater determinant constructed with single-particle wavefunctions $\phi_0, ..., \phi_{N_1-1}$. To construct the exact solution, we further note that the behavior of the many-body wavefunction under exchange between two Fermions belonging to different components is not fixed by symmetry, and requires to be fixed by additional conditions. Hence, we consider a general solution of the form [17, 23]

$$ \Psi(x_1, \dots, x_N) = \sum_{P \in S_N} a_P \chi(x_{P(1)} < \dots < x_{P(N)}) \Psi_A(x_1, \dots, x_N), \quad (3) $$

where $\chi(x_1 < \dots < x_N)$ is the indicator function of the sector $\{x_1 < \dots < x_N\} \subset \mathbb{R}^N$, i.e. it is 1 within the sector and 0 everywhere else and $S_N$ the permutation group of N elements. Here, both $\Psi$ and $\Psi_A$ are assumed to have unit normalization, and the choice of the coefficients $a_P$ will be detailed below. Note that we need only to determine $S = \frac{N!}{N_1! \dots N_r!}$ coefficients $a_P$: the wavefunction is antisymmetric under exchange of Fermions belonging to the same component, and this property is already encoded in $\Psi_A$. This observation allows us to restrict ourself to the so-called snippet basis [23, 24], i.e. consider only the global permutations modulo the permutations of particles belonging to the same species.

### General solution

In order to determine the $a_P$ coefficients for the ground state manifold, denoted $a_1, ..., a_S$ in the snippet basis, we use the same method as in [17]. It consists in a perturbative expansion of the energy to first order in $1/g \to 0$, i.e. we write $E = E_A - (\hbar^4/m^2)K/g$, where $E_A$ is the energy associated with $\Psi_A$ and $K = -(m^2/\hbar^4)(\partial E/\partial g^{-1})$ is proportional to the interaction energy, related to the Tan's contact coefficient in the two-component case [25]. We then use the Hellmann–Feynman theorem and the cusp condition to write

$$ K = \frac{1}{4 \langle \Psi | \Psi \rangle} \sum_{i<j} \int dx_1 \dots dx_N \delta(x_i - x_j) \left[ \left( \frac{\partial \Psi}{\partial x_i} - \frac{\partial \Psi}{\partial x_j} \right)_{x_i - x_j \to 0^+} - \left( \frac{\partial \Psi}{\partial x_i} - \frac{\partial \Psi}{\partial x_j} \right)_{x_i - x_j \to 0^-} \right]^2, \quad (4) $$

where we have used the natural units, namely the harmonic oscillator length $a_{ho} = \sqrt{\hbar/m\omega}$ as unit length and the harmonic oscillator energy $\hbar\omega$ as unit energy. Separating the integral over the different sectors $\{x_{P(1)} < \dots < x_{P(n)}\} \subset \mathbb{R}^N$ and recalling that $\Psi_A$ is normalized to one, we finally obtain
$$K = \frac{\sum_{P,Q \in S_N} (a_p - a_Q)^2 \alpha_{P,Q}}{\frac{1}{N!} \sum_{P \in S_N} a_p^2}, \qquad (5)$$

where the $\alpha_{P,Q}$ coefficients are non-zero if the sectors P and Q differ only by transposing two adjacent coordinates, and in this case we have, using permutational symmetry

$$\alpha_{P,Q} = \alpha_k = \int_{x_1 < \dots < x_N} dx_1 \dots dx_N \delta(x_k - x_{k+1}) \left[ \frac{\partial \Psi_A}{\partial x_k} \right]^2, \qquad (6)$$

for $k \in \{1, \dots, N-1\}$. Intuitively, $\alpha_k$ can be seen as the energy cost of an exchange between particles of different species at positions k and k + 1. Note that, thanks to the parity invariance of this problem, we also have $\alpha_k = \alpha_{N-k}$, so that we have only $[N/2]$ coefficients to compute.

In order to find the wavefunction which corresponds to the ground state at finite, large interactions, the next step is to find the solutions that minimize the energy, i.e. that maximize K. To do so, we impose that $(\partial K / \partial a_i) = 0$ for all $a_i$. This turns out to be equivalent to the diagonalization problem

$$V_A = K_A, \qquad (7)$$

with $\vec{A} = (a_1, \dots, a_S)^T$ and V is a $S \times S$ matrix depending only on the $\alpha_k$ coefficients. More precisely, the V matrix is defined in the snippet basis by

$$V_{ij} = \begin{cases} -\alpha_{i,j} & \text{if } i \neq j, \\ \sum_{k \neq i} \alpha_{i,k} & \text{if } i = j, \end{cases} \qquad (8)$$

where the $\alpha$ coefficients are defined as in equation (6) (see also [26, 27]).

In order to compute the $\alpha_k$ coefficients, we use the following expression for $\Psi_A$, based on a Vandermonde determinant result [28, 29]

$$\Psi_A(x_1, \dots, x_N) = \frac{1}{\sqrt{N!}} C_N \prod_{k=1}^{N} e^{-x_k^2 / 2} \prod_{1 \le i < j \le N} (x_i - x_j), \qquad (9)$$

where $C_N = \frac{1}{\pi^{N/4} \prod_{k=0}^{N-1} \sqrt{2^k k!}}$. We then have after some algebra

$$\left[ \frac{\partial \Psi_A}{\partial x_k} \right]_{x_k = x_{k+1}} = \frac{1}{\sqrt{N!}} C_N e^{-x_k^2 / 2} \prod_{i=1}^{N} e^{-x_i^2 / 2} (x_i - x_k)^2 \prod_{\substack{j, l \neq k, k+1 \\ 1 \le j < l \le N}} (x_j - x_l) \qquad (10)$$

and thus, using again the Vandermonde formula

$$\left[ \frac{\partial \Psi_A}{\partial x_k} \right]_{x_k = x_{k+1}} = \frac{2^{2N-3} e^{-2x_k^2}}{\pi N! (N-1)! (N-2)!} \left[ \det[\phi_{i-1}(x_j)] \right]_{N-2 \times N-2}^2 \prod_{\substack{i=1 \\ i \neq k, k+1}}^{N} (x_i - x_k)^4. \qquad (11)$$

Finally, permutation and parity invariances yield

$$\alpha_k = \frac{2^{2N-3}}{\pi N! (N-1)! (N-2)! (k-1)! (N-k-1)!} \int_{-\infty}^{+\infty} dx_k e^{-2x_k^2} \sum_{P,Q \in S_{N-2}} \epsilon(P) \epsilon(Q) \times \prod_{\substack{i=1 \\ i \neq k, k+1}}^{N} \int_{L_k}^{U_k} dx_i (x_i - x_k)^4 \phi_{P(i)-1}(x_i) \phi_{Q(i)-1}(x_i), \qquad (12)$$

where the integration limits are $(L_k, U_k) = (-\infty, x_k)$ if $i < k$ and $(x_k, +\infty)$ otherwise. An alternative derivation of the coefficients $\alpha_k$ can be found in [26].

## 2. Density profiles

As the first application, using equation (3) together with the solution of equation (7) for the coefficients $a_p$ and equation (12) for the weights $\alpha_k$ we determine the exact density profile of each component of the mixture, according to the definition

$$n_{\nu}(x) = N_{\nu} \int \prod_{j=1}^{N} dx_j \delta(x - x_{\nu}) \Psi^{*}(x_1, \dots, x_N) \Psi(x_1, \dots, x_N), \qquad (13)$$

where we have indicated by $x_\nu$ one of the coordinates corresponding to a Fermion belonging to the $\nu$th component of the mixture. In the following, for each given mixture, we focus on the ground state and first many-body excited state with a symmetry different from the ground state one.
<figure><img src="image_5.png" /><figcaption>**Figure 1.** Density profiles for the ground state (left panel) and for the first many-body excited state with a symmetry different than the ground state (right panel) for three balanced mixtures (i.e. with the same number of particles in each species $N_1 = \cdots = N_r$) of strongly interacting Fermi gases, with different numbers of components $r = 2, 3, 6$ and total particle number $N = 6$ (from top to bottom: $N_\nu = 3, 2, 1$). The density profiles are the same for each component of the mixture, i.e. $n_1(x) = n_2(x) = \cdots = n_r(x)$. The inset shows the corresponding ground state density profiles for the case of the corresponding mixtures of noninteracting Fermions. The corresponding Young tableaux (see text) are also shown in the panels near each profile.</figcaption></figure>

Quite generally, for all the ground state density profiles we find that the total density $n(x) = \sum_{\nu=1}^{r} n_\nu(x)$ coincides with the one of a noninteracting N-particle Fermi gas under external confinement, as previously reported [30]. For equal populations in the various species, the ground state density profile is the same for all the species and coincides apart to a normalization factor with the one of a noninteracting Fermi gas with N particles in harmonic confinement, which is characterized by N peaks [31], as shown in figure 1. The excited states display a variety of profiles, where in particular the two-component mixture does not change and the six-component one shows small deviations with respect to the corresponding ground-state density profiles, while the three-component mixture displays a different, two-peak structure. These differences may be accounted for by considering the different symmetry of the excited states in the three cases, see section 3.

In the case of an imbalanced mixture, the partial density profiles display a rich structure, as presented in figure 2. We first consider the case of an imbalanced two-component mixture. In the polaron case, where $N_1 = N - 1$ and $N_2 = 1$ (top panels in figure 2), for the ground-state density we observe a spatial separation of the polaron density profile at the center of the trap and the majority component at the wings of the trap. By comparing with the results for the corresponding noninteracting gas, one clearly sees the effect of repulsive interactions: the mutual repulsion among the two components push the majority component to a larger region of the trap and a hole is created around the polaron. For the excited state, we find that the density profiles are proportional to the ones of a noninteracting gas. We understand this as being related to the symmetry properties of this particular excited state—as it will be discussed in detail in section 3 below, we find that it has the same symmetry as the noninteracting gas.

In the case of more than one particle in the minority component (central panels in figure 2) for the ground state density profiles, we observe a partial demixing through a more complex structure, with the majority component occupying both the inner core and the external wings of the density profiles. This may be viewed as a mesoscopic realization of an antiferromagnetic configuration [32]. The excited state density profiles have instead a two-peak (majority component) or one-peak (minority component) structure which recalls the ground state of the polaron. Also in this case, the analysis of the symmetry of the profiles brings an explanation, since we find that this excited state has the same symmetry as the ground state of the polaron.

The three component imbalanced mixture (bottom panels in figure 2) displays an even more complex spatial-separated shell structure in the ground state profiles, generalizing the two component case: the minority component occupies the inner of the trap, and the majority component the outer shells, with the intermediate component placed spatially between the other two. For the excited state, we find again a mixed state whose density profiles are proportional to the ones of a noninteracting gas. As it will be discussed in detail in section 3 below, this is in agreement with the fact that this state has the same symmetry as the ground-state of the two-component balanced mixture.
<figure><img src="image_6.png" /><figcaption>**Figure 2.** Density profiles for the ground state (left panels) and for the first many-body excited state with a symmetry different than the ground state (right panels) for three imbalanced mixtures of strongly interacting Fermi gases, all having the same total number of particles $N = 6$. Top panels: a two-component mixture with $N_1 = 5$ (turquoise) and $N_2 = 1$ (dark blue). Central panels: a two-component mixture with $N_1 = 4$ (turquoise) and $N_2 = 2$ (dark blue). Bottom panels: a three-component mixture with $N_1 = 3$ (turquoise), $N_2 = 2$ (dark blue), $N_3 = 1$ (cyan). The insets show the corresponding ground-state density profiles for the same mixture of noninteracting gases. The corresponding Young tableau (see text) is also shown in each panel.</figcaption></figure>

### 3. Symmetry characterization

We now analyze the symmetry properties of the quantum many body states under exchange of particles [33–36]. Obviously, single component free Fermions are described by a totally anti-symmetric wavefunction. In the presence of several components interacting among each other, the wavefunction must be totally antisymmetric under exchange of particles belonging to each component but the total wavefunction has a more involved symmetry under exchange of an arbitrary pair of particles. The actual symmetry of the ground state and the excited states can be deduced from the properties of the permutation group of $N$ elements, $S_N$ [37]. More precisely, one can expand the eigenstates of the quantum system over the different irreducible representations of this group. We follow this route below, and demonstrate that the ground state (and any excited state) has a well defined symmetry, i.e. can be described by a single irreducible representation, identified by a Young tableau.

The previous diagonalization process allows us to obtain a set of $S$ values of $K$, $K_1 \le \cdots \le K_S$, with $K_1 = 0$ and $K_S = K_{max}$, and a set of associated eigenvectors $\vec{A}_1, \dots, \vec{A}_S$ that correspond to decreasing-energy solutions at finite interactions. In order to completely characterize the symmetry of the various states ($K_\ell, \vec{A}_\ell$) of the degenerate manifold at $g = \infty$, we determine to which irreducible representation of $S_N$, i.e. to which Young
<table><caption>**Table 1.** Symmetries and corresponding K eigenvalue (see equation (7)) for the multi-component Fermionic mixtures of figures 1 and 2.</caption><thead><tr><th>System</th><th>A<sub>max</sub></th><th>K<sub>max</sub></th><th>A<sub>excited</sub></th><th>K<sub>excited</sub></th></tr></thead><tbody><tr><td>r = 2, N<sub>1</sub> = N<sub>2</sub> = 3</td><td>Y<sub>-3</sub></td><td>24.97</td><td>Y<sub>-5</sub></td><td>18.91</td></tr><tr><td>r = 3, N<sub>1</sub> = N<sub>2</sub> = N<sub>3</sub> = 2</td><td>Y<sub>3</sub></td><td>30.63</td><td>Y<sub>0</sub></td><td>28.96</td></tr><tr><td>r = 6, N<sub>1</sub> = ... = N<sub>6</sub> = 1</td><td>Y<sub>15</sub></td><td>34.33</td><td>Y<sub>9</sub></td><td>33.35</td></tr><tr><td>r = 2, N<sub>1</sub> = 5, N<sub>2</sub> = 1</td><td>Y<sub>-9</sub></td><td>14.60</td><td>Y<sub>-15</sub></td><td>0</td></tr><tr><td>r = 2, N<sub>1</sub> = 4, N<sub>2</sub> = 2</td><td>Y<sub>-5</sub></td><td>22.70</td><td>Y<sub>-9</sub></td><td>14.60</td></tr><tr><td>r = 3, N<sub>1</sub> = 3, N<sub>2</sub> = 2, N<sub>3</sub> = 1</td><td>Y<sub>0</sub></td><td>28.96</td><td>Y<sub>-3</sub></td><td>24.97</td></tr></tbody></table>

tableau, the solution $\vec{A}_\ell$ corresponds. The Young tableaux are defined in the standard fashion [37]: elements belonging to the same line (column) are symmetric (antisymmetric) under exchange. Thus, for example, for six particles divided in six components, a completely antisymmetric state $\vec{A}_\ell$ will correspond to $\Box\Box\Box\Box\Box\Box$ whereas a completely symmetric state will correspond to $\Box\Box\Box\Box\Box\Box$. Note that these two states are associated with a single-component non-interacting Fermi gas and a single component Tonks gas, respectively.

In order to classify those states according to their symmetry, we use the k-cycle class sums operators $\Gamma^{(k)}$ [20, 21], defined by $\Gamma^{(k)} = \sum_{i_1<...<i_k} (i_1...i_k)$, where $(i_1...i_k)$ is the cyclic permutation of particles $i_1, ..., i_k$. In our system, on the basis of the coordinate sectors, $\Gamma^{(k)}$ is a $N! \times N!$ operator whose elements are $\Gamma_{ij}^{(k)} = (-1)^{(k-1)}$ if going from sector $i$ to sector $j$ only exchange k Fermions in a cyclic fashion, and $\Gamma_{ij}^{(k)} = 0$ otherwise. Note that, also in this case, we can reduce the dimension of such operator to an $S \times S$ matrix, by summing over the contributions of the sectors associated with one snippet.

The spectral decomposition of $\Gamma^{(k)}$ allows to associate each Young tableau Y with a given eigenvalue $\gamma_k$ [20, 21, 24]. In the following, we shall use in particular the class-sum operator $\Gamma^{(2)}$. Our method consists first in computing and diagonalizing the transposition class sum $\Gamma^{(2)}$ for a given system. Its eigenvalues $\gamma_2$ can be linked to the Young tableaux according to the expression

$$ \gamma_2 = \frac{1}{2} \sum_i \lambda_i (\lambda_i - 2i + 1), \qquad (14) $$

where $i$ and $\lambda_i$ refer respectively to the line and number of boxes in this line of Young tableau. Thus, projecting a given solution $\vec{A}_\ell$ over the eigenbasis of $\Gamma^{(2)}$ allows to characterize its symmetry and to analyze it in terms of Young tableaux. In table 1 we summarize some of our results for the ground states ($\vec{A}_{\max}$, $K_{\max}$) and the first excited states with a different symmetry ($\vec{A}_{\text{excited}}$, $K_{\text{excited}}$) of different systems for $N = 6$ particles.

We set [37]

$$ Y_{15} = \Box\Box\Box\Box\Box\Box, \quad Y_{-9} = \Box\Box\Box\Box\Box\Box, \quad Y_{-5} = \Box\Box\Box\Box\Box\Box, \qquad (15) $$

$$ Y_{-3} = \Box\Box\Box, \quad Y_0 = \Box\Box\Box\Box, $$

and the $Y_\gamma$ are obtained by taking the symmetric of $Y_\gamma$ with respect to the main diagonal.

Table 1 shows that the ground and excited states constructed with equation (3) have a definite symmetry which can be readily extracted from the associated Young tableau. An exception is provided by the ground state of the case $N_1 = N_2 = 3$ as well as the excited state of $N_1 = 3, N_2 = 2, N_3 = 1$, where the transposition class-sum operator $\Gamma^{(2)}$ does not allow to uniquely associate a Young tableau to the wavefunction, since equation (14) gives two Young tableaux corresponding to the eigenvalue $-3$, $\Box\Box$ and $\Box\Box\Box$, and also to the eigenvalue $3$, $\Box\Box\Box$ and $\Box\Box\Box\Box$. In this case, the ambiguity is lift off with the help of the 3-cycle class sum $\Gamma^{(3)}$, since these two Young tableaux correspond to two different eigenvalues of $\Gamma^{(3)}$ [20].

Furthermore, for all the mixture considered we note that the ground state, corresponding to $K = K_{\max}$, is the most symmetric one [38, 39], and the excited state is obtained by decreasing the symmetry of the state by taking out one cell from the top row and putting it in the first available lower row of the Young tableau, thus
making it more antisymmetric. The opposite case of $K = 0$ is associated with the most antisymmetric Young tableau. The above result for the ground state supports the observation of [40], where an ansatz for the ground state wavefunction was suggested, and is also in agreement with a general demonstration provided in [41]. We notice also that when the number of particles coincides with the number of components, i.e. $r = N$, the ground state is fully symmetric, and has the same symmetry as the bosonic Tonks–Girardeau gas.

This analysis provides a verification and an example of the generalization of the Lieb–Mattis theorem [22] for the case of multicomponent Fermionic mixtures. The theorem states that for *N* electrons in one-dimension, interacting by a symmetric potential, the energy of a state with total spin *S* (*S*′) is such that *E*(*S*) ≤ *E*(*S*′) if *S* < *S*′. It follows that the ground state has the smallest possible value for the spin *S*, which is realized by an antisymmetric spinor. Hence, the spatial wavefunction has the most symmetric configuration. In the case of more than two spin components, as also discussed in [38], the same feature occurs, which is displayed by our results. Since *K* tends to be maximized when the wavefunction is more symmetric, we can see *K* as an energetic indicator of the symmetry.

## 4. Summary and conclusions

In this work we have considered a multicomponent strongly correlated Fermionic mixture with up to six components. Using a generalization of the pioneering solution due to Girardeau for the Bose gas and of the recently developed solution for the two-component Fermi gas we have determined the exact many-body wavefunctions for the degenerate manifold at infinite interactions. We have identified the one which corresponds to the ground state at finite interactions as the one which has the maximum value for the parameter *K*, related to the interaction energy. We have then obtained the density profiles for the mixture under harmonic confinement. For an imbalanced mixture we have found a partial phase separation among the components, which is an effect of the strong repulsive interactions. Furthermore, we have characterized the symmetry properties of the ground and of some excited state wavefunctions of the manifold by introducing suitable class-sum operators, and we have shown that the ground-state wavefunctions have a definite symmetry, corresponding to the most symmetric (or less antisymmetric) configuration compatible with the imbalance among the components. Our exact solution for the inhomogeneous multicomponent mixture provides an important benchmark for numerical simulations of strongly correlated multicomponent Fermi gases, in a regime where the presence of the quasi-degenerate manifold challenges the convergence of the calculations, as well as for quantum simulators.

## Acknowledgments

This paper is dedicated to the memory of Marvin Girardeau, of whom AM recalls his enthusiastic guidance. We thank M Gattobigio, J Levinsen, P Massignan, M Parish and T Roscilde for fruitful discussions. AM acknowledges the hospitality of the Centre de Recherches Mathématiques de Montréal where part of this work was performed and financial support from ANR projects Mathostaq (ANR-13-JS01-0005-01) and SuperRing (ANR-15-CE30-0012-02).

## References

[1] Pagano G *et al* 2014 A one-dimensional liquid of Fermions with tunable spin *Nat. Phys.* **10** 198–201

[2] McGuire J B 1965 Interacting Fermions in one-dimension: I. Repulsive potential *J. Math. Phys.* **6** 432

[3] Gaudin M 1967 *Un système à une-dimension de Fermions en interaction* *Phys. Lett.* **A** 24 55–6

[4] Yang C N 1967 Some exact results for the many-body problem in one-dimension with repulsive delta-function interaction *Phys. Rev. Lett.* **19** 1312–5

[5] Sutherland B 1968 Further results for the many-body problem in one-dimension *Phys. Rev. Lett.* **20** 98–100

[6] Matveev K A and Glazman L I 1993 Coulomb blockade of tunneling into a quasi-one-dimensional wire *Phys. Rev. Lett.* **70** 990–3

[7] Giamarchi T 2003 *Quantum Physics in One Dimension* (Oxford: Clarendon)

[8] Guan X-W, Batchelor M T and Lee C 2013 Fermi gases in one-dimension: from Bethe Ansatz to experiments *Rev. Mod. Phys.* **85** 1633–91

[9] Orignac E and Citro R 2012 Response functions in multicomponent Luttinger liquids *J. Stat. Mech.: Theory Exp.* **2012** P12020

[10] Frahm H and Schadschneider A 1993 Critical exponents of the degenerate Hubbard model *J. Phys. A: Math. Gen.* **26** 1463

[11] He P, Lee J Y, Guan X, Batchelor M T and Wang Y 2011 Thermodynamics and spin-charge separation of one-dimensional strongly repulsive three-component Fermions *J. Phys. A: Math. Theor.* **44** 405005

[12] Zürn G, Serwane F, Lompe T, Wenz A N, Ries M G, Bohn J E and Jochim S 2012 Fermionization of two distinguishable Fermions *Phys. Rev. Lett.* **108** 075303

[13] Girardeau M D 1960 Relationship between systems of impenetrable Bosons and Fermions in one-dimension *J. Math. Phys.* **1** 516

[14] Girardeau M D and Minguzzi A 2007 Soluble models of strongly interacting ultracold gas mixtures in tight waveguides *Phys. Rev. Lett.* **99** 230402
[15] Gharashi S E and Blume D 2013 Correlations of the upper branch of 1D harmonically trapped two-component Fermi gases *Phys. Rev. Lett.* **111** 045302

[16] Cui X and Ho T-L 2014a Ground-state ferromagnetic transition in strongly repulsive one-dimensional Fermi gases *Phys. Rev. A* **89** 023611

[17] Volosniev A G, Fedorov D V, Jensen A S, Zinner N T and Valiente M 2014 Strongly interacting confined quantum systems in one-dimension *Nat. Commun.* **5** 5300

[18] Lindgren E J, Rotureau J, Forssén C, Volosniev A G and Zinner N T 2014 Fermionization of two-component few-Fermion systems in a one-dimensional harmonic trap *New J. Phys.* **16** 063003

[19] Pęcak D, Gajda M and Sowiński T 2016 Two-flavor mixture of a few Fermions of different mass in a one-dimensional harmonic trap *New J. Phys.* **18** 013030

[20] Katriel J 1993 Representation-free evaluation of the eigenvalues of the class-sums of the symmetric group *J. Phys. A: Math. Gen.* **26** 135

[21] Novoselsky A and Katriel J 1994 Hyperspherical functions with arbitrary permutational symmetry *Phys. Rev. A* **49** 833

[22] Lieb E and Mattis D 1962 Theory of ferromagnetism and the ordering of electronic energy levels *Phys. Rev.* **125** 164–72

[23] Deuretbacher F, Fredenhagen K, Becker D, Bongs K, Sengstock K and Pfannkuche D 2008 Exact solution of strongly interacting quasi-one-dimensional spinor Bose gases *Phys. Rev. Lett.* **100** 160405

[24] Fang B, Vignolo P, Gattobigio M, Miniatura C and Minguzzi A 2011 Exact solution for the degenerate ground-state manifold of a strongly interacting one-dimensional Bose–Fermi mixture *Phys. Rev. A* **84** 023626

[25] Barth M and Zwerger W 2011 Tan relations in one-dimension *Ann. Phys.* **326** 2544

[26] Deuretbacher F, Becker D, Bjerlin J, Reimann S and Santos L 2014 Quantum magnetism without lattices in strongly interacting one-dimensional spinor gases *Phys. Rev. A* **90** 013611

[27] Loft N J S, Kristensen L B, Thomsen A E and Zinner N T 2016 Comparing models for the ground state energy of a trapped one-dimensional Fermi gas with a single impurity arXiv:1508.05917

[28] Girardeau M D, Wright E M and Triscari J M 2001 Ground-state properties of a one-dimensional system of hard-core bosons in a harmonic trap *Phys. Rev. A* **63** 033601

[29] Forrester P J, Frankel N E, Garoni T M and Witte N S 2003 Finite one-dimensional impenetrable Bose systems: occupation numbers *Phys. Rev. A* **67** 043607

[30] Grining T, Tomza M, Lesiuk M, Przybytek M, Musia M, Massignan P, Lewenstein M and Moszynski R 2015 Many interacting Fermions in a one-dimensional harmonic trap: a quantum-chemical treatment *New J. Phys.* **17** 115001

[31] Vignolo P, Minguzzi A and Tosi M P 2000 Exact particle and kinetic-energy densities for one-dimensional confined gases of noninteracting Fermions *Phys. Rev. Lett.* **85** 2850–3

[32] Murmann S, Deuretbacher F, Zürn G, Bjerlin J, Reimann S M, Santos L, Lompe T and Jochim S 2015 Antiferromagnetic Heisenberg spin chain of a few cold atoms in a one-dimensional trap *Phys. Rev. Lett.* **115** 215301

[33] Harshman N L 2014 Spectroscopy for a few atoms harmonically trapped in one-dimension *Phys. Rev. A* **89** 033633

[34] Harshman N L 2016 One-dimensional, traps two-body interactions, few-body symmetries: I. One, two, and three particles *Few-Body Syst.* **57** 11

[35] Harshman N L 2016 One-dimensional traps, two-body interactions, few-body symmetries: II. N particles *Few-Body Syst.* **57** 45

[36] Guan L, Chen S, Wang Y and Ma Z-Q 2009 Exact solution for infinitely strongly interacting fermi gases in tight waveguides *Phys. Rev. Lett.* **102** 160402

[37] Hamermesh M 1989 *Group Theory and Its Applications to Physical Problems* (New York: Dover)

[38] Deuretbacher F, Becker D, Bjerlin J, Reimann S M and Santos L 2015 Quantum magnetism without lattices in strongly interacting one-dimensional spinor gases *Phys. Rev. A* **90** 013611

[39] Deuretbacher F, Becker D and Santos L 2016 Momentum distributions and numerical methods for strongly interacting one-dimensional spinor gases arXiv:1602.06816

[40] Levinsen J, Massignan P, Bruun G M and Parish M M 2015 Strong-coupling Ansatz for the one-dimensional Fermi gas in a harmonic potential *Sci. Adv.* **1** e1500197

[41] Hakobyan T 2010 Ordering of energy levels for extended SU(N) Hubbard chain *SIGMA* **6** 24