Quantum computing-based solver for interacting power grids

기술 요약: 상호작용하는 전력망을 위한 양자 컴퓨팅 기반 솔버

문제 정의
다단 터미널 송전망에서 전력 전자 장치와 재생 에너지원의 통합이 증가함에 따라 복잡한 동적 불안정성과 고조파 왜곡이 도입되었습니다. 공진 모드 분석(Resonance Mode Analysis, RMA)은 시스템 어드미턴스 행렬을 분해하여 임계 공진 모드를 식별하는 표준 기술이지만, 대규모 네트워크에 적용할 때 심각한 계산 병목 현상에 직면합니다. 이러한 시스템의 어드미턴스 행렬은 비에르미트(non-Hermitian) 성질을 가지며(활성 및 무효 전력의 존재로 인해), 고차원적입니다. 비에르미트 란초스(Lanczos) 알고리즘과 같은 고전적 반복 기법은 거대한 크릴로프 부공간(Krylov subspace) 벡터를 저장해야 하며, $n \times n$ 행렬에 대해 $O(n^3)$으로 스케일링되는 대각화 과정을 수행해야 합니다. 그리드가 수만 개의 버스로 성장함에 따라, 이러한 메모리 및 프로세싱 요구 사항은 고성능 컴퓨팅(HPC)의 한계를 초과하게 됩니다. 또한, 변분 양자 고유값 솔버(Variational Quantum Eigensolver, VQE)와 같은 표준 양자 고유값 솔버는 에르미트 연산자를 위해 설계되었으므로, 전력망 안정성 분석에 내재된 비에르미트 고유값 문제를 직접 해결할 수 없습니다.

방법론
저자들은 클래식의 메모리 제약과 표준 양자 알고리즘의 한계를 극복하기 위해 실수 분산 기반 변분 양자 고유값 솔버(Real Variance-based Variational Quantum Eigensolver, RVVQE)를 활용하는 새로운 양자-고전 하이브리드 프레임워크를 제안합니다. 방법론은 다음과 같이 진행됩니다:

1. 시스템 정식화: 그리드는 주파수 스캔 분석을 사용하여 모델링되며, 노드 전압과 주입 전류 사이의 관계는 어드미턴스 행렬 $Y_{B,f}$에 의해 정의됩니다. 공진은 $Y_{B,f}$의 고유값이 0에 접근하여 모달 임피던스가 피크를 이룰 때 식별됩니다.
2. 행렬 분해: 비에르미트 어드미턴스 행렬 $Y_{B,f}$를 양자 하드웨어와 호환되도록 두 개의 순수 에르미트 행렬 $H_f$와 $K_f$로 수학적으로 분해합니다 ($Y_{B,f} = H_f + iK_f$).
3. 양자 매핑: 이 에르미트 성분들은 $N$-큐비트 양자 처리 장치(QPU)에서 실행하기 위해 파울리 문자열(Pauli strings)로 매핑됩니다. 여기서 $N = \lceil \log_2 n \rceil$입니다. 이러한 로그 스케일링을 통해 수천 개의 버스를 가진 시스템을 상대적으로 적은 수의 큐비트(예: 70,000개 버스에 대해 17개 큐비트)로 표현할 수 있습니다.
4. RVVQE 알고리즘: 매개변수화된 양자 회로(ansatz)가 시도 상태(trial state) $|\psi(\theta)\rangle$를 준비합니다. 에너지를 최소화하는 대신, 알고리즘은 다음과 같이 정의된 실수 값 분산 비용 함수를 최소화합니다:
   $$C_{var}(\theta) = \langle H_f^2 \rangle_\theta - \langle H_f \rangle_\theta^2 + \langle K_f^2 \rangle_\theta - \langle K_f \rangle_\theta^2$$
   고전적 최적화 도구는 매개변수 $\theta$를 반복적으로 업데이트하여 이 비용 함수를 0으로 유도합니다.
