ENTCALC: Toolkit for calculating geometric entanglement in multipartite quantum systems

Resumen Técnico de "ENTCALC: Herramienta para calcular el entrelazamiento geométrico en sistemas cuánticos multipartitos"

Planteamiento del Problema
Cuantificar el entrelazamiento cuántico en sistemas multipartitos es una tarea computacionalmente exigente. Aunque existen diversas medidas de entrelazamiento (por ejemplo, la entropía de entrelazamiento, la negatividad), estas suelen enfrentar limitaciones: algunas se aplican solo a estados puros, otras no logran detectar ciertos estados entrelazados (como los estados con Transposición Parcial Positiva o PPT) y determinar la separabilidad es un problema NP-duro. Específicamente, calcular la medida geométrica del entrelazamiento para estados multipartitos mixtos generales es difícil porque implica una optimización sobre el conjunto de estados separables. Los enfoques existentes a menudo restringen el análisis a estados puros, se centran únicamente en subsistemas de pares (ignorando la estructura multipartita global) o dependen de testigos de entrelazamiento que no proporcionan una medida cuantitativa. Existe la necesidad de una herramienta que pueda proporcionar estimaciones de alta precisión del entrelazamiento geométrico tanto para estados puros como mixtos multipartitos, particularmente en casos donde las soluciones analíticas no están disponibles.

Metodología
Los autores presentan ENTCALC, un paquete de software implementado tanto en Python (entcalcpy) como en MATLAB (entcalc). La herramienta estima el entrelazamiento geométrico ($E_G$) mediante el cálculo de tanto límites inferiores como superiores, identificando así un intervalo dentro del cual se encuentra el valor real.

1. Definición de Entrelazamiento Geométrico:

   • Para estados puros $|\psi\rangle$, $E_G$ se define como $1 - \max_{|\phi\rangle \in S} |\langle\phi|\psi\rangle|^2$, donde $S$ es el conjunto de estados totalmente separables.
   • Para estados mixtos $\rho$, se define mediante la extensión de techo convexo: $E_G(\rho) = \min_{\{p_i, |\psi_i\rangle\}} \sum p_i E_G(|\psi_i\rangle)$.
   • Una formulación equivalente utiliza la fidelidad cuántica: $E_G(\rho) = 1 - \max_{\sigma \in S} F(\rho, \sigma)$.

2. Algoritmo para Estados Puros:

   • Enfoque SDP: El paquete utiliza una formulación de Programación Semidefinida (SDP) basada en relajaciones PPT (Transposición Parcial Positiva). Para un estado tripartito, el límite inferior se deriva de la maximización del solapamiento con estados PPT. Simultáneamente, se calcula un límite superior con una estimación de error ($4\sqrt{\epsilon}$), donde $\epsilon$ se relaciona con el autovalor más grande del estado PPT óptimo.
   • Descenso de Gradiente: Se utiliza un algoritmo iterativo para encontrar el estado separable más cercano, proporcionando un límite superior. Este método también se emplea como subrutina para estados mixtos.

3. Algoritmo para Estados Mixtos:

   • Límite Superior: Se calcula mediante un algoritmo de descenso de gradiente que refina iterativamente una descomposición separable del estado objetivo para maximizar la fidelidad.
   • Límites Inferiores: El paquete implementa cuatro métodos distintos de límites inferiores, permitiendo a los usuarios equilibrar la precisión frente al costo computacional:
     • Límite 1 (SDP): Basado en la maximización de la fidelidad sobre el conjunto de estados PPT. Es eficiente pero arroja cero para estados entrelazados PPT.
     • Límite 2 (Extensiones k-Simétricas): Mejora la precisión incorporando extensiones $k$-simétricas, creando una jerarquía de relajaciones más ajustadas que convergen al conjunto de estados separables cuando $k \to \infty$.
     • Límite 3 (Complementariedad Pureza-Entrelazamiento): Derivado de la relación entre la termodinámica y el entrelazamiento, utilizando una purificación del estado mixto y optimizando sobre variables auxiliares.
     • Límite 4 (Pureza-Entrelazamiento Global): Similar al Límite 3, pero opera sobre todo el sistema purificado, ofreciendo límites más ajustados a costa de una mayor dimensión computacional.

