Ferroelectric quantum critical point in superconducting hydrides: The case of H$_3$S

En employant la dynamique moléculaire par intégrale de chemin avec un potentiel appris par apprentissage automatique, cette étude révèle que la supraconductivité à haute température dans le H$_3$S se produit dans une région paraélectrique dominée par de larges fluctuations quantiques nucléaires au-dessus d'un point critique quantique ferroélectrique à environ 134 GPa, lequel appartient à la classe d'universalité d'Ising 4D.

L'essentiel

Imaginez que vous avez une minuscule bille ultra-dense composée d'atomes de soufre et d'hydrogène. Sous une pression extrême, cette bille devient un supraconducteur — un matériau qui conduit l'électricité avec une résistance nulle. Pendant longtemps, les scientifiques ont été perplexes quant à la raison pour laquelle cela se produit dans un matériau spécifique appelé H3S. Ils savaient que cela fonctionnait mieux à une certaine pression (environ 155 GPa), mais la carte de la façon dont les atomes se comportent manquait à l'appel.

Ce papier est comme le dessin de cette carte manquante. Les chercheurs ont utilisé une simulation informatique puissante pour suivre la danse des atomes, non pas comme des billes solides, mais comme des « nuages flous » de probabilité (un effet quantique). Voici ce qu'ils ont découvert, expliqué simplement :

1. Les atomes « flous » et la pression magique

Dans le monde des atomes minuscules, les choses ne sont pas solides ; elles oscillent et tressaillent. Les chercheurs ont découvert qu'à une pression spécifique (environ 134 GPa), les atomes d'hydrogène de H3S atteignent un « point de bascule ».

• L'analogie : Imaginez une balle posée dans un bol. Si le bol est profond, la balle reste au milieu. Si vous secouez le bol (chaleur) ou si vous le pressez (pression), la balle peut commencer à rouler.
• La découverte : À ce point de bascule, appelé Point Critique Quantique (PCQ), les atomes sont dans un état de confusion maximale. Ils ne sont pas installés en un point précis, mais ils ne sont pas non plus totalement aléatoires. Ils « fluctuent » sauvagement, comme une foule de gens essayant de décider dans quelle direction se tourner.

2. Le changement de phase : de « symétrique » à « asymétrique »

Le matériau peut exister sous deux formes principales (phases) :

• La phase « parfaitement équilibrée » (Paraélectrique) : Les atomes d'hydrogène sont situés pile au milieu entre les atomes de soufre. C'est symétrique, comme une balançoire parfaitement équilibrée.
• La phase « asymétrique » (Ferroélectrique) : Les atomes d'hydrogène sont poussés d'un côté. La balançoire bascule.

L'article montre que la transition de l'état « équilibré » à l'état « asymétrique » ne se produit pas exactement là où la supraconductivité est la plus forte. Au contraire, le pic de supraconductivité se trouve dans la zone « équilibrée », mais juste à côté du point de bascule où les atomes s'agitent le plus.

3. Le « point idéal » de la supraconductivité

Voici la grande surprise :

• L'ancienne idée : Les scientifiques pensaient que le pic de supraconductivité se produisait parce que le matériau passait de l'état équilibré à l'état asymétrique.
• La nouvelle découverte : L'article montre que le pic se produit en réalité dans la zone équilibrée, mais juste à côté du chaos.
• L'analogie : Pensez à un surfeur. Les meilleures vagues ne sont ni l'eau calme et plate, ni le ressac chaotique et déferlant. Les meilleures vagues sont situées précisément là où l'océan commence à s'agiter. La « rugosité » (les fluctuations quantiques) aide les électrons à s'apparier et à circuler sans résistance. L'article suggère que l'agitation sauvage des atomes d'hydrogène près du point de bascule agit comme un catalyseur pour la supraconductivité.

4. Le code de règles « Ising 4D »

Les chercheurs ont analysé la mathématique derrière ce point de bascule et ont découvert qu'il suit un code de règles très spécifique connu sous le nom de classe d'universalité d'Ising 4D.