5. 고유값 추출: 수렴되면, 원래의 비에르미트 시스템의 복소 고유값은 에르미트 성분의 기댓값을 합산하여 재구성됩니다: $\lambda = \langle H_f \rangle + i\langle K_f \rangle$. 이후 임계 공진 모달 임피던스는 $Z_{cm,f} = 1/\lambda$로 계산됩니다.
6. 축퇴 제거(Deflation): 여러 고유값을 추출하기 위해, 직교성을 보장하기 위해 이전에 발견된 고유상태와의 중첩에 대한 페널티를 비용 함수에 추가하는 순차적 축퇴 제거 기법이 사용됩니다.

주요 기여

• 알고리즘 적응: 본 논문은 비에르미트 연산자를 위해 원래 설계된 RVVQE 알고리즘을 전력 시스템 공진 분석이라는 특정 맥락에 맞게 조정함으로써, 소산성(dissipative) 전기 그리드와 개방형 양자 시스템 사이의 간극을 메웁니다.
• 확장 가능한 아키텍처: 그리드 상태를 양자 메모리에 로그 단위로 인코딩함으로써, 제안된 방법은 고전적 RAM의 $O(n^2)$ 메모리 제한을 우회하며, 대륙 규모의 네트워크를 분석하기 위한 이론적으로 확장 가능한 경로를 제공합니다.
• 하이브리드 워크플로우: 저자들은 주파수 스윕, 행렬 분해, 그리고 하이브리드 양자-고전 루프를 이용한 안정성 평가를 처리하는 완전한 "그리드-투-양자(grid-to-quantum)" 워크플로우(알고리즘 1)를 정형화했습니다.

결과
프레임워크는 표준 5-버스 송전 시스템을 사용하여 검증되었습니다. 비에르미트 어드미턴스 행렬은 100 Hz에서 3000 Hz 사이의 주파수 범위에 대해 생성되었습니다. 5x5 행렬은 3-큐비트 양자 시스템에 맞게 8x8로 패딩되었습니다.

• 정확도: RVVQE 알고리즘을 통해 계산된 임계 공진 모달 임피던스($|Z_{cm}|$)는 전체 주파수 스펙트럼에 걸쳐 표준 행렬 대각화(SciPy 라이브러리)로부터 도출된 정확한 고전적 해와 거의 완벽하게 일치했습니다.
• 하드웨어 제약: 시뮬레이션은 양자 오류 완화(QEM) 없이 수행되었습니다. 오류 수정이 없음에도 불구하고 결과는 고전적 참조 값에 가깝게 유지되었으며, 이는 알고리즘이 현재의 노이즈가 있는 하드웨어에서도 견고함을 시사합니다.
• 안정성 분석: 이 방법은 공진 주파수와 임계 모드를 성공적으로 식별하였으며, 소신호 안정성(small-signal stability) 및 참여 계수(participation factors)를 평가하는 능력을 입증했습니다.

의의 및 주장
본 논문은 이 연구가 현재 고전적 계산 경계를 넘어서는 거대 네트워크의 공진 불안정성을 진단할 수 있는 "심도 있게 확장 가능한 아키텍처"를 제공한다고 주장합니다. 저자들은 현재 하드웨어가 제한적이지만, 큐비트의 로그 스케일링이 향후 대규모 전력 시스템 분석을 위한 실행 가능한 경로를 제공한다는 점을 강조합니다. 비에르미트 행렬을 위해 RVVQE를 구현한 것은 복잡한 전력 시스템 안정성 연구(향후 더 큰 IEEE 테스트 네트워크(예: 118-bus, 300-bus) 및 ACTIVSg70k와 같은 합성 그리드 포함)에 양자 컴퓨팅을 적용하기 위한 기초적인 단계로 제시됩니다. 이 연구는 산업 규모의 문제를 실제 하드웨어로 해결했다고 주장하는 것이 아니라, 벤치마크 시스템에 대한 접근 방식의 수학적 타당성과 정확성을 입증하는 데 목적이 있습니다.