On Trace Averages and Sharp Asymptotics for Powers of Haar-Unitary Matrices

기술 요약: 하르-유니터리 행렬의 거듭제곱에 대한 트레이스 평균 및 정밀 점근법에 관하여

문제 정의
본 논문은 하르 분포를 따르는 유니터리 행렬 $U \in U(N)$의 거듭제곱의 통계적 거동을 조사한다. 구체적으로, $N \to \infty$일 때 발생하는 두 가지 뚜렷한 수렴 체제(regime)를 다룬다:

1. 전역 통계(Global Statistics): 정규화된 트레이스 $X_{N,m} = \text{Tr}(U^m)/\sqrt{m \wedge N}$의 분포와 표준 복소 가우시안 분포 사이의 총 변동 거리($d_{TV}$).
2. 국소 통계(Local Statistics): 고정된 구간 위에서 $U^m$의 언폴딩된 고윳값 스펙트럼의 분포, 특히 단위율 포아송 과정(unit-rate Poisson process)과의 거리.

이 논문의 핵심적인 동기는 Rains의 거듭제곱 맵 정리(power-map theorem)이다. 이 정리는 $U^m$의 스펙트럼이 독립적인 저차원 순환 유니터리 앙상블(CUE) 스펙트럼들의 중첩으로 분해됨을 보여준다. 본 논문은 이러한 스펙트럼 분해가 어떻게 트레이스 모멘트와 점 프로세스(point processes)에 대한 정밀한 점근 법칙으로 전이되는지를 정량화하고자 하며, 특히 거듭제곱 $m$이 $N$과 함께 스케일링되는 경우($m$이 고정되거나, 성장하거나, $N$을 초과하는 경우)를 다룬다.

방법론
저자는 표현론, 해석적 조합론, 그리고 랜덤 행렬 이론 기법의 합성을 활용한다:

