Avoided crossings in spin-boson systems: consequences for adiabatic state preparation

기술 요약: 스핀-보존 시스템에서의 회피 교차(Avoided Crossings)

문제 정의
상호작로하는 양자 시스템에서 상관된 평형 상태(correlated equilibrium states)를 준비하는 것은 근본적인 과제이며, 특히 큐비트(quсbit)가 보존 환경에 결합된 형태를 설명하는 스핀-보존 모델(spin-boson model)에서 그러하다. 양자 컴퓨팅과 같은 실제 응용 분야에서 시스템은 종종 상호작용이 없는 상태로 초기화된 후 "퀜치(quench)" 또는 급격한 상호작용 스위치 온(switch-on) 과정을 거친다. 이러한 방식은 분리하기 어려운 들뜸(excitations)과 비물리적인 과도 현상(transients)을 생성한다. 진정한 상관된 바닥 상태를 준비하기 위한 표준적인 방법은 결합 강도를 점진적으로 증가시키는 단열 스위칭 프로토콜(adiabatic switching protocols)이지만, 그 효율성은 다체 스펙트럼(many-body spectrum)의 순간적인 특성에 의해 제한된다. 구체적으로, 에너지 스펙트럼의 "회피 교차(avoided crossings)" 영역은 작은 들뜸 간극(excitation gaps)을 형성하며, 이는 단열성을 유지하기 위해 매우 엄격한 스위칭 속도 조건을 부과한다. 유한한 시스템에서 이러한 스펙트럼 특성, 결합 강도, 그리고 필요한 스위칭 시간 사이의 정확한 관계는 여전히 중요한 연구 분야이다.

방법론
저자들은 하나 또는 세 개의 보존 모드를 포함하는 유한 스핀-보존 시스템을 조사한다. 본 연구는 정적 스펙트럼 분석과 동적 시간 진화 시뮬레이션을 결합한 이중 접근 방식을 채택한다:

1. 정확한 대각화(Exact Diagonalization): 저자들은 다체 해밀토니안을 수치적으로 대각화하여 정확한 바닥 상태와 전체 에너지 스펙트럼을 얻는다. 수치적 정확성을 보장하기 위해, 보존 점유수에 기초하여 힐베르트 공간을 절단(truncation)하며, 절단 크기는 $10^{-6}$의 허용 오차에 대해 검증된다. 절단 크기는 더 강한 큐비트-배스 드레싱(dressing)을 고려하여 유효 결합 강도($\gamma$)에 따라 단조롭게 증가시킨다.
2. 단열 스위칭 프로토콜: 시스템은 상호작용이 없는 해밀토니안($g_k=0$)의 바닥 상태에서 시작된다. 상호작용은 스위치 온 시간 $t_0$ 동안 사인 함수 형태의 램프 함수 $f(t) = \sin^2(\pi t / 2t_0)$를 사용하여 천천히 켜진다.
3. 단열 조건 분석: 단열성을 위한 이론적 요구 사항은 $\max_s R_{m0}(s) \ll t_0$ 조건을 사용하여 정량화되며, 여기서 $R_{m0}(s)$는 상호작용 해밀토니안의 행렬 요소와 바닥 상태와 들뜸 상태 사이의 에너지 간격의 역제곱에 의존한다.
4. 벤치마킹: 시간 진화된 상태는 대각화를 통해 얻은 정확한 바닥 상태와 비교된다. "스위치 온 시간"($t_0$)은 특정 관측량(바닥 상태 에너지 및 파울리 행렬 기댓값)이 지정된 허용 오차(에너지는 $10^{-6}$, 진동 진폭은 $10^{-4}$) 내에서 정확한 결과에 수렴하는 데 필요한 최소 시간으로 결정된다.

주요 기여 및 결과

