A near-term quantum simulation of the transverse field Ising model hints at Glassy Dynamics

기술 요약: 가스형 역학(Glassy Dynamics)을 암시하는 횡장 이싱 모델의 근미래 양자 시뮬레이션

문제 정의
본 논문은 양자 다체 물리(quantum many-body physics)에서 나타나는 이색적인 물질 상태와 가스형 전이(glassy transitions)를 이해하는 데 필수적인 대규모 스핀 글래스(spin glass) 시스템을 시뮬레이션하는 데 따르는 계산적 과제를 다룹다. 분자 역학 및 몬테카를로 시뮬레이션과 같은 고전적 계산 접근 방식은 다전자 시스템에서의 전자 파동 함수의 지수적 성장과 무질서한 시스템을 장시간 동안 모델링하는 데 필요한 높은 계산 복잡도로 인해 어려움을 겪는다. 특히, 가스형 시스템은 수많은 국소 최솟값(local minima)과 매우 긴 이완 시간(relaxation times)을 가진 험준한 에너지 경관(rugged energy landscapes)을 나타내므로, 고전적 고성능 컴퓨팅(H용 HPC)을 사용하여 정확하게 시뮬레이션하기 어렵다. 저자들은 디지털 양자 시뮬레이션(DQS) 기술, 구체적으로 변분 양자 고유치 솔버(Variational Quantum Eigensolvers, VQE)가 근미래 양자 하드웨어에서 이러한 무질서한 스핀 구성을 효과적으로 탐사할 수 있는지 조사하고자 한다.

방법론
저자들은 종방향 자기장(longitudinal field)이 추가된 횡장 이싱 모델(Transverse Field Ising Model, TFIM)을 시뮬레이션하기 위해 VQE 알고리즘을 사용하는 하이브리드 양자-고전 접근 방식(hybrid quantum-classical approach)을 채택하였다. 연구는 다음 두 가지 특정 기하 구조에 집중하였다:

1. 1D 이싱 체인: 길이 $n=25$ (25 큐비트)인 체인.
2. 2D 이싱 격자: $6 \times 6$ 정사각형 격자 (36 큐비트).

주요 기술 구성 요소:

• 해밀토니안 구성: 시스템 해밀토니안은 z-축을 따른 인접 이웃 상호작용($J \sum \sigma_i^z \sigma_{i+1}^z$), x-방향의 횡장($-B_x \sum \sigma_i^x$), 그리고 z-방향의 종방향 자기장($-B_z \sum \sigma_i^z$)을 포함한다. 대각 상호작용(Diagonal interactions)은 제외되었다.
• 안사츠(Ansatz) 선택: 저자들은 $R_x$ 및 $R_y$ 회전 블록으로 구성된 "TwoLocal" 안사츠를 사용하였다. 결정적으로, 이 안사츠는 큐비트 간의 얽힘(entanglement) 없이 설계되었다. 이러한 선택은 회로 깊이를 줄이고 근미래 장치의 노이즈 민감도를 완화하는 동시에, 문제 구조에 필요한 양자 요동(quantum fluctuations)을 포착하기 위함이었다.
• 최적화: 해밀토니안의 기댓값을 최소화하기 위해 고전적 최적화 도구인 COBYLA가 사용되었으며, 최대 25,000회의 반복(대부분의 실행은 더 일찍 수렴함)이 수행되었다.
• 소프트웨어 구현: 저자들은 1D 및 2D 격자에 대한 해밀토니안을 모듈식으로 구성하기 위해 커스텀 파이썬 패키지인 Qiskit Glassy Dynamics(이후 QuGlassyIsing으로 지칭됨)를 개발하였다. 이 패키지는 상호작용 항($\sigma^z\sigma^z$)과 필드 항($\sigma^x, \sigma^z$)의 생성을 효율적으로 처리한다.
• 시뮬레이션 환경: 시뮬레이션은 주로 Qiskit Aer GPU 백엔드를 사용하여 수행되었다. 개념 증명 시연은 축소된 7-스핀 모델을 사용하여 실제 7-큐비트 IBM Q Oslo 장치에서 수행되었다.

