Quantum-Accelerated Nonlinear Seismic Simulation for Megathrust-Prone Urban Environments: A Three-Order-of-Magnitude Speedup via HHL-Embedded Newton–Raphson Integration

Resumen Técnico: Simulación Sísmica No Lineal Acelerada por Computación Cuántica

Planteamiento del Problema
El Análisis de Historial Temporal No Lineal (NLTHA, por sus siglas en inglés) es el estándar riguroso para evaluar el desempeño estructural sísmico, exigido por normativas como la ASCE 7-22 y la NTE E.030 de Perú para estructuras irregulares y esenciales. Sin embargo, su aplicación está severamente limitada por el costo computacional. Para un sistema con $N$ grados de libertad (DOF), la solución clásica requiere resolver un sistema lineal en cada iteración de Newton–Raphson, lo que resulta en una complejidad temporal de $O(N^3)$ por paso de tiempo. Esta escala restringe el NLTHA al estudio de edificios individuales, impidiendo evaluaciones de riesgo a escala urbana y en tiempo real en regiones de alto peligro como los entornos urbanos de la zona de subducción andina, propensos a megatrusts.

Metodología
El estudio propone un marco de integración de Newton–Raphson Mejorado por Computación Cuántica (QNR) diseñado para reducir el tiempo de cómputo de tiempo de pared (wall-clock) en tres órdenes de magnitud. La metodología integra el resolvedor lineal cuántico de Harrow–Hassidim–Lloyd (HHL) en el esquema de pasos de tiempo clásico de Hilbert–Hughes–Taylor (HHT-$\alpha$).

1. Formulación Clásica: La ecuación de movimiento de un sistema no lineal de $N$ DOF se discretiza utilizando HHT-$\alpha$. En cada paso de tiempo, el residuo no lineal se resuelve mediante iteraciones de Newton–Raphson, lo que requiere la solución de $K_T^{(i)} \Delta u^{(i)} = R^{(i)}$, donde $K_T$ es la matriz de rigidez tangente efectiva.
2. Mapeo Cuántico:
   • Codificación de Estado: El vector de fuerza residual $R$ se codifica en amplitud en un estado cuántico $|b\rangle$ utilizando $n = \lceil \log_2 N \rceil$ qubits.
   • Descomposición de Hamiltoniano: La matriz de rigidez tangente normalizada $A = K_T / \|K_T\|$ se descompone en una suma de productos tensoriales de operadores de Pauli para formar un Hamiltoniano cuántico.
   • Ejecución de HHL: El algoritmo HHL realiza la estimación de fase cuántica y la rotación controlada para obtener un estado $|x\rangle \propto A^{-1}|b\rangle$, resolviendo efectivamente el sistema lineal.
   • Lectura: Los incrementos de desplazamiento se extraen mediante tomografía de sombra clásica, requiriendo solo $O(\log N)$ mediciones para observables dispersos.
3. Protocolo de Validación: El marco fue validado mediante un estudio de Monte Carlo de 1,000 ejecuciones en un modelo de edificio de cortante no lineal de 64 DOF. Los movimientos del suelo se escalaron desde la base de datos PEER NGA-West2 para coincidir con el espectro de diseño peruano para Lima (tipo de suelo S3). La significancia estadística se evaluó mediante pruebas t de dos muestras y pruebas de normalidad de Shapiro–Wilk.

Contribuciones Clave

• Marco QNR: Una formulación matemáticamente rigurosa que incrusta HHL dentro del bucle de Newton–Raphson de HHT-$\alpha$ para la dinámica estructural.
• Formalización Física: El uso de Diagramas de Cuerpo Libre (FBD) para formalizar el mapeo de los principios físicos de SDOF y MDOF a estados cuánticos.
• Benchmarking Cuantitativo: Una matriz de comparación multimétrica que demuestra la transición de la complejidad temporal clásica de $O(N^3)$ y la complejidad espacial de $O(N^2)$ a una complejidad temporal cuántica de $O(s^2\kappa^2 \log N)$ y una complejidad espacial de $O(\log N)$.
• Robustez Estadística: Una validación exhaustiva que muestra una aceleración media de 1,002.5× con un intervalo de confianza del 95% de [980.2, 1,024.8] y un Error Cuadrático Medio (MSE) de $1.2 \times 10^{-4}$ frente a una referencia clásica de precisión completa.

Resultos
El marco QNR logró un factor de aceleración media ($\mu$) de 1,002.5× en comparación con la línea base clásica de HHT-$\alpha$.

• Significancia Estadística: El resultado es estadísticamente significativo ($t = 2,814.5$, $p < 0.0001$), con el intervalo de confianza del 95% excluyendo la hipótesis nula de ausencia de aceleración. La distribución de la aceleración fue confirmada como normal (Shapiro–Wilk $W = 0.998$, $p = 0.31$).
• Precisión: El MSE de $1.2 \times 10^{-4}$ está muy por debajo del umbral de tolerancia de ingeniería del 1% para el diseño estructural basado en deriva.
• Escalabilidad: Teóricamente, el marco soporta modelos de hasta $N \approx 10^6$ DOF, comparado con el límite práctico clásico de $N \approx 10^3$.
• Eficiencia: El estudio estima una reducción de energía relativa de aproximadamente 100× por solución.

Significancia y Reivindicaciones
El artículo sostiene que este marco transiciona el NLTHA de ser una herramienta computacionalmente prohibitiva para edificios individuales a una plataforma viable para la evaluación de riesgo sísmico a escala urbana y en tiempo real. Específicamente, para Lima, Perú —una metrópolis situada sobre cuencas de suelos blandos dentro del sistema de subducción Nazca–Sudamericana— el marco QNR permite la evaluación no lineal simultánea de miles de estructuras. Esta capacidad respalda directamente el cumplimiento de la NTE E.030 y los mandatos de reducción del riesgo de desastres de CENEPRED.

Los autores identifican tres aplicaciones transformadoras para el contexto peruano:

1. Simulación a Escala Urbana en Tiempo Real: Generar mapas de estado de daño inmediatos para la respuesta de emergencia (ej. CENEPRED) mediante la simulación simultánea de $\sim 10,000$ marcos de concreto reforzado.
2. Análisis de Peligro Sísmico Probabilístico (PSHA): Ejecutar millones de corridas para producir curvas de fragilidad no lineales que incorporen la interacción suelo–estructura específica de la geología de Lima.
3. Integración de Alerta Temprana Sísmica: Acoplar el NLTHA cuántico con el sistema SWAVE para predecir la respuesta no lineal de la infraestructura crítica (hospitales, puentes, represas) antes de la llegada de las ondas S destructivas.

Limitaciones y Trabajo Futuro
Los autores reconocen que los resultados actuales dependen de modelos teóricos de conteo de puertas para procesadores tolerantes a fallos, ya que los dispositivos actuales de Escala Intermedia con Ruido (NISQ) tienen tasas de error de puerta ($\sim 10^{-3}$) que limitan la profundidad del circuito. Además, la sobrecarga de preparación de estado para cargar matrices puede escalar como $O(N^2)$ en los peores escenarios, y el crecimiento del número de condición con el refinamiento de la malla podría erosionar la ventaja cuántica sin precondicionamiento. Se propone como trabajo futuro implementar Resolvedores Lineales Cuánticos Variacionales (VQLS) para hardware de corto plazo, realizar benchmarks en procesadores cuánticos existentes (IBM Quantum, IonQ) e incorporar modelos de suelo no lineales calibrados al perfil geotécnico de Lima.