Preserving Mass Shell Condition in the Stochastic Optimal Control Derivation of the Dirac Equation

Resumen Técnico: Preservación de la Condición de la Cáscara de Masa en la Derivación del Control Óptimo Estocástico de la Ecuación de Dirac

Planteamiento del Problema
La ecuación de Dirac es la piedra angular de la mecánica cuántica relativista para partículas de espín-1/2, asegurando que cada componente del espínor satisfaga la ecuación de Klein-Gordon y, por ende, la condición de la cáscara de masa relativista ($E^2 = p^2c^2 + m^2c^4$). Si bien la mecánica estocástica ha recuperado con éxito las ecuaciones de Schrödinger y de Klein-Gordon, la derivación de la ecuación de Dirac mediante el Control Óptimo Estocástico (SOC, por sus siglas en inglés) ha sido históricamente elusiva. Intentos previos, incluyendo el trabajo previo de los autores [10], dependieron de la linealización del término cinético relativista ($-mc\sqrt{u^\mu u_\mu}$) y de la ecuación de Hamilton-Jacobi-Bellman (HJB). Aunque esta linealización facilitó la emergencia de una estructura de tipo Dirac, no logró preservar la condición de la cáscara de masa; específicamente, la ecuación factorizada resultante omitió el término necesario $+mc$ en el operador conjugado, impidiendo la recuperación de la ecuación de Klein-Gordon para los componentes individuales.

Metodología
Este trabajo reformula la derivación por SOC de la ecuación de Dirac reteniendo el término cinético de raíz cuadrada relativista estándar en su forma no lineal original, en lugar de linealizarlo. La metodología procede de la siguiente manera:

1. Formulación Lagrangiana Covariante: Los autores definen un Lagrangiano de Lorentz-covariante para una única partícula cargada con espín-1/2 en un campo electromagnético externo. Este Lagrangiano incluye:

   • El término cinético no lineal estándar: $-mc\sqrt{w^\mu w_\mu}$.
   • Acoplamiento electromagnético mínimo: $-\epsilon e A_\mu w^\mu$.
   • Un acoplamiento espín-campo covariante (término de Pauli): $-\frac{e\hbar}{2m}\sigma_{\mu\nu}F^{\mu\nu}$.

2. Escalarización para HJB Escalar: Dado que el término de espín-campo $\sigma_{\mu\nu}F^{\mu\nu}$ es de valor matricial, no puede insertarse directamente en la ecuación HJB escalar. Para resolver esto, los autores introducen un "espínor de escalarización" $\chi$, un espínor propio constante e independiente del espacio-tiempo de la matriz de espín-campo ($\sigma_{\mu\nu}F^{\mu\nu}\chi = f\chi$). Definen un Lagrangiano escalarizado $L_\chi$ mediante la proyección bilineal $P_\chi[M] = \frac{\chi^\dagger S M \chi}{\chi^\dagger S \chi}$. Esto restringe la derivación a fondos electromagnéticos que admiten tal espínor propio (por ejemplo, campos magnéticos uniformes) y fija el sistema a un sector de espín específico.

3. SOC Complejo y Derivación de HJB: El problema se plantea en la teoría de SOC Complejo, donde la partícula se mueve en un espacio-tiempo complejo de cuatro dimensiones. La función de valor $J$ satisface la ecuación HJB. Al aplicar la "condición débil" $w^\mu w_\mu = c^2$ solo después de la diferenciación (en lugar de linealizar el término cinético), los autores derivan una política de control óptimo $w_\mu$ dependiente del gradiente de la función de valor.