• 정확한 트레이스 항등식: 논문은 프로베니우스의 캐릭터 공식(Frobenius's character formula)과 슈어 직교성(Schur orthogonality)을 사용하여 혼합 트레이스 모멘트를 유한 캐릭터 합으로 표현한다. "안정 범위(stable range, 모멘트가 가우시안 모멘트와 일치하는 구간)"를 벗어난 유한-$N$ 보정을 처리하기 위해, 저자는 **찰리-이웬스 합성(Charlier–Ewens convolution)**을 도입한다. 이는 안정적인 위크 항(Wick term)을 "보정 밴드(correction bands)"로 인덱싱된 유한 보정 급수와 분리한다. 또한, 혼합 모멘트를 치환 행렬의 행렬식의 유한 합으로 표현하는 변위 행렬식(displacement determinant) 공식을 유도하여, 캐릭터 합의 복잡성을 우회한다.
• 해석적 도구:
  • 베셀 및 프레드홀름 표현: 트레이스 모멘트의 생성 함수는 베셀 행렬식과 동일함이 확인되었다. 이들은 보로딘-오쿤코프-게로니모-케이스(Borodin–Okounkov–Geronimo–Case) 항등식을 사용하여 헹켈 연산자(Hankel operators)를 포함하는 프레드홀름 행렬식으로 더욱 분해된다.
  • 푸리에-라게르 전이(Fourier–Laguerre Transfer): 특성 함수로부터 확률 밀도로 이동하기 위해, 밀도 보정의 기저로서 라게르 다항식이 역할을 하는 전이 메커니즘을 사용한다. 이를 통해 $L^1$ 노름 계산에 필요한 정밀한 계산이 가능해지며, 이는 총 변동 경계(total variation bounds)를 산출하는 데 쓰인다.
  • 포아송-찰리(Poisson–Charlier) 전개: 독립적인 CUE 블록들의 중첩은 희소한 베르누이 과정(sparse Bernoulli processes)의 합으로 근사된다. 오차는 포아송-찰리 전개와 점 프로세스 간의 총 변동 거리 상한을 사용하여 제어된다.

주요 기여 및 결과

1. 정확한 혼합 트레이스 모멘트 및 유한-$N$ 보정
본 논문은 임의의 지수에 대해 정확한 혼합 트레이스 평균 $C_N(\mathbf{p}, \mathbf{q})$의 평가를 제공한다.

• 찰리-이웬스 보정: 정리 6은 가우시안 위크 항으로부터의 편차가 유한한 "보정 밴드"의 합임을 입증한다. 밴드의 개수는 전체 차수 $D$의 분할(partition)에 대한 기하학적 구조에 의해 결정된다.
• 변위 행렬식: 정리 7은 임의의 $N$에 대해 유효한, 치환 행렬의 행렬식의 합으로서 이러한 모멘트들을 계산하는 효율적인 공식을 제공한다.

2. 정밀한 전역 가우시안 점근법
본 논문은 정규화된 변수 $X_{N,m}$과 표준 복소 가우시안 $Z$ 사이의 총 변동 거리 $d_{TV}(X_{N,m}, Z)$에 대한 날카로운 주 차수 점근법을 결정한다. 결과는 $m$과 $N$의 관계에 의해 정의된 네 가지 체제에 대해 균일하다:

• 고정 거듭제곱 ($m=p$): 정리 1은 $e^{2/p}$를 포함하는 일관된 합산 인자를 가진 정확한 리딩 항을 제공한다. 이 인자는 첫 번째 보정 밴드에서 라게르 모드들의 일관된 합산으로부터 발생하며, 이는 단지 첫 번째 비가우시안 모멘트만을 분석해서는 놓칠 수 있는 세부 사항이다.
• 결합 성장 ($m \to \infty, N/m \to \infty$): 정리 3은 리딩 항이 고정 거듭제곱 케이스와 동일한 형태를 가짐을 보여주며, 여기서 $p \to \infty$임에 따라 $e^{2/p}$ 인자가 1로 수렴한다.
• 고정 블록 차원 ($N/m \to q$): 정리 4는 블록 차원이 유계(bounded)로 유지되는 체제를 분석한다. 리딩 보정은 두 개의 인접한 라게르 모드($L_{q+1}$ 및 $L_{q+2}$)의 선형 결합을 포함한다.
• 독립 단계 체제 ($m \ge N$): 정리 2는 $m \ge N$일 때 고윳값들이 독립적인 균등 점들이 됨을 증명한다. 본 논문은 이전의 경계보다 정밀하게 $O(N^{-1})$ 수렴 속도에 대한 명시적인 상수를 유도한다.

3. 정밀한 국소 포아송 근사
정리 5는 길이 $L$인 고정 구간 $B$ 위에서의 $U^m$의 국소 스펙트럼 통계를 다룬다.

• 포아송 극한: $N, m \to \infty$일 때, 언폴딩된 스펙트럼은 포아송 과정으로 수렴한다.
• 리딩 오차 항: 본 논문은 고윳값 개수와 포아송 변수 사이의 총 변동 거리에 대한 정확한 리딩 항을 식별한다. 오차는 블록 강도(block intensities)의 제곱 합인 $\beta_{N,m}$에 비례한다.
• 과소 분산(Under-dispersion): 두 번째 찰리 다항식을 통해 유도된 1차 보정은 국소 고윳값 개수의 "잔여 과소 분산"을 정량화하며, 대규모 $N, m$에 대해 분산이 포아송 분산보다 엄격히 작음을 보여준다.

의의 및 주장
본 논문은 단순한 차수(order-of-magnitude) 경계를 넘어 명시적인 상수와 정밀한 함수 형태를 제시함으로써, 이러한 문제들에 대한 **최초의 정밀한 점근법(sharp asymptotics)**을 제공한다고 주장한다.

• 정밀도: $O(N^{-1})$의 수렴 속도나 명시되지 않은 상수를 가진 지수적 수렴을 확립했던 이전 연구들과 달리, 본 논문은 정확한 리딩 상수(예: $e^{2/p}$ 및 특정 라게르 다항식의 적분을 포함하는 상수)를 도출한다.
• 통합: 본 논문은 모든 스케일링 체제에 걸쳐 트레이스 통계의 거동을 통합하며, 독립적인 블록들로의 스펙트럼 분해가 어떻게 전역 가우시안 변동과 국소 포아송 극한을 모두 지배하는지 보여준다.
• 보정 메커니즘: 핵심적인 통찰은 리딩 총 변동 상수가 단순히 첫 번째 비가우시안 모멘트에 의해 결정되는 것이 아니라, 라게르 보정의 **전체 밴드의 일관된 합산(coherent summation)**에 의해 결정된다는 점이다. 이는 왜 단순한 모멘트 매칭 논리가 정밀한 수렴 속도를 포착하는 데 실패하는지를 설명한다.

본 연구는 유한-$N$의 정확한 항과 정밀한 무한-$N$ 극한 사이의 간극을 메우기 위해 엄밀한 해석적 추정(베셀/프레드홀름)과 조합적 항등식(찰리/이웬스)에 의존한다.