FEM과 이동 비대칭법을 이용한 단일 스핀 최적 제어 가속화

초록

본 논문은 시간 축을 공간 좌표처럼 다루어 리우빌-폰노이만 방정식을 유한 요소법(FEM)으로 선형 시스템화하고, 이를 이동 비대칭법(MMA)과 결합해 단일 스핀의 RF 펄스 최적화를 가속한다. FEM은 스핀 궤적과 제어 그래디언트를 동시에 계산하도록 설계되었으며, adjoint 방식을 이용해 그래디언트 비용을 상태 전파와 동등하게 만든다. Helmholtz 필터와 하이퍼볼릭 탄젠트 스케일링으로 펄스 형태를 부드럽게 제약하고, MMA는 제약식(전력 제한)을 만족하면서 목표 충실도 0.995에 빠르게 수렴한다. 기존 L‑BFGS·Newton‑Raphson 대비 전체 실행 시간이 크게 단축됨을 실험적으로 입증한다.

상세 분석

이 연구는 양자 최적 제어 분야에서 두 가지 핵심 병목을 동시에 해결한다. 첫 번째는 그래디언트 계산 비용이다. 전통적인 GRAPE은 시간 구간마다 전방·후방 전파를 수행하고, 각 구간에 대해 유한 차분으로 미분을 구한다. 저자들은 시간을 1차원 메쉬로 이산화하고, 라그랑주 형태의 선형 형태함수를 이용해 리우빌‑폰노이만 방정식을 Galerkin 방법으로 투영한다. 결과적으로 전역 강성 행렬 K와 하중 벡터 f를 구성해 K·α = f 형태의 선형 시스템을 얻으며, α는 각 노드에서의 밀도 행렬 값을 담는다. 이 시스템을 한 번 풀면 스핀 궤적 전체와 동시에 K⁻¹·∂f/∂x – (∂K/∂x)·α 같은 그래디언트 항을 구할 수 있다. adjoint 방정식 Kᵀλ = (∂η/∂α)ᵀ을 추가로 푸는 것만으로 그래디언트 비용을 O(N)에서 O(1) 수준으로 낮춘다. 두 번째는 최적화 수렴 속도이다. MMA는 목적함수와 제약함수의 현재 값·그래디언트를 이용해 각 반복마다 볼록한 하위문제를 생성하고, 이 하위문제는 해석적으로 최적해를 갖는다. 따라서 전력 제한 같은 비선형 구속을 자연스럽게 포함하면서도 L‑BFGS가 겪는 비선형성에 의한 진동을 회피한다. 또한, Helmholtz 필터(K_h·x_s = f_h)를 별도 선형 시스템으로 풀어 제어 파형을 부드럽게 만들고, 하이퍼볼릭 탄젠트 스케일링(x = (1−e^{−κx_s})/(1+e^{−κx_s}))으로 진폭을