Contribuciones Clave

• Estimación de Doble Límite: A diferencia de métodos anteriores que a menudo proporcionan solo un límite o un testigo, ENTCALC calcula simultáneamente límites inferiores y superiores. Esto permite la identificación de un intervalo estrecho para el entrelazamiento geométrico real, aumentando significativamente la confianza en la estimación.
• Manejo de Estados Mixtos y PPT: La herramienta es capaz de cuantificar el entrelazamiento en estados mixtos, incluyendo estados entrelazados PPT (entrelazamiento ligado), que son indetectables por medidas basadas en la negatividad.
• Nuevas Extensiones Teóricas: Los autores extienden marcos teóricos existentes para derivar nuevos límites inferiores para el escenario multipartito basados en relaciones de complementariedad pureza-entrelazamiento.
• Implementación de Software: El lanzamiento de paquetes de código abierto tanto en Python como en MATLAB, aprovechando librerías científicas estándar (NumPy, QuTiP, CVXPY, SciPy, CVX, QETLAB), hace que el cálculo de entrelazamiento de alta precisión sea accesible.

Resultados y Aplicaciones
Los autores aplicaron ENTCALC a varios ejemplos representativos, observando que la brecha entre los límites inferior y superior era consistentemente muy pequeña (a menudo $< 10^{-5}$), lo que indica una alta precisión:

1. Estados Entrelazados PPT: El paquete cuantificó con éxito el entrelazamiento geométrico de estados entrelazados PPT de Horodecki $3 \otimes 3$, donde la diferencia entre los límites fue $\approx 5.15 \times 10^{-6}$.
2. Estados de Bajo Rango: Se cuantificó el entrelazamiento para mezclas de estados GHZ y W en sistemas de 3, 4 y 5 cúbits.
3. Entrelazamiento Térmico en Cadenas de Espín:
   • Modelo XX: El paquete detectó entrelazamiento multipartito en estados térmicos de una cadena de Heisenberg XX de tres cúbits, incluso cuando el entrelazamiento de pares era cero. Identificó regiones de entrelazamiento ligado (entrelazamiento PPT) tanto para acoplamientos ferromagnéticos ($J>0$) como antiferromagnéticos ($J<0$).
   • Modelo XXX con Campo Magnético: La herramienta observó firmas de transiciones de fase cuántica mediante discontinuidades en el entrelazamiento geométrico como función de la intensidad del campo magnético y la temperatura inversa.
   • Cadena de Espín Hexagonal: Los autores demostraron que el ajuste de un campo magnético externo puede activar el entrelazamiento entre sitios no adyacentes (por ejemplo, los espines B, D y F en un anillo de 6 cúbits) que no interactúan directamente, un fenómeno no capturado por el análisis de pares.
4. Canales Ruidosos: El paquete cuantificó la robustez de los estados GHZ y W frente al ruido de amortiguamiento de amplitud y de depolarización. Los resultados confirmaron que los estados W son generalmente más robustos contra el amortiguamiento de amplitud en niveles altos de ruido, mientras que los estados GHZ exhiben inicialmente un entrelazamiento más alto a bajos niveles de ruido.

Significancia y Reivindicaciones
El artículo afirma que ENTCALC constituye el primer paquete que permite el cálculo de alta precisión del entrelazamiento geométrico para estados mixtos multipartitos. Al proporcionar intervalos estrechos entre los límites inferiores y superiores, la herramienta ofrece un método numérico fiable para estudios teóricos donde las soluciones analíticas son intratables. Los autores enfatizan que la pequeña brecha entre los límites en sus casos de prueba confirma la precisión de los límites empleados. Posicionan el paquete como un recurso valioso para investigar fenómenos complejos como el entrelazamiento ligado en cadenas de espín, transiciones de fase cuánticas y la dinámica del entrelazamiento en entornos ruidosos, señalando que puede revelar características (como el entrelazamiento multipartito no local) que las medidas de pares pasan por alto. El trabajo no pretende resolver el problema de la separabilidad NP-duro de forma general, sino que proporciona una aproximación práctica y de alta precisión para estados específicos y físicamente relevantes.