Fermionization of a Few-Body Bose System Immersed into a Bose-Einstein Condensate

Tim Keller,Thomás Fogarty,Thomas Busch

DOI: https://doi.org/10.21468/SciPostPhys.15.3.095

2023-06-30

Abstract:We study the recently introduced self-pinning transition [Phys. Rev. Lett. 128, 053401 (2022)] in a quasi-one-dimensional two-component quantum gas in the case where the component immersed into the Bose-Einstein condensate has a finite intraspecies interaction strength. As a result of the matter-wave backaction, the fermionization in the limit of infinite intraspecies repulsion occurs via a first-order phase transition to the self-pinned state, which is in contrast to the asymptotic behavior in static trapping potentials. The system also exhibits an additional superfluid state for the immersed component if the interspecies interaction is able to overcome the intraspecies repulsion. We approximate the superfluid state in an analytical model and derive an expression for the phase transition line that coincides with well-known phase separation criteria in binary Bose systems. The full phase diagram of the system is mapped out numerically for the case of two and three atoms in the immersed component.

Quantum Gases,Quantum Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the phenomenon of self - pinning transition in a quasi - one - dimensional two - component quantum gas system when one of the components is immersed in a Bose - Einstein Condensate (BEC) and has a finite intra - species interaction strength. Specifically, the research focuses on the process in which the system reaches the self - pinning state through a first - order phase transition due to matter - wave backaction in the infinite intra - species repulsion limit, which is in contrast to the asymptotic behavior in a static trapping potential. In addition, if the inter - species interaction can overcome the intra - species repulsion, the system will also exhibit an additional superfluid state. The paper approximates the superfluid state through an analytical model and derives an expression for the phase - transition line, which is consistent with the known phase - separation criteria in binary Bose systems. Finally, the complete phase diagrams of the system with two and three atoms in the immersed component are numerically plotted.

### Core problems of the paper
1. **Self - pinning transition**: Research on the process in which the system reaches the self - pinning state through a first - order phase transition due to matter - wave backaction in the infinite intra - species repulsion limit.
2. **Existence conditions of superfluid state**: Explore whether the system will exhibit a superfluid state if the inter - species interaction can overcome the intra - species repulsion under a finite intra - species interaction strength.
3. **Analytical expression of phase - transition line**: Derive an expression for the phase - transition line and verify its consistency with the known phase - separation criteria in binary Bose systems.
4. **System phase diagram**: Numerically plot the phase diagrams of the system under different parameters, especially for the immersed components containing two and three atoms.

### Main methods and results
- **Model description**: Use the coupled Schrödinger equations to describe the interaction between the two components.
- **Numerical simulation**: Solve the coupled equations by the self - consistent imaginary - time evolution method to calculate the energy of the system.
- **Analytical model**: Based on the Thomas - Fermi approximation, derive the density distributions and energy expressions of the superfluid state and the self - pinning state.
- **Phase - diagram plotting**: Plot the phase diagram of the system by calculating the coherence (i.e., the condensate fraction) of the system, distinguishing between the superfluid state and the self - pinning state.

### Key formulas
- **Superfluid - state density**:
\[
\rho_{\text{sf}}(x)=\frac{N a_{\text{sf}}}{2}\frac{1}{\cosh^{2}(a_{\text{sf}}x)}
\]
where
\[
a_{\text{sf}}(g = 0)=N a_{0}\sqrt{\frac{1 + 2N^{2}\epsilon^{-1}}{N^{2}\epsilon}}
\]
\[
a_{0}=\frac{m g_{m}^{2}}{2g_{c}\hbar^{2}}
\]
\[
\epsilon=\frac{6a_{0}^{2}\hbar^{2}}{5m\tilde{\mu}},\quad\tilde{\mu}=\mu_{0}\left(1+\frac{g_{m}N}{g_{c}N_{c}}\right)
\]

- **Total system energy**:
\[
E_{\text{sf}}=N\left[\frac{\hbar^{2}g_{m}^{2}}{30m\tilde{\mu}g_{c}}N a_{\text{sf}}^{3}+\frac{\hbar^{2}}{m}\frac{a_{\text{sf}}^{2}}{6}-\frac{g_{m}^{2}}{6g_{c}}N a_{\text{sf}}+\frac{(N - 1)g a_{\text{sf}}}{6}\right]+\frac{\tilde{\mu}^{2}L_{c}}{2g_{c}}
\]

- **Self - pinning - state density**:
\[
\rho_{\text{pin}}(x)=\frac{a_{\text{pin}}}{2}\sum_{n = 1}