Restoring thermalization in long-range quantum magnets with staggered magnetic fields

문제 정의
상호작용 강도가 $1/r^\alpha$로 감소하며 $\alpha < d$인 강한 장거리(long-range, LR) 상호작용을 가진 양자계는 이론적으로 열화(thermalization)에 저항할 것으로 예상된다. 이러한 저항은 두 가지 주요 메커니즘에서 기인한다: 열역학적 극한에서 에너지 스펙트럼의 이산성을 생성하여 시스템 크기에 따라 스케일링되는 수명(lifetime)을 가진 준안정 상태를 보호하는 것, 그리고 치환 대칭성(permutational symmetry)의 지속이 고유 상태에 구조를 유도하고 고유 상태 열화 가설(ETH)을 위반하는 것이다. 기존 연구들은 상호작용 범위 $\alpha$가 열화에 미치는 영향에 집중해 왔으나, 상호작용 범위와 독립적으로 해밀토니언 자체의 대칭성을 감소시킴으로써 열화를 복구할 수 있는지 여부는 여전히 미해결 과제로 남아 있다.

방법론
저자들은 교차 자기장(staggered magnetic field)이 가해진 1차원 장거리 하이젠베르크 반강자성체를 조사한다. 이 시스템은 멱함수(power-law) 상호작용을 통해 결합된 두 개의 부격자(sublattice) $A$와 $B$로 구성된다. 해밀토니안은 반강자성 교환 항과 부격자 간의 대칭성을 깨뜨리고 하이젠베르크 상호작용의 회전 대칭성을 깨뜨리는 교차 자기장 항을 포함한다.

본 연구는 다음과 같은 이중 접근 방식을 채택한다:

1. 자기 일관적 평균장 이론(Self-consistent Mean-Field Theory): 저자들은 슈윙거 보존(Schwinger boson) 사상과 조던-위그너(Jordan-Wigner) 변환을 사용하여 스핀 시스템을 페르미온으로 사상한다. LR 시스템에서의 국소적 요동 억제를 근거로 스트링 연산자(string operators)에 대한 평균장 디커플링을 적용하여, 유효 이차 해밀토니안을 유도한다. 이를 통해 스펙트럼의 해석적 대각화와 집합적 관측량의 계산이 가능하다.
2. 정확한 대각화(Exact Diagonalization, ED): 분석적 예측을 검증하기 위해 유한한 시스템에 대한 수치 시뮬레이션을 수행한다. 해밀토니안 (3)은 최대 수백 개의 스핀(예: 그림 1의 $\alpha=0$인 경우 $N=500$)에 대해 대각화되며, 유한한 $\alpha$ 역학에 대한 유한 크기 스케일링은 더 작은 시스템 크기(예: 보충 자료 그림 2의 $N=18$)를 활용한다. ED는 특정 자화 섹터($S_z$) 내에서 수행되며, 힐베르트 공간의 차원을 줄이기 위해 병진 불변성과 같은 대칭성을 활용한다.

역학은 가우시안 분포를 가진 방향을 가진 순수 넬(Néel) 상태들의 통계적 앙상블인 혼합 넬 상태를 통해 탐사된다. 이를 통해 시스템 크기와는 구별되는 초기 상태의 요동($\sigma$)을 독립적으로 조절할 수 있다.

주요 기여 및 결과

• 대칭성 파괴를 통한 열화의 복구: 교차 자기장을 가하는 것이 시스템의 완전한 치환 대칭성을 깨뜨려 이분법적(bipartite) 대칭성으로 감소시킨다는 것이 핵심 발견이다. 이러한 대칭성 파기은 에너지 스펙트럼의 퇴화(degeneracy)를 해소한다. 스펙트럼은 여전히 이산적인 부분 공간들로 구성되지만, 이 부분 공간의 수는 $O(N^2)$만큼 스케일링된다(치환 불변인 LR 시스템의 $O(N)$과 비교하여). 결과적으로, 집합적 스펙트럼은 열역학적 극한($N \to \infty$)에서 조밀해지며, 평균 에너지 간격은 $O(1/N)$으로 스케일링된다. 이 조밀한 스펙트럼은 유한한 요동을 가진 초기 상태에 대해 열화를 가능하게 한다.

