Optimización de la Transmisión de Estados Cuánticos en Cadenas de Qubits usando Deep Reinforcement Learning y Algoritmos Genéticos

Comparación de estrategias de control dinámico para transferencia cuántica

¿Qué estrategia de optimización conviene más para controlar un sistema cuántico mediante pulsos magnéticos: Deep Reinforcement Learning o algoritmos genéticos? Este trabajo, publicado en Physica Scripta, compara ambos enfoques para la transmisión de estados cuánticos en cadenas de qubits, midiendo fidelidad y costo computacional en el régimen de control ideal (sin ruido). El resultado es contundente: el algoritmo genético (GA) gana con claridad, alcanzando fidelidades más altas y corriendo entre 2 y 3 órdenes de magnitud más rápido que un DQN. ¿Y si el control no es perfecto? Esa pregunta se aborda en el trabajo complementario sobre robustez frente a ruido dinámico.

En términos simples, la transferencia de estado cuántico (QST) consiste en mover información cuántica de un extremo a otro de una cadena de qubits sin perderla en el camino, algo necesario para conectar distintas partes de un futuro procesador cuántico. Por eso la QST a través de cadenas de espín homogéneas desempeña un papel crucial en la construcción de hardware cuántico escalable. Un protocolo básico de transmisión de estados cuánticos prepara un estado en un qubit y lo transfiere a otro a través de un canal, buscando minimizar el tiempo y evitar pérdida de información. La fidelidad del proceso —qué tan parecido es el estado recibido al que se quería enviar— se mide mediante funciones proporcionales a la probabilidad de transición entre ambos estados.

Tomando como punto de partida el trabajo de Zhang et al. [1], que propone el uso de deep reinforcement learning para controlar la transferencia de estados cuánticos, abordamos este problema de optimización mediante pulsos magnéticos constantes y dos estrategias complementarias: aprendizaje por refuerzo profundo (DRL), donde un agente aprende secuencias de pulsos mediante recompensas, y algoritmos genéticos (GA), que desarrollan soluciones candidatas mediante selección y mutación. Analizamos la eficiencia de ambos métodos, su capacidad para incorporar restricciones físicas y, sobre todo, cómo se comportan cuando el control deja de ser perfecto.

Control Dinámico

En términos concretos, controlar el sistema significa aplicar sobre la cadena una secuencia de pulsos magnéticos, uno después del otro, eligiendo su forma para guiar el estado cuántico de un extremo al otro. El sistema consiste en una cadena de N qubits con Hamiltoniano de intercambio XX que evoluciona bajo Hamiltonianos independientes del tiempo por intervalos. El control se implementa mediante pulsos de campo magnético aplicados en subintervalos de tiempo discretos:

  • En cada subintervalo se aplica una configuración Bk,j de campos elegida de un conjunto de configuraciones posibles
  • Se busca encontrar una secuencia óptima de 5N compuertas de un qubit que transmitan el estado desde el primer al último qubit
  • Las configuraciones posibles de campos definen el espacio de acciones

En este trabajo exploramos dos espacios de acciones diferentes:

  • Espacio de acciones original de Zhang et al.: Un conjunto de 16 acciones que aplican campos magnéticos predefinidos sobre sitios específicos de la cadena, tal como fue propuesto en el trabajo original.
  • Una acción por sitio: Un espacio de acciones simplificado donde cada acción consiste en encender un campo magnético sobre un único sitio de la cadena (N acciones posibles), permitiendo mayor granularidad en el control.

Estrategias de Control

Algoritmos Genéticos (GA)

El algoritmo genético se implementó usando PyGAD con evaluación paralela de poblaciones mediante PyTorch:

  • La población inicial consiste en secuencias aleatorias de 5N acciones
  • Se simula la aplicación de cada secuencia y se les asigna una fitness proporcional a la calidad de transmisión, evaluada sobre toda la evolución temporal
  • Se seleccionan las secuencias de mayor fidelidad como padres (selección steady-state) y se genera una nueva población mediante cruzamiento uniforme y mutación por intercambio (swap)
  • La estructura y simpleza del algoritmo genético permiten que sea fácil de paralelizar en GPU
Hiperparámetro del GA Valor
Máximo de generaciones1000
Tamaño de la población4096
Saturación (generaciones sin mejora)30
Tipo de mutaciónswap
Probabilidad de mutación0.99
Genes mutadosN
Tipo de selección de padres'sss' (steady-state)
Padres apareados / elitismo409 / 409
Tipo y probabilidad de cruzamientouniforme, 0.8

Tabla I. Hiperparámetros del algoritmo genético utilizados para obtener los resultados presentados en el paper.

Un barrido sistemático de los parámetros físicos del control —el campo magnético h y el paso temporal Δt— sobre una cadena de N=64 qubits mostró que el GA encuentra secuencias con probabilidad de transmisión P(t)>0.99 en una región amplia de ese espacio de parámetros, sin necesidad de un ajuste fino.

