Spectral Asymptotics of a Singular Harmonic Oscillator with Aharonov--Bohm Flux

Resumen Técnico: Asintótica Espectral de un Oscilador Armónico Singular con Flujo de Aharonov–Bohm

Planteamiento del Problema
El artículo investiga la asintótica espectral del operador oscilador armónico singular $T_h$ definido en el intervalo $(0, 1)$:
$$T_h = -h^2 \frac{d^2}{dx^2} + x^2 + \frac{h^2(\nu^2 - 1/4)}{x^2},$$
sujeto a condiciones de contorno de Dirichlet en $x=0$ y $x=1$. Aquí, $h > 0$ es un parámetro semiclásico y $\nu \ge 0$ es un parámetro fijo.

Este operador surge naturalmente de dos contextos físicos:

1. Teoría de Control: Como la transformada de Fourier parcial del operador de Baouendi-Grushin.
2. Mecánica Cuántica Magnética: Como la reducción radial del operador de Schrödinger magnético en el disco unitario en presencia de un campo magnético constante $B$ y un flujo de Aharonov–Bohm $\alpha$. En este contexto, $h = 2/B$ y $\nu = |m + \alpha|$, donde $m$ es el momento angular. El límite $h \to 0^+$ corresponde al límite de campo magnético fuerte ($B \to +\infty$).

El problema de autovalores para $T_h$ es equivalente a hallar los ceros de la función hipergeométrica de Kummer $M(a, b, z)$ con respecto a su primer parámetro $a$, específicamente la condición $M\left(-\frac{\lambda + 2(1+\nu)h}{4h}, 1+\nu, \frac{1}{h}\right) = 0$.

Metodología
Los autores emplean una combinación de análisis semiclásico y la teoría asintótica uniforme de funciones especiales. La metodología se divide en tres regímenes espectrales distintos basados en el nivel de energía relativo a la barrera de potencial en el contorno $x=1$:

1. Regimen de Tunelización (Estados de Baja Energía a Semi-Excitados): Para niveles de energía donde el punto de retorno clásico se encuentra estrictamente dentro del dominio $(0, 1)$, los autovalores están exponencialmente cerca de aquellos del operador en la semirrecta $\mathbb{R}^+$. Los autores utilizan:

   • Estimaciones de Agmon para establecer el decaimiento exponencial de las autofunciones en la región prohibida.
   • Asintóticas de Plancherel–Rotach para polinomios de Laguerre asociados (para bajas energías) y asintóticas de funciones de Whittaker (para energías más altas) desarrolladas por Dunster, Gabutti y Gatteschi.
   • Construcción de estados de prueba mediante una técnica de reducción de orden para derivar expansiones asintóticas precisas para el desdoblamiento de autovalores $\lambda_k(h) - e_k(h)$.

2. Regimen Supercrítico (Alta Energía): Para energías $E > 1$ donde todo el intervalo $(0, 1)$ es clásicamente permitido, el efecto de tunelización desaparece. El análisis se basa en:

   • Expansiones asintóticas uniformes de las funciones de Whittaker $M_{\kappa, \mu}(z)$.
   • Inversión de fórmulas asintóticas para los ceros positivos de estas funciones para derivar una regla de cuantización de Bohr–Sommerfeld desplazada.

3. Regimen de Transición: Cerca de la energía crítica $E=1$, donde el punto de retorno clásico toca el contorno, las expansiones estándar fallan. Los autores utilizan:

   • Asintóticas de Airy uniformes (Dunster) para describir los autovalores en una ventana de tamaño $O(h^{2/3})$ alrededor del umbral crítico.
   • Argumentos de emparejamiento (matching) para demostrar que la expansión de Airy conecta consistentemente las asintóticas de tunelización (subcrítico) y la regla de Bohr–Sommerfeld (supercrítico).

Contribuciones Clave y Resultados

