Full Quantum Work Statistics for Non-Homogeneous Many-Body Systems

Este artículo establece un marco de primeros principios utilizando la teoría del funcional de la densidad dependiente del tiempo térmico para calcular las estadísticas completas del trabajo cuántico y los momentos del trabajo disipado en sistemas de muchos cuerpos con interacción, demostrando su poder predictivo al analizar el cruce de aislante de Mott a aislante de banda dentro del modelo de Hubbard.

Lo esencial

La visión general: Midiendo el "costo" de sacudir un sistema cuántico

Imagina que tienes una pista de baile muy compleja y concurrida (un sistema cuántico) donde todos se toman de las manos y se mueven al unísono. Esto es un "sistema de muchos cuerpos". Ahora, imagina que de repente empujas la pared de la sala o cambias el tempo de la música (una fuerza externa). Los bailarines tropezarán, chocarán entre sí y, finalmente, se asentarán en un nuevo ritmo.

La energía perdida durante este tropiezo —la "fricción" de la pista de baile— se llama trabajo disipado. En el mundo cuántico, esto no es solo un deslizamiento suave; es un evento caótico y agitado lleno de fluctuaciones aleatorias.

Este artículo presenta un nuevo mapa de alta precisión (un marco matemático) para predecir exactamente cuánta energía se pierde y qué tan caótica será esa pérdida, sin necesidad de simular a cada bailarín individualmente.

El problema: La "caja negra" del caos cuántico

Durante mucho tiempo, los científicos tuvieron dos formas de estudiar estos sistemas:

1. La forma fenomenológica: Adivinaban cómo reaccionaría el sistema basándose en reglas generales, como decir: "Normalmente se calienta cuando lo empujas". Esto es como adivinar el clima mirando el cielo sin un termómetro. Es útil, pero no muy preciso.
2. La forma exacta: Intentaban calcular el movimiento de cada partícula individual. Para un sistema con miles de millones de partículas, esto es como intentar contar cada grano de arena en una playa mientras sopla un huracán. Es computacionalmente imposible.

Los autores buscaban una solución "Goldilocks" (punto medio ideal): un método que fuera lo suficientemente preciso para ver los detalles, pero lo suficientemente simple como para poder ejecutarse realmente en una computadora.

La solución: El truque de las "sombras chinescas"

Los autores utilizaron una técnica llamada Teoría del Funcional de la Densidad dependiente del tiempo térmico (thTDDFT).

Piensa en el sistema cuántico real y complejo como un espectáculo de sombras chino gigante e intrincado con miles de marionetas interactuando. Es demasiado difícil rastrear cada cuerda y cada articulación.

• El truque: En lugar de rastrear las marionetas reales, crean un espectáculo de sombras mucho más simple (es un sistema de partículas que no interactúan), pero que está diseñado matemáticamente para proyectar la misma sombra exacta (densidad) en la pared que el sistema real y complejo.
• El beneficio: Al estudiar la sombra simple, pueden descubrir exactamente qué está haciendo el sistema complejo. No necesitan conocer los secretos de cada una de las interacciones; solo necesitan saber cómo se mueve la "sombra".

El descubrimiento clave: Dividir la "fricción"

El artículo hace una distinción inteligente entre dos tipos de "fricción" o pérdida de energía:

1. La parte "adiabática" (El estiramiento lento): Imagina estirar lentamente una banda elástica. Incluso si lo haces perfectamente lento, la banda ofrece resistencia porque su forma está cambiando. Esta es la pérdida de energía debida al cambio de forma del sistema, no debido al caos.
2. La parte "no adiabática" (El chasquido repentino): Imagina que esa banda elástica se rompe de repente. La pérdida de energía aquí proviene de los sacudones y transiciones repentinos y caóticos.

Los autores desarrollaron una forma de separar estas dos. Demostraron que la parte "caótica" (no adiabática) está directamente vinculada a cómo responde el sistema a un toque rápido (una "función de relajación"). Al usar su método de "sombras chinescas", pueden calcular esta función de respuesta a partir de los primeros principios (leyes básicas de la física) en lugar de adivinar.

La prueba: La pista de baile del "Modelo de Hubbard"

Para demostrar que su mapa funciona, lo probaron en un modelo teórico famoso llamado modelo de Hubbard.

• La configuración: Imagina una línea de bailarines (electrones) en una cuadrícula. Pueden saltar al siguiente lugar, pero si dos bailarines intentan ocupar el mismo lugar, reciben un "choque" (repulsión).
• El experimento: Aplicaron un empuje "escalonado" (empujando a los bailarines en números impares hacia un lado y a los de números pares hacia el otro).
• El resultado: A medida que cambiaban la fuerza del empuje y la temperatura, el sistema cambiaba entre diferentes "estados de la materia":
  • Aislante de Mott: Los bailarines están atrapados en su lugar porque tienen miedo de chocar con sus vecinos.
  • Aislante de Banda: Los bailarines están atrapados porque el suelo mismo está inclinado.
  • Aislante de Orden de Enlace: Un punto medio extraño donde los bailarines se emparejan en un patrón específico.

