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

기술 요약: 양자 가속 비선형 지진 시뮬레이션

문제 정의
비선형 이력 해석(NLTHA)은 불규칙 구조물 및 필수 구조물에 대해 ASCE 7-22 및 페루의 NTE E.030과 같은 규정에서 요구하는 지진 성능 평가의 엄격한 표준이다. 그러나 이의 적용은 계산 비용으로 인해 심각하게 제한된다. $N$개의 자유도(DOF)를 가진 시스템의 경우, 고전적 해법은 매 뉴턴-랩슨(Newton–Raphson) 반복 단계마다 선형 시스템을 풀어야 하며, 결과적으로 시간 복잡도는 시간 단계당 $O(N^3)$이 된다. 이러한 스케일링은 NLTHA를 단일 건물 연구로 국한시키며, 안데스 섭입대의 거대 역단층(megathrust) 영향권에 있는 고위험 도시 지역과 같은 도시 규모의 실시간 위험 평가를 불가능하게 만든다.

방법론
본 연구는 벽시계 계산 시간(wall-clock computation time)을 3개 차수(orders of magnitude)만큼 줄이기 위해 설계된 양자 강화 뉴턴-랩슨(QNR) 통합 프레임워크를 제안한다. 이 방법론은 해로우-해시디즘-로이드(HHL) 양자 선형 솔버를 고전적인 힐버-휴즈-테일러(HHT-$\alpha$) 시간 단계 기법에 통합한다.

1. 고전적 정식화: $N$-DOF 비선형 운동 방정식은 HHT-$\alpha$를 사용하여 이산화된다. 각 시간 단계에서 비선형 잔차(nonlinear residual)는 $K_T^{(i)} \Delta u^{(i)} = R^{(i)}$를 해결하는 과정을 통해 해결되며, 여기서 $K_T$는 유효 접선 강성 행렬이다.
2. 양자 매핑:
   • 상태 인코딩: 잔류력 벡터 $R$은 $n = \lceil \log_2 N \rceil$개의 큐비트를 사용하여 양자 상태 $|b\rangle$로 진폭 인코딩된다.
   • 해밀토니안 분해: 정규화된 접선 강성 행렬 $A = K_T / \|K_T\|$는 양자 해밀토니안을 형성하기 위해 파울리 연산자의 텐서 곱 합으로 분해된다.
   • HHL 실행: HHL 알고리즘은 양자 위상 추정 및 제어 회전을 수행하여 $|x\rangle \propto A^{-1}|b\rangle$인 상태를 생성하며, 이는 효과적으로 선형 시스템을 해결한다.
   • 판독(Readout): 변위 증분은 고전적 섀도 토모그래피(classical shadow tomography)를 통해 추출되며, 희소 관측량(sparse observables)에 대해 $O(\log N)$번의 측정만을 요구한다.
3. 검증 프로토콜: 본 프레임워크는 64-DOF 비선형 전단 건물 모델에 대해 1,000회 실행된 몬테카를로 연구를 통해 검증되었다. 지진동은 리마(Lima)의 설계 스펙트럼(토양 유형 S3)에 부합하도록 PEER NGA-West2 데이터베이스에서 스케일링되었다. 통계적 유의성은 두 표본 t-검정(two-sample t-tests) 및 샤피로-윌크 정규성 검정(Shapiro–Wilk normality tests)을 사용하여 평가되었다.

주요 기여

• QNR 프레임워크: 구조 역학을 위해 HHL을 HHT-$\alpha$ 뉴턴-랩슨 루프 내에 내장한 수학적으로 엄격한 정식화.
• 물리적 정식화: SDOF 및 MDOF의 물리적 원리를 양자 상태로 매핑하기 위한 자유물체도(FBD)의 사용.
• 정량적 벤치마킹: 고전적 $O(N^3)$ 시간 복잡도 및 $O(N^2)$ 공간 복잡도에서 양자 $O(s^2\kappa^2 \log N)$ 시간 복잡도 및 $O(\log N)$ 공간 복잡도로의 전환을 입증하는 다중 지표 비교 매트릭스.
• 통계적 견고성: 평균 1,002.5배의 속도 향상(95% 신뢰 구간 [980.2, 1,024.8])과 고정밀 고전 참조 모델 대비 $1.2 \times 10^{-4}$의 평균 제곱 오차(MSE)를 보여주는 종합적인 검증.

