Exceptional-Point-Anchored Variational Quantum Eigensolver for Non-Hermitian Many-Body Phase Diagrams: Bridging Skin-Effect Topology and Entanglement Criticality on NISQ Hardware

Resumen Técnico: VQE Basado en Puntos Excepcionales para Diagramas de Fase de Muchos Cuerpos No Hermíticos

Planteamiento del Problema
El estudio de los sistemas cuánticos no hermíticos se ha expandido desde las curiosidades matemáticas hacia una frontera de la física, que abarca fenómenos como los puntos excepcionales (EPs) y el efecto de piel no hermítico (NHSE). Sin embargo, simular sistemas de muchos cuerpos no hermíticos en hardware de Escala Intermedia con Ruido (NISQ) presenta tres desafíos fundamentales:

1. Autovalores Complejos: Los VQE estándar dependen de funciones de costo de valores reales derivadas de Hamiltonianos hermíticos; los autovalores complejos imposibilitan un límite inferior variacional natural.
2. Biorortogonalidad: Los sistemas no hermíticos requieren la preparación y medición independiente de los autoestados izquierdo ($\langle \tilde{\psi}_L|$) y derecho ($|\psi_R\rangle$), ya que el producto interno estándar es reemplazado por una métrica biorortogonal.
3. Sobrecarga de Recursos: Simular la dinámica no hermítica en computadoras cuánticas unitarias requiere típicamente dilatación basada en ancillas o postselección, lo que conduce a una sobrecarga exponencial que hace que las simulaciones sean inviables más allá de ~8 qubits.

Las propuestas existentes han abordado aspectos parciales (por ejemplo, evolución en tiempo imaginario o dilatación de la matriz de densidad), pero carecen de un marco unificado capaz de acceder simultáneamente a los estados fundamentales biorortogonales, detectar EPs y mapear diagramas de fase de muchos cuerpos.

Metodología: El Marco B-VQE
Los autores introducen el Algoritmo Variacional de Autovalores Biorortogonales (B-VQE), un algoritmo de circuito dual diseñado específicamente para Hamiltonianos no hermíticos en hardware NISQ.

• Ansatz de Circuito Dual: A diferencia del VQE estándar, el B-VQE emplea dos unitarias parametrizadas independientes, $U_R(\theta)$ y $U_L(\phi)$, para aproximar el autoestado derecho $|\psi_R^\theta\rangle$ y el autoestado izquierdo $|\tilde{\psi}_L^\phi\rangle$, respectivamente.
• Función de Costo Biorortogonal: El algoritmo optimiza una función de costo $L_{B-VQE}$ definida como:
  $$L_{B-VQE}(\theta, \phi) = \text{Re}[E_{bio}(\theta, \phi)] + \lambda [\text{Im}(E_{bio}(\theta, \phi))]^2$$
  donde $E_{bio}$ es el cociente de Rayleigh biorortogonal. El segundo término penaliza las desviaciones de energía imaginaria, impulsando al sistema hacia autovalores reales en la fase de simetría Paridad-Tiempo (PT) no rota.
• Detector de Puntos Excepcionales (EPD): Un módulo operacional identifica los EPs mediante la medición de la métrica de coalescencia $C_\lambda = 1 - |\langle \tilde{\psi}_L^\phi | \psi_R^\theta \rangle|^2$. En un EP, los autoestados izquierdo y derecho coalescen ($C_\lambda \to 0$), proporcionando una firma nativa del hardware.
• Tensor Geométrico Cuántico No Hermítico (NH-QGT): El marco extiende el tensor geométrico cuántico al entorno biorortogonal. Esto permite la caracterización separada de la topología del estado (vía la curvatura de Berry biorortogonal) y la topología de banda, resolviendo discrepancias específicas de los sistemas no hermíticos.
• Mitigación por Muestreo de Importancia (IS): Para superar la sobrecarga exponencial de la postselección, los autores proponen un esquema de reponderación clásica. Al muestrear cadenas de bits de los circuitos cuánticos y aplicar pesos clásicos basados en la relación de las amplitudes derecha e izquierda, el método recupera una sobrecarga polinomial sin requerir postselección de ancillas.