• 스펙트럼 구조와 회피 교차: 본 연구는 저에너지 다체 스펙트럼이 복잡한 회피 교차 구조를 가짐을 확인한다. 이러한 교차의 위치와 간극 크기는 터널링 진폭($\epsilon$)과 결합 강도($\gamma$)에 따라 달라진다.
• 스위칭 시간에 대한 제약: 다체 에너지 준위 간격이 작은 영역(회피 교차와 관련된)은 단열적 준비에 가장 강력한 제약을 가한다.
  • $\epsilon = 1$ (터널링이 큐비트 에너지와 동일할 때), 요구되는 스위치 온 시간은 단조 증가하는 adiabatic ratio $R_{m0}(s)$와 일치하게 $\gamma$에 따라 단조 증가한다.
  • 더 작은 터널링 진폭($\epsilon = 0.5, 0.1$)의 경우, 스위치 온 시간은 비단조적인 거동을 보인다. 요구되는 시간의 국소적 감소는 스펙트럼의 특정 회피 교차 및 준위 번칭(level bunching)과 일치하며, 이는 단열 조건에 대한 지배적인 기여가 서로 다른 들뜸 상태 사이에서 이동함을 나타낸다.
  • $\epsilon = 0.1$인 경우, 강한 비단열 진동보다는 포화(saturation)에 도달하는 속도가 느리기 때문에 약 3배 더 긴 스위치 온 시간이 필요하다.
• 관측량 의존적 수렴: 전역적(global) 관측량과 국소적(local) 관측량 사이의 불일치는 중요한 발견이다. 바닥 상태 에너지는 전반적인 스펙트럼 추세를 따르며 $t_0$가 증가함에 따라 단조롭게 수렴하는 반면, 국소적 관측량(특히 파울리 행렬 기댓값 $\langle \sigma_z \rangle$ 및 $\langle \sigma_x \rangle$)은 비단조적인 스위칭 시간 거동을 보인다.
  • 단일 모드 케이스에서 $\epsilon = 1$일 때, 파울리 행렬을 위한 스위치 온 시간은 임계 결합($\gamma \approx 3.3$)에서 정점을 찍은 후 감소한다. 이 정점은 근접한 퇴화 상태(near-degenerate excited states) 사이의 코히어런트 진동이 지속되는 근공명(near-resonant) 영역에 대응한다.
  • 이러한 국소적 관측량들은 좁은 간극 근처에서 코히어런트 진동을 유지하며, 바닥 상태 에너지가 이미 수렴했음에도 불구하고 진동 진폭을 수렴 허용 오차 아래로 평균화하기 위해 더 긴 스위치 온 시간을 요구한다.

의의 및 주장
본 논문은 열린 양자 시스템의 이론적 처리와 양자 시뮬레이션 및 큐비트 제어의 실제 구현을 위한 구체적인 벤치마크를 제공한다고 주장한다. 유한한 환경 내에서의 고유한 양자 역학을 분리함으로써, 본 연구는 연속체 배스(continuum baths)에 사용되는 근사법 없이 상관된 평형 상태가 어떻게 출현하는지를 명확히 한다.

저자들은 단열적 상태 준비의 효율성이 단순히 전역 에너지 간격에 의해 결정되는 것이 아니라, 측정되는 특정 관측량에 매우 민다는 점을 강조한다. 유한한 크기 효과와 회피 교차에 의해 유지되는 코히어런트 진동으로 인해, 국소적 양은 전역적인 에너지 척사와 상당히 다른 시간 척도에서 평형에 도달할 수 있다. 따라서 큐비트 제어 및 양자 정보 처리 프로토콜은 과도 현상을 피하고 시스템이 의도된 상관된 평형 상태에 도달하도록 하기 위해 이러한 국소적 동적 특성을 반드시 고려해야 한다. 이러한 결과는 상관된 평형 상태로부터 시작하는 것이 외부 구동(external driving)에 대한 시스템의 반응을 조사하는 통제된 경로를 제공함을 시사한다.