Los autores descubrieron que su método podía ver claramente las "firmas" de estas diferentes fases en la pérdida de energía. Por ejemplo, justo en el límite donde el sistema cambia de una fase a otra, la "fricción" (pérdida de energía) aumentaba drásticamente. Esto confirmó que su método puede detectar cambios sutiles en el mundo cuántico simplemente midiendo cuánta energía se desperdicia.

Por qué esto es importante

Este artículo no inventa una nueva batería o un nuevo chip de computadora. En su lugar, proporciona una nueva herramienta de medición.

• Antes: Los científicos tenían que adivinar cómo se comportarían los sistemas cuánticos cuando se les presionaba, o tenían que esperar a que las supercomputadoras fallaran al intentar calcularlo.
• Ahora: Tienen una receta fiable, basada en "primeros principios", para calcular exactamente cuánta energía se pierde y cómo fluctúa el sistema, incluso en sistemas cuánticos complejos y concurridos.

Cierra la brecha entre la realidad desordenada de las partículas que interactúan y las matemáticas limpias y resolubles de los sistemas de "sombras", permitiendo a los científicos predecir el costo termodinámico de los procesos cuánticos con alta precisión.

Resumen técnico

Resumen Técnico: Estadísticas Completas del Trabajo Cuántico para Sistemas de Muchos Cuerpos No Homogéneos

Planteamiento del Proble de Investigación
El artículo aborda el desafío de caracterizar la termodinámica fuera del equilibrio de sistemas cuánticos de muchos cuerpos interactuantes, centrándose específicamente en las estadísticas del trabajo disipado durante protocolos de conducción de tiempo finito. Si bien la teoría de respuesta lineal generalizada (GLRT) ha establecido que la distribución completa del trabajo disipado puede expresarse mediante una función de relajación, la aplicación práctica de este marco se ha visto limitada. En ausencia de datos microscópicos, la función de relajación suele modelarse de forma fenomenológica, sacrificando la precisión predictiva. Además, los enfoques existentes tienen dificultades para proporcionar una ruta desde los primeros principios para reconstruir esta función en sistemas no homogéneos y correlacionados, donde la coherencia cuántica y las correlaciones de muchos cuerpos juegan un papel crítico. Los autores pretenden cerrar la brecha entre la teoría del funcional de la densidad térmica (thDFT) y las estadísticas del trabajo fuera del equilibrio para proporcionar un marco predictivo y microscópico para estos sistemas.

Metodología
Los autores desarrollan un formalismo basado en la Teoría del Funcional de la Densidad Dependiente del Tiempo Térmica de Respuesta Lineal (LR-thTDDFT). El núcleo del enfoque consiste en derivar la función de relajación, $\psi(t)$, directamente de los primeros principios microscópicos en lugar de recurrir al modelado fenomenológico.

1. Marco Teórico:

   • Las estadísticas del trabajo se analizan dentro del marco de la GLRT, donde todos los cumulantes de la distribución del trabajo disipado son funcionales bilineales de la derivada del protocolo de conducción, ponderados por la transformada de Fourier de la función de relajación, $\tilde{\psi}(\omega)$.
   • La función de relajación se descompone en componentes adiabáticos ($\psi_{ad}$) y no adiabáticos ($\psi_{na}$). La parte adiabática contabiliza la irreversibilidad debida a la deformación temporal del espectro de energía (independiente de la tasa de conducción), mientras que la parte no adiabática captura transiciones disipativas genuinas entre niveles de energía.
   • Utilizando los teoremas de Runge-Gross y van Leeuwen extendidos a temperaturas finitas, los autores establecen una correspondencia biunívoca entre el potencial externo dependiente del tiempo y la densidad del sistema.
   • El componente no adiabático de la función de relajación está vinculado explícitamente a la parte imaginaria de la función de respuesta densidad-densidad, $\chi_{ij}^\beta(\omega)$, del sistema interactuante.
   • La función de respuesta interactuante se obtiene de la función de respuesta de Kohn-Sham (KS) no interactuante mediante el formalismo de Gross-Kohn, que incorpora el kernel de Hartree-intercambio-correlación (Hxc) a temperatura finita, $f_{Hxc}^\beta(\omega)$.

2. Aproximaciones:

   • Para hacer el problema tratable, los autores emplean la Aproximación de Densidad Local Adiabática Térmica (thALDA). En esta aproximación, el kernel Hxc es instantáneo en el tiempo (independiente de la frecuencia) y local en el espacio.
   • Esta elección descuida los efectos de memoria y las características excitónicas más allá de las excitaciones de partícula-hueco individuales, pero retiene la simplicidad y la eficiencia computacional de los funcionales locales.