• L'analogie : Imaginez différents jeux (comme les échecs, les dames ou le Go). Même s'ils se ressemblent différemment, ils peuvent tous suivre la même logique sous-jacente sur la façon dont les pièces se déplacent. Les chercheurs ont découvert que la façon dont ces atomes se comportent suit la même « logique » qu'un modèle mathématique complexe utilisé pour décrire comment les choses changent d'état en quatre dimensions. Cela confirme que leur découverte est une loi fondamentale de la physique, et non un simple hasard.

5. Pourquoi cela compte pour la « carte »

Avant cette étude, la carte de H3S était floue. Les chercheurs ont utilisé un nouveau type de « cerveau d'IA » (un potentiel d'apprentissage automatique) pour exécuter des simulations trop coûteuses pour les anciennes méthodes.

• Ils ont découvert que si l'on ignore le caractère « flou » quantique des atomes (en les traitant comme des billes de billard solides), on obtient une carte erronée. Les oscillations quantiques déplacent la pression de transition d'une quantité énorme (environ 50 GPa).
• En incluant ces oscillations, ils ont enfin localisé le « point de bascule » (PCQ) et montré que le pic de supraconductivité se situe dans une région de fortes fluctuations quantiques, juste au-dessus du point de bascule.

Résumé

L'article révèle que la magie de la supraconductivité dans H3S n'est pas causée par le simple changement de forme du matériau. Elle se produit parce que le matériau plane juste à côté d'un « point de bascule quantique » où les atomes vibrent sauvagement. Ces vibrations sauvages agissent comme un catalyseur, aidant l'électricité à circuler parfaitement. Les chercheurs ont désormais cartographié exactement où cela se produit et ont prouvé que cela suit une règle mathématique spécifique et universelle.

Résumé technique

Énoncé du Problème
L'hydrure de soufre (H$_3$S) est un supraconducteur à haute température critique ($T_c$) proéminent, présentant une $T_c$ de 203 K à environ 155 GPa. Bien que le dôme supraconducteur présente un maximum similaire à celui des cuprates et du SrTiO$_3$ dopé, la nature précise du diagramme de phase, en particulier la relation entre le maximum supraconducteur et la transition de phase structurelle, reste non résolue. Des études théoriques antérieures suggéraient une transition de phase ferroélectrique de la phase cubique centrée sur les corps $Im\bar{3}m$ vers la phase trigonale $R3m$, pilotée par le déplacement de l'hydrogène. Cependant, ces théories plaçaient la pression de transition à un niveau nettement inférieur (de plus de 40 GPa) au pic de $T_c$ expérimental. De plus, une compréhension globale du diagramme de phase de l'H$_3$S sur une plage étendue de températures et de pressions manquait jusqu'à présent.