• 3단계 평형 역학: 넬 파라미터 $\langle n_z \rangle$의 시간 진화는 세 가지 뚜렷한 역학적 단계를 보인다:

  1. 단기(Short-time): 주파수 $\omega_p = \sqrt{Jh}$를 갖는 고전적 진자와 유사한 진동적 평균장 역학을 보인다.
  2. 중기(Intermediate-time): 진동의 주기적 회복(revival)과 감쇠가 나타난다. 감쇠 시간 $T_{dec}$은 넬 파라미터의 초기 상태 요동에 반비례($T_{dec} \propto 1/\sigma_{n_z}$)하지만, 시스템 크기 $N$에는 무관하다.
  3. 장기(Long-time): 평균적으로 마이크로캐노니컬 앙상블로 평형을 이룬다. 평형 시간 $T_{eq}$는 균일 자화의 요동에 따라 $1/\sigma_{sz}^2$로 스케일링된다. 결정적으로, 유한한 요동을 가진 상태의 경우, 평형 시간은 시스템 크기에 관계없이 유한한 시간 내에 발생하는데, 이는 평형 시간이 $N$에 따라 발산하는 치환 불변 LR 시스템의 메타스테이빌리티(metastability)와 대조된다.

• 세퍼라트릭스(Separatrix)에서의 에르고드성 파괴: 저에너지 및 중에너지 상태는 열화되지만(장기 평균이 마이크로캐노니컬 앙상블과 일치함), 세퍼트릭스 에너지(결합된 궤적과 자유로운 궤적 사이의 경계) 근처의 상태들은 에르고드성 파괴를 보인다. 이 영역에서는 장기 평균이 초기 상태에 의존한다. 이러한 거동은 불안정한 고정점(수직 진자)에 국소화된 "스카(scar)-와 유사한" 고유 상태에 기인한다. 이 상태들은 낮은 얽힘 엔트로피와 특정 초기 상태와의 비정상적으로 큰 중첩을 보이며, 열화를 방지한다.

• 유한한 $\alpha$ 효과: 유한한 상호작용 범위($0 < \alpha < 1$)에 대해, 부격자 $A$와 $B$의 집합적 스핀 길이($S_A^2, S_B^2$)는 더 이상 보존되지 않으며 시간에 따라 감쇠한다. 그러나 이 집합적 스핀의 수명은 $N^{3/2}$로 스케일링되며, 열역학적 극한에서 발산한다. 따라서 열역학적 극한에서 $0 < \alpha < 1$의 역학은 넬 파라미터는 평형을 이루는 동안 집합적 스핀 길이는 동결되는 완전 연결 극한($\alpha = 0$)의 역학과 실질적으로 일치한다.

의의
본 논문은 교차 자기장을 도입함으로써 강한 장거리 양자 자석에서 광범위한 초기 상태에 대해 열화를 복구할 수 있음을 입증한다. 이 메커니즘은 치환 대칭성을 깨뜨려 조밀한 에너지 스펙트럼을 생성함으로써, LR 시스템의 전형적인 이산 스펙트럼과 메타스테이빌리티를 우회한다. 결과는 스펙트럼이 조밀하다면 평형 시간이 시스템 크기가 아닌 초기 상태의 요동에 의존함을 보여준다. 또한, 본 연구는 스펙트럼 중앙에서 불안정한 고전적 지점에 국소화된 스카 유사 상태를 통한 에르고드성 파괴의 특정 메커니즘을 식별하였다. 저자들은 이 제안이 국소적인 교차 자기장과 장거리 상호작용만을 필요로 하여, 트랩된 이온, 리드베리 원자 배열, 공동 결합 원자 등의 플랫폼에서 실험적으로 구현 가능하다는 점을 언급하였다.