Tracking Adiabaticity in Non-Equilibrium Many-Body Systems: The Hard Case of the X-ray Photoemission in Metals

Resumen Técnico: Seguimiento de la Adiabaticidad en Sistemas de Muchos Cuerpos Fuera del Equilibrio

Planteamiento del Problema
El artículo aborda el desafío de rastrear de manera fiable la adiabaticidad en sistemas cuánticos complejos de muchos cuerpos fuera del equilibrio. Aunque el Criterio Adiabático Cuántico (QAC) es la herramienta estándar para evaluar la evolución adiabática, enfrenta limitaciones significativas: depende de suposiciones perturbativas, tiene dificultades con los espectros de energía continuos (comunes en metales) y no logra contabilizar los efectos de memoria o la dinámica no markoviana. Los autores proponen probar métodos basados en métricas recientemente desarrollados —específicamente aquellos que utilizan la distancia de Bures, la distancia de traza y la distancia de Densidad Local Natural (NLD)— en un escenario de "caso difícil": la fotoemisión de rayos X en metales. Este fenómeno se caracteriza por un carácter fuertemente fuera del equilibrio, un espectro de energía continuo y la catástrofe de ortogonalidad de Anderson, donde el estado fundamental del sistema se vuelve ortogonal al estado inicial debido a un potencial de dispersión localizado.

Metodología
Los autores investigan la evolución temporal de un gas de Fermi sujeto a un potencial de dispersión localizado dependiente del tiempo $K(t)$, modelando la interacción de un fotón de rayos X con los electrones de la capa interna en un metal. El sistema se describe mediante un modelo de red de enlaces (tight-binding) 1D con una banda semillena.

1. Marco Teórico:

   • Los autores utilizan el enfoque de espacio métrico de la mecánica cuántica. Definen la adiabaticidad de un sistema que evoluciona desde un estado fundamental inicial $|\phi_0(0)\rangle$ hacia un estado dependiente del tiempo $|\Psi(t)\rangle$ midiendo la distancia entre $|\Psi(t)\rangle$ y el estado fundamental instantáneo $|\phi_0(t)\rangle$.
   • Derivan un límite superior analítico explícito para la distancia NLD ($D_n$) en términos de la distancia de traza ($D_T$) y la probabilidad de solapamiento $|c_0(t)|^2 = |\langle \Psi(t) | \phi_0(t) \rangle|^2$. Específicamente, demuestran que $D_n \leq 2 D_T$, estableciendo una conexión rigurosa entre la métrica de densidad local, computacionalmente más económica, y las métricas del estado cuántico.
   • Se deriva una nueva solución analítica para la evolución temporal de los coeficientes del sistema. Esta solución asume que, para un $K(t)$ continuo, la dinámica está dominada por excitaciones de partícula-hueco simples, omitiendo excitaciones de orden superior. Esto permite un cálculo computacionalmente eficiente de $|c_0(t)|^2$ en un amplio rango de parámetros.

2. Implementación Numérica:

   • Para validar el enfoque analítico y manejar la complejidad total de muchos cuerpos, los autores emplean el método de Grupo de Renormalización Numérica en Espacio Real (eNRG). Esto implica la discretización de la banda metálica y el uso de un procedimiento de suavizado (integración sobre el parámetro de desfase $\theta$) para eliminar oscilaciones no físicas.
   • La evolución temporal se computa utilizando el método de Crank-Nicolson.

3. Parámetros:

   • El estudio varía el potencial de dispersión máximo $\bar{K}$, el tiempo de rampa $T$ y la brecha de energía $\Delta\varepsilon$ (controlada por el tamaño del sistema $N$).
   • Los autores comparan los resultados contra el QAC estándar, el cual es aproximado utilizando únicamente los estados fundamental y de primera excitación.

Resultos Clave

• Validez de los Métodos Basados en Métricas: El estudio confirma que los métodos basados en métricas (distancias de Bures, de traza y NLD) siguen siendo válidos y efectivos para rastrear la adiabaticidad incluso en presencia de la catástrofe de ortogonalidad de Anderson y espectros continuos.
• Precisión de la Solución Analítica: La solución analítica derivada, que considera únicamente excitaciones de partícula-hueco simples, muestra un excelente acuerdo con los resultados numéricos de eNRG. Se mantiene precisa incluso cuando el sistema está lejos de la adiabaticidad (hasta $|c_0(t)|^2 \approx 0.6$). Las desviaciones ocurren solo cuando las excitaciones de dos partículas-hueco se vuelven significativas.
• Comportamiento Universal: Los autores identifican un comportamiento universal para la ocupación del estado fundamental $|c_0(T)|^2$ al final de la rampa del potencial. Este comportamiento depende de la relación entre el tiempo de rampa y la escala de tiempo de la brecha ($T/T_{\Delta\varepsilon}$) y la fuerza del potencial de dispersión.
• Fallo del QAC: El Criterio Adiabático Cuántico estándar falla al predecir o rastrear la adiabaticidad para este sistema. El QAC sugiere incorrectamente que la adiabaticidad mejora con el tiempo o depende linealmente del producto $\bar{K} \cdot T$, mientras que el análisis basado en métricas revela una dependencia más compleja de $\bar{K}$ y de la relación $T/T_{\Delta\varepsilon}$. El QAC subestima la región no adiabática para rampas rápidas y sobreestima la adiabaticidad para rampas lentas en el límite continuo.
• Superioridad de la Distancia de Densidad Local: La distancia NLD no solo rastrea la adiabaticidad, sino que también captura información dinámica que las distancias de Bures y de traza omiten. Específicamente, para $t > T$ (después de que el potencial deja de cambiar), la distancia de Bures permanece constante (ya que la proyección sobre el estado fundamental se conserva en un sistema cerrado), mientras que la distancia NLD disminuye. Esta disminución refleja la relajación del sistema mediante oscilaciones de Friedel y la redistribución de carga hacia la configuración del estado fundamental instantáneo, una dinámica fuera del equilibrio invisible para las métricas de proyección de estado.

Significancia y Reivindicaciones
El artículo afirma que la distancia de Densidad Local Natural es una herramienta robusta, computacionalmente eficiente y experimentalmente accesible para caracterizar la adiabaticidad en sistemas complejos de muchos cuerpos. Al establecer un límite superior explícito que relaciona la distancia NLD con la distancia de traza, los autores proporcionan una justificación teórica para el uso de la densidad local (que es más fácil de computar, por ejemplo, mediante la Teoría del Funcional de la Densidad) como un sustituto (proxy) de la evolución completa del estado cuántico.

El estudio demuestra que, para la fotoemisión de rayos X, la métrica de densidad local proporciona una imagen más completa de la dinámica del sistema que las distidades de estado cuántico tradicionales, particularmente al capturar los procesos de relajación post-quench. Además, el trabajo destaca la insuficiencia del QAC estándar para sistemas con espectros continuos y correlaciones fuertes, abogando por enfoques basados en métricas que no dependan de expansiones perturbativas o condiciones de brecha discreta. La solución analítica derivada ofrece un método rápido y fiable para determinar los regímenes adiabáticos en gases de Fermi con potenciales localizados, válido mucho más allá del límite estrictamente adiabático.