2025-11-24T06:13:17.730458

Nonequilibrium Statistics of Biased Kondo Resonance

Han
Numerical renormalization group (NRG) is formulated for nonequilibrium steady-state by converting finite-lattice many-body eigenstates into scattering states. Extension of the full-density-matrix NRG for a biased Anderson impurity model, simplified by formulating with the original orbital basis as the Hamiltonian, enables detailed studies of the sub-Kondo spectral evolution in the zero-temperature limit, confirming the double-resonance structure at bias of the Kondo energy scale $T_K$. The distribution shows distinct multi-scale spectral features at energy $ω$ below the Kondo scale ($ω\lesssim T_K$) and near the bias ($ω\gtrsim V$), leading to the nonequilibrium temperature $T_{\rm loc}$ local to the Kondo dot scaling as $k_BT_{\rm loc}\approx V$ for $V\gg T_K$. The current-voltage relation in the low-temperature limit ($T\ll T_K$) deviates from the unitary limit as the bias exceeds the Kondo scale ($V/2\gtrsim T_K$) and reaches the current saturation regime.
academic

Statistiques Hors Équilibre de la Résonance Kondo Biaisée

Informations Fondamentales

  • ID de l'article: 2503.14400
  • Titre: Nonequilibrium Statistics of Biased Kondo Resonance
  • Auteur: Jong E. Han (State University of New York at Buffalo)
  • Classification: cond-mat.str-el (Systèmes électroniques fortement corrélés)
  • Date de publication: 14 octobre 2025
  • Lien de l'article: https://arxiv.org/abs/2503.14400

Résumé

Cet article formule une méthode de groupe de renormalisation numérique (NRG) pour les états stationnaires hors équilibre en convertissant les états propres à réseau fini en états de diffusion. En utilisant la base orbitale primitive comme formulation du hamiltonien, l'auteur étend la NRG de matrice densité complète au modèle d'impureté Anderson biaisé, permettant une étude détaillée de l'évolution spectrale sous-Kondo à température zéro. Cela confirme la structure de double résonance sous le biais à l'échelle Kondo TKT_K. La distribution révèle des caractéristiques spectrales multi-échelles uniques sous l'échelle Kondo (ωTK\omega \lesssim T_K) et près du biais (ωV\omega \gtrsim V), conduisant à une température hors équilibre locale TlocT_{\rm loc} du point quantique Kondo qui s'échelonne comme kBTlocVk_BT_{\rm loc} \approx V lorsque VTKV \gg T_K. À la limite basse température (TTKT \ll T_K), la relation courant-tension s'écarte de la limite unitaire et atteint une région de saturation du courant lorsque le biais dépasse l'échelle Kondo (V/2TKV/2 \gtrsim T_K).

Contexte de Recherche et Motivation

Description du Problème

  1. Défi fondamental: La résolution des problèmes de mécanique statistique hors équilibre dans les systèmes quantiques biaisés a représenté un défi majeur au cours des dernières décennies, particulièrement pour les modèles d'impureté quantique fortement corrélés.
  2. Importance:
    • Le problème Kondo sert de banc d'essai pour les techniques hors équilibre de systèmes fortement corrélés, essentiel pour comprendre le comportement hors équilibre des systèmes quantiques à plusieurs corps
    • Les modèles d'impureté quantique jouent un rôle clé dans la résolution de la limite de réseau de matière condensée via la théorie du champ moyen dynamique
  3. Limitations des méthodes existantes:
    • La rigueur théorique et l'efficacité numérique des solveurs d'impureté hors équilibre n'ont pas atteint un niveau satisfaisant
    • Les difficultés de calcul limitent la compréhension de la façon dont le singulet Kondo se décompose en régions sous-Kondo bimodales sous biais
    • La compréhension statistique hors équilibre du spectre est insuffisante par rapport aux caractéristiques spectrales du modèle
  4. Motivation de la recherche:
    • Construire une extension hors équilibre d'état stationnaire accessible basée sur la méthode NRG mature
    • Fournir une compréhension claire de l'évolution de la résonance Kondo sous biais
    • Développer une méthode générale applicable aux interactions d'impureté complexes

