페르미온–보손 매핑: Jordan–Wigner 변환

Jordan–Wigner 변환은 페르미온 생성·소멸 연산자를 스핀(큐비트) 연산자로 정확히 변환하는 대수적 사상(mapping)이다. 이 변환을 통해 전자 구조 계산이나 허바드 모델 같은 페르미온 해밀토니안을 양자 컴퓨터의 파울리 연산자 체계로 표현할 수 있다. 변환 과정에서 비국소적 스트링 연산자가 등장하는 것이 핵심 비용이다.

개념 소개

양자화학·응집물질 물리의 중요한 모델들은 대부분 페르미온 입자(전자, 쿼크 등)로 기술된다. 페르미온은 반교환(anti-commutation) 관계를 만족한다.

${c_{i}, c_{j}^{†}} = δ_{ij}, {c_{i}, c_{j}} = 0$

반면 양자 컴퓨터의 기본 소자인 큐비트는 스핀-1/2 대수, 즉 파울리 연산자 ${X, Y, Z}$ 로 표현된다. 두 대수 체계 사이를 잇는 것이 Jordan–Wigner(JW) 변환이다. 1928년 Pascual Jordan과 Eugene Wigner가 1차원 페르미온 시스템을 스핀 체인으로 변환하기 위해 처음 도입하였다.

핵심 원리

점유수 표현과 큐비트 대응

모드 $i$ 의 점유 여부를 $∣ n_{i} ⟩ \in {∣0 ⟩, ∣1 ⟩}$ 으로 나타낸다. 단일 모드에서는 다음 대응이 자연스럽다.

$c_{i} \leftrightarrow \frac{X _{i} + i Y _{i}}{2} = σ_{i}^{-}, c_{i}^{†} \leftrightarrow \frac{X _{i} - i Y _{i}}{2} = σ_{i}^{+}$

그러나 다중 모드에서 페르미온의 반교환 관계는 서로 다른 모드 간 부호를 요구한다. 스핀 연산자는 서로 다른 자리(site)에서 교환하므로 이 부호를 직접 재현하지 못한다.

Jordan–Wigner 스트링

이를 해결하기 위해 JW 변환은 사이트 $i$ 앞에 놓인 모든 모드의 점유수 패리티를 누적하는 스트링 연산자를 도입한다.

$c_{j} = (i < j \prod Z_{i}) σ_{j}^{-}$

$c_{j}^{†} = (i < j \prod Z_{i}) σ_{j}^{+}$

여기서 $\prod_{i < j} Z_{i}$ 가 Jordan–Wigner 스트링이다. $Z_{i}$ 는 해당 모드가 점유( $∣1 ⟩$ )되면 $- 1$ , 비어( $∣0 ⟩$ )있으면 $+ 1$ 을 곱하므로 정확히 반교환 부호를 생성한다.

반교환 관계 검증

$j < k$ 일 때:

${c_{j}, c_{k}^{†}} = 0 (j \neq = k)$

을 위 정의로 직접 계산하면, $j$ 위치의 $σ^{-}$ 와 $k$ 위치 스트링 내의 $Z_{j}$ 사이에서 ${Z, σ^{-}} = 0$ 관계 덕분에 전체 반교환이 성립함을 확인할 수 있다.

해밀토니안 변환: 1차원 허바드 호핑 항

$H = - t \sum_{j} (c_{j}^{†} c_{j + 1} + h.c.)$ 에 JW를 적용하면:

$c_{j}^{†} c_{j + 1} = σ_{j}^{+} Z_{j} \cdot Z_{j} σ_{j + 1}^{-} = σ_{j}^{+} σ_{j + 1}^{-}$

인접 항은 스트링이 상쇄되어 국소적 파울리 표현이 된다. 그러나 장거리 호핑 $c_{j}^{†} c_{k}$ ( $∣ j - k ∣ > 1$ )에서는 $Z_{j + 1} \dots Z_{k - 1}$ 이 남아 $O (N)$ 길이 파울리 스트링이 발생한다.

예시·응용

수소 분자( $H_{2}$ ) 양자 시뮬레이션

최소 기저(STO-3G)에서 $H_{2}$ 는 4개의 스핀 궤도를 갖는다. JW 변환 후 해밀토니안은 다음과 같은 형태가 된다.

$H = g_{0} I + g_{1} Z_{0} + g_{2} Z_{1} + g_{3} Z_{0} Z_{1} + g_{4} Y_{0} Y_{1} + g_{5} X_{0} X_{1}$

계수 $g_{i}$ 는 단전자·이전자 적분으로 결정되며, 이 파울리 표현에 VQE(변분 양자 고유값 분해) 알고리즘을 직접 적용한다.

# OpenFermion + Qiskit을 이용한 H2 JW 변환 예시 (개념 코드)
from openfermion import FermionOperator, jordan_wigner
from openfermion.utils import count_qubits

