Kibble--Zurek Scaling in the Dicke Model at Mesoscopic Scales

기술 요약: 메조스코픽 규모에서의 디케 모델 내 키블-주렉(Kibble–Zurek) 스케일링

문제 정의
디케(Dicke) 모델은 집단적 빛-물질 물리학 및 초방사 상전이를 위한 전형적인 프레임워크 역할을 한다. 근본적인 해밀토니안은 단순하지만, 실험적으로 접근 가능한 메조스코픽 영역(유한한 입자 수 $N$)에서 임계 지수를 정량적으로 추출하는 것은 상당한 도전 과제이다. 두 가지 주요 장애물이 점근적 스케일링 거동을 가로막는다:

1. 상관 시간의 느린 발산: 모든 입자가 결합된(all-to-all coupled) 시스템에서는 임계 상관 시간이 $N$에 따라 느리게 성장하여, 유한 크기 보정(finite-size corrections)이 지배하게 되고 동역학적 지수(dynamical exponent) 추정치를 편향시킨다.
2. 소산적 교차(Dissipative Crossover): 피할 수 없는 광자 손실은 유의미한 섭동으로 작용하여 시스템을 별도의 소산적 보편성 클래스로 유도한다. 메조스코픽 규모에서 코히어런트(닫힌 계)와 소산적(열린 계) 극한 사이의 교차는 폭이 넓으며, 이는 겉보기 지수가 두 점근적 고정점 중 어느 곳과도 일치하지 않게 만든다.
   또한, 기존의 모든-입자 결합 모델에 대한 동역학적 스케일링 법칙은 점근적인 대규모 $N$ 극한을 위해 공식화되었으므로, 강한 유한 크기 드리프트와 얽혀 있는 교차 효과(크기, 소산, 램프 속도)가 복잡하게 작용하는 유한 시스템에서 제어된 추출 방법을 위한 간극을 남겨두고 있다.

방법론
저자들은 닫힌 디케 모델과 열린 디케 모델 모두에 대해 정적, 스펙트럼 및 구동(램핑) 동역학을 동일한 층위에서 다루는 통합된 메조스코픽 스케일링 프레임워크를 개발하였다.

1. 대규모 $N$ 장론(Large-$N$ Field Theory): 저자들은 임계 지수를 분석하기 위해 통합된 켈디시(Keldysh) 장론을 활용한다. 선도적인 $1/N$ 보정을 명시적으로 유지함으로써, 다음과 같은 뚜렷한 고정점을 식별한다:

   • 닫힌 경우 ($\kappa=0$): 시스템은 모든-입자 결합 이징(Ising) 보편성 클래스에 속하며, 지수는 $z=1/3$ 및 $\nu=3/2$이다.
   • 열린 경우 ($\kappa \neq 0$): 시스템은 Model-A 보편성 클래스로 흐르며, 지수는 $z=1/2$ 및 $\nu=2$이다.
   • 이 이론은 키블-주렉(KZ) 스케일링 관계식 $1/\mu = 1/\nu + z$를 확립하며, 여기서 $\mu$는 램핑 속도를 특징짓는다.

2. 통합 스케일링 프로토콜: 메조스코픽 데이터로부터 지수를 추출하기 위해, 저자들은 선도적인 무관한 보정(irrelevant correction, 지수 $\omega$)을 명시적으로 고려하는 2단계 프로토콜을 제안한다:

   • 1단계 ( $\omega$ 고정): 무차원 관측량(예: 바인더 비 $U$)을 임계점에서 분석한다. 유한-$N$ 드리프트는 $N^{-\omega}$의 거듭제곱에 대한 절단된 급수로 피팅되어 무관한 지수를 결정한다.
   • 2단계 ($z, \nu, \mu$ 추출): 고정된 $\omega$를 사용하여, 다음으로부터 나머지 지수들을 독립적으로 추출한다:
     • $z$: 임계 갭(열린 시스템의 경우 리우빌리안 갭)의 스케일링.
     • $\nu$: 튜닝 파라미터에 대한 무차원 관측량 미분값의 크기 의존성.
     • $\mu$: 램프 동역학에서의 주요 속도 보정(예: $\partial_s U$ 또는 $\partial^2_s U$)의 크기 의존성.

주요 결과
프레임워크를 검증하기 위해 수치 시뮬레이션과 대규모 $N$ 전개가 수행되었다:

• 닫힌 디케 모델:

  • 무관한 지수는 $\omega \approx 0.372$로 결정되었다.
  • 추출된 지수들은 이론적인 대규모 $N$ 값을 회복했다: $z \approx 0.33$ (이론값 $1/3$), $\nu \approx 1.52$ (이론값 $3/2$), $\mu \approx 0.99$ (이론값 $1$).
  • KZ 관계식 $1/\mu = 1/\nu + z$가 검증되었다.
  • 램프 동역학은 시간 역전 대칭성에 부합하는 이차 속도 의존성(선형 보정의 부재)을 보였다.

• 열린 디케 모델:

  • $\kappa \in [0.1, 3]$ 범위의 소산율에 대한 공동 피팅을 통해 $\omega \approx 0.363$을 얻었다.
  • $\kappa=1$인 경우, 정확한 대각화(exact diagonalization)는 $z \approx 0.48$을 나타냈다. 그러나 일반적인 $\kappa$에 대해, 정확한 대각화가 접근 가능한 작은 크기($N \le 24$)는 모드 교차(mode crossings)와 재배열을 보여 점근적 스케일링을 가렸다.
  • 이를 해결하기 위해 $N \in [10^5, 10^7]$에 대해 수치적 대규모 $N$ 전개를 수행하였으며, 이를 통해 추출된 지수가 서로 다른 $\kappa$에 대해 이론적 값 $z=1/2$에 접근함을 확인하였다.
  • 지수 $\nu$는 $\nu \approx 2.02$ ($\kappa=1$일 때)로 추출되었으며, 서로 다른 $\kappa$에 대해 기대값 주변에서 변동하였다.
  • 지수 $\mu$는 $\kappa=1$일 때 $\approx 0.96$, 서로 다른 $\kappa$에 대해 $\approx 1.02(5)$로 나타났다.
  • 열린 경우에도 KZ 관계식이 성립하였다.
  • 램프 동역학은 시간 역전 대칭성 파괴로 인한 선형 속도 보정을 보였다.

• 교차 동역학:

  • 본 연구는 유한 크기($N$), 소산($\kappa$), 램프 속도($s$) 사이의 경쟁을 명확히 하였다.
  • 교차 스케일링 변수 $\kappa N^{1/3}$이 식별되었으며, 이는 작은 규모/약한 소산 시스템에서의 "닫힌 계와 유사한" 거동과 "열린 계와 유사한" 거동을 구분한다.
  • KZ 시간 척도($t_{KZ} \sim s^{-1/3}$)와 소산 시간 척도($t_\kappa \sim \kappa^{-1}$) 사이의 상호작용은 왜 가공되지 않은 램프 데이터가 복잡해 보일 수 있는지를 설명하며, 이는 이러한 효과들을 분리하기 위해 통합된 스케일링 처리가 필요함을 시사한다.

의의 및 주장
본 논문은 닫힌 계와 열린 계 모두에서 장거리 상호작용을 갖는 메조스코픽 시스템의 정적 및 동적 임계 스케일링을 해결하는 통합 프레임워크를 구축했다고 주장한다. 주요 기여는 다음과 같다:

1. 메조스코픽 모호성 해결: 선도적인 무관한 보정을 포함함으로써, 표준적인 유한 크기 스케일링이 실패하거나 편향된 추정치를 낼 수 있는 메조스코픽 데이터로부터 올바른 열역학적 지수를 회복하는 프로토콜을 제공한다.
2. KZ 스케일링 검증: 공간적 도메인 형성이 없는 완전 연결 모델에서 키블-주렉 스케일링($1/\mu = 1/\nu + z$)을 처음으로 제어된 방식으로 검증하였으며, 이를 코히어런트 및 소산 영역 모두에서 확인하였다.
3. 일반적 적용 가능성: 저자들은 이 프레임워크가 캐비티 QED에 국한되지 않고, 트랩된 이온(trapped ions)이나 Sachdev–Ye–Kitaev 유형의 모델과 같이 장거리 상호작용을 갖는 다른 메조스코픽 시스템에도 널리 적용될 수 있다고 밝히며, 특히 제어된 대규모 $N$ 이론을 사용할 수 없고 직접적인 데이터 추출이 필수적인 경우에 유용하다.

이 연구는 열역학적 극한이 원칙적으로는 명확하지만, 실험적 접근성은 종종 메조스코픽 규모로 제한된다는 점을 강조하며, 이 프로토콜이 해당 영역에서 보편성 클래스를 탐구할 수 있는 실질적인 경로를 제공함을 보여준다.