Pressure and Size Dependence of Roton Emission and Vortex Creation by Moving Objects in He~II in $T \to 0$ Limit: Generalized Nonlocal Gross-Pitaevskii Model

Cet article présente une étude numérique utilisant un modèle généralisé non local de Gross-Pitaevskii pour analyser comment la pression et la taille de l'obstacle influencent les vitesses critiques pour l'émission de rotons et la nucléation de vortex par des objets en mouvement dans l'hélium-4 superfluide à température nulle, marquant le premier cadre théorique à traiter simultanément ces mécanismes dépendants de la pression.

L'essentiel

Imaginez un monde où un liquide peut s'écouler sans aucune friction, tel un fantôme glissant entre vos doigts. Il s'agit de l'Hélium Superfluide (He II), un état particulier de la matière qui n'existe que lorsque l'hélium est refroidi à des températures proches du zéro absolu.

Ce papier est comparable à un laboratoire de simulation haute technologie où les auteurs soumettent ce liquide fantomatique à l'épreuve. Ils voulaient comprendre ce qui se produit lorsque l'on pousse un petit objet (comme une minuscule bulle ou une particule chargée) à travers ce superfluide. Plus précisément, ils souhaitaient savoir : À quelle vitesse pouvez-vous le pousser avant que le liquide ne « brise » son écoulement parfait ?

Voici l'histoire de leur découverte, décomposée en concepts simples :

1. Les Deux Façons dont le Liquide Se Brise

Lorsque vous poussez un objet à travers l'hélium superfluide, le liquide ne subit pas simplement une « traînée » comme l'eau. Au lieu de cela, il réagit de deux manières distinctes une fois atteint une certaine limite de vitesse (appelée vitesse critique) :