# 간단한 2-모드 호핑 항 생성
hop = FermionOperator('0^ 1', -1.0) + FermionOperator('1^ 0', -1.0)

# JW 변환
jw_op = jordan_wigner(hop)
print(jw_op)
# 출력 예: -0.5 [X0 X1] + -0.5 [Y0 Y1]

스케일링과 한계

항목	Jordan–Wigner	Bravyi–Kitaev
스트링 길이	$O (N)$	$O (lo g N)$
큐비트 수	$N$	$N$
구현 복잡도	낮음	중간

장거리 상호작용이 많을수록 JW 스트링의 깊이 비용이 커져, 대형 분자에서는 Bravyi–Kitaev 변환이나 페르미온-프리 부분 공간 기법이 대안으로 사용된다.

정리

Jordan–Wigner 변환은 페르미온 대수의 반교환 관계를 파울리 Z 스트링으로 인코딩하는 정확한 동형(isomorphism)이다. 1차원 국소 모델에서는 스트링이 상쇄되어 효율적이지만, 비국소 항이 많은 경우 회로 깊이가 $O (N)$ 으로 증가한다. 양자화학 해밀토니안의 양자 컴퓨터 구현, VQE, 위상 추정 알고리즘의 전처리 단계로서 핵심적인 역할을 한다.

연습문제

Q1.단일 모드 페르미온 연산자 $c$에 대해 $\{c, c^\dagger\} = 1$이 $\sigma^-, \sigma^+$ 표현에서도 성립함을 직접 계산으로 보여라.

힌트 보기

$\sigma^+ = |0\rangle\langle 1|$, $\sigma^- = |1\rangle\langle 0|$로 놓고 $\sigma^-\sigma^+ + \sigma^+\sigma^-$를 계산한다.

해설 보기

$\sigma^-\sigma^+ = |1\rangle\langle 0|\cdot|0\rangle\langle 1| = |1\rangle\langle 1|$, $\sigma^+\sigma^- = |0\rangle\langle 1|\cdot|1\rangle\langle 0| = |0\rangle\langle 0|$. 합하면 $|0\rangle\langle 0| + |1\rangle\langle 1| = I$. 따라서 $\{c, c^\dagger\} = 1$ 성립.

Q2.2-모드 페르미온 호핑 $c_0^\dagger c_1$을 JW 변환으로 파울리 연산자 합으로 표현하라.

힌트 보기

$c_0^\dagger = \sigma_0^+$, $c_1 = Z_0 \sigma_1^-$ 임을 이용하고, $\sigma^+ = \frac{X-iY}{2}$, $\sigma^- = \frac{X+iY}{2}$로 전개한다.

해설 보기

$c_0^\dagger c_1 = \sigma_0^+ \cdot Z_0 \sigma_1^-$. $Z_0 \sigma_0^+ = Z_0 \frac{X_0-iY_0}{2}$에서 $ZX = iY$, $ZY = -iX$이므로 $Z_0\sigma_0^+ = -\sigma_0^+$. 따라서 $c_0^\dagger c_1 = -\sigma_0^+\sigma_1^- = -\frac{(X_0-iY_0)}{2}\cdot\frac{(X_1+iY_1)}{2} = -\frac{1}{4}(X_0X_1 + iX_0Y_1 - iY_0X_1 + Y_0Y_1)$. h.c.를 더하면 호핑 항은 $-\frac{1}{2}(X_0X_1 + Y_0Y_1)$이 된다.

Q3.3-모드 시스템에서 $c_0^\dagger c_2$ 항을 JW 변환하면 어떤 파울리 스트링이 나타나는가? 국소성 측면에서 어떤 의미를 갖는가?

해설 보기

$c_0^\dagger c_2 = \sigma_0^+ \cdot (Z_0 Z_1) \sigma_2^-$. $Z_0\sigma_0^+ = -\sigma_0^+$이므로 $c_0^\dagger c_2 = -\sigma_0^+ Z_1 \sigma_2^-$. 파울리 표현으로 전개하면 $-\frac{1}{2}(X_0 Z_1 X_2 + Y_0 Z_1 Y_2)$ 형태의 3-체 스트링이 된다. 중간 사이트 $Z_1$이 비국소 연산으로 남아 회로 구현 시 CNOT 게이트가 추가로 필요하다. 이것이 JW 변환의 $O(N)$ 비용 문제의 직접적 예시다.

개념 소개

핵심 원리

점유수 표현과 큐비트 대응

Jordan–Wigner 스트링

반교환 관계 검증

해밀토니안 변환: 1차원 허바드 호핑 항

예시·응용

수소 분자(H2​) 양자 시뮬레이션

스케일링과 한계

정리

연습문제

관련 용어

수소 분자( $H_{2}$ ) 양자 시뮬레이션