트리디아곤 모델의 β 다이슨 브라운 운동 한계

초록

본 논문은 GβE 행렬 과정에 하우스홀더 트리디아곤화 알고리즘을 적용한 뒤, 그 결과로 얻어지는 트리디아곤 행렬 과정의 고정된 k×k 부분행렬이 n→∞에서 어떻게 수렴하는지를 연구한다. β=1에 대해 엄밀히 증명하고, β=2,4에 대해서는 수치 실험을 통해 동일한 형태의 수렴을 제시한다. 이를 바탕으로 동적 β-스톡캐스틱 에어리 연산자의 후보 형태를 제안한다.

상세 분석

논문은 먼저 GβE(β‑Hermite) 행렬을 하우스홀더 변환으로 삼대각(tridiagonal) 형태로 바꾸는 고전적 결과를 상기한다. 이때 얻어지는 정방행렬은 대각 원소가 정규분포, 비대각 원소가 χ‑분포를 따르는 독립적인 확률변수들로 구성된다. 저자들은 이 정적 모델을 시간에 따라 진화하는 GβE 과정에 적용함으로써, 각 시점마다 β‑Dyson Brownian Motion(DBM)으로 기술되는 고유값 동역학을 갖는 트리디아곤 행렬 과정 T(t)를 정의한다. 핵심 질문은 “n이 무한대로 커질 때, 고정된 k×k 왼쪽 위 블록의 엔트리들은 어떤 확률 과정으로 수렴하는가?”이다.

주요 정리(Theorem 3.1)는 β=1인 경우, 대각 원소 a_j(t)와 비대각 원소 b_j(t) (j≤k)가 각각 a_j(t) → √2 OU(2j‑1)(t), (1/√β n)(b_j(t)^2‑β n) → √2 OU(2j)(t) 로 수렴함을 보인다. 여기서 OU(m) 은 평균 0, 변동계수 √(2m)인 정규화된 정상 오르엔슈틴 과정이며, 서로 독립이다. 비대각 원소의 스케일링은 χ‑분포의 중심극한을 이용해 OU 과정으로 변환한다는 점이 핵심이다. 증명은 하우스홀더 변환 단계마다 발생하는 x_k(t) 벡터가 고차원 OU 과정임을 이용하고, 차원 n이 충분히 크면 (n‑k‑1) 차원에서의 구형 대칭성으로 인해 정규화된 반경 OU(R OU) 과정으로 수렴함을 보인다. 저자들은 이 과정을 정밀히 다루기 위해 Chebyshev 다항식, 비교단 파티션, 자유 확률 이론 등을 활용한다.

β=2,4에 대해서는 동일한 구조가 기대되지만, 현재는 수치 실험에 의존한다. 시뮬레이션 결과는 k=1,2,3에 대해 (1) 대각 원소가 √2 OU(odd)와 일치하고, (2) 비대각 원소가 적절히 중심·스케일링된 후 √2 OU(even)와 일치함을 보여준다. 그러나 저자들은 전체 n×n 행렬의 고유값(특히 가장 큰 k개)과의 직접적인 연결 고리를 아직 확립하지 못했으며, 제안된 동적 β‑스톡캐스틱 에어리 연산자(Section 5)가 실제로는 기대와 다르게 동작한다는 부정적 증거도 제시한다. 이는 제한된 k에 대해서는 수렴이 성립하지만, k가 n에 비례적으로 커질 때는 복잡한 상호작용이 발생한다는 점을 시사한다.

결과적으로, 이 논문은 트리디아곤화된 DBM 과정의 로컬(고정 k) 동역학을 명확히 규명함으로써, 기존의 정적 β‑트리디아곤 모델을 동적 버전으로 확장하는 첫 걸음을 제공한다. 동시에, 전체 스펙트럼 수준에서의 연산자 수렴 문제는 아직 미해결이며, 향후 연구 과제로 남는다.