Meshfree versus grid-based Schrödinger solvers for modeling the interactions between free-electron wave packets and light

Resumen Técnico: Solucionadores de Schrödinger sin Malla frente a Solucionadores Basados en Rejilla para el Modelado de Interacciones de Paquetes de Ondas de Electrones Libres con la Luz

Planteamiento del Problema
La interacción entre paquetes de ondas de electrones libres y campos electromagnéticos es fundamental para la microscopía electrónica ultrarrápida y la nanofotónica cuántica, lo que permite la ingeniería de paquetes de ondas de electrones para aplicaciones tales como peines de energía y superposiciones de estados de momento. Las descripciones teóricas existentes a menudo dependen de aproximaciones eicónicas o de no retroceso (no-recoil), las cuales fallan al no capturar la evolución dinámica completa de los paquetes de ondas de electrones en regímenes que involucran el retroceso del electrón, efectos de campo fuerte o dinámicas más allá de la aproximación de rotación de onda.

Se requieren enfoques numéricos basados en el Hamiltoniano de acoplamiento mínimo para un tratamiento completo. Sin embargo, los métodos estándar basados en rejillas (grid-based) que acoplan la ecuación de Schrödinger dependiente del tiempo con las ecuaciones de Maxwell enfrentan cuellos de botella computacionales significativos. Esto surge de la sustancial disparidad entre las escalas de discretización espacial y temporal requeridas para la ecuación de Schrödinger frente a las ecuaciones de Maxwell, haciendo que el solucionador de Schrödinger sea el factor limitante en las simulaciones híbridas.

Metodología
Los autores presentan un marco numérico sin malla (mesh-free) que resuelve directamente la ecuación de Schrödinger de partícula única dependiente del tiempo para potenciales electromagnéticos arbitrarios. El método se basa en la aproximación de Paquete de Ondas Gaussiano Descongelado (TGWP, por sus siglas en inglés).

• Formulación Matemática: La función de onda $\psi(t, r)$ se aproxima como una suma de $N$ paquetes de ondas gaussianos:
  $$\psi(t, r) \approx \sum_{j=1}^N c_j e^{iS_j(t)/\hbar} g_j(t, r)$$
  donde $c_j$ son coeficientes complejos independientes del tiempo, y $g_j(t, r)$ son gaussianas dependientes del tiempo definidas por sus centros en posición ($q_j$) y momento ($p_j$), matrices de ancho ($C_j$), y factores de normalización/fase ($\gamma_j$).
• Inicialización: El paquete de ondas gaussiano inicial se representa como un conjunto sobrecompleto de gaussianas. Los coeficientes de expansión se determinan mediante una cuadratura de cuasi-Monte Carlo con puntos de Sobol' ponderados por una función gaussiana.
• Evolución: Debido a la linealidad de la ecuación de Schrödinger, cada gaussiana se evoluciona semiclásicamente e independientemente. Los parámetros ($q_j, p_j, Q_j, P_j, S_j$) siguen las ecuaciones de movimiento clásicas derivadas de la aproximación armónica local del potencial $V(t, r)$ alrededor del centro del paquete. Las matrices de ancho se parametrizan utilizando la parametrización de Hagedorn ($C_j = P_j Q_j^{-1}$). Las ecuaciones de movimiento se integran utilizando un integrador de segundo orden de Störmer-Verlet.
• Evaluación Comparativa (Benchmarking): El método TGWP se evalúa frente a un solucionador de referencia basado en rejilla utilizando el método de Fourier de pasos divididos (split-step Fourier) con división de Strang. Ambos métodos se aplican a dos escenarios físicos distintos:
  1. Campo cercano cuasi-estático: Interacción con el campo cercano óptico de una nanobarra de oro (configuración PINEM).
  2. Dispersión Compton Estimulada: Interacción con dos pulsos láser superpuestos e inclinados en el espacio libre.

Resultos Clave
El estudio compara la aproximación TGWP contra el solucionador basado en rejilla a través de varios conjuntos de parámetros, incluyendo diferentes intensidades de campo eléctrico, energías cinéticas de electrones y longitudes de onda de láser.

• Precisión:
  • Campo cercano (PINEM): El método TGWP reproduce la deformación transitoria del paquete de ondas y la distribución final de momento con alta fidelidad. Aunque el método produce oscilaciones artificiales en regiones de intensidad de función de onda insignificante (contribuyendo al error $L^2$), la envolvente espectral, las posiciones de los bandas laterales y la redistribución de intensidad coinciden estrechamente con la solución basada en rejilla. Los errores de altura de pico se mantuvieron por debajo del 1.2% y los errores de ancho por debajo del 2.3%.
  • Dispersión Compton Estimulada: El acuerdo fue excelente, con el método TGWP capturando con precisión la redistribución coherente de la función de onda en el espacio de momento y las estructuras de interferencia oscilatorias. El error regularizado se mantuvo por debajo del 1% durante toda la propagación.
  • Robustez: El método mantuvo su precisión incluso cuando los parámetros se desviaron de las condiciones ideales de concordancia de fase, reproduciendo correctamente las tendencias físicas en regímenes fuera de resonancia.
• Desempeño:
  • El enfoque sin malla demostró aceleraciones computacionales significativas. En el caso del campo cercano, se observó una aceleración de aproximadamente 100x en hardware comparable.
  • En el caso de la dispersión Compton estimulada, la aceleración alcanzó aproximadamente 800x (el TGWP tomó ~1 segundo frente a ~14 minutos para el método de pasos divididos por cada conjunto de parámetros).
  • El método escala eficientemente al eliminar la necesidad de una malla espacial, que es el principal cuello de botella en el acoplamiento con solucionadores de Maxwell basados en rejilla.

Significancia y Reivindicaciones
El artículo establece el método TGWP como una herramienta computacional versátil, eficiente y escalable para la óptica electrónica cuántica. Sus principales contribaciones son:

1. Más allá de las aproximaciones: Permite la simulación de interacciones electrón-luz más allá de las aproximaciones eicónicas y de no retroceso, capturando la dinámica completa de acoplamiento mínimo.
2. Eficiencia: Al eliminar el requisito de la malla espacial, el método supera la disparidad de escala de discretización que afecta a las simulaciones híbridas de Schrödinger-Maxwell, ofreciendo aceleraciones de hasta dos órdenes de magnitud.
3. Validación: El método es rigurosamente validado contra un benchmark basado en rejilla tanto para interacciones de campo cercano localizadas como para la dispersión en espacio libre, mostrando un excelente acuerdo cuantitativo.

Los autores concluyen que este marco proporciona una base sólida para modelar entornos electromagnéticos complejos dependientes del tiempo. Sugieren que, combinado con algoritmos de optimización y aprendizaje automático, este enfoque podría facilitar el diseño inverso de entornos electromagnéticos para el moldeado bajo demanda y el control coherente de paquetes de ondas de electrones. El trabajo también posiciona al TGWP como un candidato prometedor para resolver la dinámica completa de Maxwell–Schrödinger, particularmente dado su aplicación exitosa previa a Hamiltonianos magnéticos (acoplamiento mínimo).