Strong-disorder expansion of the root-averaged density of states for the Anderson model on the Bethe lattice

Cet article démontre que pour le modèle d'Anderson sur le réseau de Bethe dans le régime de fort désordre avec des distributions de site unique à support compact et localement analytiques, la densité d'états moyennée à la racine est absolument continue et admet un développement réel analytique d'ordre fini dont tous les coefficients impairs s'annulent et dont les termes d'ordre supérieur sont déterminés par les marches fermées courtes sur l'arbre.

L'essentiel

Imaginez que vous vous tenez au milieu d'une forêt infinie et parfaitement symétrique. Chaque arbre de cette forêt a exactement le même nombre de voisins (disons $q+1$). Il s'agit du réseau de Bethe, une forme mathématique qui ressemble à un arbre mais qui s'étend à l'infini sans aucun cycle.

Maintenant, imaginez que chaque arbre de cette forêt possède un « poids » caché et aléatoire qui lui est attaché. Certains sont lourds, d'autres légers, et les poids sont choisis aléatoirement selon une règle spécifique. Il s'agit du modèle d'Anderson.

Les physiciens et les mathématiciens veulent savoir : « Si j'envoie une onde d'énergie à travers cette forêt, comment se propage-t-elle ? À quoi ressemble la « densité » de ces ondes d'énergie ? » C'est ce qu'on appelle la densité d'états.

Habituellement, le calcul de cette grandeur est incroyablement difficile car l'aléatoire des poids fait rebondir les ondes de manière chaotique et imprévisible. Cependant, cet article se concentre sur un scénario spécifique : le désordre fort. Cela signifie que les poids aléatoires sur les arbres sont si lourds et variés qu'ils dominent le système. Le « saut » entre les arbres (la connexion) devient une perturbation minuscule, presque négligeable, par rapport aux masses énormes des poids.

Voici une décomposition simple de ce que l'auteur, Masahiro Kaminaga, a découvert :

1. La vue « zoomée »

