밀도 상태 계산의 계산 복잡도

초록

이 논문은 지역 해밀토니안의 바닥 상태 퇴화도와 에너지 밀도 상태를 계산하는 문제의 복잡도를 조사한다. 저자들은 이를 양자 복잡도 클래스 QMA의 카운팅 버전인 #BQP로 정확히 규정하고, #BQP가 고전적인 카운팅 클래스 #P보다 더 어렵지 않음을 증명한다. 따라서 고전 해밀토니안의 경우에도 양자 해밀토니안과 동일한 난이도를 가진다.

상세 분석

논문은 먼저 “밀도 상태(Density of States, DOS)”와 “바닥 상태 퇴화도(Ground State Degeneracy, GSD)”라는 두 물리량을 정확히 정의한다. DOS는 특정 에너지 구간에 존재하는 고유값의 개수를, GSD는 최저 에너지 고유값의 중복도를 의미한다. 두 문제 모두 지역 해밀토니안(H)이라는 입력을 받아, H의 스펙트럼 정보를 추출해야 하는데, 이는 일반적인 양자역학 시뮬레이션에서 가장 어려운 작업 중 하나이다.

복잡도 이론적 관점에서 저자들은 기존에 널리 사용되는 QMA(Quantum Merlin‑Arthur)와 그 카운팅 버전인 #QMA를 검토한다. QMA는 “양자 증명”을 검증하는 문제를 다루며, #QMA는 해당 증명의 수를 세는 문제이다. 그러나 #QMA는 정의상 비결정론적 양자 회로의 출력 개수를 세는 것이므로, 실제 물리적 해밀토니안 문제와 직접 연결하기엔 부적절했다. 이를 보완하기 위해 저자들은 #BQP라는 새로운 클래스를 도입한다. #BQP는 “양자 다항시간(Quantum Polynomial time) 알고리즘이 출력하는 확률 분포의 가중 합을 세는” 문제를 의미한다. 구체적으로, 양자 회로 C가 입력 x에 대해 accept 확률 p(x)를 반환할 때, #BQP는 모든 x에 대해 p(x)의 합을 구하는 작업이다.

핵심 정리는 두 단계로 전개된다. 첫 번째 단계에서는 GSD와 DOS 문제를 #BQP‑완전(complete)임을 보인다. 이를 위해 저자들은 임의의 #BQP 인스턴스를 지역 해밀토니안 H에 효율적으로 매핑하는 다항시간 변환을 설계한다. 변환 과정에서 “에너지 창(energy window)”을 도입해, 특정 에너지 구간에 해당하는 고유값만을 선택적으로 포함하도록 H를 구성한다. 이렇게 하면 #BQP 인스턴스의 카운팅 값이 바로 H의 바닥 상태 퇴화도 혹은 해당 구간 내 DOS와 일치한다. 반대로, GSD와 DOS를 #BQP 문제로 환원함으로써 두 물리 문제의 하드니스가 #BQP와 동등함을 증명한다.

두 번째 단계에서는 #BQP와 고전적인 카운팅 클래스 #P 사이의 관계를 분석한다. 저자들은 양자 회로를 고전적인 비결정론적 회로로 시뮬레이션하는 “양자‑고전 변환” 기법을 사용한다. 구체적으로, 양자 회로의 각 게이트를 고전적인 논리 게이트와 추가적인 확률 변수로 분해하고, 전체 회로의 accept 확률을 다항시간 내에 #P‑함수로 표현한다. 이를 통해 #BQP ⊆ #P임을 보이며, 역방향 포함 관계는 이미 알려진 #P ⊆ #BQP (양자 회로가 고전 회로를 포함한다는 사실)와 결합해 #BQP = #P임을 결론짓는다.

이 결과는 물리학적 의미가 크다. 고전적인 스핀 시스템(예: 이징 모델)과 양자 스핀 시스템(예: 하이젠베르크 모델) 모두에 대해, 바닥 상태 퇴화도와 DOS를 정확히 계산하는 문제는 동일한 복잡도 난이도를 가진다. 즉, 양자 상호작용이 추가된다고 해서 계산 난이도가 급격히 상승하지 않으며, 기존의 #P‑완전 문제와 동등한 수준임을 보여준다. 이는 양자 물질의 스펙트럼 분석이 이론적으로는 고전적인 카운팅 문제와 같은 어려움을 가진다는 중요한 통찰을 제공한다.

또한, 저자들은 실용적인 함의도 논한다. 현재 양자 시뮬레이션 알고리즘은 근사적인 스펙트럼 추정에 초점을 맞추고 있지만, 정확한 DOS 계산은 #P‑완전 문제와 동등하므로, 효율적인 근사 알고리즘이 존재하더라도 정확한 해답을 얻는 것은 일반적인 다항시간 알고리즘으로는 불가능할 가능성이 높다. 따라서 연구자들은 문제의 난이도를 인식하고, 근사도와 오류 한계를 명시적으로 관리해야 함을 강조한다.