Dependencia de la probabilidad de transmisión del GA respecto de los parámetros físicos del control

Probabilidad de transmisión del GA según los parámetros físicos del control (h, Δt), para N=64.

Deep Reinforcement Learning (DRL)

El agente DRL se implementó como una Deep Q-Network (DQN) con replay de experiencia priorizado usando TensorFlow, adaptado de la implementación de Zhang et al.:

  • El espacio de estados codifica la parte real e imaginaria de cada coeficiente del vector de estado de la cadena, s(t) = (Re{ψ₁},...,Re{ψN},Im{ψ₁},...,Im{ψN})
  • Se entrena el agente con una política ε-greedy: explora acciones al azar con probabilidad ε y explota la mejor acción conocida el resto del tiempo, minimizando el error cuadrático entre la red Q y una red objetivo periódicamente actualizada
  • La recompensa es una función escalonada de la probabilidad de transmisión: crece suavemente para P(t)≤0.8, se acelera hasta P(t)=1-χ, y salta a un valor muy alto (2500) apenas se supera ese umbral, empujando al agente a buscar transmisiones casi perfectas
  • Inyección de ruido durante el entrenamiento episódico (domain randomization, es decir, variar aleatoriamente las condiciones de entrenamiento para que el agente generalice mejor) para mejorar la robustez frente a control imperfecto
  • El desempeño de DRL depende fuertemente de la elección de hiperparámetros y condiciones de entrenamiento
