🔬 mesoscale physics

Quasiparticle-resolved variational theory of Andreev spin qubits

Cet article développe une théorie variationnelle en continuum pour les qubits d'Andreev qui révèle comment les configurations cohérentes de deux quasiparticules dominent les couplages Josephson dépendants du spin et fournit des expressions sous forme fermée pour les paramètres clés, faisant le pont entre les modèles de points quantiques fortement corrélés et les hamiltoniens de circuits effectifs.

Auteurs originaux : Teodor Iličin, Rok Žitko

Publié 2026-10-07

📖 5 min de lecture🧠 Analyse approfondie

CC BY 4.0

Auteurs originaux : Teodor Iličin, Rok Žitko

Article original sous licence CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Ceci est une explication générée par l'IA de l'article ci-dessous. Elle n'a pas été rédigée ni approuvée par les auteurs. Pour une précision technique, consultez l'article original. Lire la clause de non-responsabilité complète

Dans le monde microscopique de l'informatique quantique, l'information est stockée dans l'état le plus fragile qui soit : le spin d'un seul électron. Ce spin agit comme une minuscule aiguille de boussole, pointant soit vers le haut, soit vers le bas pour représenter un zéro ou un un. Pour protéger cette information délicate du bruit de l'environnement, les scientifiques piègent souvent l'électron à l'intérieur d'un circuit supraconducteur, un matériau qui conduit l'électricité avec une résistance nulle. Lorsqu'un électron est confiné dans un tel circuit, il forme un type spécial d'état lié connu sous le nom d'état lié d'Andreev. Ces états sont le cœur d'un nouveau type de bit quantique, ou qubit, appelé qubit de spin d'Andreev. Le défi pour les chercheurs est de comprendre exactement comment cet électron piégé interagit avec la mer de supraconducteurs d'électrons qui l'entoure. Bien que l'électron porte l'information, il n'existe pas de manière isolée ; il échange constamment de l'énergie et des particules virtuelles avec le supraconducteur, créant une danse complexe de mécanique quantique qui détermine comment le qubit se comporte et comment il peut être contrôlé.

Une équipe de physiciens de l'Institut Jožef Stefan en Slovénie a développé une nouvelle façon de cartographier cette interaction avec une haute précision. Ils se sont concentrés sur une configuration spécifique où un point quantique, agissant comme une minuscule île pour les électrons, est connecté à deux broches supraconductrices. Dans ce système, l'électron sur l'île peut interagir avec les broches supraconductrices par deux chemins différents : un chemin direct à travers l'île elle-même, et un chemin de fond qui contourne l'île. Les chercheurs voulaient savoir comment ces interactions créent les différences d'énergie entre les états de spin haut et de spin bas du qubit, une propriété essentielle pour lire et écrire l'information quantique. Pour ce faire, ils ont créé un modèle mathématique détaillé qui suit le comportement de l'électron et des « quasiparticules » — des excitations dans le supraconducteur qui agissent comme des particules — qui l'entourent.

Le cœur de leur travail repose sur une méthode appelée théorie variationnelle, qui est essentiellement une façon de deviner la forme la plus probable de la fonction d'onde quantique, puis d'affiner cette supposition jusqu'à ce qu'elle corresponde le plus étroitement possible aux lois de la physique. Dans ce cas, les chercheurs ont construit un modèle qui inclut la possibilité que le système possède zéro, une ou deux quasiparticules dans les broches supraconductrices à n'importe quel moment donné. Bien que la probabilité de trouver deux quasiparticules soit très faible, l'équipe a découvert que ces configurations rares sont étonnamment importantes. Même si elles apparaissent peu fréquemment, elles interfèrent avec l'état plus commun à zéro quasiparticule d'une manière qui modifie considérablement l'énergie du qubit. Plus précisément, ces configurations à deux quasiparticules contribuent exactement à la moitié de l'effet principal qui fait que les états de spin haut et de spin bas présentent des énergies différentes, même lorsque la répulsion électron-électron est faible. À mesure que la répulsion entre les électrons s'intensifie, ces états à deux quasiparticules deviennent le facteur dominant déterminant le comportement du qubit.

Les chercheurs ont testé leur modèle par rapport à de puissantes simulations informatiques connues sous le nom de calculs de groupe de renormalisation numérique, qui sont considérées comme la référence pour résoudre ce type de problèmes. Leurs résultats concordent très étroitement avec les simulations, confirmant que leur modèle simplifié capture l'essentiel de la physique. Ils ont découvert que la force de l'interaction entre le qubit et le supraconducteur dépend fortement de la répulsion de Coulomb, qui est la force qui repousse les électrons. Lorsque cette répulsion est forte, le courant de super-courant conventionnel qui conduit l'énergie est supprimé beaucoup plus rapidement que l'énergie dépendante du spin. Cela signifie que dans les systèmes à forte répulsion électronique, les effets dépendants du spin deviennent relativement plus importants, même si le couplage global avec le supraconducteur s'affaiblit. Cette découverte aide à expliquer pourquoi certaines configurations expérimentales présentent des signaux dépendants du spin forts malgré des interactions électroniques fortes.

