Stationary perturbation theory without sums over intermediate states: Supersymmetric Expansion Algorithm

Résumé Technique : Théorie de la Perturbation Stationnaire sans Sommes sur les États Intermédiaires : Algorithme d'Expansion Supersymétrique

Énoncé du Problème
L'article traite d'une limitation fondamentale de la théorie de la perturbation de Rayleigh-Schrödinger (RSPT) conventionnelle. Bien que la RSPT fournisse un cadre systématique pour calculer les corrections des valeurs propres d'énergie et des états propres, elle repose sur des sommes sur les états intermédiaires. La correction de premier ordre des états propres nécessite déjà une « somme embarrassante », et les corrections d'ordres supérieurs nécessitent des sommes de plus en plus complexes sur les états intermédiaires à chaque ordre. Ce formalisme devient calculatoirement encombrant et souvent impraticable au-delà des ordres les plus bas. Des alternatives antérieures, telles que les règles de somme ou la formulation logarithmique de l'équation de Schrödinger proposée par Aharonov et Au, ont réussi à réduire les calculs de l'état fondamental à des formes de quadrature (intégrales) pour les systèmes unidimensionnels, mais restent moins pratiques pour les états excités en raison de la présence de nœuds.

Méthodologie : L'Algorithme d'Expansion Supersymétrique (SEA)
Les auteurs proposent l'Algorithme d'Expansion Supersymétrique (SEA) pour obtenir les corrections perturbatives d'ordre arbitraire pour tous les états propres (fondamentaux et excités) sans sommer sur les états intermédiaires. La méthode exploite la supersymétrie cachée inhérente aux problèmes de mécanique quantique non relativiste.

La procédure centrale comprend les étapes suivantes :

1. Formulation Logarithmique : L'équation de Schrödinger est formulée sous une forme logarithmique non linéaire impliquant un superpotentiel $W(x, \lambda) = -\frac{d}{dx} \ln u(x, \lambda)$, où $u$ est la fonction d'onde.
2. Expansion Perturbative : Le superpotentiel, le potentiel et l'énergie sont développés en puissances d'un petit paramètre $\lambda$. Cela transforme l'équation différentielle non linéaire en une cascade hiérarchique d'équations différentielles linéaires du premier ordre couplées pour les coefficients d'expansion.
3. États « Limites » sans Nœuds : L'algorithme résout d'abord l'état propre sans nœud (l'état limite) de l'Hamiltonien à un niveau supersymétrique donné. Pour les systèmes unidimensionnels, il s'agit de l'état fondamental ; pour les systèmes multidimensionnels réductibles à des potentiels unidimensionnels effectifs (comme les potentiels centraux), ce sont les états ayant le moment angulaire $l$ le plus élevé pour un nombre quantique principal $n$ donné.
4. Solution Récursive via la Supersymétrie :
   • La solution sans nœud pour le $r$-ième Hamiltonien partenaire $H_r$ est résolue en utilisant l'expansion logarithmique.
   • Le terme d'ordre de grandeur pour la solution sans nœud de $H_r$ est dérivé des solutions non perturbées de l'Hamiltonien original $H_0$ en appliquant de manière récursive l'inverse des opérateurs d'élévation supersymétriques ($a^\dagger$).
   • Une fois la solution sans nœud de $H_r$ établie, les corrections perturbatives pour son énergie et son superpotentiel sont calculées à l'aide d'intégrales pondérées par la densité de probabilité non perturbée de cet état limite spécifique.
   • Les états propres de l'Hamiltonien original $H_0$ possédant $r$ nœuds sont ensuite reconstruits en appliquant la chaîne d'opérateurs d'élévation ($a^\dagger_0, a^\dagger_1, \dots$) à la solution sans nœud de $H_r$.

Contributions Clés et Résultats

• Élimination des Sommes sur les États Intermédiaires : La principale contribution est la réduction des corrections perturbatives à des « formes de quadrature ». Les corrections d'énergie et de fonction d'onde sont exprimées comme des intégrales pondérées par les densités de probabilité des états limites pertinents, évitant totalement la sommation sur les états intermédiaires requise dans la RSPT.
• Ordre Arbitraire et Applicabilité Générale : Le formalisme est applicable à des ordres arbitraires de la théorie de la perturbation et s'étend à tous les états propres, y compris les états excités avec des nœuds. Il fonctionne aussi bien pour les potentiels polynomiaux que non polynomiaux.
• Solutions Analytiques pour des Potentiels Spécifiques : Les auteurs démontrent l'efficacité de la méthode en retrouvant des solutions analytiques connues pour les potentiels de Yukawa, de Hulthén, anharmonique et de Cornell, où les coefficients du superpotentiel sont des polynômes.
• Exemple Non Polynomial : L'article applique la méthode à une particule dans une boîte avec une perturbation d'oscillateur harmonique. Ce cas implique des superpotentiels non polynomiaux. Les auteurs ont dérivé avec succès les corrections d'énergie jusqu'à l'ordre $\lambda^3$ pour l'état fondamental et le premier état excité, fournissant des expressions analytiques explicites pour les valeurs propres d'énergie.

Signification et Revendications
L'article affirme que le SEA fournit une « solution analytique complète » aux problèmes de perturbation sous la forme d'une série de puissances de le paramètre d'expansion. En exploitant la supersymétrie cachée des systèmes quantiques, la méthode offre une alternative transparente et systématique à la RSPT. Les auteurs soulignent que, bien que les méthodes antérieures (comme celles d'Aharonov et Au) aient été limitées dans le traitement des états excités, le SEA connecte naturellement la solution des états excités à la hiérarchie supersymétrique, rendant le calcul des corrections pour tous les états simple. Le travail affirme que cette approche est particulièrement puissante pour les systèmes où les coefficients de l'expansion du potentiel ne sont pas des polynômes, car elle ne repose pas sur la réduction du problème à des relations de récurrence algébriques, mais résout directement la hiérarchie différentielle.