[논문 리뷰] Monte Carlo Search for Very Hard KSAT Realizations for Use in Quantum Annealing

Neuhaus Thomas|arXiv (Cornell University)|2014. 12. 17.

Quantum Computing Algorithms and Architecture참고 문헌 3인용 수 5

한 줄 요약

이 논문은 고유의 만족 가능한 할당을 가진 K-만족성 문제(KSAT)에서 가장 어려운 인스턴스를 최소한의 절차-스핀 비율(α ≈ 1)에서 식별하기 위해 다중칸노니컬 엔세임블 기법을 사용한 몬테카를로 검색 방법을 제안한다. 이는 에너지 E=1에서의 상태 밀도 Ω₁이 Ω₁ ∝ exp[rN] 형태로 지수적으로 발산하며, 비율 상수 r은 K=2일 때 0.348에서 K=6일 때 0.680으로 증가하여 ln 2에 수렴함을 규명한다. 이는 양자 애너일링 효율성 테스트를 위한 기준 엔세임블을 제공한다.

ABSTRACT

Using powerful Multicanonical Ensemble Monte Carlo methods from statistical physics we explore the realization space of random K satisfiability (KSAT) in search for computational hard problems, most likely the 'hardest problems'. We search for realizations with unique satisfying assignments (USA) at ratio of clause to spin number $α=M/N$ that is minimal. USA realizations are found for $α$-values that approach $α=1$ from above with increasing number of spins $N$. We consider small spin numbers in $2 \le N \le 18$. The ensemble mean exhibits very special properties. We find that the density of states of the first excited state with energy one $Ω_1=g(E=1)$ is consistent with an exponential divergence in $N$: $Ω_1 \propto { m exp} [+rN]$. The rate constants for $K=2,3,4,5$ and $K=6$ of KSAT with USA realizations at $α=1$ are determined numerically to be in the interval $r=0.348$ at $K=2$ and $r=0.680$ at $K=6$. These approach the unstructured search value ${ m ln}2$ with increasing $K$. Our ensemble of hard problems is expected to provide a test bed for studies of quantum searches with Hamiltonians that have the form of general Ising models.

연구 동기 및 목표

• 무작위 K-만족성 문제(KSAT)에서 고유의 만족 가능한 할당(USA)을 가진 가장 계산적으로 어려운 인스턴스를 식별하는 것.
• 절차-스핀 비율 α = M/N가 1보다 약간 큰 방향에서 최소 비율에 도달할 때 KSAT의 실현 공간을 탐색하는 것.
• 첫 번째(excited) 에너지 수준(E=1)에서의 상태 밀도 Ω₁ = g(E=1)를 특성화하고, 시스템 크기 N에 따른 지수적 발산을 정량화하는 것.
• 양자 어드아바틱 계산 및 양자 애너일링 효율성 테스트를 위한 어려운 문제 엔세임블을 생성하는 것.
• 알려진 복잡도 스케일링을 갖는 물리적으로 타당한, 계산적으로 도전적인 양자 검색 알고리즘 테스트베드를 제공하는 것.

제안 방법

• 통계역학의 다중칸노니컬 엔세임블 몬테카를로 방법을 사용하여 고유의 만족 가능한 할당(USA)을 가진 희귀한, 어려운 KSAT 실현을 샘플링한다.
• 스핀 매핑 sᵢ = 2Xᵢ − 1을 사용하여 불리안 변수를 이징 스핀으로 변환하고, E=0에서 지구 상태를 갖는 이징 해밀토니안 H_KSAT를 구성한다.
• 제약 조건을 설정: (1) 고유의 지구 상태(g(E=0)=1), (2) 최소 α = M/N > 1, (3) E=1에서의 위상공간 부피 최대화(Ω₁).
• 작은 시스템(2 ≤ N ≤ 18)에 대해 수치 시뮬레이션을 수행하고 상태 밀도 g(E)를 계산하여 Ω₁ = g(E=1)를 추출한다.
• Ω₁의 지수적 스케일링 Ω₁ ∝ exp[rN]을 분석하고, K=2에서 6까지의 비율 상수 r(K)를 결정한다.
• 어드아바틱 양자 계산 프레임워크에서 양자 분할 함수 Z_Q를 사용하여 양자 갭-상관 길이 Δξ_GAP와 그 스케일링을 추정한다.

실험 결과

연구 질문

• RQ1 무작위 KSAT 문제에서 고유의 만족 가능한 할당(USA)이 존재하는 최소 절차-스핀 비율 α = M/N는 얼마이며, 이는 시스템 크기 N에 따라 어떻게 스케일링되는가?
• RQ2 고유의 만족 가능한 할당을 가진 KSAT 문제에서 첫 번째(excited) 에너지 수준(E=1)에서의 상태 밀도 Ω₁ = g(E=1)는 시스템 크기 N에 따라 어떻게 스케일링되는가?
• RQ3 K=2에서 6까지의 지수적 스케일링 Ω₁ ∝ exp[rN]에서 비율 상수 r(K)는 무엇이며, 이는 어떻게 비정제된 검색 한계 ln 2에 수렴하는가?
• RQ4 어려운 KSAT 엔세임블에서 양자 갭-상관 길이 Δξ_GAP는 시스템 크기 N에 따라 어떻게 스케일링되는가? 이는 양자 애너일링 실행 시간에 어떤 의미를 갖는가?
• RQ5 E=1에서 Ω₁를 최대화한 어려운 KSAT 문제 엔세임블은 양자 검색 알고리즘에 대해 강력한 기준 테스트로 활용될 수 있는가?

주요 결과

• 에너지 E=1에서의 상태 밀도 Ω₁ = g(E=1)는 시스템 크기 N에 따라 Ω₁ ∝ exp[rN] 형태로 지수적으로 발산하며, 비율 상수 r은 K=2일 때 0.348에서 K=6일 때 0.680으로 증가한다.
• K=2에서 6까지의 비율 상수 r(K)는 비정제된 검색 한계 ln 2 ≈ 0.693에 수렴하며, 이는 큰 K에 대해 거의 모든 스핀 구성이 E=1에 위치해 있고 유일한 지구 상태(E=0)를 제외한 나머지는 모두 E=1에 위치함을 시사한다.
• K=3일 때, 양자 갭-상관 길이 Δξ_GAP는 Δξ_GAP ∝ exp[r_Q N] 형태로 스케일링되며, r_Q ≈ 0.60(1)로 계산되며, 이는 고전적 비율 r(K=3)=0.569(1)과 유사하여 강한 양자 복잡도 스케일링을 시사한다.
• 최소 α ≈ 1에서 고유의 만족 가능한 할당을 가진 어려운 KSAT 문제 엔세임블은 랜덤 KSAT 엔세임블에서 지수적으로 희귀하므로, 양자 검색 알고리즘 테스트에 이상적이다.
• Ω₁의 지수적 특이성은 고전적 검색 시간이 τ_classical ∝ exp[rN]로 스케일링되며, 양자 애너일링 실행 시간이 Δξ_GAP의 제곱에 비례할 수 있음을 시사하여 기하급수적 저하를 초래할 수 있다.
• 높은 순서의 스핀 상호작용을 포함한 이 어려운 문제 엔세임블은, 두 스핀 상호작용 및 자외장 상호작용만을 포함하는 최대 독립 집합(MIS) 형태로 변환될 수 있으며, 이 과정에서 경직성이 유지된다.