Engineered Dissipation for Thermal-State Tracking Across Quantum Criticality

Resumen Técnico: Disipación Diseñada para el Seguimiento de Estados Térmicos a través de la Criticalidad Cuántica

Planteamiento del Problema
La preparación de estados de equilibrio en sistemas de muchos cuerpos cuánticos es un desafío central para las tecnologías cuánticas, que abarcan desde el recocido cuántico (quantum annealing) hasta la computación adiabática. Los protocolos estándar dependen del accionamiento adiabático, donde un parámetro de control se varía lentamente para mantener al sistema en el autoestado instantáneo del Hamiltoniano. Sin embargo, cerca de un punto crítico cuántico (QCP), la brecha de energía se cierra, lo que provoca que la escala de tiempo adiabática requerida diverja. En consecuencia, el accionamiento en tiempo finito genera inevitablemente excitaciones no adiabáticas y defectos, un fenómeno descrito por el mecanismo de Kibble–Zurek (KZ).

Si bien el acoplamiento con el entorno (disipación) suele verse como una fuente de decoherencia, también puede ser un recurso para la ingeniería de estados. Los enfoques existentes se han centrado mayoritariamente en cómo un entorno natural y fijo modifica la dinámica crítica o en estrategias de control coherente (por ejemplo, atajos hacia la adiabaticidad) para suprimir excitaciones. Existe una brecha significativa en el diseño sistemático de tasas de transición disipativas markovianas para dirigir activamente un sistema a lo largo de una trayectoria térmica prescrita de tiempo finito durante el cruce de un QCP. Específicamente, los métodos actuales no abordan plenamente el requisito dual de controlar las excitaciones no adiabáticas y, simultáneamente, generar la redistribución de poblaciones necesaria para seguir una trayectoria térmica cambiante.

Metodología
Los autores desarrollan un marco disipativo markoviano diseñado para seguir un estado de Gibbs instantáneo, $\rho^\star(t)$, durante un rampado de tiempo finito del Hamiltoniano del sistema $H_S(t)$.

1. Marco Teórico:

   • El estado objetivo se define como el estado de Gibbs instantáneo $\rho^\star(t) = e^{-\beta H_S(t)} / Z(t)$.
   • La evolución temporal de este estado objetivo requiere tanto una rotación unitaria de la base de autoestados (gestionada mediante el accionamiento contra-diabático) como una redistribución de las probabilidades de ocupación (poblaciones).
   • Dado que la evolución unitaria no puede cambiar los autovalores, los autores construyen un término no unitario disipativo, $D[\rho]$, dentro de la ecuación maestra de Gorini-Kossakowski-Sudarshan-Lindblad (GKSL).
   • La construcción central implica la ingeniería de corrientes de transición entre autoestados de energía instantáneos. La corriente de población neta está restringida para coincidir con la derivada temporal de las poblaciones objetivo, $\dot{\lambda}_n(t)$.
   • Los autores derivan una familia de tasas de transición no negativas, $W_{ij}(t)$, que satisfacen esta restricción. Esta familia incluye un término de corriente común arbitrario, $Q(t) \geq 0$, lo que permite un control adicional sobre la amortiguación de la coherencia sin alterar la transferencia neta de población. La realización mínima ($Q=0$) activa solo la corriente direccional necesaria (excitación o relajación) dictada por la trayectoria objetivo.

2. Aplicación al Modelo de Ising de Campo Transverso (TFIM):

   • El marco se aplica al TFIM unidimensional impulsado linealmente a través de sus puntos críticos cuánticos ($g/J = \pm 1$).
   • Utilizando las transformaciones de Jordan-Wigner y de Fourier, el modelo se descompone en sectores de momento de Bogoliubov-de Gennes (BdG) independientes, cada uno actuando como un sistema efectivo de dos niveles.
   • El disipador diseñado se aplica modo a modo. Las tasas de transición se calculan basándose en la brecha de energía de cuasipartícula instantánea y en las poblaciones térmicas prescritas para cada modo $k$.

3. Análisis Comparativo:

   • El rendimiento de la dinámica diseñada se evalúa frente a otros dos escenarios: (1) evolución puramente unitaria (accionamiento coherente) y (2) acoplamiento a un reservorio térmico natural (baño Ohmic con temperatura fija).
   • Las métricas incluyen la dinámica de la población de estados excitados, la desviación de la excitación térmica (distancia del estado de Gibbs objetivo), la supresión de la coherencia generada dinámicamente y la escala de los errores residuales con la tasa de accionamiento.

