BCₙ형 폴리크로노코스‑프라흠 스핀 체인의 통계와 베리‑타보르 예측 위반

본 논문은 su(m) 대칭을 갖는 Polychronakos‑Frahm(PF) 스핀 체인의 BCₙ형 버전을 대상으로, 그 파티션 함수와 스펙트럼 통계적 특성을 심도 있게 분석한다. 연구는 크게 네 부분으로 구성된다. 첫 번째 부분에서는 “freezing trick”이라 불리는 방법을 도입한다. 이는 Calogero‑Sutherland 모델의 강 결합 한계에서 입자 좌표가 잠재적 최소점 ξ_i 로 고정되고, 스핀 자유도만 남는 과정을 의미한다. Hamiltonian (1)은 두 종류의 장거리 상호작용 1/(ξ_i−ξ_j)²와 1/(ξ_i+ξ_j)², 그리고 경계 항 β·S_i·ξ_i⁻² 로 이루어져 있다. 여기서 ε=±1 은 스핀 반전 대칭을 나타낸다. 최소점 ξ_i 는 일반화 라게르 다항식 L_{N}^{β−1} 의 영점이며, N≫β 일 때 그 분포는 원형법칙 ρ_N(x)=1/(2π)√{8N−x²} 로 근사된다. 이는 Aₙ형 PF 체인과 유사한 밀도이며, 이후 레벨 밀도 분석에 핵심적인 역할을 한다. 두 번째 부분에서는 스칼라와 스핀 Hamiltonian을 동일한 형태의 연산자 H′ 로 변환한다. Dunkl 연산자 J_i 를 이용해 H′ 를 상삼각 행렬 형태로 표현하고, 고유값은 E′_n = a|n|+E₀ (a는 결합 상수, |n|는 다항식 차수) 로 단순히 차수에 비례한다. 스핀 연산자 S_{ij}, S_i 를 좌표 교환 연산자 K_{ij}, K_i 로 치환함으로써 스핀 체인의 스펙트럼을 정확히 얻는다. 세 번째 부분에서는 파티션 함수를 구한다. 스칼라 파티션 Z_sc 은 차이 p_i 로 전개해 ∏_{i=1}^N (1−q^{2i})^{-1} 형태가 된다. 스핀 파티션 Z_ε 은 스핀 자유도에 대한 제약(조건 i‑iii)으로 인해 다중 파티션 ν∈P_N 와 그에 대응하는 중복도 d(ν)를 도입한다. 짝수 m 경우 d(ν)=∏_{i=1}^r C(m/2, ν_i) 로 간단히 표현되고, 최종 파티션 함수는 Z(T)=q^{-N}∏_{i=1}^N(1−q^{2i}) ∑_{ν∈P_N} d(ν) q^{∑_{j=1}^{ℓ(ν)}N_j}/∏_{j=1}^{ℓ(ν)}(1−q^{N_j}) 형태를 갖는다. 여기서 N_j=j∑_{i=1}^{ν_i}이며 ℓ(ν)는 파티션의 파트 수이다. 특히 m=2(스핀 ½) 에서는 Z(T)=q^{N(N−1)/2}∏_{i=1}^N(1+q^i) 로 단순화되고, 에너지 스펙트럼은 E_j= N(N−1)/2 + j (j=0,…,N(N+1)/2) 로 등간격을 이룬다. 각 레벨의 퇴화도는 정수 j 를 N 이하의 서로 다른 부분으로 분할하는 방법 수 Q_N(j) 로 주어지며, 이는 Ramanujan의 5차 모크 시타 함수와도 연결된다. 홀수 m 경우에는 짝·홀 차수를 별도 정렬하고, 스핀 중복도 d_ε^s(ν) 를 두 개의 조합식으로 나눠 계산한다. 최종 파티션 함수는 (29)‑(30) 식에 따라 추가 인자 q^{−(N+N_s)} 가 등장한다. 네 번째 부분에서는 파티션 함수를 이용해 스펙트럼 통계량을 조사한다. 레벨 밀도 ρ(E)를 수치적으로 계산하면, N≥30 정도에서 평균 μ와 분산 σ² 를 갖는 가우시안 ρ_G(E)와 거의 일치한다. 이는 Haldane‑Shastry(HS) 체인에서도 관찰된 동일한 현상이며, 장거리 상호작용 적분 가능한 모델의 일반적 특성으로 해석된다. 다음으로 인접 레벨 간격 s_i=(E_{i+1}−E_i)/ΔE(ΔE는 평균 간격)를 정규화하고 누적 분포 P(s)=Prob(S≤s)를 구한다. 수치 결과는 포아송(P_P(s)=1−e^{−s})도, Wigner(P_W(s)=1−e^{−πs²/4})도 맞지 않으며, “거의 선형 구간 + 급격한 급락” 형태의 독특한 분포를 보인다. 이를 설명하기 위해 저자들은 (i) 레벨 밀도가 가우시안, (ii) 에너지 스펙트럼이 등간격, (iii) 스펙트럼이 대칭적이라는 세 가지 사실을 활용한다. 이들로부터 누적 간격 분포를 P(s)≈1−½ erfc