Contributions Principales

  1. Innovation méthodologique: Propose une nouvelle méthode NRG hors équilibre qui traite les problèmes hors équilibre d'état stationnaire en convertissant les états propres à réseau fini en états de diffusion
  2. Simplification technique: Utilise la formulation hamiltonienne en base orbitale primitive, simplifiant la mise en œuvre et la rendant plus adaptée à l'extension vers des problèmes généralisés
  3. Découvertes physiques:
    • Confirme le comportement de division du pic Kondo à V/2TKV/2 \approx T_K
    • Découvre les caractéristiques de distribution multi-échelles et le phénomène d'inversion de population
    • Révèle la relation d'échelle de la température hors équilibre TlocVT_{\rm loc} \approx V
  4. Propriétés de transport: Vérifie le phénomène de saturation du courant à basse température dans la région TKV/23TKT_K \lesssim V/2 \lesssim 3T_K

Détails de la Méthode

Cadre Théorique

Modèle d'Impureté Anderson

Le hamiltonien est défini comme: H^=σϵdndσ+U(nd12)(nd12)+αkσϵkcαkσcαkσtΩαkσ(dσcαkσ+cαkσdσ)\hat{H} = \sum_\sigma \epsilon_d n_{d\sigma} + U\left(n_{d\uparrow} - \frac{1}{2}\right)\left(n_{d\downarrow} - \frac{1}{2}\right) + \sum_{\alpha k\sigma} \epsilon_k c^\dagger_{\alpha k\sigma} c_{\alpha k\sigma} - \frac{t}{\sqrt{\Omega}} \sum_{\alpha k\sigma} (d^\dagger_\sigma c_{\alpha k\sigma} + c^\dagger_{\alpha k\sigma} d_\sigma)

α=L,R\alpha = L,R désigne les réservoirs source et drain, UU est le paramètre de répulsion Coulombienne, et ϵd\epsilon_d est le niveau d'impureté.

Matrice Densité Hors Équilibre

La matrice densité hors équilibre d'état stationnaire est construite comme: ρ^neq=1Zneqeβ(H^VY^)\hat{\rho}_{\rm neq} = \frac{1}{Z_{\rm neq}} e^{-\beta(\hat{H} - V\hat{Y})}

où l'opérateur Y de Hershfield est défini comme: Y^=12kσ(ψLkσψLkσψRkσψRkσ)\hat{Y} = \frac{1}{2}\sum_{k\sigma}(\psi^\dagger_{Lk\sigma}\psi_{Lk\sigma} - \psi^\dagger_{Rk\sigma}\psi_{Rk\sigma})

Construction des États de Diffusion

L'innovation clé réside dans la construction de l'opérateur Y via la relation: Y^=iηL+iηY^0\hat{Y} = \frac{i\eta}{-\mathcal{L} + i\eta}\hat{Y}_0

L\mathcal{L} est l'opérateur Liouvillien et Y^0\hat{Y}_0 est l'opérateur Y initial avec l'impureté découplée. Dans la représentation des états propres énergétiques: αY^β=iηEα+Eβ+iηαY^0β\langle\alpha|\hat{Y}|\beta\rangle = \frac{i\eta}{-E_\alpha + E_\beta + i\eta}\langle\alpha|\hat{Y}_0|\beta\rangle

Implémentation Numérique

Représentation en Chaîne de Wilson

  • Utilise la chaîne de Wilson (WC) pour simuler les réservoirs, longueur de chaîne NWC=100N_{\rm WC} = 100, paramètre de renormalisation Λ=1.2\Lambda = 1.2
  • Traitement des orbitales d'environnement: considère les changements d'énergie libre des orbitales d'environnement sous biais
  • Le poids de la matrice densité atteint un pic à la longueur de chaîne de Wilson n=Mn = M, satisfaisant DMV/2D_M \approx V/2

NRG de Matrice Densité Complète (FDM-NRG)

  • Nombre d'états conservés Nkept=2000N_{\rm kept} = 2000
  • Moyenne z utilisant 8 valeurs z pour réduire la discrétisation de la chaîne de Wilson
  • Dimension de l'espace de Hilbert à chaque itération NRG: 16×Nkept16 \times N_{\rm kept}

Configuration Expérimentale

