A Quantum Spectral Solver for Periodic Incompressible Stokes Flow

Resumen Técnico: Un Solucionador Espectral Cuántico para el Flujo de Stokes Incompresible Periódico

Planteamiento del Problema
El artículo aborda el desafío de resolver las ecuaciones de Stokes incompresibles, estacionarias y bidimensionales en un dominio periódico utilizando computación cuántica. Dado que los métodos espectrales son altamente efectivos para problemas periódicos —transformando operadores diferenciales en multiplicadores algebraicos diagonales en el espacio de Fourier—, representar la restricción de incompresibilidad (la condición de divergencia nula) y el acoplamiento presión-velocidad asociado dentro de un marco unitario cuántico sigue siendo un obstáculo significativo. Los autores buscan desarrollar un solucionador espectral cuántico que gestione explícitamente la naturaleza vectorial del campo de velocidades, imponga la restricción de incompresibilidad y recupere tanto los campos de velocidad como de presión, además de observables homogeneizados, sin necesidad de reconstruir el campo completo.

Metodología
El método propuesto aprovecha la Transformada de Fourier Cuántica (QFT) como un cambio de base coherente para moverse entre el espacio físico y el espacio espectral. El núcleo del algoritmo se basa en una factorización espectral del operador de Stokes que separa el problema en operaciones cuánticas manejables:

1. Factorización Espectral: En el espacio de Fourier, las ecuaciones de Stokes se desacoplan modo a modo. La solución implica una proyección de Helmholtz para imponer la incompresibilidad (eliminando la componente longitudinal de la fuerza) seguida de un escalado de Laplaciano inverso escalar.
2. Proyección Basada en Rotación: Para el caso 2D, la proyección de Helmholtz se realiza mediante una rotación dependiente del modo que transforma las componentes cartesianas de la velocidad $(u_0, u_1)$ en coordenadas longitudinales ($\parallel$) y transversales ($\perp$) relativas al vector de onda $k$. La restricción de incompresibilidad se impone entonces anulando la componente longitudinal y escalando la componente transversal por el factor del Laplaciano inverso ($1/(\mu|k|^2)$).
3. Arquitectura del Circuito Cuántico:
   • Registros: El circuito utiliza registros de modo ($k_0, k_1$) para los modos de Fourier, un registro de componentes ($c$) para las direcciones cartesianas, un registro objetivo ($t$) para la codificación de amplitud y un registro de trabajo ($w$) para la aritmética.
   • Preparación del Estado: El campo de forzamiento de entrada se codifica en un estado cuántico normalizado.
   • QFT: Una QFT 2D transforma el estado al espacio de Fourier.
   • Rotación Angular ($U_R$): Se aplica una rotación dependiente del modo $R(k)$ al registro de componentes para alinear con la base longitudinal-transversal. El ángulo de rotación se aproxima mediante polinomios por partes codificados mediante rotaciones controladas.
   • Escalado Escalar ($U_\Gamma, U_\Lambda$): La componente transversal se escala por la función de Green $\Gamma(k) = 1/(\mu|k|^2)$, y la componente longitudinal se escala por $\Lambda(k) = 1/|k|$ (para la recuperación de la presión). Estos escalados se implementan mediante la codificación de amplitud polinómica en el registro objetivo.
   • Post-selección y Reconstrucción: Se recupera la solución deseada post-seleccionando ramas específicas de los registros auxiliares (por ejemplo, $c=1$ para la velocidad transversal, $c=0$ para la presión). La QFT inversa reconstruye el campo físico.
4. Extracción de Observables: Los autores demuestran que las cantidades homogeneizadas, como la energía cinética, pueden estimarse directamente desde el estado espectral antes de la QFT inversa, evitando el costo de la reconstrucción completa del campo.

Contribuciones Clave

1. Primer Solucionador de Stokes Cuántico Basado en QFT: Los autores presentan lo que afirman es el primer solucionador espectral cuántico para las ecuaciones de Stokes incompresibles, estacionarias y periódicas en 2D.
2. Factorización Espectral Explícita: Derivan una factorización compatible con la computación cuántica que separa la proyección de Helmholtz, el escalado del Laplaciano inverso y la recuperación de la presión, utilizando una representación basada en rotaciones para el proyector 2D.
3. Diseño de Circuito Modular: El artículo proporciona una arquitectura de circuito explícita que incluye rotaciones dependientes del modo y bloques espectrales polinómicos condicionados por la componente. Los parámetros de aproximación (grado polinómico y teselado) se hacen explícitos.
4. Recuperación de Observables Homogeneizados: El método se valida mediante un benchmark inspirado en un RVE, mostrando que la energía cinética homogeneizada puede recuperarse sin reconstruir el campo de velocidad completo, alineándose con los flujos de trabajo de la mecánica computacional multiescala.

Resultados Numéricos
El método fue validado utilizando tres benchmarks en un simulador de estado de vector sin ruido:

• Vórtice de Taylor-Green: Una solución suave manufacturada utilizada para verificar la consistencia de los circuitos ensamblados (convenciones de Fourier, normalización y proyección).
• Dipolo de Fuerza Regularizado: Un benchmark periódico con forzamiento localizado y soporte de Fourier amplio. Esta prueba distinguió entre errores de truncamiento por resolución finita y errores introducidos por la aproximación polinómica de los operadores espectrales. Los resultados mostraron que, mientras que el solucionador discreto exacto converge con el refinamiento de la malla, la aproximación polinómica introduce un plateau de error para los modos de alta frecuencia si el teselado es fijo.
• Fluctuación Inspirada en RVE: Un benchmark que utiliza un espectro de fluctuación tipo Kolmogorov ($k^{-5/3}$). Los resultados demostraron que la energía cinética homogeneizada converge hacia la solución de referencia a medida que se refina la malla, siempre que la aproximación polinómica se ajuste al espectro ponderado por la energía. El estudio destacó que el conteo de puertas está dominado por la codificación polinómica de la rotación y los operadores de la función de Green, no por la propia QFT.

Significancia y Reivindicaciones
Los autores posicionan este trabajo no como un reemplazo directo de la DNS clásica o de un modelo de turbulencia, sino como un bloque espectral de Stokes para la mecánica computacional cuántica. Su principal significancia radica en:

• Compatibilidad Multiescala: El solucionador está diseñado para ser embebido en arquitecturas cuántico-clásicas de Elemento de Volumen Representativo (RVE) más amplias, donde el paralelismo cuántico puede actualizar simultáneamente todos los estados de los RVE asociados con los puntos de cuadratura macroscópicos.
• Eficiencia en la Estimación de Observables: Demuestra la capacidad de extraer cantidades promediadas (como la energía cinética) directamente del estado espectral, lo cual es crucial para los flujos multiescala donde la reconstrucción completa del campo es innecesaria y costosa.
• Escalabilidad: Bajo supuestos estándar de preparación de estado eficiente y estimación de observables, el circuito exhibe una dependencia polilogarítmica de la resolución de la malla ($O(\text{polylog}(N))$), con el grado polinómico y el conteo de teselas sirviendo como parámetros de aproximación explícitos.

Los autores reconocen que los resultados actuales están limitados a simulaciones sin ruido y que una evaluación completa de recursos debe considerar los costos de preparación de estado, las probabilidades de éxito de la post-selección y los gastos generales de tolerancia a fallos. Señalan que extender este enfoque a 3D o dominios no periódicos probablemente requeriría formulaciones de Combinación de Unitarias (LCU) para el proyector vectorial en lugar del enfoque basado en rotaciones utilizado aquí.