Relativistic Quantum Simulation under Periodic and Dirichlet Boundary Conditions: A First-Quantised Framework for Near-Term Devices

Resumen Técnico: Simulación Cuántica Relativista bajo Condiciones de Contorno Periódicas y de Dirichlet

Planteamiento del Problema
La simulación de sistemas cuánticos es computacionalmente exigente debido al crecimiento exponencial de la complejidad con el tamaño del sistema, un desafío que se exacerba cuando deben incorporarse efectos relativistas para describir fenómenos físicos con precisión. Mientras que la teoría de perturbaciones es un enfoque estándar para estudiar efectos relativistas más allá del marco no relativista, la implementación de estas simulaciones en dispositivos cuánticos de corto plazo sigue siendo difícil. Los métodos existentes suelen depender de la segunda cuantización, que es el estándar para la química cuántica, pero puede carecer de un manejo intuitivo para sistemas con no linealidades. Existe la necesidad de un marco práctico de primera cuantización capaz de estimar las energías del estado fundamental relativistas para sistemas bajo diferentes condiciones de contorno en hardware de corto plazo.

Metodología
Los autores proponen un marco de primera cuantización para la simulación cuántica relativista que discretiza la función de onda en una rejilla finita representada por qubits del sistema. La metodología central involucra los siguientes componentes:

1. Discretización y Operadores de Traslación:
   La función de onda continua 1D se mapea a una rejilla discreta de $2^L$ puntos utilizando $L$ qubits, donde la posición $x$ se escala al dominio $0 \le x < 1$. El operador de momento al cuadrado $\hat{p}^2$ se aproxima mediante el método de diferencias finitas basado en operadores de traslación cuántica. El operador de traslación $\hat{A}_x$ (sumador cuántico) desplaza el estado de posición $|j\rangle$ a $|j+\delta x\rangle$, mientras que su conjugado $\hat{A}_x^\dagger$ (sustractor cuántico) desplaza en la dirección opuesta. El momento al cuadrado se expresa como $\hat{p}^2 = -(mc L_m)^2 (\hat{A}^{(1)} - 2\hat{1})$, donde $L_m = \lambda_m / \delta x$ es un parámetro adimensional que involucra la longitud de onda de Compton reducida.

2. Expansión Perturbativa de la Energía Cinética:
   El Hamiltoniano de energía cinética relativista $\hat{H}_{re}^K$ se deriva de la forma de raíz cuadrada de la relación energía-momento relativista. Se expande en una serie perturbativa de términos de momento de orden superior ($\hat{p}^2, \hat{p}^4, \hat{p}^6, \dots$) válida en el régimen donde $\langle \hat{p}^2 \rangle \ll m^2 c^2$. La energía cinética relativista total se aproxima truncando esta serie, incorporando coeficientes derivados de la ecuación de Klein-Gordon.

3. Manejo de Condiciones de Contorno:
   El marco aborda dos condiciones de contorno distintas:

   • Condiciones de Contorno Periódicas (PBC): Los operadores de traslación preservan la unitariedad envolviendo la rejilla ($\hat{A}_x |2^L-1\rangle = |0\rangle$). Los valores de expectativa de las potencias pares del momento se expresan como combinaciones de $\langle \hat{A}^{(l)} \rangle$.
   • Condiciones de Contorno de Dirichlet (DBC): Aquí, la función de onda no se envuelve ($\psi(0) \neq \psi(1)$). Los autores introducen un operador de corrección $\hat{E}_0 = |2^L-1\rangle\langle 0| + |0\rangle\langle 2^L-1|$ para contabilizar el contorno. Los valores de expectativa para DBC se formulan como los términos de PBC más términos adicionales que involucran a $\hat{E}_0$ y operadores relacionados ($\hat{E}_1, \hat{E}_2$), los cuales dependen de las amplitudes de probabilidad en los bordes de la rejilla.

4. Algoritmo de Simulación Cuántica:
   Los autores diseñan circuitos cuánticos para estimar los valores de expectativa requeridos:

   • Energía Cinética: Se aplican puertas de traslación controladas de forma iterativa a un qubit de control. La medición del qubit de control en la base Z produce la parte real de $\langle (\hat{A}_x^\dagger)^l \rangle$, lo que corresponde a los términos de momento necesarios.
   • Energía Potencial: Se utiliza una matriz de densidad $\hat{\rho}_V$ que representa el potencial en un circuito para estimar $\langle \hat{V} \rangle$ a través de las probabilidades del estado del sistema.
   • Optimización Variacional: La energía total se estima minimizando la suma de los valores de expectativa de la energía cinética y potencial sobre un conjunto de parámetros variacionales $\{\bar{\kappa}\}$ en un estado ansatz $|\psi(\{\bar{\kappa}\}, x)\rangle$.

Contribuciones Clave

• Receta de Primera Cuantización: El artículo proporciona una receta específica para la simulación relativista utilizando la primera cuantización, ofreciendo una alternativa a los métodos de segunda cuantización ampliamente utilizados, particularmente beneficiosos para sistemas con no linealidades.
• Marco de Contorno Unificado: Presenta un enfoque unificado para manejar tanto las condiciones de contorno Periódicas como las de Dirichlet. Crucialmente, demuestra que los datos obtenidos para PBC pueden reutilizarse para estimar las energías de DBC añadiendo términos de corrección de contorno específicos, reduciendo la carga computacional.
• Compatibilidad con Dispositivos de Corto Plazo: El método utiliza primitivas de circuito cuántico estándar (sumadores controlados, puertas Hadamard y mediciones) y optimización variacional, lo que lo hace adecuado para su implementación en dispositivos cuánticos de corto plazo.
• Diseños de Circuitos Explícitos: Los autores proporcionan diseños de circuitos detallados para computar los valores de expectativa de los operadores de traslación y los términos de corrección de contorno requeridos para la expansión perturbativa.

Resultados y Reclamaciones
El artículo afirma ofrecer una "vía práctica" para simular efectos relativistas en dispositivos cuánticos de corto plazo. Los resultados teóricos muestran que:

• Los términos de energía cinética relativista hasta un orden específico de perturbación pueden descomponerse en combinaciones lineales de valores de expectativa de operadores de traslación ($\langle \hat{A}^{(l)} \rangle$) y operadores de contorno ($\langle \hat{E}_d \rangle$).
• Los coeficientes para estas expansiones pueden precalcularse basándose en parámetros físicos ($m, c, \hbar, \delta x$).
• El método permite la optimización variacional de estados ansatz para estimar energías del estado fundamental tanto no relativistas como relativistas.

Significancia
Los autores posicionan este trabajo como un paso fundacional hacia el apoyo de desarrollos futuros en física y química cuántica. Al permitir la simulación de efectos relativistas en hardware de corto plazo, el marco aborda una brecha en las capacidades actuales de simulación cuántica. El artículo enfatiza que, si bien aumentar el número de qubits mejora la resolución de la rejilla, se debe tener cuidado para asegurar que el parámetro adimensional $L_m$ permanezca lo suficientemente pequeño para que la expansión perturbativa siga siendo válida. El trabajo no pretende resolver la ecuación de Dirac completa ni manejar límites ultra-relativistas, sino que se enfoca específicamente en el régimen donde las correcciones relativistas al marco no relativista son relevantes. Se sugiere como trabajo futuro explorar operadores de Laplaciano de orden superior y métodos de construcción de ansatz más precisos.