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

Resumen Técnico: Teoría de Perturbación Estacionaria sin Sumas sobre Estados Intermedios: Algoritmo de Expansión Supersimétrica

Planteamiento del Problema
El artículo aborda una limitación fundamental en la teoría de perturbación de Rayleigh-Schrödinger (RSPT) convencional. Si bien la RSPT proporciona un marco sistemático para calcular correcciones a los autovalores y autoestados de energía, depende de sumas sobre estados intermedios. La corrección de primer orden para los autoestados ya requiere una "suma incómoda", y las correcciones de órdenes superiores requieren sumas cada vez más complejas sobre estados intermedios en cada orden. Este formalismo se vuelve computacionalmente engorroso y a menudo impráctico más allá de los órdenes más bajos. Alternativas previas, como las reglas de suma o la formulación logarítmica de la ecuación de Schrödinger propuesta por Aharonov y Au, han logrado reducir los cálculos del estado fundamental a formas de cuadratura (integrales) para sistemas unidimensionales, pero siguen siendo menos convenientes para estados excitados debido a la presencia de nodos.

Metodología: El Algoritmo de Expansión Supersimétrica (SEA)
Los autores proponen el Algoritmo de Expansión Supersimétrica (SEA) para obtener correcciones perturbativas de orden arbitrario para todos los autoestados (fundamentales y excitados) sin sumar sobre estados intermedios. El método aprovecha la supersimetría oculta inherente a los problemas de mecánica cuántica no relativista.

El procedimiento central consta de los siguientes pasos:

1. Formulación Logarítmica: La ecuación de Schrödinger se expresa en una forma logarítmica no lineal que involucra un superpotencial $W(x, \lambda) = -\frac{d}{dx} \ln u(x, \lambda)$, donde $u$ es la función de onda.
2. Expansión Perturbativa: El superpotencial, el potencial y la energía se expanden en potencias de un pequeño parámetro $\lambda$. Esto transforma la ecuación diferencial no lineal en una cascada jerárquica de ecuaciones diferenciales lineales de primer orden acopladas para los coeficientes de expansión.
3. Estados "Borde" sin Nodos: El algoritmo primero resuelve el autoestado sin nodos (el "estado borde") del Hamiltoniano en un nivel supersimétrico dado. Para sistemas unidimensionales, este es el estado fundamental; para sistemas de mayores dimensiones reducibles a potenciales unidimensionales efectivos (como potenciales centrales), estos son los estados con el mayor momento angular $l$ para un número cuántico principal $n$ dado.
4. Solución Recursiva vía Supersimetría:
   • Se resuelve la solución sin nodos para el $r$-ésimo Hamiltoniano compañero $H_r$ utilizando la expansión logarítmica.
   • El término de orden principal para la solución sin nodos de $H_r$ se deriva de las soluciones no perturbadas del Hamiltoniano original $H_0$ aplicando recursivamente la inversa de los operadores de elevación supersimétricos ($a^\dagger$).
   • Una vez establecida la solución sin nodos para $H_r$, las correcciones perturbativas para su energía y superpotencial se calculan mediante integrales ponderadas por la densidad de probabilidad no perturbada de ese estado sin nodos específico.
   • Los autoestados del Hamiltoniano original $H_0$ con $r$ nodos se reconstruyen aplicando la cadena de operadores de elevación ($a^\dagger_0, a^\dagger_1, \dots$) a la solución sin nodos de $H_r$.

Contribuciones Clave y Resultados

• Eliminación de Sumas sobre Estados Intermedios: La contribución principal es la reducción de las correcciones perturbativas a "formas de cuadratura". Las correcciones de energía y de la función de onda se expresan como integrales ponderadas por las densidades de probabilidad de los estados borde relevantes, evitando completamente la sumatoria sobre estados intermedios requerida en la RSPT.
• Orden Arbitrario y Aplicabilidad General: El formalismo es aplicable a órdenes arbitrarios de la teoría de perturbación y se extiende a todos los autoestados, incluyendo estados excitados con nodos. Funciona tanto para potenciales polinómicos como no polinómicos.
• Soluciones Analíticas para Potenciales Específicos: Los autores demuestran la eficacia del método recuperando soluciones analíticas conocidas para los potenciales Yukawa, Hulthén, anharmónico y Cornell, donde los coeficientes del superpotencial son polinomios.
• Ejemplo No Polinómico: El artículo aplica el método a una partícula en una caja con una perturbación de oscilador armónico. Este caso involucra superpotenciales no polinómicos. Los autores derivan con éxito las correcciones de energía hasta el orden $\lambda^3$ tanto para el estado fundamental como para el primer estado excitado, proporcionando expresiones analíticas explícitas para los autovalores de energía.

Significancia y Reivindicaciones
El artículo afirma que el SEA proporciona una "solución analítica completa" a los problemas de perturbación en forma de una serie de potencias en el parámetro de expansión. Al aprovechar la supersimetría oculta de los sistemas cuánticos, el método ofrece una alternativa transparente y sistemática a la RSPT. Los autores enfatizan que, mientras que los métodos anteriores (como el de Aharonov y Au) estaban limitados en el manejo de estados excitados, el SEA conecta naturalmente la solución de los estados excitados con la jerarquía supersimétrica, haciendo que el cálculo de las correcciones para todos los estados sea directo. El trabajo sostiene que este enfoque es particularmente poderoso para sistemas donde los coeficientes de la expansión del potencial no son polinomios, ya que no depende de la reducción del problema a relaciones de recurrencia algebraicas, sino que resuelve directamente la jerarquía diferencial.