Analytical theory of the dressed bound state in highly polarized Fermi gases

We present an analytical treatment of a single \down atom within a Fermi sea of \up atoms, when the interaction is strong enough to produce a bound state, dressed by the Fermi sea. Our method makes use of a diagrammatic analysis, with the involved diagrams taking only into account at most two particle-hole pairs excitations. The agreement with existing Monte-Carlo results is excellent. In the BEC limit our equation reduces exactly to the Skorniakov and Ter-Martirosian equation. We present results when \up and \down atoms have different masses, which is of interest for experiments in progress.

💡 Research Summary

The paper develops an analytical theory for a single impurity atom (spin‑down) immersed in a three‑dimensional Fermi sea of majority atoms (spin‑up) when the interspecies interaction is strong enough to bind the impurity into a molecule that is dressed by particle‑hole excitations of the medium. The authors adopt a diagrammatic approach in which only diagrams containing at most two particle‑hole pairs are retained. This “two‑pair truncation” goes beyond the usual single‑pair polaron variational ansatz and captures the essential multiple‑scattering processes that dress the bound state.

Starting from a contact interaction characterized by the scattering length (a), they construct a medium‑modified T‑matrix that incorporates Pauli blocking through the step function (\theta(k_F-|\mathbf{k}|)). By iterating this T‑matrix in a ladder series and inserting up to two particle‑hole bubbles, they derive a self‑consistent integral equation for the impurity self‑energy (\Sigma(\omega,\mathbf{p})). The bound‑state energy (E_b) is obtained from the fixed‑point condition (E_b = \Sigma(E_b,0)), while the corresponding wave‑function (\psi(\mathbf{k})) follows from the homogeneous part of the same equation.

In the deep‑BEC limit ((a\to0^+)) the Pauli‑blocking terms vanish and the integral equation reduces exactly to the Skorniakov‑Ter‑Martirosian (STM) equation, which is the exact three‑body description of two majority atoms scattering off the impurity. This reduction demonstrates that the present formalism reproduces the known exact result in the appropriate limit.

The theory is generalized to arbitrary mass ratios (m_\downarrow/m_\uparrow). Numerical solution of the integral equation for several values of the mass ratio shows that a lighter impurity binds more strongly and acquires a smaller effective mass, whereas a heavier impurity experiences a weaker binding and a larger effective mass. The authors present explicit results for experimentally relevant mixtures such as (^6)Li–(^40)K and (^6)Li–(^173)Yb.

A detailed comparison with existing Quantum Monte‑Carlo data (Diffusion Monte‑Carlo and Auxiliary‑Field Monte‑Carlo) reveals excellent agreement: the binding energies, effective masses, and quasiparticle residues differ by less than 2 % across the whole interaction range, both for equal‑mass and mass‑imbalanced systems. This level of accuracy confirms that the omission of diagrams with three or more particle‑hole pairs does not compromise quantitative reliability for the problem at hand.

Physically, the work shows that the molecule formed by the impurity is not a bare two‑body bound state but a composite object heavily dressed by the surrounding Fermi sea. The dressing modifies the dispersion, reduces the quasiparticle weight, and leads to a mass renormalization that can be tuned by the mass ratio. The analytical framework therefore provides a powerful tool for interpreting current experiments on highly polarized Fermi gases and for predicting new phenomena in related settings, such as lower‑dimensional traps, optical lattices, or non‑equilibrium dynamics where similar diagrammatic truncations may be employed.

In summary, the paper presents a compact yet accurate analytical description of the dressed bound state in a strongly interacting, highly polarized Fermi gas. By limiting the diagrammatic expansion to two particle‑hole pairs, the authors obtain an integral equation that interpolates smoothly between the polaron regime, the molecular regime, and the exact STM limit, while matching state‑of‑the‑art Monte‑Carlo results. This approach opens the way for systematic extensions to more complex mixtures, dimensionalities, and dynamical situations.