4. Transformación a la Ecuación de Dirac:

   • La ecuación HJB se transforma utilizando el mapa de Hopf-Cole ($J = i\epsilon\hbar \ln \phi$), convirtiendo la HJB escalar no lineal en una ecuación lineal de segundo orden para el campo escalar $\phi$.
   • Se demuestra que esta ecuación lineal es equivalente a una ecuación "cuadrada-Dirac": $(i\hbar\gamma^\nu\partial_\nu - e\gamma^\nu A_\nu - mc)(i\hbar\gamma^\mu\partial_\mu - e\gamma^\mu A_\mu + mc)\phi = 0$.
   • El espínor de Dirac $\psi$ se construye como $\psi = (i\hbar\gamma^\mu\partial_\mu - e\gamma^\mu A_\mu + mc)\phi \chi$.
   • En el límite estacionario ($\partial_\tau \to 0$), esto produce la ecuación de Dirac estándar: $(i\hbar\gamma^\mu\partial_\mu - e\gamma^\mu A_\mu - mc)\psi = 0$.

5. Análisis de la Cáscara de Masa: Los autores analizan la condición de la cáscara de masa tanto en el límite clásico ($\hbar \to 0$) como en el cuántico. Demuestran que la relación de la cáscara de masa corregida cuánticamente incluye contribuciones de la covariación cuadrática de Itô (potencial cuántico) y del término de Pauli. Crucialmente, muestran que el observable bilineal físico $\langle \pi_\mu \pi^\mu \rangle_{\text{phys}}$ recupera la condición de la cáscara de masa exactamente cuando se aplica la corrección de la covariación de Itô, cancelando el término del potencial cuántico y dejando el invariante de Pauli.

Resultados Clave

• Restauración de la Condición de la Cáscara de Masa: A diferencia de las derivaciones linealizadas previas, este enfoque preserva la condición de la cáscara de masa relativista. La ecuación derivada se factoriza naturalmente en el operador de Dirac y su conjugado, asegurando que las soluciones satisfagan la ecuación de Klein-Gordon.
• Restricción al Sector de Espín: La derivación es válida para fondos donde existe un espínor propio del espacio-tiempo independiente de un $\sigma_{\mu\nu}F^{\mu\nu}$. Esto permite tratar la interacción de espín de valor matricial como un escalar dentro de un sector de espín fijo.
• Validación Numérica (Problema de Dirac-Landau): La teoría fue probada numéricamente utilizando un electrón en un campo magnético uniforme (problema de Dirac-Landau).
  • Minimización de la Acción: Las simulaciones estocásticas confirmaron que la acción estocástica promedio alcanza un mínimo local en la deriva óptima analíticamente derivada $w^\star$.
  • Acuerdo Cuantitativo: Los componentes promedio de la acción estocástica (cinética, electromagnética y de espín) computados bajo la deriva óptima coincidieron con los valores analíticos de Dirac-Landau dentro de intervalos de confianza del 95%.
  • Optimalidad Local: El análisis de perturbación mostró que las desviaciones de la deriva óptima resultan en un aumento cuadrático de la acción esperada, confirmando la minimidad local.
• Interpretación Geométrica: Los autores visualizaron el espacio de momento complejo, mostrando que la cáscara de masa cuántica es una deformación de la cáscara de masa clásica débil ($\pi_\mu \pi^\mu = m^2c^2$) por pequeñas correcciones cuánticas, que permanecen acotadas y oscilatorias.

Significancia y Reivindicaciones
El artículo afirma proporcionar una "derivación de SOC relativísticamente consistente" de la ecuación de Dirac que resuelve la inconsistencia de los enfoques linealizados anteriores respecto a la condición de la cáscara de masa. Al retener el término cinético no lineal y utilizar la escalarización para el acoplamiento espín-campo, los autores cierran la brecha entre el control óptimo estocástico y la mecánica cuántica relativista sin sacrificar la relación fundamental de energía-momento.

El trabajo destaca que la ecuación de Dirac puede derivarse de los primeros principios de la mecánica estocástica (específicamente del SOC Complejo), donde el límite $\hbar \to 0$ recupera la ecuación de Hamilton-Jacobi relativista de tiempo propio. Los autores posicionan esto como un marco robusto donde la recuperación de espectros y funciones de onda se enmarca como un problema de optimización de alta dimensión sobre campos de deriva, potencialmente aplicable a futuros algoritmos cuánticos variacionales (VQE, QAOA). El artículo no propone nuevos montajes experimentales, sino que ofrece una validación teórica y numérica del marco de SOC para partículas de espín-1/2 relativistas.