On the electronic path integral normal modes of the Meyer-Miller-Stock-Thoss representation of nonadiabatic dynamics

Resumen Técnico: Sobre los modos normales de la integral de trayectoria electrónica en la representación de Meyer–Miller–Stock–Thoss de la dinámica no adiabática

Planteamiento del Problema
La simulación precisa y eficiente de la dinámica no adiabática es esencial para comprender la transferencia de carga y energía en sistemas complejos. Un criterio crítico para la validez de tales métodos es la conservación de la Distribución de Boltzmann Cuántica (QBD, por sus siglas en inglés). Mientras que la dinámica de Matsubara logra conservar la QBD para la dinámica de núcleos de una sola superficie mediante la truncación de los modos normales superiores de la integral de trayectoria en el tiempo imaginario, su extensión a sistemas no adiabáticos (NA-Mats) presenta un desafío. El método NA-Mats existente aplica la truncación de Matsubara a las variables nucleares pero retiene la representación de Meyer–Miller–Stock–Thoss (MMST) para las variables electrónicas sin truncación de modos normales. En consecuencia, este método no logra conservar la QBD para sistemas generales. Esto plantea una pregunta fundamental: ¿el truncamiento de los modos normales superiores de la integral de trayectoria de las variables electrónicas dentro de la representación MMST conduce a la conservación de la QBD, o son estos modos fundamentalmente diferentes de sus contrapartes nucleares?

Metodología
Los autores investigan las propiedades de los modos normales electrónicos dentro del marco de mapeo MMST utilizando un sistema exclusivamente electrónico para aislar el comportamiento de los grados de libertad electrónicos.

• Marco Teórico: El estudio utiliza el Hamiltoniano MMST, que mapea un sistema cuántico de múltiples niveles en osciladores armónicos acoplados con variables de posición ($q$) y momento ($p$) electrónicas. Los autores definen los modos normales electrónicos ($\check{q}_k, \check{p}_k$) utilizando la misma matriz de transformación ($T$) empleada para los modos de Matsubara nucleares, efectivamente diagonalizando la matriz de la constante de resorte del polímero de anillo libre.
• Análisis del Liouvilliano: Se deriva el Liouvilliano electrónico exacto en la base de los modos normales. Para un sistema con grados de libertad nucleares y electrónicos desacoplados, se muestra que las ecuaciones de evolución para los modos normales son separables e idénticas en forma a las ecuaciones de evolución de las cuentas (beads).
• Funciones de Correlación: Los autores formulan las funciones de correlación de la Transformada de Kubo Generalizada (KT) en términos tanto de las variables de las cuentas como de las variables de los modos normales. Específicamente, examinan si los observables (como las poblaciones de estado) pueden expresarse como funciones de un número finito de los modos normales más bajos.
• Protocolo de Truncación: Las pruebas computacionales consisten en truncar la integral de trayectoria de dos maneras:
  1. Truncación de Cuentas (Bead Truncation): Propagar solo las primeras $j$ cuentas mientras el resto permanece en $t=0$.
  2. Truncación de Modos Normales: Propagar solo los $k$ modos normales más bajos mientras los modos superiores permanecen en $t=0$.
• Modelo del Sistema: Se simula un sistema electrónico de dos niveles con una matriz de potencial ($V$) asimétrica utilizando unidades reducidas ($\beta=1$) y un polímero de anillo de 8 cuentas.

Contribuciones Clave y Resultados
El estudio presenta el primer análisis de los modos normales electrónicos en la representación MMST, arrojando varios hallazgos contraintuitivos en comparación con la dinámica de Matsubara nuclear convencional:

1. Observables No Finitos: A diferencia de los observables nucleares, que a menudo dependen de un número finito de modos de Matsubara de baja frecuencia, los observables electrónicos en la representación MMST son funciones de todos los modos normales. La matriz de transformación no puede eliminarse al realizar la truncación, lo que significa que el observable depende del conjunto completo de modos.
2. Falta de Estrechamiento de la Distribución: En la dinámica de Matsubara nuclear, la distribución de los modos normales superiores se estrecha, lo que permite integrarlos. En la representación electrónica MMST, los modos normales superiores no están restringidos por la distribución; su distribución no se estrecha, lo que impide la integración de los modos de alta frecuencia sin pérdida de precisión.
3. Fallo en la Conservación de la QBD para Trayectorias Individuales: El análisis algebraico y los resultados computacionales demuestran que el truncamiento de la dinámica en los modos normales MMST no conserva la QBD para una sola trayectoria. La conservación de la QBD en el caso de múltiples cuentas requiere la completitud de todo el bucle de la integral de trayectoria; la truncación rompe esta condición.
4. Artefacto de Promedio de Ensamble: Aunque la QBD parece conservarse al promediar sobre un ensamble de trayectorias (debido a un efecto de promedio), no se conserva para trayectorias individuales. Esto distingue al método de una dinámica de verdadera conservación de la QBD.
5. Funciones de Correlación Inexactas: El truncamiento de la dinámica en los modos normales MMST conduce a funciones de correlación inexactas, ya que se requieren todos los modos para obtener una dinámica razonable.

Significancia y Reivindicaciones
El artículo concluye que los modos normales de la integral de trayectoria MMST no son óptimos para desarrollar un método de dinámica no adiabática preciso que conserve la QBD. Los autores argumentan que las propiedades específicas requeridas para que la dinámica de Matsubara conserve la QBD —específicamente, que los observables sean funciones de un número finito de modos y el estrechamiento de las distribuciones de los modos superiores— no se cumplen para las variables electrónicas MMST. Por consiguiente, la simple aplicación de la truncación de Matsubara a las variables electrónicas en el marco NA-Mats no produce un método que conserve la QBD para una sola trayectoria. El trabajo sugiere que pueden ser necesarias representaciones de mapeo o métricas alternativas para lograr un método de dinámica no adiabática que satisfaga simultáneamente los criterios de conservación de la QBD, reproducción de oscilaciones de Rabi y escalamiento clásico.