주요 결과

• 에너지 경관: 시뮬레이션 결과, 종방향 자기장($B_z$)을 증가시키면 스핀을 정렬시켜 시스템 에너지를 낮추는 반면, 횡장($B_x$)을 증가시키면 에너지를 크게 낮추지 않으면서도 양자 요동과 스핀 플립을 유도한다는 것을 밝혀냈다.
• 무질서 상(Disordered Phases): 연구는 무질서한 스핀 구성이 나타나는 특정 영역을 식별하였다. 이 "무질서 상"은 강자성(ferromagnetic)과 반강자성(anti-ferromagnetic) 도메인의 혼합으로 특징지어진다.
  • 종방향 자기장이 없는 경우($B_z=0$), 시스템은 완벽한 반강자성 질서를 보였다.
  • 무질서 상태는 특정 $B_x$와 $B_z$의 조합(예: $B_z=1, B_x=2$)에서만 관찰되었는데, 이는 두 자기장 사이의 경쟁이 균일한 질서를 방해하여 반강자성 배경 내에 국소적인 강자성 모멘트를 생성하기 때문이다.
  • 저자들은 2D 시뮬레이션에서 가장 확률이 높은 세 가지 스핀 구성이 거의 동일하다는 점에 주목하였으며, 이는 무질서 상 내에서도 일관된 패턴화된 질서가 존재함을 시사한다.
• 하드웨어 시연: IBM Q Oslo 장치(7 큐비트)에서의 실행은 상당한 변동을 동반한 에너지 수렴 프로파일을 보여주었으며, 최종 에너지 값($-4.0139$)은 이론적 최솟값보다 높았다. 이 결과는 NISQ(노이즈가 있는 중간 규모 양자) 장치 고유의 게이트 오류, 판독 오류 및 큐비트 결맞음(decoherence)에 기인한다.

의의 및 주장
본 논문은 근미래 양자 시뮬레이터가 고전적 방법으로는 접근하기 어려운 이싱 모델의 무질서 상을 탐지할 수 있다는 "개념 증명"을 제공한다고 주장한다. 저자들은 다음과 같은 점을 강조한다:

• VQE의 효용성 검증: 얽히지 않은 안사츠를 사용하더라도 VQE가 복잡한 파동 함수를 다항 시간 내에 모델링할 수 있음을 입증하여, 고전적 방법과 비교하여 계산 비용을 절감할 수 있음을 보여준다.
• 가스형 역학 규명: 결과는 횡장과 종방향 자기장의 상호작용이 험준한 에너지 경관을 만들어내며, 이는 무질서가 비효율적인 반응 메커니즘을 초래하는 촉매 과정이나 스핀 글래스의 거동을 반영한다는 것을 시사한다.
• 도구 개발: QuGlassyIsing 패키지의 도입은 가스형 역학을 위한 해밀토니안 구성을 위한 모듈식의 사용자 친화적인 인터페이스를 제공하여, 복잡한 양자 시스템에 대한 향후 연구를 용이하게 한다.
• 미래 전망: 저자들은 현재의 NISQ 한계(결맞음 시간, 게이트 충실도)가 더 크고 복잡한 시스템으로의 확장을 제한하고 있지만, 본 연구가 미래의 결함 허용 양자 컴퓨팅(FTQC) 응용을 위한 토대를 마련했다고 겸허히 밝힌다. 그들은 양자 오류 수정이 개선됨에 따라, 이러한 방법론이 3D 시스템 및 더 복잡한 기하 구조의 가스형 거동 시뮬레이션을 가능하게 하여 새로운 자기 저장 물질 및 촉매 시스템 설계에 정보를 제공할 수 있다고 제안한다.

결론적으로, 본 논문은 현재의 구현이 하드웨어 제약으로 인해 제한적이지만, 가스형 시스템의 역학과 구조에 대한 미묘한 이해를 가져오기 위한 양자 정보 처리 도구의 잠재력을 성공적으로 부각시켰다고 밝히고 있다.