Au-delà du simple calcul des énergies, l'équipe a cartographié l'endroit où les quasiparticules virtuelles résident réellement dans l'espace. Ils ont découvert que ces particules ne sont pas confinées au minuscule point quantique mais sont réparties sur une distance comparable à la longueur de cohérence du supraconducteur, qui est une mesure de la distance jusqu'à laquelle l'ordre supraconducteur s'étend. De plus, ils ont découvert que le spin de l'électron n'est pas entièrement localisé sur le point ; une partie de son spin est partagée avec les broches supraconductrices. Ce partage est influencé par une propriété appelée couplage spin-orbite, qui lie le spin de l'électron à son mouvement. Lorsqu'un champ magnétique est appliqué, ce partage de spin crée un effet unique : un supercourant qui circule même lorsqu'il n'y a pas de différence de tension à travers la jonction. Ce courant est le même pour les états de spin haut et de spin bas, ce qui signifie qu'il ne dépend pas de l'information du qubit mais plutôt de la redistribution physique du spin lui-même.

L'étude offre un pont clair entre la réalité microscopique complexe des électrons en interaction et les modèles effectifs plus simples utilisés pour concevoir des circuits quantiques. En montrant exactement comment la structure de la fonction d'onde détermine les couplages, les chercheurs ont offert un outil aux ingénieurs pour prédire comment ces qubits se comporteront dans des dispositifs réels. Ils ont démontré qu'ignorer les événements rares à deux quasiparticules conduirait à des erreurs significatives dans la prédiction des performances du qubit. Ce travail confirme que pour véritablement comprendre et contrôler ces bits quantiques, il faut tenir compte de l'interaction subtile et cohérente entre l'électron piégé et les excitations virtuelles du supraconducteur, une relation qui est bien plus complexe et influente qu'on ne l'avait supposé auparavant.

Résumé technique : Théorie variationnelle résolue en quasi-particules des qubits de spin d'Andreev

Énoncé du problème
Les qubits de spin d'Andreev (ASQ) encodent l'information quantique dans le spin d'une quasi-particule unique confinée dans un état lié d'Andreev au sein d'un contact ponctuel supraconducteur ou d'une hétérostructure hybride. Bien que les Hamiltoniens phénoménologiques décrivent avec succès les ASQ dans de nombreux régimes, une compréhension microscopique du couplage Josephson dépendant du spin ( $E_S$ ) et du rôle des configurations de quasi-particules (QP) virtuelles reste incomplète. Plus précisément, l'interaction entre les corrélations de Coulomb, le couplage spin-orbite (SOC) et le continuum d'états supraconducteurs nécessite une description qui dépasse les simples modèles effectifs ou la théorie des perturbations d'ordre faible. Les questions ouvertes incluent l'étendue spatiale de l'état lié porteur de spin, l'apparition d'un écrantage de type Kondo partiel, et l'origine microscopique précise de la levée de dégénérescence du spin en présence d'états de continuum.

Méthodologie
Les auteurs développent une description variationnelle en continuum pour le doublet de parité impaire dans une jonction de point quantique modélisée par un modèle d'impureté d'Anderson (AIM) couplé à deux réservoirs supraconducteurs. Le système comprend :

Hamiltonien : Une orbitale impure avec interaction (répulsion de Coulomb $U$ , décalage $\delta$ ) couplée à des réservoirs BCS via un saut (hopping) conservant le spin ( $V$ ) et un terme de tunnel direct entre les réservoirs ( $W_\sigma$ ) qui incorpore des amplitudes dépendantes du spin et le SOC.
Ansatz variationnel : La fonction d'onde est construite comme une superposition de configurations contenant jusqu'à deux quasi-particules de Bogoliubov dans les réservoirs. Cela inclut l'état de vide (0-QP), les états à 1-QP (single-QP), et les états à 2-QP (double-QP), ces derniers étant de plus décomposés en paires de bain singulet et triplet.
Symétrie et canaux : Les auteurs utilisent une symétrie unitaire ( $\eta$ ) et introduisent des canaux de QP symétriques/antisymétriques ( $\gamma_{s,k,\sigma}$ ) pour simplifier les équations variationnelles.
Schémas d'approximation : Les équations intégrales couplées dérivées du principe variationnel sont résolues en utilisant trois niveaux de troncature de noyau (Niveau 1 : 1-QP uniquement ; Niveau 2 : 2-QP sans rétroaction ; Niveau 3 : 2-QP avec rétroaction de noyau mixte).
Validation : Les résultats variationnels sont comparés aux calculs de groupe de renormalisation numérique (NRG) et de groupe de renormalisation de matrice de densité (DMRG) utilisant le solveur open-source qdjj solver.