3. Sistema Modelo:

   • El marco se aplica al modelo de Hubbard unidimensional sujeto a un potencial externo estacado. Este sistema se elige por su rico diagrama de fases, que incluye fases de aislante de Mott (MI), aislante de banda (BI) y aislante de orden de enlace (BOI).
   • El estudio analiza el sistema a través del cruce de aislante de Mott a aislante de banda, variando la fuerza de interacción ($U$) y la amplitud del potencial estacado ($v_0$).

Contribuciones Clave

• Derivación desde Primeros Principios: El artículo establece una ruta directa para calcular la función de relajación a partir de densidades microscópicas y potenciales de Kohn-Sham, eliminando la necesidad de un modelado fenomenológico del kernel de disipación.
• Descomposición Adiabática/No Adiabática: Los autores proporcionan una interpretación física clara de la función de relajación al separar las contribuciones adiabáticas de las no adiabáticas, vinculando la primera a la susceptibilidad estática y la segunda a la respuesta dinámica.
• Estadísticas de Trabajo Microscópicas: El trabajo demuestra cómo calcular la distribución completa del trabajo cuántico (específicamente sus primeros pocos cumulantes) para sistemas interactuantes de muchos cuerpos utilizando thTDDFT, cerrando la brecha entre la DFT de equilibrio y la termodinámica fuera del equilibrio.
• Benchmarking: El enfoque se compara con la diagonalización exacta para el dímero de Hubbard (modelo de Hubbard de dos sitios), mostrando una buena concordancia cualitativa y cuantitativa incluso en regímenes fuertemente correlacionados.

Resultados

• Firmas de Fase: En el modelo de Hubbard, el primer momento del trabajo disipado (producción de entropía irreversible) exhibe un pico pronunciado que traza el cruce entre los regímenes MI y BI, proporcionando una firma termodinámica de la región BOI. Esta mejora se atribuye a la mayor sensibilidad de las densidades térmicas a las variaciones de los parámetros de control cerca de fases competidoras.
• Comportamiento de los Cumulantes: Se analizaron los tres primeros cumulantes de la distribución del trabajo disipado como funciones de la duración de la perturbación o quench ($\tau$).
  • El tercer cumulante se anula en el límite adiabático ($\tau \to \infty$), lo que es consistente con que la distribución se aproxime a una forma gaussiana donde solo permanecen los dos primeros momentos.
  • El primer y segundo cumulante retienen firmas de la transición de fase (BOI a BI) incluso en el régimen adiabático, impulsados por la contribución adiabática al trabajo irreversible.
  • El sistema en el régimen de aislante de banda alcanza la adiabaticidad a tasas de quench más rápidas en comparación con las fases de Mott o BOI.
• Relación de Incertidumbre Termodinámica (TUR): El estudio confirma que el límite inferior de la TUR ($F_W \geq 2/\beta$) no se satura en presencia de fluctuaciones cuánticas, incluso en el régimen de respuesta lineal, distinguiendo el comportamiento cuántico de las expectativas clásicas.
• Precisión: Los resultados de LR-thTDDFT para el dímero de Hubbard muestran errores relativos medios de aproximadamente $1.36 \times 10^{-2}$, $2.76 \times 10^{-2}$ y $1.76 \times 10^{-2}$ para el primer, segundo y tercer cumulante, respectivamente, en comparación con la diagonalización exacta.

Significancia y Reivindicaciones
El artículo afirma proporcionar un marco microscópico y transferible para la termodinámica cuántica en sistemas correlacionados. Al conectar las estadísticas del trabajo fuera del equilibrio directamente con las densidades térmicas y los potenciales de Kohn-Sham, el enfoque ofrece una herramienta predictiva que evita las suposiciones fenomenológicas. Los autores enfatizan que, si bien la implementación actual utiliza la thALDA (que tiene limitaciones conocidas en regímenes fuertemente correlacionados dominados por interacciones no locales), el formalismo en sí es exacto dentro de la respuesta lineal y permite mejoras sistemáticas mediante kernels de intercambio-correlación más refinados, no locales o dependientes de la frecuencia.

El trabajo destaca la utilidad de los métodos basados en DFT para investigar la termodinámica fuera del equilibrio en sistemas grandes (por ejemplo, gases atómicos ultrafríos, materiales de baja dimensionalidad) donde las técnicas de muchos cuerpos computacionalmente exigentes (como DMRG o QMC) pueden ser prohibitivas. Los autores señalan modestamente que, aunque el estudio actual captura las firmas de fase y las tendencias generales, no intenta extraer leyes de escala universales (como la escala de Kibble-Zurek) cerca de los puntos críticos, lo que requeriría kernels más avanzados y un análisis dedicado. La contribución principal es el establecimiento del vínculo teórico y la demostración de su viabilidad para extraer la irreversibilidad termodinámica a partir de modelos microscópicos.