만약 당신이 Rust 개발자라면, 양자 컴퓨팅(quantum computing)은 대부분 다른 사람의 언어로 이루어져 왔을 것입니다. 사람들이 연구용으로 사용하는 시뮬레이터들은 Python이나 C++로 구축되어 있습니다. Rust 생태계를 떠나지 않고 양자 회로(quantum circuits)를 탐구하고 싶었다면, 선택지는 매우 적었습니다. 외부 런타임(runtimes)에 대한 바인딩(Bindings), 버려진 크레이트(crates), 혹은 stable 버전에서 컴파일되지 않는 학술용 코드들 말입니다.

Sirraya QuTub은 Rust로 구축된 완전한 양자 시뮬레이션 스택입니다. 상태 벡터(State vectors), 밀도 행렬(Density matrices), 실제 하드웨어에 맞춰 교정된 노이즈 채널(Noise channels), 문헌에 따라 구현된 얽힘 측정치(Entanglement measures)까지 모두 하나의 크레이트에 담겨 있습니다. cargo add sirraya-qutub 명령어로 추가하기만 하면 끝납니다.

하지만 질문은 또 다른 양자 시뮬레이터가 존재하는지 여부가 아닙니다. 질문은 당신이 자신의 언어를 떠나지 않고도 소스 코드를 읽고, 각 레이어가 무엇을 하는지 이해하며, 스스로 확장할 수 있느냐 하는 것입니다.

이 포스트에서는 단 하나의 Rust 프로그램을 통해 다음 과정을 살펴봅: Bell state 생성, 측정, 밀도 행렬 검토, 얽힘(entanglement) 정량화, 그리고 큐비트(qubit)를 잃었을 때 어떤 일이 발생하는지 관찰하기. 물리학적 배경 지식은 필요하지 않습니다. 출력되는 모든 숫자에 대해 설명을 제공합니다.

프로그램

Bell state는 가장 단순한 얽힌 양자 상태(entangled quantum state)입니다: 두 개의 큐비트가 완벽하게 상관되어 있는 상태를 말합니다. 하나를 측정하면 즉시 다른 하나의 상태를 알 수 있습니다. 이는 양자 컴퓨팅의 "hello world"와 같습니다.

전체 프로그램은 다음과 같습니다:

전체 프로그램은 다음과 같습니다:

use sirraya_qutub::core::{create_bell_state, DensityMatrix};
fn main() -> Result<(), String> {
let bell = create_bell_state()?;
bell.print_state();
let distribution = bell.get_probability_distribution();
for (state, prob) in &distribution {
println!(" |{}⟩: {:.4}", state, prob);
}
let density = DensityMatrix::from_state_vector(bell.get_state_vector())?;
density.print_density_matrix();
println!(" Purity: {:.4} (1.0 = pure)", density.purity());
println!(" Concurrence: {:.4} (1.0 = max entangled)", density.concurrence()?);
println!(" Negativity: {:.4} (>0 certifies entanglement)", density.negativity(&[1])?);
let reduced = density.partial_trace(&[0])?;
println!(" Reduced purity: {:.4} (0.5 = maximally mixed)", reduced.purity());
println!(" Entropy: {:.4} bits (1.0 = one ebit)", reduced.von_neumann_entropy());
Ok(())
}

이것으로 끝입니다. 임포트와 출력문까지 포함하여 50줄입니다. 이것이 무엇을 생성하는지, 그리고 각 숫자가 무엇을 의미하는지 설명합니다.

출력 결과 (Line by Line)

상태 벡터:

|00⟩: 0.707107 (probability: 0.500000)
|11⟩: 0.707107 (probability: 0.500000)

0.707107은 $1/
oot{2}{}$입니다. 이 상태는 $|00
angle$와 $|11
angle$의 동일한 중첩(superposition) 상태입니다. 두 결과 모두 동등하게 발생할 확률을 가지며, $|01
angle$나 $|10
angle$은 절대 볼 수 없습니다. 두 큐비트가 연결되어 있습니다.

밀도 행렬 (Density matrix)

0.500000 0.000000 0.000000 0.500000  
0.000000 0.000000 0.000000 0.000000  
0.000000 0.000000 0.000000 0.000000  
0.500000 0.000000 0.000000 0.500000

상태 벡터는 진폭(amplitudes)을 알려줍니다. 밀도 행렬은 모든 것을 알려줍니다: 대각선 요소에는 확률이, 비대각선 요소에는 양자 결맞음(quantum coherence)이 담겨 있습니다. 구석에 있는 0.5 값들은 일종의 결맞은 중첩 상태를 나타냅니다—즉, 이 상태는 두 결과 사이의 고전적인 동전 던지기가 아니라 $|00
angle$와 $|11
angle$가 동시에 존재하는 것입니다.

얽힘 측정값 (Entanglement measures):

• Purity: 1.0000 (1.0 = pure)
• Concurrence: 1.0000 (1.0 = maximally entangled)
• Negativity: 0.5000 (>0 certifies entanglement)

