양자 상태 단층 촬영: 밀도 행렬의 완전한 재구성

개념 소개

고전적 CT 스캔이 여러 각도에서 X선을 조사하여 3차원 내부 구조를 재현하듯, 양자 상태 단층 촬영은 동일하게 준비된 앙상블을 다양한 측정 기저로 반복 관측하여 밀도 행렬 $\rho$를 결정한다. 단일 측정은 양자 상태를 붕괴시켜 부분 정보만 얻으므로, 동일 상태를 $N$회 재현·측정하는 앙상블 접근이 필수적이다.

QST는 양자 컴퓨팅 회로 검증, 얽힘 상태 품질 평가, 양자 오류 정정 코드워드 진단 등 실험적·공학적 장면에서 가장 기본적인 진단 도구로 쓰인다.

핵심 원리

밀도 행렬과 자유도

$n$-큐비트 계의 상태는 $2^n\times 2^n$ 밀도 행렬 $\rho$로 완전히 기술된다. $\rho$는 다음 세 조건을 만족해야 한다.

$$\rho ={\rho }^{†},\rho \ge 0,\mathrm{Tr}(\rho )=1$$

이를 만족하는 실수 자유도는 $4^n-1$개이다. 단일 큐비트($n=1$)에서는 자유도 3개로 블로흐 벡터 $(r_x,r_y,r_z)$가 상태를 결정한다.

$$\rho =\frac{1}{2}\left(I+r_{x}X+r_{y}Y+r_{z}Z\right)$$

순수 상태는 $|\mathbf{r}|=1$, 혼합 상태는 $|\mathbf{r}|<1$이다.

정보 완전 측정

측정 연산자 집합 $\{M_k\}$는 밀도 행렬 공간을 완전히 생성(informationally complete)해야 한다. 단일 큐비트에서 파울리 기저 $\{X,Y,Z\}$ 방향 측정이 이를 충족한다.

$$⟨P⟩=\mathrm{Tr}(\rho P),P\in \{X,Y,Z\}$$

각 기대값은 반복 측정의 통계 평균으로 추정한다.

선형 반전과 최대 가능도 추정

선형 반전(Linear Inversion): 측정 결과 벡터 $\mathbf{f}$와 측정 행렬 $A$로부터

$$\hat{\rho }=A^{-1}\mathbf{f}$$

를 직접 계산한다. 구현이 단순하지만, 유한 표본 통계 오차로 인해 $\hat{\rho }\ge 0$ 조건이 위반될 수 있다.

최대 가능도 추정(MLE): 측정 결과 $\{n_k\}$에 대한 로그 가능도 함수를 물리적 구속 조건 하에서 최대화한다.

$${\hat{\rho }}_{\mathrm{MLE}}=\underset{\rho \ge 0,\mathrm{Tr}\rho =1}{\mathrm{argmax}}\sum _k{n}_k\ln \mathrm{Tr}(M_k\rho )$$

MLE는 항상 물리적으로 유효한 밀도 행렬을 반환하며, 유한 샘플에서도 안정적이다.

$n$-큐비트 확장

$n$-큐비트 시스템에서는 파울리 텐서곱 $\{I,X,Y,Z{\}}^{\otimes n}$의 $4^n$개 원소를 기저로 사용한다. 필요한 독립 측정 설정 수가 $3^n$에 비례하여 지수적으로 증가하므로, 대규모 시스템에는 압축 센싱(compressed sensing) 또는 행렬 곱 상태(MPS) 단층 촬영 같은 효율적 기법이 요구된다.

예시·응용

단일 큐비트 QST (Python)

import numpy as np

I = np.eye(2)
X = np.array([[0,1],[1,0]], dtype=complex)
Y = np.array([[0,-1j],[1j,0]])
Z = np.array([[1,0],[0,-1]], dtype=complex)

def reconstruct_1qubit(rx, ry, rz):
    """파울리 기대값으로부터 밀도 행렬 선형 반전"""
    return 0.5 * (I + rx*X + ry*Y + rz*Z)

def fidelity(rho, sigma):
    """충실도 F(ρ, σ) 계산 (순수 상태 기준)"""
    return np.real(np.trace(rho @ sigma))

# |+⟩ 상태: ⟨X⟩=1, ⟨Y⟩=0, ⟨Z⟩=0
rho_est = reconstruct_1qubit(1.0, 0.0, 0.0)
rho_true = np.array([[0.5, 0.5],[0.5, 0.5]])

print("재구성 밀도 행렬:\n", rho_est)
print(f"충실도: {fidelity(rho_true, rho_est):.4f}")

충실도와 품질 지표

재구성된 $\hat{\rho }$와 목표 상태 ${\rho }_{\mathrm{target}}$ 사이의 충실도는

$$F(\rho ,\hat{\rho })={\left(\mathrm{Tr}\sqrt{\sqrt{\rho }\hat{\rho }\sqrt{\rho }}\right)}^2$$

로 정의되며, $F=1$이면 완전 일치이다. 이온 트랩·초전도 큐비트 플랫폼에서 단일 큐비트 QST 충실도 99% 이상이 보고되고 있다.

주요 응용 분야

• 게이트 검증: 양자 회로 출력 상태가 이상적 설계와 얼마나 일치하는지 정량화
• 얽힘 증인: 재구성된 $\rho$로부터 부분 전치(partial transpose) 기법 등으로 얽힘 검출
• 양자 오류 정정: 코드워드 상태의 순도(purity) $\mathrm{Tr}({\rho }^2)$ 추적
• 양자 공정 단층 촬영(QPT): 채널 $\mathcal{E}$를 QST를 기초로 확장하여 완전한 과정 행렬 재구성

정리

QST는 밀도 행렬을 실험적으로 결정하는 가장 근본적인 양자 진단 프로토콜이다. 단일 큐비트에서는 파울리 3축 측정만으로 세 실수 파라미터를 복원하며, $n$-큐비트에서는 $4^n-1$개의 파라미터가 필요하다는 지수적 복잡도가 핵심 도전 과제다. 선형 반전은 단순하지만 물리적 유효성을 보장하지 못하고, MLE는 이를 해결하는 표준 접근이다. 압축 센싱과 구조화된 앙사츠 기법은 대규모 시스템에서 이 한계를 극복하는 현재진행형 연구 방향이다.