Resumen Técnico: Un Algoritmo Cuántico Sin Penalizaciones para Encontrar Estados Propios de Energía

Enunciado del Problema

Encontrar los valores propios y los estados propios de Hamiltonianos de muchos cuerpos es un desafío fundamental en física y ciencia computacional. Los métodos clásicos enfrentan una escalabilidad exponencial (almacenamiento $O(2^{2n})$ y tiempo $O(2^{3n})$) a medida que aumenta el número de qubits $n$. Aunque los algoritmos cuánticos ofrecen un camino a seguir, los enfoques existentes tienen limitaciones significativas:

• Eigensolver Cuántico Variacional (VQE): Un enfoque híbrido cuántico-clásico que a menudo sufre de "mesetas estériles" (gradientes que se desvanecen) durante la optimización. Encontrar estados excitados típicamente requiere términos de penalización para suprimir las contribuciones del estado fundamental, lo que aumenta la complejidad y la posibilidad de error.
• Recocido Cuántico: Requiere tiempos de evolución largos para garantizar la adiabaticidad y lucha con pequeños huecos de energía entre el estado fundamental y los estados excitados.
• Evolución en Tiempo Imaginario (ITE) y QITE: Aunque prometedoras, la ITE estándar enfrenta problemas de escalabilidad para realizar operadores no unitarios. La ITE Cuántica (QITE) a menudo depende de la optimización variacional o de términos de penalización para estados excitados, reintroduciendo los desafíos de las mesetas estériles y requiriendo ansatzes complejos.
• Otros Métodos: Técnicas como la Diagonalización por Filtro Cuántico (QFD) y la Estimación de Fase Cuántica Filtrada (FQPE) dependen fuertemente de estados de referencia específicos o funciones de filtro optimizadas.

Los autores identifican la necesidad de un algoritmo totalmente cuántico que pueda encontrar tanto estados fundamentales como excitados sin funciones de penalización, pasos variacionales o bucles híbridos clásico-cuánticos, evitando así las mesetas estériles.

Metodología

El algoritmo propuesto es un método totalmente cuántico y sin penalizaciones basado en el muestreo estocástico de una matriz de densidad de estado mixto combinado con una simulación basada en estados.

1. Descomposición del Hamiltoniano y Mapeo de la Matriz de Densidad

El algoritmo asume que el Hamiltoniano objetivo $H$ puede descomponerse en dos componentes resolubles, $H = h_1 + h_2$, donde los valores propios ($\lambda^j_i$) y los estados propios ($|\psi^j_i\rangle$) de $h_1$ y $h_2$ son conocidos y fáciles de preparar.

• Los autores desplazan y reescalan los valores propios de $h_1$ y $h_2$ para crear una distribución de probabilidad válida.
• Construyen una matriz de densidad de estado mixto $\rho$ equivalente al Hamiltoniano desplazado:
  $$\rho = \frac{1}{C} (H + (c_1 + c_2)I) = \sum_{i,j} p^j_i |\psi^j_i\rangle\langle\psi^j_i|$$
  donde $p^j_i$ son probabilidades derivadas de los valores propios desplazados. Los estados propios de $\rho$ son idénticos a los de $H$.

2. Evolución Estocástica en Tiempo Imaginario (ITE) para el Estado Fundamental

Para encontrar el estado fundamental, el algoritmo realiza ITE sobre $\rho$ utilizando muestreo estocástico:

• Muestreo: En cada paso de tiempo $\eta$, un proyector $P^j_i = |\psi^j_i\rangle\langle\psi^j_i|$ se muestrea de la distribución $\{p^j_i\}$.
• Operador de Evolución: En lugar de implementar el no unitario completo $e^{-\eta \rho}$, el algoritmo implementa el operador $I - \eta P^j_i$ para el proyector muestreado.
• Implementación del Circuito: Esto se logra mediante un circuito unitario con post-selección que involucra un qubit auxiliar. La unitaria $U$ es una puerta controlada-SU(2) de $n$ qubits. Si el sistema está en el estado correspondiente al proyector muestreado, el auxiliar experimenta una rotación; de lo contrario, permanece sin cambios. La post-selección del auxiliar en el estado $|0\rangle$ produce la operación deseada $I - \eta P^j_i$.
• Convergencia: Repetir este proceso impulsa al sistema hacia el estado fundamental de $\rho$ (y por tanto de $H$) sin optimización variacional.

3. Simulación Basada en Estados para Estados Excitados