Resultos Clave
El marco fue validado en tres modelos paradigmáticos utilizando simulaciones de Qiskit-Aer calibradas con modelos de ruido de hardware IBM Heron r2:

1. Cadena de Hubbard No Hermítica: El B-VQE mapeó con éxito un diagrama de fases que contiene cuatro fases distintas: ergódica, localización de muchos cuerpos no hermítica (NH-MBL), ergódica con PT rota y localización de piel. El algoritmo logró errores de energía relativa $\epsilon_E < 5 \times 10^{-3}$ e identificó correctamente la transición MBL mediante razones de espaciamiento de niveles y escalamiento de la entropía de entrelazamiento.
2. Cadena de Espín XXZ No Hermítica: El método detectó puntos excepcionales con una precisión de $\delta\lambda < 0.02t$. También caracterizó las "cicatrices cuánticas de muchos cuerpos mejoradas por puntos excepcionales", observando que los tiempos de coherencia de las cicatrices se prolongan exponencialmente cerca de los EPs. La entropía de entrelazamiento biorortogonal reveló un escalamiento de teoría de campo conforme (CFT) con una carga central efectiva $c_{eff} \approx 0.85$ en el punto crítico del EP.
3. Modelo t-J 2D No Hermítico: El algoritmo resolvió una piel de Fermi de muchos cuerpos, donde los estados se acumulan en el límite. Identificó con éxito el número de giro (winding number) topológico, demostrando una transición aguda de fases triviales ($W=0$) a topológicas ($W=1$) coincidiendo con el inicio de la acumulación de piel macroscópica.

Desempeño y Escalabilidad

• Precisión: En todos los modelos, el B-VQE mantuvo errores de energía por debajo de $5 \times 10^{-3}$ en simuladores sin ruido y se mantuvo robusto bajo ruido de hardware calibrado (ej. IBM Kingston, Fez, Marrakesh).
• Detección de EP: El módulo EPD localizó puntos excepcionales dentro de $\delta\lambda < 0.02t$ en simuladores sin ruido y $\delta\lambda < 0.08t$ bajo condiciones de ruido realista.
• Eficiencia de Recursos: Comparado con la Trotterización, el B-VQE requirió aproximadamente 2 veces menos circuitos profundos para alcanzar una fidelidad de estado equivalente. El esquema de mitigación IS redujo la escala de la sobrecarga de exponencial (postselección estándar) a polinomial ($O(n^{1.55})$), permitiendo simulaciones de hasta $n=16$ qubits con probabilidades de éxito $>0.4$.

Significación y Reivindicaciones
El artículo posiciona al B-VQE como una metodología de NISQ completa y escalable para la física de muchos cuerpos no hermíticos. Sus contribuciones primarias son:

1. Resolución de la Brecha Algorítmica: Proporciona el primer marco variacional que accede simultáneamente a los estados fundamentales biorortogonales, detecta sistemáticamente los puntos excepcionales y mapea diagramas de fase de muchos cuerpos.
2. Resolución Topológica: Al utilizar el NH-QGT, resuelve la discrepancia entre la topología del estado y la topología de banda inherente a los sistemas no hermíticos, una característica inaccesible para el VQE hermítico estándar o los enfoques de solo estructura de bandas.
3. Viabilidad Práctica: El esquema de mitigación de muestreo de importancia elimina la barrera de la postselección exponencial, haciendo que la simulación de sistemas interactuantes no hermíticos sea factible en el hardware actual y futuro cercano.

Los autores concluyen que el B-VQE establece una vía para la exploración nativa de hardware de la detección cuántica mejorada por EP, el transporte basado en el efecto de piel y la criticidad cuántica no unitaria, tendiendo un puente efectivo entre la física teórica de muchos cuerpos no hermíticos y la simulación cuántica experimental.