Comme le désordre est si fort, l'auteur suggère de « zoomer » ou de rééchelonner notre point de vue. Au lieu d'examiner les nombres d'énergie bruts, nous les examinons par rapport à la force du désordre ($\lambda$). C'est comme regarder une chaîne de montagnes à travers une lunette astronomique ; les rochers individuels (les poids aléatoires) deviennent le trait principal, et les petits sentiers qui les relient (les connexions de l'arbre) deviennent des détails secondaires.

2. La magie de la forme « arbre »

La forêt n'a pas n'importe quelle forme ; c'est un arbre. Dans un arbre, si vous partez de la racine et marchez un certain nombre de pas, vous ne pouvez revenir au départ qu'en faisant un nombre pair de pas. Si vous faites un nombre impair de pas, vous êtes garanti d'être ailleurs.

L'auteur utilise ce fait simple pour prouver quelque chose de surprenant : Toutes les corrections « impaires » à la densité d'énergie s'annulent.

• Considérez le calcul comme une recette. Vous avez un ingrédient principal (les poids aléatoires).
• Vous ajoutez des ingrédients « correctifs » pour tenir compte des connexions de l'arbre.
• L'auteur prouve que les ingrédients correctifs 1er, 3e, 5e, etc., sont exactement nuls. Vous n'avez besoin de vous soucier que des 2e, 4e, 6e, etc.

3. L'analogie de la « marche »

Pour déterminer exactement à quoi ressemble la densité d'énergie, l'auteur imagine un « marcheur aléatoire » se déplaçant à travers la forêt.

• Le marcheur commence à la racine, fait quelques pas et doit revenir à la racine.
• L'auteur calcule combien de façons différentes le marcheur peut le faire et à quelle fréquence il visite des arbres spécifiques.
• Parce que la forêt est un arbre, ces « marches » sont très structurées. Elles ne restent pas coincées dans des cycles (car il n'y a pas de cycles).
• La formule finale de la densité d'énergie est une somme de ces motifs de marche spécifiques.

4. Le résultat : une courbe lisse et prévisible

Même si les poids sont aléatoires, l'auteur prouve que si vous regardez la densité d'énergie « moyenne » sur une plage spécifique, elle est lisse et prévisible.

• Le terme dominant : La partie la plus importante de la réponse est simplement la distribution des poids aléatoires eux-mêmes. Si les poids sont uniformément distribués (comme une ligne plate), la densité d'énergie commence par une ligne plate.
• Les corrections : Les connexions de l'arbre ajoutent de petites ondulations à cette ligne. L'auteur fournit une formule précise pour ces ondulations.
  • La première ondulation (la correction d'ordre 2) dépend du nombre de voisins de chaque arbre ($q$) et de la forme de la distribution des poids aléatoires.
  • L'auteur calcule explicitement cette première ondulation pour le cas où les poids sont uniformément distribués.

5. Pourquoi cela compte (selon l'article)

Avant cet article, nous savions que la densité d'énergie existait, mais nous n'avions pas de recette précise, étape par étape, pour la calculer dans le cas d'un désordre fort.

• L'article fournit un développement d'ordre fini. Cela signifie que vous pouvez calculer la réponse avec la précision souhaitée en ajoutant plus de termes à la recette.
• Il prouve que la réponse est analytique, ce qui signifie qu'il s'agit d'une courbe très lisse, sans rupture brutale ni arêtes vives dans la région étudiée.
• Il relie directement les mathématiques complexes des « marches aléatoires sur les arbres » à la propriété physique de « la répartition de l'énergie ».

Analogie récapitulative

Imaginez que vous essayez de prédire la taille moyenne d'une foule de personnes debout sur un sol bosselé et irrégulier (les poids aléatoires).

• Ancienne méthode : Vous essayez de mesurer chaque personne et chaque bosse, ce qui est impossible.
• Méthode de cet article : Vous réalisez que le sol est si bosselé que les propres tailles des personnes comptent le plus. Les bosses entre elles (les connexions de l'arbre) ne provoquent que de minuscules ajustements spécifiques.
• La découverte : Parce que le sol a la forme d'un arbre, les « oscillations » causées par les connexions s'annulent d'une manière très spécifique (les termes impairs disparaissent). L'auteur vous donne une formule pour calculer exactement comment la forme du sol ajuste la taille moyenne, terme par terme.

En bref, l'article prend un système chaotique et aléatoire et montre que, sous un désordre fort, il se comporte d'une manière étonnamment ordonnée, calculable et lisse, grâce à la géométrie unique de la forêt en forme d'arbre.

Résumé technique

1. Énoncé du problème

L'article étudie le modèle d'Anderson sur l'arbre infini $(q+1)$-régulier (réseau de Bethe), défini par l'hamiltonien :
$$H_{\lambda, \omega} = A + \lambda V_\omega$$
où $A$ est l'opérateur d'adjacence, $\lambda > 0$ est l'intensité du désordre, et $V_\omega$ est un opérateur de multiplication par des potentiels aléatoires indépendants et identiquement distribués (i.i.d.) $\{\omega_x\}$, suivant une loi commune $\mu$.

L'objectif principal est d'analyser la densité d'états (DOS) moyennée à la racine dans le régime de désordre fort ($\lambda \to \infty$). Plus précisément, l'article se concentre sur la fenêtre d'énergie redimensionnée $E = \lambda \xi$, où $\xi$ appartient au support de la distribution mono-site $\mu$.

Distinction clé : Contrairement aux graphes amenables (par exemple $\mathbb{Z}^d$) où la DOS est définie via des limites de comptage d'é valeurs dans des volumes finis, le réseau de Bethe est non amenable. Le volume de la frontière est comparable au volume du volume interne. Par conséquent, l'auteur définit la mesure DOS $N_\lambda$ strictement comme la mesure spectrale moyennée à la racine :
$$N_\lambda(B) = \mathbb{E}[\langle \delta_0, E_{H_{\lambda, \omega}}(B) \delta_0 \rangle]$$
où $\delta_0$ est le vecteur de base à la racine. L'objectif est de prouver l'absolue continuité de cette mesure et de dériver un développement asymptotique explicite de sa densité $n_\lambda(E)$ en puissances de $\lambda^{-1}$.

2. Méthodologie

L'auteur utilise une combinaison de développement en marche aléatoire et d'analyse complexe pour surmonter la difficulté de contrôler la résolvante près de l'axe réel (où la DOS est définie).

A. Développement en marche aléatoire

La résolvante $G_\lambda(0,0;z) = \langle \delta_0, (H_{\lambda, \omega} - z)^{-1} \delta_0 \rangle$ est développée à l'aide de la série de Neumann dans la région où $\text{Im}(z)$ est grand. En traitant le terme de saut $A$ comme une perturbation du terme diagonal $\lambda V_\omega - z$, la résolvante est exprimée comme une somme sur les marches fermées $\gamma$ sur l'arbre :
$$G_\lambda(0,0;z) = \sum_{n=0}^\infty (-1)^n \sum_{\gamma \in \Gamma_n(0,0)} \prod_{j=0}^n \frac{1}{\lambda \omega_{x_j} - z}$$
où $\Gamma_n(0,0)$ est l'ensemble des marches fermées de longueur $n$ commençant et se terminant à la racine. En prenant l'espérance $\mathbb{E}$ et en redimensionnant $z = \lambda \zeta$, on obtient un développement impliquant des produits de transformées de Stieltjes mono-site $s_k(\zeta) = \int \frac{d\mu(t)}{(t-\zeta)^k}$.

B. Prolongement analytique complexe

Un obstacle majeur est que la série de Neumann standard ne converge que pour $\text{Im}(z) > q+1$, ce qui est loin de l'axe réel où la DOS est nécessaire.

• Hypothèse : La distribution mono-site $\mu$ a un support compact et une densité $\rho$ qui est localement analytique sur un intervalle $I^\sharp$ contenant l'intervalle cible $I$.
• Technique : L'auteur utilise un argument d'intégrale de Cauchy pour montrer que les transformées de Stieltjes mono-site $s_k(\zeta)$ admettent un prolongement holomorphe du demi-plan supérieur vers un voisinage complexe $\Omega_\delta(I)$ de l'intervalle réel $I$.
• Extension : Puisque les coefficients du développement en marche aléatoire sont des sommes finies de produits de ces $s_k(\zeta)$, la résolvante moyennée complète $m_\lambda(\lambda \zeta)$ hérite de ce prolongement holomorphe.

C. Exploitation de la structure de l'arbre

La géométrie spécifique du réseau de Bethe est cruciale :

1. Nature bipartite : L'arbre est biparti, ce qui signifie que toute marche fermée commençant et se terminant à la racine doit avoir une longueur paire. Par conséquent, tous les termes d'ordre impair dans le développement s'annulent ($M_n \equiv 0$ pour $n$ impair).
2. Profils d'occupation : Puisque le graphe est un arbre (pas de boucles autres que le retour en arrière), la contribution d'une marche fermée est entièrement déterminée par le profil d'occupation (combien de fois chaque sommet est visité) et le sous-arbre fini visité. Cela permet de calculer les coefficients comme des sommes combinatoires finies.

3. Contributions et résultats clés

A. Théorème principal (Développement de la résolvante)

Sous l'hypothèse d'analyticité locale sur $\mu$, l'auteur prouve que pour $\lambda$ suffisamment grand, la résolvante diagonale moyennée redimensionnée $m_\lambda(\lambda \zeta)$ admet un prolongement holomorphe vers un voisinage complexe de $I$. Elle possède un développement asymptotique uniforme :
$$m_\lambda(\lambda \zeta) = \sum_{n=0}^N \lambda^{-n-1} M_n(\zeta) + R_{N,\lambda}(\zeta)$$
où :

• $M_n(\zeta)$ sont des fonctions holomorphes données par des sommes finies sur les marches fermées de longueur $n$.
• Le reste $R_{N,\lambda}$ est borné par $O(\lambda^{-N-2})$.
• Les coefficients impairs s'annulent : $M_n \equiv 0$ pour tout $n$ impair.

B. Développement de la densité d'états

En appliquant la formule d'inversion de Stieltjes au prolongement holomorphe, l'article établit que la mesure DOS moyennée à la racine est absolument continue sur la fenêtre redimensionnée $\lambda I$. Sa densité $n_\lambda(E)$ admet le développement :
$$n_\lambda(\lambda \xi) = \sum_{n=0}^N \lambda^{-n-1} a_n(\xi) + r_{N,\lambda}(\xi)$$
où $a_n(\xi) = \pi^{-1} \text{Im} M_n(\xi)$.

Propriétés spécifiques des coefficients :

1. Terme dominant ($n=0$) : $a_0(\xi) = \rho(\xi)$, la densité de la distribution mono-site. Cela confirme l'intuition heuristique que dans la limite de désordre fort, la DOS est dominée par le potentiel local.
2. Termes impairs nuls : $a_n(\xi) = 0$ pour tout $n$ impair.
3. Termes d'ordre supérieur : Les coefficients $a_n$ pour $n$ pair $\geq 2$ sont des sommes finies déterminées par les profils d'occupation des marches fermées courtes sur l'arbre. Ils dépendent explicitement du nombre de branches $q$ et des dérivées des transformées de Stieltjes de $\mu$.

C. Calcul explicite pour la distribution uniforme

L'article calcule explicitement le premier terme de correction non nul pour le cas où $\mu$ est la distribution uniforme sur $[-a, a]$.

• $a_0(\xi) = \frac{1}{2a}$
• $a_2(\xi) = -\frac{q+1}{2a(a^2 - \xi^2)}$
• Le développement devient :
  $$n_\lambda(\lambda \xi) = \frac{1}{2a}\lambda^{-1} - \frac{q+1}{2a(a^2 - \xi^2)}\lambda^{-3} + O(\lambda^{-5})$$
  Ce résultat met en évidence que le terme de correction diverge lorsque $\xi$ approche les bords du support ($\pm a$), ce qui est cohérent avec l'exigence que le développement soit valable strictement à l'intérieur du support.

4. Importance

• Premier développement explicite au niveau des coefficients : Avant ce travail, les développements en désordre fort pour la DOS sur le réseau de Bethe étaient soit qualitatifs, soit dépourvus de formules explicites pour les coefficients. Cet article fournit la première dérivation rigoureuse de développements explicites d'ordre fini où les coefficients sont des sommes combinatoires sur les marches de l'arbre.
• Justification rigoureuse des heuristiques : Il valide mathématiquement l'intuition physique selon laquelle un désordre fort conduit à une DOS dominée par la distribution mono-site, les termes de saut fournissant des corrections systématiques et calculables.
• Avancement méthodologique : La combinaison des développements en marche aléatoire avec le prolongement analytique complexe des transformées de Stieltjes offre un cadre robuste pour étudier les propriétés spectrales des opérateurs aléatoires sur des graphes non amenables, contournant les limites des approximations de volume fini.
• Insight structurel : La preuve utilise élégamment la structure de l'arbre (bipartition et absence de cycles) pour simplifier le développement, montrant que la « nature arborescente » conduit à l'annulation des termes d'ordre impair et à la calculabilité des termes d'ordre supérieur.

En résumé, Kaminaga fournit une fondation mathématique rigoureuse pour le comportement en désordre fort du modèle d'Anderson sur le réseau de Bethe, offrant une série asymptotique précise pour la densité d'états qui relie les statistiques locales du désordre aux propriétés spectrales globales à travers la géométrie de l'arbre sous-jacent.