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

Resumo Técnico: Teoria de Perturbação Estacionária sem Somas sobre Estados Intermediários: Algoritmo de Expansão Supersimétrica

Definição do Problema
O artigo aborda uma limitação fundamental da convencional teoria de perturbação de Rayleigh-Schrödinger (RSPT). Embora a RSPT forneça um arcabouço sistemático para o cálculo de correções aos autovalores e autoestados de energia, ela depende de somas sobre estados intermediários. A correção de primeira ordem para os autoestados já exige uma "soma incômoda", e as correções de ordens superiores necessitam de somas cada vez mais complexas sobre estados intermediários em cada ordem. Esse formalismo torna-se computacionalmente pesado e frequentemente impraticável além das ordens mais baixas. Alternativas anteriores, como as regras de soma ou a formulação logarítmica da equação de Schrödinger proposta por Aharonov e Au, conseguiram reduzir os cálculos do estado fundamental a formas de quadratura (integrais) para sistemas unidimensionais, mas permanecem menos convenientes para estados excitados devido à presença de nós.

Metodologia: O Algoritmo de Expansão Supersimétrica (SEA)
Os autores propõem o Algoritmo de Expansão Supersimétrica (SEA) para obter correções perturbativas de ordem arbitrária para todos os autoestados (fundamentais e excitados) sem somar sobre estados intermediários. O método explora a supersimetria oculta inerente aos problemas de mecânica quântica não relativística.

O procedimento central envolve as seguintes etapas:

1. Formulação Logarítmica: A equação de Schrödinger é posta em uma forma logarítmica não linear envolvendo um superpotencial $W(x, \lambda) = -\frac{d}{dx} \ln u(x, \lambda)$, onde $u$ é a função de onda.
2. Expansão Perturbativa: O superpotencial, o potencial e a energia são expandidos em potências de um pequeno parâmetro $\lambda$. Isso transforma a equação diferencial não linear em uma cascata hierárquica de equações diferenciais lineares de primeira ordem acopladas para os coeficientes de expansão.
3. Estados "Borda" Sem Nós: O algoritmo primeiro resolve o autoestado sem nós (edge state) do Hamiltoniano no nível supersimétrico dado. Para sistemas unidimensionais, este é o estado fundamental; para sistemas de dimensões superiores redutíveis a potenciais unidimensionais efetivos (como potenciais centrais), estes são os estados com o maior momento angular $l$ para um dado número principal $n$.
4. Solução Recursiva via Supersimetria:
   • A solução sem nós para o $r$-ésimo Hamiltoniano parceiro supersimétrico $H_r$ é resolvida usando a expansão logarítmica.
   • O termo de ordem líder para a solução sem nós de $H_r$ é derivado das soluções não perturbadas do Hamiltoniano original $H_0$ aplicando recursivamente o inverso dos operadores de levantamento supersimétricos ($a^\dagger$).
   • Uma vez estabelecida a solução sem nós para $H_r$, as correções perturbativas para sua energia e superpotencial são calculadas usando integrais ponderadas pela densidade de probabilidade não perturbada desse estado sem nós específico.
   • Os autoestados do Hamiltoniano original $H_0$ com $r$ nós são então reconstruídos aplicando a cadeia de operadores de levantamento ($a^\dagger_0, a^\dagger_1, \dots$) à solução sem nós de $H_r$.

Principais Contribuições e Resultados

• Eliminação de Somas de Estados Intermediários: A principal contribuição é a redução das correções perturbativas para "formas de quadratura". As correções de energia e de função de onda são expressas como integrais ponderados pelas densidades de probabilidade dos estados de borda relevantes, evitando completamente a soma sobre estados intermediários exigida na RSPT.
• Ordem Arbitrária e Aplicabilidade Geral: O formalismo é aplicável a ordens arbitrárias de teoria de perturbação e estende-se a todos os autoestados, incluindo estados excitados com nós. Funciona tanto para potenciais polinomiais quanto não polinomiais.
• Soluções Analíticas para Potenciais Específicos: Os autores demonstram a eficácia do método recuperando soluções analíticas conhecidas para os potenciais Yukawa, Hulthén, anarmônico e Cornell, onde os coeficientes do superpotencial são polinômios.
• Exemplo Não Polinomial: O artigo aplica o método a uma partícula em uma caixa com uma perturbação de oscilador harmônico. Este caso envolve super potenciais não polinomiais. Os autores derivam com sucesso as correções de energia até a ordem $\lambda^3$ tanto para o estado fundamental quanto para o primeiro estado excitado, fornecendo expressões analíticas explícitas para os autovalores de energia.

Significância e Alegações
O artigo afirma que o SEA fornece uma "solução analítica completa" para problemas de perturbação na forma de uma série de potências no parâmetro de expansão. Ao alavancar a supersimetria oculta dos sistemas quânticos, o método oferece uma alternativa transparente e sistemática à RSPT. Os autores enfatizam que, enquanto métodos anteriores (como o de Aharonov e Au) eram limitados no tratamento de estados excitados, o SEA conecta naturalmente a solução de estados excitados à hierarquia supersimétrica, tornando o cálculo das correções para todos os estados direto. O trabalho afirma que esta abordagem é particularmente poderosa para sistemas onde os coeficientes da expansão do potencial não são polinômios, pois não depende da redução do problema a relações de recorrência algébricas, mas sim da resolução direta da hierarquia diferencial.