Méthodologie
Pour combler ces lacunes, les auteurs ont employé la Dynamique Moléculaire par Intégrale de Chemin (PIMD) afin de traiter exactement les effets quantiques nucléaires (NQE). Étant donné le coût computationnel prohibitif de la combinaison de la PIMD avec les méthodes de structure électronique ab initio (telles que la Théorie de la Fonctionnelle de la Densité, DFT) pour les tailles de système et les températures requises, les auteurs ont utilisé un Potentiel Interatomique par Apprentissage Automatique (MLIP). Plus précisément, ils ont entraîné un potentiel MACE (Réseau de Passage de Messages Équivariant d'Ordre Supérieur) sur un ensemble de données de configurations DFT-BLYP couvrant des pressions de 40 à 220 GPa à 100 K et 200 K.

Les simulations ont été réalisées à l'aide du package i-PI avec des ensembles NVT pour des tailles de système de cellules primitives $L \times L \times L$ ($L \in \{2, 3, 4\}$) afin de permettre une analyse de mise à l'échelle de taille finie (finite-size scaling). L'étude a couvert des températures allant de 50 K à 300 K. Les observables clés comprenaient :

• Paramètres d'Ordre : Déplacement global des protons ($\Delta$) et sa valeur absolue ($\Delta_{abs}$) pour distinguer les phases paraélectrique et ferroélectrique.
• Observables Locaux : Distribution des déplacements locaux des protons et des moments dipolaires locaux.
• Propriétés Dynamiques : Fonctions de Green de phonons dans le temps imaginaire ($G_{ij}(\tau)$) pour extraire les fréquences de phonons anharmoniques et étudier les effets de retard.
• Analyse de Mise à l'Échelle : Mise à l'échelle de taille finie du paramètre d'ordre près du point critique pour déterminer la classe d'universalité.

Résultats Clés

1. Diagramme de Phase et Point Critique Quantique (QCP) : L'étude révèle une ligne de transition ferroélectrique qui se déplace vers des pressions plus basses lorsque la température diminue. En extrapolant à la limite thermodynamique et à $T=0$ K, les auteurs identifient un Point Critique Quantique (QCP) ferroélectrique à $p_{QCP} \approx 134 \pm 2$ GPa.
2. Classe d'Universalité : L'analyse de mise à l'échelle de taille finie du paramètre d'ordre près du QCP produit un exposant critique $\nu \approx 0,483 \pm 0,014$. Cette valeur est cohérente avec la classe d'universalité d'Ising 4D (où $\nu = 1/2$), en supposant l'invariance de Lorentz ($z=1$).
3. Discrépance avec la $T_c$ Expérimentale : La ligne de transition ferroélectrique calculée (par exemple, à 124 GPa pour 200 K) est nettement inférieure au pic supraconducteur expérimental de 155 GPa. Par conséquent, le maximum de $T_c$ expérimental se situe dans la phase paraélectrique $Im\bar{3}m$, et non dans la phase ferroélectrique.
4. Rôle des Fluctuations Quantiques : La région entourant le QCP et s'étendant dans la phase paraélectrique est caractérisée par de fortes fluctuations quantiques nucléaires. Les auteurs quantifient cela en utilisant le rapport entre la variance de la fonction de Green locale au temps imaginaire $\beta/2$ et celle à $\tau=0$ ($\sigma^2_g(\beta/2)/\sigma^2_g(0)$). Ce rapport culmine dans la phase paraélectrique près du QCP, indiquant de forts effets de retard et des moments dipolaires fluctuants, qui sont supprimés dans la phase ferroélectrique où les protons sont gelés.
5. Comparaison avec la Dynamique Classique : Les simulations de Dynamique Moléculaire (MD) classique montrent que la transition ferroélectrique se produit à des pressions beaucoup plus élevées (décalées d'environ 50 GPa à 200 K) par rapport au cas quantique, soulignant le rôle critique des NQE pour stabiliser la phase paraélectrique et abaisser la pression de transition.
6. Adoucissement des Phonons : Les modes de phonons optiques pilotant la transition présentent un adoucissement à mesure que le système approche de la transition du côté paraélectrique, suivi d'un saut de fréquence brusque lors de l'entrée dans la phase ferroélectrique, ce qui est cohérent avec une transition de phase de second ordre.

Signification et Revendications
Cet article établit que la supraconductivité à haute $T_c$ dans l'H$_3$S ne se produit pas à la transition de phase ferroélectrique elle-même, mais plutôt dans la région paraélectrique dominée par de fortes fluctuations quantiques et des effets de retard. Les auteurs suggèrent que ces fluctuations, potentiellement médiées par la proximité du QCP ferroélectrique, pourraient améliorer la supraconductivité. Ce scénario implique que la théorie standard de Migdal-Eliashberg, qui néglige souvent la dépendance complète en fréquence et les corrections de vertex, pourrait être insuffisante pour décrire l'H$_3$S. La présence d'un QCP ferroélectrique et l'adoucissement marqué des phonons optiques soutiennent la possibilité de mécanismes d'appariement d'électrons au-delà de l'ordre linéaire, similaires à ceux proposés pour les métaux dilués comme le SrTiO$_3$. Ce travail fournit un diagramme de phase détaillé et de précision quantique qui réconcilie la localisation du dôme supraconducteur avec les instabilités structurelles sous-jacentes, soulignant la nécessité d'inclure les effets quantiques nucléaires et des évaluations de couplage électron-phonon non perturbatives dans les futurs traitements théoriques des hydrures supraconducteurs.