Locally Purified Maximally Mixed States At Scale: Entanglement Pruning and Symmetries

Resumen Técnico: Estados Máximamente Mezclados Purificados Localmente a Escala: Poda de Entrelazamiento y Simetrías

Planteamiento del Problema
Los Operadores de Densidad Localmente Purificados (LPDO, por sus siglas en inglés) son un ansatz de redes de tensores líder para representar estados cuánticos mezclados a escala, ofreciendo una reducción cuadrática en la dimensión virtual respecto a los MPO (Operadores de Producto de Matrices) mientras imponen explícitamente la positividad. Sin embargo, los autores identifican una ineficiencia crítica: la representación del estado máximamente mezclado (el límite de temperatura infinita, $\rho_1$) suele ser subóptima en computaciones prácticas.

Aunque el estado máximamente mezclado es físicamente separable (requiriendo una dimensión de enlace virtual $\chi=1$), los protocolos estándar que involucran canales cuánticos ruidosos (por ejemplo, desfasamiento y bitflip) aplicados a estados entrelazados resultan en un LPDO donde la dimensión de enlace virtual $\chi$ permanece inflada desde el estado puro inicial. Esto ocurre porque los canales de ruido aumentan el índice de mezcla $\kappa$ (que codifica correlaciones clásicas) pero dejan la dimensión de enlace coherente $\chi$ sin cambios. En consecuencia, la representación resultante de un estado separable conlleva "entrelazamiento" espurio en el espacio virtual, lo que conduce a una asignación de memoria innecesaria y sobrecarga computacional. El artículo aborda el desafío de mapear estas representaciones subóptimas de vuelta a su forma de complejidad mínima y óptima ($\chi=1$).

Metodología
Los autores emplean una combinación de herramientas numéricas y analíticas para resolver la suboptimalidad de las representaciones de los LPDO, centrándose específicamente en la "poda de entrelazamiento" de la dimensión de enlace virtual $\chi$.

1. Truncamiento Preservador de la Fidelidad:
   Los autores generalizan el truncamiento estándar de la Descomposición en Valores Singulares (SVD) utilizado en los MPS (Estados de Producto de Matrices). A diferencia de los MPS, donde el error de truncamiento escala linealmente con el corte $\Lambda$, los autores demuestan que para el estado máximamente mezclado, el error de aproximación (medido como infidelidad) es más permisivo. Proponen un protocolo donde los tensores de vecinos cercanos se contraen y se realiza una SVD con un corte $\Lambda$ significativamente más alto (por ejemplo, $0.1$a$0.5$) que los utilizados típicamente en MPS. Esto es seguido por una renormalización tanto en el subespacio de mezcla ($\kappa$) como en el coherente ($\chi$). Al realizar barridos a través de la cadena múltiples veces, esta poda de fuerza bruta reduce exitosamente $\chi$ a 1 manteniendo una fidelidad unitaria con el estado ideal máximamente mezclado.

2. Optimización de Riemann:
   Para optimizar aún más la representación, los autores utilizan la libertad de gauge isométrica inherente al espacio de mezcla virtual del LPDO (subespacio $\kappa$). Encuadran la poda de las correlaciones espurias de $\chi$ como un problema de optimización en la variedad de Stiefel.

   • Funciones Objetivo: Se minimizan dos métricas basadas en la entropía: la segunda entropía de Rényi ($S_{sr}$) de la matriz de densidad reducida de dos sitios y la entropía de Von-Neumann ($S_{vn}$) de los valores singulares en el enlace.
   • Protocolo: El algoritmo busca isometrías que actúan sobre el subespacio $\kappa$ para minimizar estos objetivos entrópicos. La optimación se realiza utilizando resolvedores de región de confianza (Manopt.jl) con diferenciación numérica (Zygote.jl).
   • Sinergia: Los autores encuentran que la optimización de Riemann es particularmente efectiva cuando se combina con el truncamiento. Aunque la optimización por sí sola puede reducir la entropía, puede aumentar inadvertidamente la cola de los valores singulares, estancando la reducción de $\chi$ si se aplica un corte estricto. Por lo tanto, el protocolo acopla la optimización con un paso de truncamiento de corte alto para podar eficazmente $\chi$ hasta 1.