• Asintótica de Tunelización Refinada (Teoremas 1.1 y 1.2):
  El artículo extiende los resultados conocidos para $k$ fijo a un régimen uniforme donde el número cuántico $k$ crece conforme $h \to 0$.

  • Para "estados semi-excitados" ($k \lesssim h^{-\gamma}$ con $\gamma < 1/2$), los autores recuperan la fórmula de Gannot con un prefactor explícito que involucra el parámetro $\nu$.
  • Para estados más excitados hasta el umbral de tunelización ($\ell(h) \le e_k(h) \le \rho^* < 1$), derivan el desdoblamiento asintótico:
    $$\lambda_k(h) = e_k(h) + \frac{2}{\pi} h \exp\left(-\frac{2S_{k,h}}{h}\right) (1 + o(1)),$$
    donde $S_{k,h}$ es la distancia de Agmon desde el contorno a la región clásica. Esto cuantifica el mecanismo de tunelización para estados excitados y exhibe explícitamente la influencia del flujo de Aharonov–Bohm $\nu$ en el prefactor.

• Regla de Bohr–Sommerfeld Desplazada (Teorema 1.4):
  Para energías $E > 1$, los autores derivan una regla de cuantización de Bohr–Sommerfeld adaptada al potencial singular:
  $$\lambda_k(h) = E_0\left(4h\left(k + \frac{\nu}{2} - \frac{1}{4}\right)\right) + O(h^2).$$
  El desplazamiento $\frac{\nu}{2} - \frac{1}{4}$ se identifica como la corrección de Bohr–Sommerfeld–Langer derivada de la singularidad de inverso al cuadrado en el origen.

• Expansión del Régimen de Transición (Teorema 1.6):
  En la ventana crítica $|4kh - 1| \le T h^{2/3}$, los autores proporcionan una expansión asintótica de tres términos que involucra la función de fase de Airy $\theta$. Esta expansión interpola suavemente entre el comportamiento de tunelización exponencial y el comportamiento algebraico de Bohr–Sommerfeld.

• Brecha Espectral Uniforme (Teorema 1.7):
  El artículo demuestra un límite inferior no asintótico para la brecha espectral:
  $$\lambda_{k+1}(h) - \lambda_k(h) > 4h,$$
  válido para todo $k \ge 1$ y $h > 0$. Este resultado se establece utilizando la teoría de Sturm–Liouville y la construcción de un operador de tipo "creación", independiente de los límites asintóticos.

Significación y Reivindicaciones
Los autores posicionan este trabajo como una descripción exhaustiva de los regímenes espectrales de un oscilador armónico singular, motivado por observaciones recientes (Vanlaere, 2026) sobre la persistencia de la tunelización hacia niveles de energía más altos.

• Unificación: El artículo unifica el tratamiento de la tunelización de baja energía, la cuantización de alta energía y la transición crítica, proporcionando una imagen asintótica completa a través de todas las escalas de energía.
• Precisión: Va más allá de los límites superiores/inferiores para proporcionar fórmulas asintóticas uniformes explícitas, incluyendo las tasas de decaimiento exponencial óptimas y los prefactores algebraicos específicos dependientes del flujo $\nu$.
• Funciones Especiales: El trabajo traduce los resultados espectrales en fórmulas asintóticas explícitas para los ceros de la función hipergeométrica de Kummer, aprovechando mejoras recientes en la teoría asintótica uniforme de las funciones de Whittaker.
• Interpretación Física: Los resultados aclaran el papel del flujo de Aharonov–Bohm en el desdoblamiento espectral y la condición de cuantización, mostrando cómo el flujo modifica el número cuántico efectivo en el límite semiclásico.

El artículo no propone nuevos montajes experimentales, sino que proporciona fundamentos matemáticos rigurosos para comprender las propiedades espectrales de sistemas modelados por este operador, con relevancia específica para la teoría de control y sistemas cuánticos magnéticos.