세 가지 독립적인 측정값, 세 가지 확인. 1.0의 순도 (Purity)는 상태를 완벽하게 알고 있음을 의미합니다. 즉, 노이즈가 없고 누락된 정보가 없다는 뜻입니다. 1.0의 Concurrence (동시성)는 2-큐비트 얽힘 (two-qubit entanglement)의 이론적 최댓값입니다. 0.5의 Negativity (부성)는 부분 전치 고유값 (partial transpose eigenvalues) {0.5, 0.5, 0.5, -0.5}에서 비롯됩니다. 음수 부분의 합은 얽힘이 존재함을 알려주며 그 양을 정량화합니다.

큐비트 1을 손실한 후:

감소된 순도 (Reduced purity): 0.5000 (0.5 = maximally mixed)
엔트로피 (Entropy): 1.0000 bits (1.0 = one ebit of entanglement)

이것이 핵심 결과입니다. 순수한 Bell 상태 (Bell state)로 시작합니다. 하나의 큐비트를 트레이스 아웃 (trace out)합니다. 물리적으로는 그 큐비트에 대한 접근 권한을 잃는 것입니다. 남은 큐비트는 이제 최대 혼합 상태 (maximally mixed)가 됩니다. 당신은 그것에 대해 아무것도 알 수 없습니다. 그것은 1.0 bits의 엔트로피를 가진 완벽한 50/50 동전 던지기와 같으며, 이는 단일 큐비트가 가질 수 있는 가능한 최대치입니다.

순수한 얽힘 상태 (pure entangled state)는 혼합된 축소 상태 (mixed reduced state)를 가집니다. 전체는 확정적이지만, 부분은 최대 불확실성을 가집니다. 정보는 어느 한 쪽 큐비트에 저장되어 있지 않습니다. 그것은 큐비트들 사이의 상관관계 (correlations)에 저장되어 있습니다. 이것이 양자 얽힘 (quantum entanglement)을 한 문장으로 정의한 것입니다.

Rust 개발자에게 이것이 중요한 이유

이 과정 중 그 어느 것도 밑바탕에 깔린 선형 대수학 (linear algebra)을 이해할 필요를 요구하지 않았습니다. API는 물리적 개념을 직접적으로 노출합니다: purity(), concurrence(), negativity(), partial_trace(), von_neumann_entropy(). 각 메서드 이름은 문헌 용어와 일치합니다. 각 메서드는 명확한 물리적 해석이 가능한 숫자를 반환합니다.

이것이 설계 철학입니다: 만약 양자 연산이 시뮬레이션되고 있다면, 그 구현은 이해 가능하고, 검사 가능하며, 재현 가능해야 합니다. 시뮬레이터는 블랙박스가 아닙니다.

Rust 개발자들에게 이는 다른 Rust 크레이트 (crate)를 탐색하는 것과 같은 방식으로 양자 컴퓨팅을 탐구할 수 있음을 의미합니다. 소스 코드를 읽으세요. 테스트를 실행하세요. 예제를 수정하세요. 게이트 (gates)를 확장하세요. 노이즈를 추가하고 얽힘에 어떤 일이 일어나는지 관찰하세요. GHZ 상태를 구축하고 이를 W 상태와 비교해 보세요. Ising 결합 게이트 (Ising coupling gates)를 사용하여 VQE ansatz를 구현해 보세요. 진입 장벽은 cargo add sirraya-qutub 명령어와 호기심뿐입니다.

향후 계획

Bell 상태 (Bell state)는 아주 작습니다. 하지만 그 이면에 담긴 아이디어는 결코 작지 않습니다.

QuTub에서 중요한 것은 단순히 또 하나의 시뮬레이터(simulator)를 갖는 것이 아닙니다. 개발자, 연구자, 그리고 학생들이 직접 검사하고 그 위에 구축할 수 있는 투명한 토대(foundation)를 갖는 것입니다. 현재 이 크레이트(crate)는 전체 상태 벡터 시뮬레이터 (full state-vector simulator), 하드웨어 교정 채널을 포함한 밀도 행렬 노이즈 시뮬레이션 (density-matrix noise simulation), 문헌에 기반한 구현을 갖춘 얽힘 측정 (entanglement measures), VQE를 위한 Pauli-string 관측량 (observables), 양자 푸리에 변환 (quantum Fourier transform), 양자 리저버 컴퓨팅 (quantum reservoir computing), 그리고 교차 엔트로피 벤치마킹 검증 (cross-entropy benchmarking validation) 기능을 제공합니다.

다음 주요 단계는 QuTub 트랜스파일러 (transpiler)입니다. 이는 추상적인 양자 게이트 (quantum gates)를 네이티브 하드웨어 게이트 세트 (native hardware gate sets)로 분해하고, 결과로 나온 회로를 최적화하며, 실제 하드웨어 충실도 (fidelity) 수치에 맞춰 교정된 의미 있는 게이트별 노이즈 시뮬레이션을 가능하게 하는 계층입니다. 해당 계층은 현재 Sirraya Labs에서 활발히 개발 중입니다.

우리는 이를 한 단계씩, 공개적으로 구축해 나가고 있습니다.

Sirraya QuTub은 crates.io에서 사용할 수 있습니다. 전체 문서는 docs.rs/sirraya-qutub에서 확인하세요. MIT 라이선스 하에 제공됩니다.