• L'Explosion « Roton » : Imaginez le liquide comme un océan calme. Si vous poussez l'objet trop vite, vous ne créez pas seulement des vagues ; vous créez soudainement un essaim de minuscules particules énergétiques appelées rotons. C'est comme si le liquide décidait soudainement de se briser en un million d'étincelles énergétiques minuscules. Cela se produit à une vitesse spécifique.
• Le Tourbillon « Vortex » : Si vous le poussez encore plus vite (ou si l'objet est suffisamment grand), le liquide commence à tourner. Il crée de minuscules tornades microscopiques appelées vortex quantiques. Ce sont comme de minuscules tourbillons qui s'accrochent à l'objet, l'entraînant vers le bas.

L'objectif principal du papier était de déterminer exactement à quelle vitesse il faut aller pour déclencher les « étincelles » (rotons) par rapport aux « tourbillons » (vortex).

2. L'Expérience du Cocotte-minute

Les auteurs n'ont pas seulement observé le liquide à une seule pression. Ils ont simulé ce qui se produit à mesure qu'ils compriment l'hélium de plus en plus fort, depuis le vide (0 bar) jusqu'au point où il se transformerait en rocher solide (environ 25 bars).

Ils ont utilisé un modèle mathématique spécial (un « modèle de Gross-Pitaevskii non local généralisé ») qui agit comme un moteur de jeu vidéo ultra-précis. Ce moteur a été programmé pour imiter le comportement réel et complexe des atomes d'hélium, y compris le comportement étrange des « rotons » que les équations physiques standards manquent généralement.

3. La Grande Découverte : La Compression Modifie les Règles

Voici ce qu'ils ont découvert, en utilisant une analogie simple :

Imaginez que vous essayez de courir à travers une foule de personnes.

• À basse pression (foule lâche) : Les gens sont espacés. Il est en fait assez difficile de déclencher une émeute (créer un roton) car ils sont loin les uns des autres. Mais si vous courez assez vite, vous pourriez faire trébucher quelqu'un et déclencher une réaction en chaîne de personnes qui tombent (vortex).
• À haute pression (foule serrée) : Les gens sont serrés épaule contre épaule. Maintenant, il est beaucoup plus facile de déclencher une émeute (rotons) car ils sont si proches les uns des autres. Cependant, il devient plus difficile de déclencher une réaction en chaîne de personnes qui tombent (vortex) car la foule est si dense et rigide qu'elle résiste à la rotation.

Les Résultats :

• Vitesse des Rotons : À mesure qu'ils comprimaient l'hélium (augmentation de la pression), la vitesse nécessaire pour créer ces « étincelles » (rotons) a diminué. Vous n'avez pas besoin de courir aussi vite pour briser l'écoulement.
• Vitesse des Vortex : À mesure qu'ils comprimaient l'hélium, la vitesse nécessaire pour créer les « tourbillons » (vortex) a augmenté. Vous devez courir beaucoup plus vite pour faire tourner le liquide.

4. Le « Point Doux » pour la Détection

Cela crée un écart fascinant. À haute pression, il existe une large gamme de vitesses où vous pouvez créer les « étincelles » (rotons) sans créer les « tourbillons » (vortex).

Par le passé, les scientifiques avaient du mal à étudier les rotons car ils étaient souvent cachés par les tourbillons désordonnés. Les auteurs suggèrent qu'en comprimant l'hélium à haute pression, nous pouvons créer un environnement « propre » où les rotons apparaissent seuls, les rendant beaucoup plus faciles à étudier.

5. La Taille Compte

Le papier a également examiné la taille de l'objet se déplaçant à travers le liquide.

• Objets minuscules (comme un ion unique) : Ils sont très sensibles. Ils atteignent d'abord la « limite des rotons ».
• Grands objets (comme un grand disque) : Ils sont moins sensibles aux rotons. Ils ont tendance à atteindre d'abord la « limite des vortex », quelle que soit la pression.

Résumé

Les auteurs ont construit un microscope numérique pour observer l'hélium superfluide sous pression. Ils ont découvert que comprimer l'hélium rend plus facile la création d'étincelles d'énergie (rotons) mais plus difficile la création de tourbillons tournants (vortex).

Cela explique pourquoi les expériences passées ont observé des comportements différents à différentes pressions et suggère que si nous voulons étudier les mystérieuses particules « rotons », nous devrions mener nos expériences sous haute pression, où le liquide est plus susceptible de révéler ses secrets sans se salir avec des tourbillons.

Note : Les auteurs admettent que leur simulation a été réalisée en deux dimensions (une tranche plate du monde) car la faire en 3D complète est trop lourde en calculs, mais ils estiment que la physique qu'ils ont découverte reste vraie pour le monde réel en 3D.

Résumé technique

Résumé technique : Dépendance en pression et en taille de l'émission de rotons et de la création de vortex par des objets en mouvement dans l'hélium II à la limite T → 0

Énoncé du problème
La rupture de la superfluidité dans l'hélium II (He II) est régie par deux mécanismes principaux lorsqu'un objet se déplace dans le fluide : l'émission de rotons (excitations élémentaires) et la nucléation de vortex quantifiés. Bien que la dépendance en pression de la vitesse critique pour l'émission de rotons ($v_{rot}^{crit}$) et la création de vortex ($v_{v}^{crit}$) ait été investiguée expérimentalement — notamment par le groupe de Lancaster utilisant des ions négatifs —, ces phénomènes ont historiquement manqué d'une description théorique unifiée. Les observations expérimentales indiquent qu'à basse pression, la nucléation de vortex limite la vitesse des ions, tandis qu'à haute pression (au-dessus de $\approx 10-12$ bar), l'émission de rotons devient le facteur limitant. De plus, la vitesse critique pour l'émission de rotons diminue avec la pression, alors que la vitesse critique pour la nucléation de vortex augmente. Les cadres théoriques existants, tels que l'équation de Gross-Pitaevskii standard (GPE), échouent à capturer l'émission de rotons car ils possèdent un spectre d'excitation strictement monotone, incapable de reproduire la courbe de dispersion de Landau avec son minimum caractéristique « roton ». À l'inverse, les théories microscopiques peinent souvent à décrire les deux mécanismes dans un cadre unique.

Méthodologie
Les auteurs emploient un modèle de Gross-Pitaevskii généralisé et non local (gnlGPE) pour simuler la dynamique de l'hélium $^4$ superfluide à température nulle ($T \to 0$). Le modèle est défini par un paramètre d'ordre complexe $\psi(x,t)$ et intègre un potentiel d'interaction interatomique non local $V(r)$ ainsi qu'un terme de correction de densité d'ordre supérieur.

1. Calibration du modèle : Pour assurer une précision physique sur une gamme de pressions (de 0 à 24 bar, approchant la pression de solidification d'environ 25 bar), les auteurs n'utilisent pas une forme analytique fixe pour le potentiel d'interaction. À la place, ils inversent la relation de dispersion de Bogoliubov généralisée en utilisant des données expérimentales de Godfrin et al. (2021). Cela permet au modèle de reproduire les courbes de dispersion exactes dépendantes de la pression, incluant les branches phonon, maxon et roton, ainsi que le plateau de Pitaevskii.
2. Configuration de la simulation : L'étude est principalement menée dans une approximation bidimensionnelle (2D). Les auteurs simulent un obstacle circulaire (disque) se déplaçant dans un He II au repos. L'obstacle est modélisé comme un potentiel externe avec un profil gaussien qui épuise complètement la densité superfluide.
3. Détermination de la vitesse critique :
   • Émission de rotons ($v_{rot}^{crit}$) : Déterminée par l'analyse du diagramme de bifurcation des solutions stationnaires. Les auteurs utilisent une méthode de Newton-Raphson couplée à une méthode de gradient conjugué bi-conjugué stabilisée pour trouver le nombre de Mach maximal ($M = U/c$) pour lequel des solutions stationnaires stables existent. Au-delà de ce point, le système devient instable et des rotons sont émis.
   • Nucléation de vortex ($v_{v}^{crit}$) : Déterminée via une approche dynamique. L'obstacle est accéléré jusqu'à atteindre une vitesse cible, et le système est évolué sur une échelle de temps d'environ $10\xi/c$ (où $\xi$ est la longueur de guérison) pour observer si des vortex quantifiés sont nucléés.
4. Paramètres : Les simulations couvrent les pressions $P = 0, 2, 5, 10,$ et $24$ bar. Les tailles d'obstacles varient de $D = 10\xi$ à $D = 1000\xi$. Une attention particulière est portée à la taille des ions négatifs, qui varie avec la pression (d'environ 34,5 Å à 0 bar à environ 21,7 Å à 24 bar).

Résultats clés

• Dépendance en pression des vitesses critiques : Les simulations confirment les tendances opposées observées expérimentalement. À mesure que la pression augmente, le nombre de Mach critique pour l'émission de rotons ($M_{rot}^{crit}$) diminue, passant d'environ 0,247 à 0 bar à environ 0,128 à 24 bar. À l'inverse, le nombre de Mach critique pour la nucléation de vortex ($M_{v}^{crit}$) augmente avec la pression.
• Dépendance en taille : Pour les petits obstacles (comparables aux tailles des ions négatifs, $D < 100$ nm), la vitesse critique est déterminée par le plus bas des deux seuils ($v_{rot}^{crit}$ ou $v_{v}^{crit}$). À basse pression, $v_{v}^{crit}$ est le facteur limitant ; à haute pression, $v_{rot}^{crit}$ est limitant. À mesure que la taille de l'obstacle augmente ($D \ge 1000\xi$), les vitesses critiques pour toutes les pressions s'effondrent sur une courbe unique régie par la nucléation de vortex, similaire aux prédictions de la GPE standard.
• Le « gap » : Une découverte significative est l'élargissement de l'écart entre $v_{rot}^{crit}$ et $v_{v}^{crit}$ à mesure que la pression augmente. À haute pression, il existe une fenêtre de vitesse où des rotons sont émis, mais où la nucléation de vortex ne s'est pas encore produite. Cela suggère que l'émission de rotons peut être isolée et observée expérimentalement à haute pression sans la complication immédiate de la génération de vortex.
• Visualisation : Les visualisations de densité confirment qu'à des nombres de Mach intermédiaires (par exemple, $M=0,28$), des rotons sont émis à toutes les pressions, tandis que la nucléation de vortex ne se produit qu'à des nombres de Mach plus élevés ($M \approx 0,3$), le seuil de nucléation se déplaçant vers des valeurs plus élevées à mesure que la pression augmente.

Portée et affirmations
Les auteurs affirment que ce travail représente la première tentative d'analyser la dépendance en pression de l'émission de rotons et de la nucléation de vortex, ainsi que de leurs vitesses critiques associées, dans un cadre théorique unique. En utilisant la GPE généralisée et non locale, l'étude comble avec succès le fossé entre le critère de Landau (émission de rotons) et la nucléation hydrodynamique de vortex.

L'article souligne que, bien que l'étude soit restreinte à la 2D en raison de la complexité computationnelle des simulations 3D (notamment concernant la détermination de la graine minimale pour la nucléation de vortex en 3D), les résultats 2D fournissent une explication qualitative et semi-quantitative de la physique sous-jacente qui s'aligne avec les observations expérimentales. Les auteurs notent que leurs valeurs numériques pour les vitesses critiques sont comparables aux données expérimentales, validant l'approche consistant à utiliser des relations de dispersion dérivées expérimentalement pour contraindre le potentiel d'interaction. Le travail offre une explication physiquement motivée de la raison pour laquelle le caractère du mouvement des ions change à la pression critique $P_{crit} \approx 10-12$ bar, l'attribuant au croisement des courbes dépendantes de la pression pour $v_{rot}^{crit}$ et $v_{v}^{crit}$.

Les auteurs restent modestes concernant les implications en 3D, reconnaissant que si les simulations 2D expliquent les tendances générales, la topologie spécifique de la nucléation de vortex en 3D (par exemple, boucles ion-vortex contre anneaux) reste une question ouverte nécessitant des investigations futures. Ils notent également que modéliser l'ion comme une sphère dure parfaite est une simplification, bien que les tendances en pression devraient rester robustes.