The odd fermion at the edge: odd-even staggering in the trapped, unitary Fermi gas

Résumé Technique : Le fermion impaire à la bordure d'un gaz de Fermi unitaire piégé

Énoncé du problème
L'article étudie l'oscillation impaire-paire de l'énergie de l'état fondamental dans un gaz de Fermi unitaire à piège harmonique avec une grande charge de nombre de particules $Q$. Dans ce système, les fermions interagissent avec une longueur de diffusion $s$-wave infinie. Pour un $Q$ pair, l'état fondamental est un superfluide apparié. Pour un $Q$ impair, un fermion unique non apparié subsiste. La question centrale est le comportement et la mise à l'échelle de l'énergie de ce quasiparticle non apparié. Des travaux antérieurs par Son [11] ont soutenu qu'à grand $Q$, ce quasiparticle ne réside pas dans le volume (où le gap d'appariement est maximal) mais se localise à la bordure du nuage (la surface de Thomas-Fermi), là où le gap local s'annule. Cette localisation implique une énergie d'excitation paramétriquement petite par rapport au gap du volume. L'article vise à fournir une base systématique pour cet argument du « quasiparticle de bordure », à dériver les lois de mise à l'échelle spécifiques et à estimer le coefficient universel régissant la division d'énergie.

Méthodologie
Les auteurs emploient une approche en trois volets pour analyser le système :

1. Approximation de la sonde (Analyse BdG) :
   Le fermion non apparié est traité comme un quasiparticle sonde se déplaçant dans un fond fixe de condensat superfluide, en négligeant la rétroaction. Les auteurs utilisent l'Hamiltonien de Bogoliubov-de Gennes (BdG) dans l'approximation de la densité locale (LDA). En développant l'Hamiltonien au voisinage du rayon classique $R_{cl}$ du nuage, ils identifient une « limite distinguée » où les termes cinétiques et potentiels évoluent de manière similaire. Cela réduit le problème microscopique à un système d'Airy couplé universel dans la direction perpendiculaire à la bordure, couplé à un mouvement sur la frontière sphérique.

2. Théorie effective des champs (EFT) :
   Pour aller au-delà de l'approximation de la sonde et des hypothèses de champ moyen, les auteurs construisent une EFT à grande charge. Ce cadre développe l'action en puissances négatives de la grande charge $Q$. Le volume est décrit par le mode de Goldstone du superfluide, tandis que la bordure est augmentée par un mode fermionique confiné à la frontière. En utilisant le théorème de co-aire, l'action de bordure est reformulée pour décrire un fermion couplé au champ de Goldstone du superfluide. Cette méthode capture les interactions entre le fermion non apparié et les phonons du volume (modes de Goldstone) qui sont gelés dans la limite de la sonde.

3. Analyse numérique :
   Les auteurs effectuent une diagonalisation numérique des équations radiales BdG complètes piégées sans supposer a priori la localisation de l'état. Ils discrétisent l'Hamiltonien sur une grille et utilisent un solveur Arnoldi pour cibler les modes de basse énergie proches de zéro. Ils analysent les valeurs propres et les fonctions propres résultantes pour vérifier la localisation à la bordure et extraire les exposants de mise à l'échelle et les coefficients.

Contributions clés et résultats

• Localisation à la bordure et mise à l'échelle :
  Les trois méthodes confirment que pour un $Q$ impair, le fermion non apparié forme un quasiparticle localisé à la bordure. L'énergie de cet état, qui domine l'oscillation impaire-paire $\Delta_{\text{stag}}(Q)$, suit la loi :
  $$\Delta_{\text{stag}}(Q) \approx \chi \xi^{1/6} (24Q)^{1/9} \hbar\omega$$
  où $\xi$ est le paramètre de Bertsch et $\chi$ est un coefficient dimensionnel sans unité universel.

• Structure spectrale :
  Il est démontré que le spectre de basse énergie consiste en des grappes d'états. L'énergie de premier ordre est déterminée par la plus petite valeur propre positive ($\chi_0$) d'un problème d'Airy transverse universel perpendiculaire à la bordure. Il existe une structure plus fine provenant des modes de moment angulaire sur la sphère, séparés par une échelle d'énergie de l'ordre de $1/R_{cl}^2$. La masse effective du fermion de bordure, $M_\parallel$, est trouvée être de l'ordre de l'unité mais distincte de la masse du fermion de volume en raison du mélange particule-trou.

• Coefficient universel $\chi$ :
  L'analyse numérique donne une valeur pour le paramètre de Son $\chi$. En utilisant des paramètres d'entrée motivés par la QMC pour le rapport de gap $y$, les auteurs trouvent :
  $$\chi = 0,72(2)$$
  Cette valeur est cohérente avec les ajustements de données existantes de diffusion de type Monte Carlo à nœuds fixes (FN-DMC).

• Couplages EFT :
  La formulation EFT identifie deux constantes de couplage (LEC) : le coefficient $e_0$ (identifié avec $\chi$) et la masse effective $M_\parallel$. L'EFT dérive également explicitement le couplage entre le courant du fermion de bordure et la vitesse du champ de Goldstone, décrivant un décalage de type Doppler pour les excitations de bordure se propageant avec ou contre l'écoulement superfluide.

Signification et revendications
L'article prétend fournir un fondement théorique systématique pour l'hypothèse du quasiparticle de bordure de Son. En combinant des arguments de mise à l'échelle analytiques (BdG et EFT) avec une vérification numérique, les auteurs démontrent que l'oscillation impaire-paire est un phénomène de bordure plutôt que de volume.

Les auteurs soulignent que leur principale contribution est la robustesse de la loi de mise à l'échelle à grand $Q$ ($Q^{1/9}$) et l'existence de l'état localisé à la bordure. Ils déclarent explicitement que la comparaison avec les données de simulation de nombreux corps existantes sert principalement de test de cohérence de l'image de la bordure et de la mise à l'échelle prédite, plutôt que d'une détermination de précision des énergies à $Q$ fini. Les énergies détaillées à $Q$ fini sont notées comme étant sensibles au choix des paramètres d'entrée du volume (tels que le paramètre de Bertsch et le rapport de gap), tandis que les exposants de mise à l'échelle et la nature d'ordre un du coefficient universel sont des prédictions robustes du cadre.

Ce travail fait le pont entre la physique nucléaire de nombreux corps et les gaz atomiques froids fortement interagissants, en utilisant le gaz de Fermi unitaire comme un laboratoire pour la dynamique conforme non relativiste et les techniques d'EFT à grande charge.