3. Simetría Analítica e Inyectividad:
   Los autores proporcionan una justificación analítica para la existencia de isometrías de desentrelazamiento. Invocan el concepto de simetría débil (invariancia bajo la conjugación $U\rho U^\dagger = \rho$) y inyectividad débil. Dado que el estado máximamente mezclado es débilmente simétrico bajo cualquier unitaria, la acción de una unitaria sobre los índices físicos (que infla artificialmente $\chi$) es equivalente a una isometría en el subespacio virtual $\kappa$.

   • Derivan una expresión de forma cerrada para el desentrelazador: para una unitaria $U$ que actúa en el espacio físico, la isometría equivalente en el subespacio $\kappa$ es $V = e^{-i\phi}U$. Así, aplicar $V^\dagger$ en el espacio virtual invierte exactamente el entrelazamiento espurio inducido por $U$, restaurando la representación óptima de $\chi=1$.

4. Truncamiento $\Delta$-Adaptativo para Estados Parcialmente Mezclados:
   Los autores extienden su análisis a estados parcialmente depolarizados ($\rho_0$) caracterizados por una infidelidad $\Delta$ respecto al estado máximamente mezclado. Introducen una estrategia de truncamiento $\Delta$-adaptativa. Demuestran que a medida que el estado se aleja del límite máximamente mezclado ($\Delta$ aumenta), el corte de truncamiento debe ajustarse. La fidelidad del estado truncado interpola linealmente entre los límites de precisión estándar de MPS (para estados puros, $\Delta=1$) y los límites de truncamiento preservador de la fidelidad (para estados máximamente mezclados, $\Delta=0$).

Resultados Clave

• Representación Óptima: Los autores demuestran con éxito que la representación subóptima del LPDO del estado máximamente mezclado (derivada de circuitos ruidosos) puede transformarse en la forma óptima ($\chi=1, \kappa=2$) utilizando truncamiento de alto corte y optimización de Riemann.
• Eficiencia de Truncamiento: Los valores de corte altos ($\Lambda \approx 0.1 - 0.5$) en el subespacio $\chi$ no degradan la fidelidad para el estado máximamente mezclado, contrario a la intuición derivada de los MPS de estados puros. Esto permite una poda agresiva de $\chi$.
• Rendimiento de Optimización: La optimización de Riemann utilizando $S_{sr}$ converge más rápido y es computacionalmente más eficiente que el uso de $S_{vn}$. Sin embargo, la minimización de $S_{vn}$ produce medidas de infidelidad ligeramente menores.
• Verificación Analítica: La derivación analítica de forma cerrada confirma que los optimizadores numéricos están encontrando efectivamente la inversa de las acciones unitarias que causaron el entrelazamiento espurio, validando el argumento de inyectividad.
• Comportamiento de Escala: Para estados parcialmente mezclados, el corte de truncamiento requerido para mantener una infidelidad específica escala con la distancia hacia el estado máximamente mezclado. El método proporciona una interpolación suave entre el comportamiento de truncamiento estándar de MPS y el comportamiento único de preservación de la fidelidad del estado máximamente mezclado.

Significado y Reivindicaciones
El artículo afirma resolver la suboptimalidad representacional de los LPDO para el estado máximamente mezclado, un límite físicamente relevante como punto de convergencia de algoritmos cuánticos de profundidad polinómica en hardware ruidoso. Al minimizar los recursos (específicamente la dimensión de enlace virtual $\chi$) requeridos para representar estos estados, los autores aseguran que la eficiencia de los algoritmos de redes de tensores para estados mezclados puede incrementarse sustancialmente.

Este trabajo se posiciona como un paso hacia el descubrimiento de los límites fundamentales de escalabilidad de las redes de tensores y su potencial latente para representar la amplia gama de estados cuánticos mezclados accesibles en dispositivos cuánticos de tecnología actual (near-term). Los autores sugieren que estos métodos pueden impactar la ingeniería de qubits, reducir los requisitos de recursos clásicos para la simulación de sistemas cuánticos abiertos (incluyendo circuitos ruidosos) y mejorar la fiabilidad de las estrategias de mitigación de errores basadas en la extrapolación de ruido. El artículo no propone nuevos montajes experimentales, sino que ofrece mejoras algorítmicas para la simulación y el análisis de modelos de computación cuántica existentes.