먼저 읽으면 좋은 용어
개념 소개
VQE는 양자 위상 추정(QPE)처럼 긴 결맞음 시간을 요구하는 대신, 얕은 파라미터화 회로와 고전 최적화기를 결합하여 해밀토니안 의 최솟값 고유값(바닥 상태 에너지 )을 구한다. 아이디어의 출발점은 변분 원리다.
임의의 정규화된 시험 상태 에 대해 항상
이 성립한다. 따라서 파라미터 로 제어되는 양자 회로가 만드는 상태 의 에너지 기댓값을 최소화하면 에 수렴할 수 있다.
핵심 원리
1. 변분 부등식의 엄밀한 유도
의 고유 에너지를 , 고유 상태를 이라 하면,
등호는 일 때 성립한다.
2. 안사츠(Ansatz) 설계
파라미터화된 유니터리 를 통해 시험 상태를 생성한다.
안사츠 유형 선택은 표현력·깊이·노이즈 민감도 사이의 트레이드오프를 결정한다.
| 유형 | 특징 |
|---|---|
| UCCSD | 결합 클러스터 이론 기반 물리 동기 안사츠, 회로 깊이가 깊음 |
| 하드웨어 효율적(HEA) | 장치 연결성에 맞춘 얕은 회로, 노이즈에 강함 |
| ADAPT-VQE | 연산자 풀에서 기울기 크기 순으로 게이트를 동적 추가 |
3. 해밀토니안의 파울리 분해
분자 해밀토니안을 파울리 연산자의 선형 결합으로 변환한다.
에너지 기댓값은 각 항을 개별 측정하여 합산한다.
4. 파라미터 이동 규칙
기울기 기반 최적화를 사용할 때, 회로 내 회전 게이트 (는 파울리 연산자)의 편미분은 추가 회로 실행만으로 해석적으로 계산된다.
고전 수치 미분 없이 하드웨어에서 그대로 실행 가능하다는 점이 핵심이다.
예시·응용
수소 분자 기저 상태 에너지 (Qiskit 예시)
from qiskit.circuit.library import TwoLocal
from qiskit.quantum_info import SparsePauliOp
from qiskit.primitives import StatevectorEstimator
from scipy.optimize import minimize
import numpy as np
# STO-3G 기저의 H2 해밀토니안 (단순화)
hamiltonian = SparsePauliOp.from_list([
("II", -1.0523732),
("ZI", 0.3979374),
("IZ", -0.3979374),
("ZZ", -0.0112801),
("XX", 0.1809270),
])
ansatz = TwoLocal(2, ["ry", "rz"], "cx", reps=1)
estimator = StatevectorEstimator()
def cost(params):
job = estimator.run([(ansatz, hamiltonian, params)])
return float(job.result()[0].data.evs)
init = np.random.uniform(0, 2 * np.pi, ansatz.num_parameters)
res = minimize(cost, init, method="COBYLA",
options={"maxiter": 500})
print(f"바닥 상태 에너지: {res.fun:.6f} Hartree")
# 참값: -1.137270 Hartree
주요 응용 분야
- 양자화학: 분자 기저 상태 에너지, 전이 상태, 분자 스펙트럼 계산
- 재료과학: Hubbard 모델, 강상관 전자 시스템 시뮬레이션
- 조합 최적화: QAOA와 함께 MaxCut, 포트폴리오 최적화에 적용
현재의 한계
- 바렌 고원(Barren Plateau): 큐비트 수 증가 시 기울기가 지수적으로 감소하여 최적화가 어려워진다.
- 측정 비용: 파울리 항 수 에 비례하는 측정 횟수가 필요하다.
- 노이즈: 실제 하드웨어에서 에너지 기댓값에 계통 오차가 발생한다.
정리
VQE는 변분 원리를 하이브리드 구조로 구현한 알고리즘으로, 양자 하드웨어의 역할을 기댓값 추정기로 한정함으로써 NISQ 장치에서 실행 가능성을 확보한다. 안사츠 선택, 파라미터 이동 규칙을 통한 기울기 계산, 파울리 분해 기반 측정이 세 가지 핵심 구성 요소다. 바렌 고원과 측정 비용 문제가 남아 있으나, 오류 완화 기법·ADAPT-VQE·양자 자연 기울기(QNG) 등의 발전과 함께 분자 시뮬레이션에서 고전 한계를 돌파하기 위한 연구가 활발히 진행 중이다.
Exercises
연습문제
Q1변분 원리의 부등식 $\langle \psi | \hat{H} | \psi \rangle \geq E_0$이 등호가 될 조건을 수학적으로 서술하고, VQE에서 이 조건이 실제로 달성되기 어려운 이유를 안사츠의 관점에서 설명하라.
힌트 보기
고유 상태 전개를 이용해 기댓값을 쓴 뒤, 등호 조건을 확인하라. 안사츠가 실제 바닥 상태를 포함할 수 있는지 생각해 보라.
해설 보기
기댓값을 $\sum_n |\langle E_n|\psi\rangle|^2 E_n$으로 전개하면, 등호는 $|\langle E_0|\psi\rangle|^2 = 1$, 즉 $|\psi\rangle = e^{i\phi}|E_0\rangle$일 때만 성립한다. 현실에서는 안사츠가 표현할 수 있는 상태 공간(안사츠 다양체)이 제한되어 있어 바닥 상태가 그 안에 없으면 등호 달성이 불가능하다. 이를 안사츠 편향(ansatz bias) 또는 표현력 부족이라 한다.
Q2파라미터 이동 규칙 $\frac{\partial E}{\partial \theta_k} = \frac{E(\theta_k + \pi/2) - E(\theta_k - \pi/2)}{2}$를 유도하라. 단, 회전 게이트가 $U(\theta_k) = e^{-i\theta_k G}$ ($G^2 = I/4$인 파울리 생성자)로 주어진다고 가정한다.
힌트 보기
$\frac{\partial}{\partial \theta_k} e^{-i\theta_k G} = -iG\, e^{-i\theta_k G}$를 사용하고, $e^{-i\alpha G}$의 행렬 지수를 $G$의 고유값으로 전개하라.
해설 보기
$G$의 고유값이 $\pm 1/2$이므로 $e^{-i\theta G} = \cos(\theta/2)I - 2i\sin(\theta/2)G$로 쓸 수 있다. 편미분을 계산하면 $\partial_{\theta_k} E = \langle \psi | (-iG)^\dagger \hat{H} U + U^\dagger \hat{H} (-iG)|\cdots\rangle$이 되고, 이를 정리하면 정확히 $[E(\theta_k+\pi/2)-E(\theta_k-\pi/2)]/2$가 도출된다. 즉 수치 미분 없이 두 번의 회로 실행만으로 정확한 기울기를 얻는다.
Q3HEA(하드웨어 효율적 안사츠)와 UCCSD 안사츠를 회로 깊이, 물리적 동기, 노이즈 민감도 측면에서 비교하라. 어떤 상황에서 각각이 더 적합한가?
해설 보기
UCCSD는 결합 클러스터 이론에서 유도된 물리 동기를 가지므로 정확도가 높지만, 회로 깊이가 큐비트 수 $n$에 대해 $O(n^4)$ 이상으로 증가해 현재 NISQ 하드웨어에서 노이즈로 인한 오차가 크다. 반면 HEA는 장치 연결성을 직접 활용하는 얕은 회로로 노이즈에 강하지만, 표현력이 제한적이어서 안사츠 편향이 발생할 수 있다. 소규모 분자나 노이즈 수준이 낮은 장치에서는 UCCSD, 큐비트 수가 많거나 노이즈가 심한 환경에서는 HEA(또는 ADAPT-VQE)가 더 실용적이다.
관련 용어