Hiperparámetro del DQN (línea base, Zhang et al.) Valor
Tamaño de minibatch32
Tasa de aprendizaje0.01
Factor de descuento (γ)0.95
Tamaño de la memoria de replay40000
Neuronas (1ª / 2ª capa oculta)120 / 40
Tasa ε-greedy1 → 0.01 (decay 0.0001)
Número de episodios50000
Período de aprendizaje (t')5 timesteps
Período de reemplazo de la red objetivo (C)200 episodios

Tabla III. Hiperparámetros de la implementación original del DQN, adaptada de Zhang et al. [1], usada como línea base de comparación.

Resultados sin Ruido: GA vs DRL

GA: acciones por sitio vs acciones de Zhang

Comparación de conjuntos de acciones del GA

(a) Probabilidad de transmisión máxima obtenida por el GA en función de la longitud de la cadena, usando el conjunto de acciones por sitio (rojo) y el conjunto propuesto por Zhang et al. (azul). Los cuadrados son el mejor valor sobre 30 poblaciones iniciales, los círculos el promedio de esos máximos. (b) Evolución temporal de P(t) para N=32 con ambos conjuntos de acciones, comparada con la evolución natural sin control (verde).

DQN original: fidelidad vs longitud de cadena

Fidelidad del DQN original

(a) Probabilidad de transmisión máxima obtenida por el DQN (hiperparámetros originales de Zhang et al.) en función de la longitud de cadena, sobre 30 realizaciones de los pesos iniciales de la red. (b) Evolución temporal de P(t) para N=32 en cada realización (negro), comparada con la evolución natural (verde).

  • Con el conjunto de acciones por sitio, el GA mantiene P(t)≳0.99 para todas las longitudes de cadena estudiadas; con el conjunto de Zhang et al., su desempeño decae con N
  • El DQN original iguala al GA solo para cadenas muy cortas (N≈8-12); a partir de ahí la brecha se agranda claramente a medida que crece la longitud de la cadena
  • Ambos métodos encuentran secuencias que llevan la transmisión a su máximo en tiempos cercanos al límite cuántico de velocidad (QSL) —el tiempo mínimo teórico que la dinámica permite para completar la transferencia—, pero solo el GA lo hace de forma consistente en calidad

Comparación de Desempeño y Costo Computacional

Ambos métodos se evaluaron en términos de fidelidad máxima, desempeño promedio y costo computacional.

Fidelidad: GA vs DRL

Comparación GA vs DRL

Comparación directa entre el GA con acciones por sitio (rojo) y el DQN original (negro): probabilidad de transmisión máxima vs longitud de cadena (a) y evolución temporal de P(t) para N=32 (b).

Tiempos de CPU

Tiempos de CPU

Tiempos de ejecución requeridos por GA (rojo) y DQN (negro) para alcanzar sus criterios de convergencia, en escala logarítmica.

Observaciones clave:

  • El algoritmo genético permite obtener secuencias que alcanzan fidelidades más altas que el DQN original, salvo para cadenas muy cortas donde ambos son comparables
  • Para cadenas de más qubits es conveniente cambiar el conjunto de acciones y usar el esquema por sitio
  • El algoritmo DQN emplea tiempos de ejecución entre 2 y 3 órdenes de magnitud mayores que los requeridos por el GA; esa diferencia de costo es la que permitió automatizar la búsqueda de hiperparámetros y probar múltiples funciones de fitness y conjuntos de acciones para el GA, algo inviable de replicar de la misma forma con el DQN

Distribución de Acciones en las Secuencias del GA

Para entender qué hace el GA en la práctica, se analizó, sobre 1000 secuencias exitosas (P(t)≳0.99) con el esquema por sitio, qué tan seguido se actúa sobre cada sitio de la cadena, para distintas longitudes:

N=8

Distribución de acciones N=8

Frecuencia de uso de cada sitio (0 = sin acción).

N=16

Distribución de acciones N=16

Frecuencia de uso de cada sitio (0 = sin acción).

N=32

Distribución de acciones N=32

Frecuencia de uso de cada sitio (0 = sin acción).

  • La acción "cero" (no actuar) se vuelve menos relevante a medida que crece la longitud de la cadena: para lograr una transferencia exitosa hace falta accionar compuertas de fase no triviales la mayor parte del tiempo
  • Los sitios centrales de la cadena se usan más que el resto, excepto por los extremos; el sitio penúltimo destaca con una frecuencia notablemente mayor a la de los demás, consistente con el "congelamiento" de la excitación en el último sitio que se observa en las curvas P(t)

GA para Cadenas Largas

GA para cadenas largas

Probabilidad de transmisión vs longitud de cadena, hasta N=128. Los cuadrados son el mejor valor sobre M poblaciones iniciales (alcanza el objetivo de 0.99 hasta N≤60); los círculos son el promedio de los máximos, que decae más rápido. Aun así, el GA mantiene probabilidades máximas mayores a 0.97 en cadenas de 128 qubits.

La dificultad de la tarea crece con la longitud de la cadena: manteniendo fijo el tamaño de la población inicial, el criterio de saturación termina deteniendo el algoritmo antes de alcanzar el objetivo para N>60. Mejorar estos resultados para cadenas muy largas requeriría, sencillamente, más individuos y más poblaciones iniciales.

Conclusiones

  • Es sorprendente lo bien que el GA entrega secuencias de control capaces de llevar la transmisión a valores muy altos sobre un rango amplio de longitudes de cadena. Aunque cadenas de cientos de qubits todavía son poco realistas, los resultados muestran que sistemas cuánticos con dimensión de Hilbert de uno o dos cientos son controlables mediante secuencias generadas por GA.
  • A pesar de explorar varias funciones de recompensa y conjuntos de acciones, el DQN no alcanza la calidad de transmisión del GA, ni siquiera en cadenas cortas, y la diferencia empeora con la longitud. Una posible explicación es que la fitness del GA evalúa el efecto de la secuencia a lo largo de toda la evolución temporal, mientras que el DQN recompensa acciones locales —una estrategia que funciona en muchos otros problemas de DRL, pero que aquí no alcanza para igualar al GA. No se encontró un patrón reconocible en las acciones usadas por ambos métodos que sugiera una forma práctica de sesgar la búsqueda hacia mejores políticas.
  • El trabajo deja una línea abierta para el GA: introducir criterios de fitness "suavizados" que limiten la variabilidad de las secuencias tipo escalón (todavía sin una propuesta convincente). Para el DQN, mejorar su calidad en este régimen ideal requeriría una exploración más efectiva del espacio de acciones junto con una política multi-paso cuya recompensa contemple que la probabilidad de transmisión permanece despreciable durante la mayor parte de la evolución, en línea con el límite cuántico de velocidad. ¿Y qué pasa cuando el control deja de ser ideal? Esa pregunta —y por qué cambia por completo el panorama entre GA y DQN— se desarrolla en el trabajo complementario sobre robustez frente a ruido dinámico.

Charla

Charla sobre este trabajo, "Aprendizaje automático aplicado a la transmisión de estados cuánticos":

Tecnologías

Python
PyGAD
PyTorch
TensorFlow
Optuna
NumPy

Estado del Trabajo

Este trabajo fue desarrollado como parte de mi investigación doctoral en el Instituto de Física Enrique Gaviola (CONICET-UNC) y la FaMAF-UNC. El artículo fue publicado el 26 de junio de 2026 en Physica Scripta (vol. 101, 265101) y representa una contribución original al campo del control cuántico mediante técnicas de optimización computacional. El estudio de robustez frente a ruido dinámico se desarrolló en colaboración y continúa en una página aparte.

Referencias

[1] X.-M. Zhang et al., Physical Review A 97, 052333 (2018). DOI: 10.1103/PhysRevA.97.052333

[2] S. Perón Santana, A. Fiuri, M. Domínguez y O. Osenda, "The pros and cons of using deep reinforcement learning or genetic algorithms to design control schemes for quantum state transfer on qubit chains", Physica Scripta 101, 265101 (2026). DOI: 10.1088/1402-4896/ae7cc1

Links