Contributions clés et résultats

Expressions analytiques pour les couplages Josephson :
Les auteurs dérivent des expressions analytiques en forme fermée pour l'énergie Josephson conventionnelle ( $E_J$ ) et l'énergie Josephson dépendante du spin ( $E_S$ ) à gap supraconducteur $\Delta$ fini.

Couplage conventionnel ( $E_J$ ) : Récupéré au second ordre en hybridation ( $\Gamma^2$ ), ce qui est cohérent avec les résultats de la théorie des perturbations au quatrième ordre pour les points quantiques multi-niveaux.
Couplage dépendant du spin ( $E_S$ ) : Dérivé à l'ordre dominant en amplitude de tunnel direct ( $W$ ) et en hybridation ( $\Gamma$ ). L'expression montre que $E_S$ est supprimé par le facteur $\Delta/(\Delta + U/2)$ par rapport à la limite atomique.

Rôle des configurations à deux quasi-particules :
Une conclusion centrale est que les configurations à 2-QP, malgré des poids de probabilité faibles, sont essentielles pour l'amplitude de la levée de dégénérescence du spin.

Cohérence intersectorielle : Le couplage dépendant du spin provient de la cohérence entre le secteur 0-QP (vide) et le secteur 2-QP, ainsi qu'entre le secteur 1-QP.
Contribution quantitative : À répulsion de Coulomb nulle ( $U=0$ ), l'interférence 0–2-QP contribue exactement à la moitié du couplage dépendant du spin dominant. À mesure que $U$ augmente, ce terme d'interférence domine la contribution 1-QP, fournissant le comportement dominant en $1/U$ du couplage.

Précision et convergence :

Comparaison (Benchmarking) : L'approximation variationnelle de Niveau 3 (incluant la rétroaction 2-QP) montre un excellent accord avec NRG et DMRG pour les énergies et les composantes de la fonction d'onde, les écarts provenant principalement de l'omission des secteurs 3-QP et supérieurs.
Erreurs de troncature : L'erreur introduite par la troncature du noyau au sein de l'espace 2-QP est jugée mineure (moins de 0,255 % pour l'énergie à $\phi=0$ ). Les écarts plus importants par rapport à NRG lors d'une forte hybridation sont attribués à l'exclusion des secteurs à plus hautes énergies de QP plutôt qu'à l'approximation du noyau elle-même.

Redistribution du spin et réponse magnétique :

Transfert induit par le SOC : L'inclusion du SOC induit un transfert de spin entre le point quantique et les réservoirs. Cette redistribution est capturée par le changement de la fraction de moment local $q_d$ .
Courant Josephson commun : Lorsque l'énergie Zeeman sur le point quantique diffère de celle dans les réservoirs ( $b_{rel} \neq 0$ ), le transfert de spin induit par le SOC génère une contribution de courant Josephson commune aux deux états de qubit. Ce courant est fini même à biais de phase nul, représentant un effet Josephson anomal ( $\phi_0$ -junction) piloté par l'interaction entre le SOC et les gradients de champ magnétique.

Structure spatiale :
Les auteurs cartographient la localisation des quasi-particules virtuelles et la densité de spin physique dans les réservoirs.

Localisation : Les quasi-particules virtuelles s'étendent sur la longueur de cohérence de BCS $\xi_0$ . La composante 2-QP est dominée par une seule orbitale naturelle par canal, avec des probabilités d'occupation dépassant 95 %.
Densité de spin : La densité de spin physique dans les réservoirs est positive, mais la fonction de corrélation spin-point quantique change de signe avec la distance. Le terme d'interférence 0–2-QP domine la redistribution spatiale du spin sans altérer la valeur intégrée du spin.

Signification
Ce travail établit un lien microscopique rigoureux entre les modèles de points quantiques interagissants et les Hamiltoniens effectifs $E_J$ – $E_S$ utilisés pour concevoir les circuits ASQ. Il clarifie que la structure de la fonction d'onde, spécifiquement la cohérence entre les différents secteurs de quasi-particules, dicte les couplages. Les résultats fournissent des contraintes quantitatives pour la conception de circuits, démontrant que la répulsion de Coulomb supprime les couplages conventionnels et dépendants du spin différemment ( $E_J \propto U^{-2}$ contre $E_S \propto U^{-1}$ dans la limite de grand $U$ ). Cela implique que si la dépendance relative au spin ( $E_S/E_J$ ) peut augmenter avec l'interaction, l'amplitude absolue du couplage dépendant du spin diminue. Le travail identifie également le mécanisme spécifique (cohérence 0–2-QP) responsable de la levée de dégénérescence du spin, offrant un cadre prédictif plus précis que les approches phénoménologiques ou de perturbation d'ordre faible précédentes.

Articles similaires