[논문 리뷰] Simulations of Many-Body Quantum Systems by a Quantum Computer

Stephen Wiesner|ArXiv.org|1996. 03. 26.

Quantum Computing Algorithms and Architecture인용 수 128

한 줄 요약

이 논문은 큐비트 얽힘을 활용하여 복잡한 양자 상태를 표현함으로써 양자 컴퓨터가 양자 다체계(예: 원자와 분자)를 효율적으로 시뮬레이션할 수 있다고 제안한다. 고전적 컴퓨터는 자원이 지수적으로 증가하는 반면, 양자 방법은 해밀토니안 분해를 통해 유니터리 진화 단계를 적용함으로써 다항식 시간 내에 시뮬레이션을 달성하여 정확하고 확장 가능한 양자 역학적 시뮬레이션을 가능하게 한다.

ABSTRACT

We suggest that quantum computers can solve quantum many-body problems that are impracticable to solve on a classical computer.

연구 동기 및 목표

• 고전적 컴퓨터가 처리할 수 없는 양자 다체계를 양자 컴퓨터가 효율적으로 시뮬레이션할 수 있음을 보여주기 위해.
• 고전적 시뮬레이션의 근본적 한계를 다루기 위해, 시스템 크기가 증가함에 따라 $2^N$차원 힐베르트 공간을 추적해야 하는 이유로 인해 지수적으로 증가하는 스케일링 문제를 해결하기 위해.
• 큐비트 기반 상태 표현과 시간 진화 연산자를 사용하여 일차 및 이차 상호작용 해밀토니안에 의해 지배되는 시스템을 위한 실용적인 양자 알고리즘을 제안하기 위해.
• 시뮬레이션 시간이 시스템 크기에 대해 다항식 스케일링됨을 보여주어 고전적 지수적 스케일링과 대비하기 위해.

제안 방법

• 각 시스템 내의 양자 입자를 $k$ 큐비트로 표현하여 공간 파동함수를 $2^k$개의 위치 빈으로 이산화하기.
• 소규모 시간 간격 $\Delta t$를 사용하여 각 해밀토니안 항(일차 및 이차 상호작용)에 대해 $I + iH\Delta t/\hbar$ 형태의 시간 진화 연산자를 순차적으로 적용하기.
• 빠른 푸리에 변환(FFTs)을 큐비트 레지스터에 적용하여 위치 및 운동량 표현 간 전환을 가능하게 하여 $P^2/2m$ 항의 효율적 진화를 실현하기.
• 위치에 따라 달라지는 위치 에너지 항 $V_{ij}$를 현재 큐비트 상태로부터 값 계산하고, $\Delta t V_{ij}/\hbar$ 비례하는 단계 이동을 적용하기.
• 시간 $T$ 동안의 역학을 시뮬레이션하기 위해 전체 시간 진화 절차를 $T/\Delta t$번 반복한 후 최종 상태를 측정하여 관측 가능량을 추출하기.
• 양자 중첩과 얽힘을 활용하여 전체 진폭 벡터를 명시적으로 저장하지 않고도 다체 파동함수를 자연스럽게 표현하고 진화시키기.

실험 결과

연구 질문

• RQ1 양자 컴퓨터는 고전적 컴퓨터보다 양자 다체계를 더 효율적으로 시뮬레이션할 수 있는가?
• RQ2 일차 및 이차 상호작용을 가진 $N$ 입자 시스템의 시간 진화를 시뮬레이션하는 데 필요한 계산 복잡도는 무엇인가?
• RQ3 큐비트 레지스터를 사용하여 상호작용하는 입자의 슈뢰딩거 방정식을 어떻게 효율적으로 구현할 수 있는가?
• RQ4 다체계의 지수적 힐베르트 공간을 다항식 시간 내에 양자 알고리즘을 사용해 시뮬레이션할 수 있는가?
• RQ5 얽힘은 양자 컴퓨터에서 양자 역학적 시뮬레이션을 효율적으로 가능하게 하는 데 어떤 역할을 하는가?

주요 결과

• 일차 및 이차 해밀토니안을 가진 $M$개의 $N$ 이중 상태 시스템에 대한 시뮬레이션은 고전적 지수적 스케일링을 피하기 위해 $N$에 대해 다항식 시간 내에 수행될 수 있다.
• 각 입자에 대해 $k$ 큐비트를 사용하고 1차원에서 $N$ 개의 입자를 다룰 경우, 총 시뮬레이션 시간은 $N$과 $k$에 대해 다항식 스케일링되며, 시간 단위당 구체적으로 $O(Nk \log k)$이다.
• FFT의 사용으로 운동량 공간에서의 진화가 효율적으로 가능해져, 각 입자당 운동 에너지 항의 비용을 $O(2^k)$에서 $O(k \log k)$로 감소시킨다.
• 이 방법은 위치 상태의 중첩을 자연스럽게 처리하고, 얽힌 구성에 기반한 단계 이동을 적용하여 양자 중첩성을 유지한다.
• 측정 결과는 시뮬레이션을 반복하고 평균을 내어 관측 가능량의 정확한 기대값을 도출할 수 있다.
• 이 접근법은 고전적 컴퓨터가 처리할 수 없는 $2^{Nk}$차원의 전체 진폭 벡터를 저장할 필요가 없어진다.