Para encontrar estados excitados, el algoritmo debe suprimir las contribuciones de los estados de menor energía ya encontrados (por ejemplo, el estado fundamental $P_g$).

• Desafío: Implementar directamente proyectores para estados fundamentales altamente entrelazados es difícil.
• Solución: Los autores utilizan un protocolo de simulación basada en estados. Esto implica un sistema dual de $n$ qubits (auxiliar) preparado en el estado fundamental conocido $P_g$, un qubit de control y el sistema principal.
• Mecanismo: Se aplica una puerta controlada-SWAP entre el sistema principal y el auxiliar, condicionada al qubit de control. Medir el qubit de control en la base $|+\rangle$ implementa efectivamente la evolución $e^{-\eta P_g}$ (aproximada como $I - \eta P_g$), proyectando fuera la componente del estado fundamental.
• Procedimiento de Elevación: Para encontrar el primer estado excitado, la probabilidad de muestrear el proyector del estado fundamental se "eleva" artificialmente (aumenta) en la distribución de la matriz de densidad. Luego, el algoritmo realiza ITE sobre esta distribución modificada. Cuando se muestrea el proyector del estado fundamental, se utiliza la simulación basada en estados para filtrarlo. Este proceso puede iterarse para encontrar estados excitados superiores.

4. Variante para Hamiltonianos Generales

Para Hamiltonianos que no se dividen fácilmente en dos partes resolubles, los autores proponen una variante donde $H$ se trata como una combinación lineal de operadores de Pauli ($H = \sum \alpha_i O_i$). El muestreo estocástico se aplica directamente a estos operadores de Pauli, con implementaciones de circuito similares para los pasos de evolución.

Resultados Clave

Los autores validaron el algoritmo utilizando simulaciones clásicas del Modelo de Ising con Campo Transversal (una cadena 1D con condiciones de frontera periódicas), que sirve como punto de referencia porque puede dividirse en componentes resolubles en la base $Z$ y en la base $X$.

• Estado Fundamental (sistema de 4 qubits): Se logró una fidelidad del 99.9% con el estado fundamental exacto.
• Estados Excitados (sistema de 4 qubits):
  • 1er estado excitado: 98.6% de fidelidad.
  • 2do estado excitado: 96.8% de fidelidad.
  • 3er estado excitado: 96.6% de fidelidad.
• Prueba de Escalabilidad (sistema de 10 qubits): El algoritmo encontró exitosamente el estado fundamental (98.6% de fidelidad) y el primer estado excitado (97.8% de fidelidad). También resolvió el tercer estado propio más bajo (90.2% de fidelidad) a pesar de la casi degeneración entre los dos primeros estados.
• Variante de Pauli: La variante de muestreo de operadores de Pauli también demostró altas fidelidades (97.1%–99.5%) para los primeros cuatro estados propios en el caso de 4 qubits.
• Escalabilidad de Recursos: El análisis teórico indica que el algoritmo escala con $O(n)$ puertas de dos qubits por paso, una mejora significativa sobre los métodos ITE basados en Trotter, que escalan como $O(n^2)$.

Significado y Afirmaciones

El artículo afirma presentar un algoritmo totalmente cuántico que aborda el problema de los estados propios sin las desventajas de los métodos híbridos o variacionales actuales:

1. Sin Penalizaciones: Encuentra estados excitados sin agregar términos de penalización a una función de pérdida, evitando la necesidad de ansatzes complejos y el riesgo de mesetas estériles.
2. Sin Pasos Híbridos: Elimina la necesidad de bucles de optimización clásica, convirtiéndolo en un procedimiento puramente cuántico.
3. Amplia Aplicabilidad: Se argumenta que la suposición de que $H$ puede dividirse en partes resolubles no es restrictiva, cubriendo muchos modelos físicos e incluso problemas completos de QMA.
4. Robustez: El método demuestra la capacidad de encontrar múltiples estados propios de baja energía secuencialmente, incluso en sistemas con niveles de energía casi degenerados.

Los autores reconocen que la tasa de éxito general depende de la post-selección, lo cual puede ser un desafío en los dispositivos cuánticos de escala intermedia ruidosos (NISQ) actuales debido a errores de puerta y decoherencia. Sin embargo, postulan que el método ofrece una vía concreta para utilizar circuitos cuánticos en la resolución de problemas de estados propios intratables clásicamente, sirviendo como una adición fundamental a la caja de herramientas de la computación cuántica.