결과
QNR 프레임워크는 고전적 HHT-$\alpha$ 베이스라인 대비 1,002.5배의 평균 속도 향상 계수($\mu$)를 달値했다.

• 통계적 유의성: 결과는 통계적으로 유의미하며($t = 2,814.5$, $p < 0.0001$), 95% 신뢰 구간은 속도 향상이 없다는 귀무가설을 배제한다. 속도 향상 분포는 정규성을 확인하였다 (Shapiro–Wilk $W = 0.998$, $p = 0.31$).
• 정확도: $1.2 \times 10^{-4}$의 MSE는 변위 기반 구조 설계의 1% 공학적 허용 오차 범위를 훨씬 밑돈다.
• 확장성: 이론적으로 본 프레임워크는 실질적인 고전적 한계인 $N \approx 10^3$에 비해 최대 $N \approx 10^6$ DOF의 모델을 지원한다.
• 효율성: 본 연구는 솔루션당 약 100배의 상대적 에너지 감소를 추정한다.

의의 및 주장
본 논문은 이 프레임워크가 NLTHA를 단일 건물을 위한 계산적으로 부담스러운 도구에서 실시간 도시 규모 위험 평가를 위한 실행 가능한 플랫폼으로 전환한다고 주장한다. 특히, 나스카-남미 섭입 시스템 내의 연약 지반 분지에 위치한 메트로폴리스인 페루 리마의 경우, QNR 프레임워크는 수천 개의 구조물을 동시에 비선형 평가하는 것을 가능하게 한다. 이러한 능력은 NTE E.030 준수 및 CENEPRED의 재난 위험 감소 명령을 직접적으로 지원한다.

저자들은 페루의 맥락에서 세 가지 변혁적인 응용 분야를 식별하였다:

1. 실시간 도시 규모 시뮬레이션: 약 10,000개의 철근 콘크리트 프레임을 동시에 시뮬레이션함으로써 긴급 대응(예: CENEPRED)을 위한 즉각적인 손상 상태 지도를 생성한다.
2. 확률론적 지진 위험 분석(PSHA): 리마의 지질 특성에 특화된 지반-구조물 상호작용을 포함하는 비선형 취약도 곡선을 생성하기 위해 수백만 번의 실행을 수행한다.
3. 지진 조기 경보 통합: 파괴적인 S파가 도착하기 전 핵심 인프라(병원, 교량, 댐)의 비선형 응답을 예측하기 위해 양자 NLTHA를 SWAVE 시스템과 결합한다.

한계 및 향래 과제
저자들은 현재의 결과가 결함 허용 프로세서를 위한 이론적 게이트 수 모델에 의존하고 있음을 인정한다. 현재의 노이즈가 있는 중간 규모 양자(NISQ) 장치는 게이트 오류율($\sim 10^{-3}$)이 높아 회로 깊이를 제한하기 때문이다. 또한, 행렬 로딩을 위한 상태 준비 오버헤드는 최악의 경우 $O(N^2)$로 스케일링될 수 있으며, 격자 세분화에 따른 조건수(condition number) 성장은 프리컨디셔닝(preconditioning) 없이 양자 이점을 저해할 수 있다. 향후 과제로 근접 미래 하드웨어를 위한 변분 양자 선형 솔버(VQLS) 구현, 기존 양자 프로세서(IBM Quantum, IonQ)에서의 벤치마킹, 그리고 리마의 지질학적 프로파일에 맞춰 보정된 비선형 토양 모델의 통합이 제안되었다.