표면 부호(Surface Code)의 구조와 안정자 형식론

개념 소개

표면 부호는 Kitaev의 토릭 부호(toric code)를 평면 경계 조건으로 수정한 안정자 부호다. 핵심 특징은 근접 이웃 결합만으로 구현 가능하다는 점이다. 대다수의 오류 정정 부호가 원거리 큐비트 사이의 게이트를 요구하는 반면, 표면 부호는 격자 위 인접 큐비트 사이의 상호작용만 사용한다. 이 덕분에 초전도 회로, 이온 트랩 등 2차원 배치가 자연스러운 하드웨어와 잘 부합하며, Google과 IBM 모두 대규모 내결함성 양자컴퓨터 로드맵에 표면 부호를 채택하고 있다.

핵심 원리

안정자 형식론

표면 부호는 안정자 부호(stabilizer code)의 일종이다. $n$개 물리 큐비트의 힐베르트 공간에서, 파울리 군 ${\mathcal{G}}_n$의 교환 가능한 부분군 $\mathcal{S}$(안정자 군)를 정의하고, 모든 원소에 대해 $+1$ 고유값을 갖는 상태들의 부분공간을 논리 부호공간으로 삼는다:

$$\mathcal{C}=\{|\psi ⟩|S|\psi ⟩=+|\psi ⟩,\forall S\in \mathcal{S}\}$$

$n$개의 독립 안정자 생성원이 있을 때 $k=n-|\text{생성원}|$개의 논리 큐비트가 인코딩된다.

회전 표면 부호의 격자 구조

실용적으로 가장 많이 쓰이는 **회전 표면 부호(rotated surface code)**는 $d\times d$ 격자 절점에 데이터 큐비트를 배치한다. 거리(code distance) $d$인 부호의 구성 요소는 다음과 같다.

구성 요소	수량
데이터 큐비트	$d^2$
X형 안정자 (보조 큐비트)	$(d^2-1)/2$
Z형 안정자 (보조 큐비트)	$(d^2-1)/2$
논리 큐비트	1

격자를 체스판으로 시각화하면, 흰 셀마다 X형 안정자, 검은 셀마다 Z형 안정자가 대응된다. 내부 셀의 안정자는 네 꼭짓점의 데이터 큐비트에 파울리 연산자를 동시에 적용하는 4체 연산자이고, 경계 셀은 2체 연산자다:

$$S_X^{(p)}=X_a{X}_b{X}_c{X}_d,S_Z^{(q)}=Z_a{Z}_b{Z}_c{Z}_d$$

임의의 X형 안정자와 Z형 안정자는 격자 위에서 공유하는 큐비트 수가 항상 짝수(0 또는 2)이므로 서로 교환 가능하다($[S_X,S_Z]=0$). 이것이 군 $\mathcal{S}$의 교환 가능성 조건을 만족시킨다.

증후군 측정과 오류 검출

보조 큐비트를 CNOT 회로로 데이터 큐비트와 결합한 뒤 보조 큐비트를 측정하면, 데이터 큐비트를 파괴하지 않고도 각 안정자의 고유값 $\pm 1$을 얻을 수 있다. 이를 **증후군 측정(syndrome measurement)**이라 한다.

• 오류 없음 → 모든 안정자 $+1$
• 데이터 큐비트 $i$에 $Z$ 오류 발생 → 큐비트 $i$를 포함하는 X형 안정자만 $-1$로 반전

$Z_i$는 $X_i$와 반교환($\{X,Z\}=0$)하므로, 해당 X형 안정자의 측정값이 뒤집히는 것이다. Z형 안정자는 $Z_i$와 교환하므로 영향받지 않는다. 반대로 $X$ 오류는 Z형 안정자만 들뜨게 한다.

디코더(decoder)는 $-1$을 반환한 안정자 집합(증후군)을 입력받아 가장 가능성 높은 오류 패턴을 추론하고 수정 연산을 적용한다. 가장 대표적인 디코딩 알고리즘은 **최소 가중치 완전 매칭(minimum-weight perfect matching, MWPM)**이다.

논리 연산자와 부호 거리

논리 파울리 연산자 $\bar{X}$와 $\bar{Z}$는 안정자 군의 원소가 아니면서 부호공간을 보존하는 연산자다. 회전 표면 부호에서 이들은 각각 격자의 한 경계에서 반대 경계까지 이어지는 파울리 체인이다.

$$\bar{X}=X_{i_1}{X}_{i_2}\cdots X_{i_d},\bar{Z}=Z_{j_1}{Z}_{j_2}\cdots Z_{j_d}$$

이 체인의 최소 길이, 즉 논리 연산자의 최소 무게가 부호 거리 $d$다. 거리 $d$ 부호는 최대 $\lfloor (d-1)/2\rfloor$개의 임의 단일 큐비트 오류를 정정할 수 있다.

예시·응용

거리-3 표면 부호 시뮬레이션

# Stim 라이브러리로 거리-3 회전 표면 부호 회로 생성
import stim

circuit = stim.Circuit.generated(
    "surface_code:rotated_memory_z",
    rounds=3,          # 증후군 측정 반복 횟수
    distance=3,
    after_clifford_depolarization=0.001,  # 물리 오류율 0.1%
)
# 샘플러로 측정 결과(증후군) 수집
sampler = circuit.compile_detector_sampler()
detection_events, observable_flips = sampler.sample(
    shots=1000, separate_observables=True
)
print(f"검출 이벤트 형상: {detection_events.shape}")

$d=3$이면 데이터 큐비트 9개, 안정자 8개로 단 하나의 논리 큐비트를 인코딩하며, 임의의 단일 큐비트 오류 1개를 정정할 수 있다.

오류 역치와 물리 큐비트 오버헤드

물리 오류율 $p$가 임계값 $p_{th}$ 이하일 때, 거리 $d$를 늘릴수록 논리 오류율이 지수적으로 감소한다:

$$p_L\approx A{\left(\frac{p}{p_{th}}\right)}^{\lfloor (d+1)/2\rfloor }$$

디폴라라이징 잡음 모형에서 수치 시뮬레이션으로 추정한 표면 부호의 임계값은 $p_{th}\approx 1\%$다. 현재 초전도 소자의 게이트 오류율이 $\sim 0.1\%$ 수준이므로 원칙적으로 역치 이하 동작이 가능하다.

그러나 오버헤드가 크다. 거리 $d$ 표면 부호는 총 $2d^2-1$개의 물리 큐비트로 논리 큐비트 1개를 구현하며, 실용적인 알고리즘 수행을 위해 $d\sim 20$$30$이 필요하다고 추산되므로 논리 큐비트 하나당 최소 수백수천 개의 물리 큐비트가 필요하다.

정리

표면 부호는 2차원 격자의 X형·Z형 안정자로 하나의 논리 큐비트를 인코딩하고, 증후군 측정을 통해 비파괴적으로 오류를 검출·정정한다. 근접 이웃 결합, 높은 오류 임계값($\sim 1\%$), MWPM 디코딩의 효율성이 결합되어 현재 내결함성 양자컴퓨터의 가장 유력한 후보 아키텍처로 자리잡았다. 부호 거리 $d$를 늘릴수록 논리 오류율은 지수적으로 감소하지만, 물리 큐비트 수는 $O(d^2)$으로 증가하는 트레이드오프가 핵심 공학 과제로 남아 있다.