Resultados Clave
Las simulaciones numéricas para el TFIM revelan lo siguiente:

• Supresión de la Desviación de la Excitación Térmica: El protocolo diseñado reduce significativamente la desviación del estado real respecto a la trayectoria térmica prescrita en comparación con la dinámica unitaria y la del baño natural. Mientras que las evoluciones unitarias y de baño natural exhiben desviaciones persistentes tras cruzar el punto crítico, la dinámica diseñada impulsa activamente las poblaciones de vuelta hacia el estado objetivo.
• Comportamiento de Escala:
  • Dinámica Unitaria: Sigue el escalamiento convencional de Kibble–Zurek en el régimen de tasa intermedia, donde la desviación de la excitación térmica residual escala como $\delta^{nth}_{ex} \sim (v/v_{ch})^{1/2}$. A tasas muy lentas, los efectos de tamaño finito modifican este escalamiento.
  • Dinámica Diseñada: La desviación residual es fuertemente suprimida y exhibe un escalamiento sustancialmente más pronunciado, aproximándose a $\delta^{nth}_{ex} \sim (v/v_{ch})^2$ sobre un amplio rango de tasas de accionamiento. Esto indica un comportamiento de accionamiento lento cuadrático en lugar del comportamiento de raíz cuadrada de KZ.
  • Baño Natural: Exhibe un comportamiento de cruce dependiente de la competencia entre la escala de tiempo del rampado y la escala de tiempo de relajación del baño fijo, fallando en adaptarse al cambio del objetivo térmico.
• Control de Coherencia: La dinámica diseñada suprime naturalmente las coherencias fuera de la diagonal generadas durante el accionamiento. Además, el parámetro libre $Q(t)$ permite una mejora sintonizable de la amortiguación de la coherencia (mediante el aumento de las tasas de transición bidireccionales) sin perturbar el seguimiento de la población prescrito.
• Robustez: El protocolo mantiene su efectividad cuando el sistema es impulsado a través de ambos puntos críticos ($g/J = 1$ y $g/J = -1$), evitando la acumulación de errores que típicamente plagan las dinámicas de baño natural en traversos multicríticos.

Significancia y Reivindicaciones
El artículo sostiene que la disipación diseñada ofrece una ruta distinta para controlar la preparación de estados térmicos en sistemas críticos cuánticos impulsados. La significancia principal reside en la capacidad de remodelar activamente el acercamiento en tiempo finito a una trayectoria térmica, en lugar de depender pasivamente de un entorno fijo o de un control puramente coherente.

Reivindicaciones clave:

1. Direccionamiento Activo: A diferencia de los baños naturales, las corrientes de transición diseñadas se construyen directamente a partir de la trayectoria térmica prescrita, permitiendo que el sistema siga activamente el estado objetivo.
2. Escalamiento Mejorado: El protocolo cambia la dependencia de la tasa finita de los errores residuales desde el comportamiento estándar de Kibble–Zurek hacia un escalamiento cuadrático más pronunciado, lo que sugiere una mayor eficiencia para la preparación de estados en tiempo finito.
3. Control Dual: El marco proporciona un control simultáneo sobre la redistribución de poblaciones (seguimiento del objetivo térmico) y la amortiguación de la coherencia (mediante el parámetro $Q(t)$).
4. Viabilidad Experimental: Los autores señalan que los avances recientes en ingeniería de reservorios (por ejemplo, en circuitos superconductores e iones atrapados) proporcionan los ingredientes necesarios, tales como el bombeo incoherente programable y selectivo de energía y tasas de disipación sintonizables, para potencialmente realizar estos protocolos. Sin embargo, reconocen que mapear las tasas resueltas por sector de momento a controles de espín locales en sistemas de muchos cuerpos presenta desafíos experimentales relacionados con la no localidad y las restricciones de ancho de banda.

El trabajo concluye que, si bien la construcción actual se formula para sectores cuadráticos independientes, establece un marco fundamental para la preparación de estados disipativos que podría extenderse a sistemas de muchos cuerpos interactuantes y optimizarse bajo restricciones específicas de recursos experimentales.