Paramètres du Modèle

  • Région Kondo forte: U/Γ=10U/\Gamma = 10, Γ=0.1\Gamma = 0.1, TK=0.0044T_K = 0.0044
  • Symétrie particule-trou: ϵd=0\epsilon_d = 0
  • Asymétrie particule-trou: ϵd=U/4\epsilon_d = U/4
  • Température: T=0.1TKT = 0.1T_K (limite basse température)

Paramètres de Calcul

  • Paramètre de renormalisation: Λ=4\Lambda = 4
  • Paramètre d'élargissement spectral: δ=0.2×max(TK,T,V/2)\delta = 0.2 \times \max(T_K, T, V/2)
  • Paramètre infinitésimal: η=0.1×max(TK,T,V/2)\eta = 0.1 \times \max(T_K, T, V/2)

Métriques d'Évaluation

  • Fonction spectrale: ρΣ(ω)\rho_\Sigma(\omega) (calculée via technique d'auto-énergie)
  • Fonction de distribution: f(ω)=Gd<(ω)/[2πiρ(ω)]f(\omega) = G^<_d(\omega)/[2\pi i\rho(\omega)]
  • Température hors équilibre: Tloc2=6π2ω[f(ω)Θ(ω)]dωT_{\rm loc}^2 = \frac{6}{\pi^2}\int_{-\infty}^{\infty} \omega[f(\omega) - \Theta(-\omega)]d\omega
  • Courant: calculé via technique graphique

Résultats Expérimentaux

Analyse de l'Évolution Spectrale

Vérification Non-Interagissante

  • Valide la méthode dans le modèle de niveau de résonance non-interagissant
  • À petit biais, le courant suit I=G0VI = G_0 V, où G0=1/(2π)G_0 = 1/(2\pi) est la conductance quantique unitaire
  • À grand biais, conforme à la relation exacte I=(Γ/π)tan1(V/2Γ)I = (\Gamma/\pi)\tan^{-1}(V/2\Gamma)

Division de la Résonance Kondo

  1. Biais zéro: La fonction spectrale satisfait la limite unitaire à ρΣ(0)1/(πΓ)\rho_\Sigma(0) \approx 1/(\pi\Gamma)
  2. Biais critique: Lorsque V/2TKV/2 \approx T_K, la résonance Kondo commence à afficher une division de pic
  3. Haut biais: Les positions de pic s'alignent à ±V/2\pm V/2, mais la largeur de ligne est bien supérieure à TKT_K, à considérer comme faiblement corrélée

Caractéristiques Statistiques Hors Équilibre

Évolution de la Fonction de Distribution

  • Bas biais: La fonction de distribution est proche de la fonction de Fermi-Dirac
  • Biais modéré: Développe des pics larges et importants à ωV|ω| \gtrsim V
  • Inversion de population: Apparition de f(0)<0f'(0) < 0 entre les pics Kondo divisés (V/2ωV/2-V/2 \lesssim \omega \lesssim V/2)

Échelle de Température Hors Équilibre

Découverte clé: à la limite de haut biais TlocVpourV/2TKT_{\rm loc} \approx V \quad \text{pour} \quad V/2 \gg T_K

Cela indique que la température électronique locale est proportionnelle au biais.

Propriétés de Transport

Relation Courant-Tension

  1. Limite basse température (T=0.1TKT = 0.1T_K):
    • À petit biais, satisfait la limite unitaire G0=(2π)1G_0 = (2\pi)^{-1}
    • Commence à s'écarter de la limite unitaire lorsque V/2TKV/2 \gtrsim T_K
    • Saturation du courant dans la région TKV/23TKT_K \lesssim V/2 \lesssim 3T_K
  2. Haute température (T=TKT = T_K):
    • Aucun phénomène de saturation du courant
    • Prouve que la saturation du courant provient des effets fortement corrélés

Travaux Connexes

Méthodes d'Impureté Quantique Hors Équilibre

L'article examine systématiquement diverses techniques hors équilibre pour traiter le problème Kondo:

  • Théories analytiques: Méthode de Bethe ansatz
  • Méthodes graphiques: Formule de Landauer, théorie des perturbations
  • Méthodes numériques: Simulation graphique en temps réel, simulation en temps imaginaire, méthodes de réseau tensoriel
  • Groupe de renormalisation: RG en temps réel, RG perturbatif, RG fonctionnel

Comparaison avec les Méthodes NRG Existantes

  • Méthode SNRG: NRG d'état stationnaire proposée par Anders utilisant des états de diffusion non-interagissants pour construire Y^0\hat{Y}_0
  • Méthode présente: Utilise la base du hamiltonien primitif, implémentation plus directe, adaptée à l'extension vers des interactions d'impureté complexes

Conclusions et Discussion

Conclusions Principales

  1. Validation de la méthode: Combine avec succès la théorie des états de diffusion de Hershfield avec la méthode NRG, fournissant un outil numérique fiable pour les systèmes fortement corrélés hors équilibre
  2. Image physique: Confirme la division du pic Kondo à V/2TKV/2 \approx T_K, révélant la structure multi-échelles complexe de la région sous-Kondo
  3. Caractéristiques hors équilibre: Découvre l'inversion de population et la relation d'échelle linéaire de la température hors équilibre, fournissant de nouvelles perspectives pour comprendre les états fortement corrélés hors équilibre
  4. Phénomènes de transport: Vérifie le phénomène de saturation du courant, confirmant les prédictions théoriques

Limitations

  1. Coût de calcul: En raison du traitement de deux chaînes de Wilson, la croissance de l'espace de Hilbert est rapide, limitant le nombre d'états conservés
  2. Dépendance paramétrique: La méthode est sensible au choix du paramètre infinitésimal η\eta, nécessitant un ajustement minutieux
  3. Limitation de température: Principalement applicable à la région basse température TTKT \ll T_K

Directions Futures

  1. Applications étendues: Applicable aux modèles d'interactions d'impureté complexes
  2. Théorie de réseau: Application aux modèles de réseau hors équilibre via la théorie du champ moyen dynamique
  3. Comparaison expérimentale: Comparaison détaillée avec les expériences de transport de points quantiques

Évaluation Approfondie

Points Forts

  1. Rigueur théorique: Basée sur le cadre théorique rigoureux de Hershfield, dérivations mathématiques claires
  2. Innovation méthodologique: Conversion astucieuse des états propres à réseau fini en états de diffusion, évitant la complexité de la construction explicite d'états de diffusion
  3. Fiabilité numérique: Vérifiée via la limite non-interagissante, assurant la correction de la méthode
  4. Perspectives physiques: Fournit une compréhension approfondie de l'évolution hors équilibre de la résonance Kondo

Insuffisances

  1. Complexité de calcul: Coût de calcul significativement augmenté par rapport à la NRG d'équilibre
  2. Ajustement paramétrique: Le choix du paramètre η\eta nécessite l'expérience, pouvant affecter l'universalité des résultats
  3. Vérification expérimentale: Manque de comparaison directe avec les expériences réelles de points quantiques

Impact

  1. Contribution méthodologique: Fournit un nouvel outil numérique pour les systèmes fortement corrélés hors équilibre
  2. Valeur théorique: Approfondit la compréhension du comportement hors équilibre de l'effet Kondo
  3. Perspectives d'application: Extensible aux modèles d'impureté plus complexes et aux systèmes de réseau

Scénarios d'Application

  • Analyse théorique des expériences de transport de points quantiques
  • Recherche hors équilibre de systèmes électroniques fortement corrélés
  • Solveurs d'impureté pour la théorie du champ moyen dynamique
  • Modélisation théorique de l'électronique moléculaire

Références

L'article cite les travaux importants du domaine, notamment:

  • Travaux originaux de Kondo et manuel classique de Hewson
  • Théorie d'état stationnaire hors équilibre de Hershfield
  • Méthode de groupe de renormalisation numérique de Wilson
  • Méthode NRG d'état stationnaire d'Anders
  • Diverses théories d'impureté quantique hors équilibre et méthodes numériques

Cet article apporte une contribution importante au domaine de la théorie quantique à plusieurs corps hors équilibre, fournissant une méthode numérique à la fois rigoureuse et pratique pour étudier les systèmes fortement corrélés hors équilibre, avec une valeur particulière pour comprendre le comportement hors équilibre de l'effet Kondo.