[논문 리뷰] Convergence Rate for Spectral Distribution of Addition of Random Matrices

Zhigang Bao, László Erdős외 1명|arXiv (Cornell University)|2016. 06. 09.

Random Matrices and Applications참고 문헌 19인용 수 3

한 줄 요약

이 논문은 $ A $와 $ B $가 결정적 헤르미트 행렬이고 $ U $가 하르 분포를 가진 유니터리 행렬일 때, 합 $ H = A + UBU^* $의 스펙트럼 분포에 대해 최적 수렴 속도 $ N^{-1+\gamma} $를 확립한다. 특별히 구성된 보조 양 $ Z_i $를 통한 새로운 변동 평균화 메커니즘을 도입함으로써, 일반적인 요소 간 독립성이 없는 하르 유니터리 행렬의 경우에도, 이론적 고유값 간격에 로그 인자까지 고려한 최적 수렴 속도를 달성한다.

ABSTRACT

Let $A$ and $B$ be two $N$ by $N$ deterministic Hermitian matrices and let $U$ be an $N$ by $N$ Haar distributed unitary matrix. It is well known that the spectral distribution of the sum $H=A+UBU^*$ converges weakly to the free additive convolution of the spectral distributions of $A$ and $B$, as $N$ tends to infinity. We establish the optimal convergence rate ${\frac{1}{N}}$ in the bulk of the spectrum.

연구 동기 및 목표

• 스펙트럼 영역에서 $ H = A + UBU^* $의 스펙트럼 분포가 자유加산 병합 $ \mu_A \boxplus \mu_B $로 수렴하는 최적 수렴 속도를 확립하는 것.
• 하르 유니터리 행렬에서 요소 간 독립성이 부족하여 기존의 변동 평균화 기법이 복잡해지는 문제를 해결하는 것.
• 고차 모멘트 계산을 피하기 위해 적분 by parts를 통한 재귀적 모멘트 추정 기법을 활용한 새로운 변동 평균화 방법을 개발하는 것.
• 녹그 함수의 정규화된 트레이스 $ \frac{1}{N}\mathrm{Tr}\,G $가 최대 $ N^{-1+\gamma} $의 정도로 변동함을 보여, 최적 수렴을 의미하는 것.
• empirical 스펙트럼 측도 $ \mu_H $가 $ N^{-1} $ 속도로 농도를 보임을 증명하는 것. 이는 고유값 간격 스케일과 일치한다.

제안 방법

• 녹그 함수의 대각성분 $ G_{ii} $와 $ (UBU^*G)_{ii} $의 선형 조합인 보조 양 $ Z_i $를 도입하여, $ N^{-1/2} $ 스케일의 변동을 대수적 상쇄 효과로 활용하는 것.
• 잔여 $ N^{-1/2} $ 스케일 변동을 흡수하기 위해 수정항 $ \Upsilon $를 정의하고, 이가 실제로는 $ N^{-1/2} $ 가 아니라 $ N^{-1+\gamma} $ 의 크기임을 보이는 것.
• 복잡한 고차 모멘트 전개를 피하기 위해 가우시안 변수에 기반한 적분 by parts를 활용한 재귀적 모멘트 추정 기법을 사용하는 것.
• 녹그 함수 성분의 종속성 구조를 제어하기 위해 하르 유니터리 행렬의 부분 무작위성 분해 성질을 활용하는 것.
• 녹그 함수 트레이스 추정치로부터 스펙트럼 측도의 수렴 속도를 도출하기 위해 헬퍼-시올스트란드 함수 해석학을 적용하는 것.
• 스케일 $ N^{-1+\gamma} $ 에서 국소 법칙을 기초로 삼는, 버블 영역에서의 분석을 수행하는 것.

실험 결과

연구 질문

• RQ1 스펙트럼 영역에서 $ H = A + UBU^* $의 경험적 스펙트럼 분포가 자유加산 병합 $ \mu_A \boxplus \mu_B $로 수렴하는 최적 수렴 속도는 무엇인가?
• RQ2 하르 유니터리 행렬에서 요소 간 독립성이 없을 경우, 변동 평균화는 어떻게 달성할 수 있는가?
• RQ3 $ N^{-2/3} $ 보다 수렴 속도를 향상시킬 수 있으며, 만약 가능하면 어떤 새로운 기법이 필요할까?
• RQ4 empirical 스펙트럼 측도 $ \mu_H $가 기대값 주위에서의 농도 속도는 무엇인가?
• RQ5 녹그 함수 트레이스 $ \frac{1}{N}\mathrm{Tr}\,G $의 행동은 $ \mu_H $의 수렴과 어떻게 관련이 있는가?

주요 결과

• 스펙트럼 영역에서 $ H = A + UBU^* $의 분포가 $ \mu_A \boxplus \mu_B $로 수렴하는 속도는 $ N^{-1+\gamma} $이며, 임의의 $ \gamma > 0 $ 에 대해 이는 $ N^\gamma $ 인자까지 고려한 최적 수렴 속도이다.
• 녹그 함수의 정규화된 트레이스 $ \frac{1}{N}\mathrm{Tr}\,G $는 최대 $ N^{-1+\gamma} $ 정도의 변동을 보이며, 이는 최적 스케일과 일치한다.
• 특히 구성된 보조 양 $ Z_i $를 통한 새로운 변동 평균화 메커니즘이 성공적으로 구현되었으며, 이는 대수적 구조를 통해 $ N^{-1/2} $ 스케일 변동을 상쇄시킨다.
• 잔여 변동을 흡수하기 위해 도입된 수정항 $ \Upsilon $는 직관적인 크기와는 달리, 실제로는 $ N^{-1/2} $ 가 아니라 $ N^{-1+\gamma} $ 의 크기임을 입증하였다.
• $ \mu_H $가 기대값 주위에서의 농도 속도는 $ N^{-1} $ 이며, 이는 버블 영역에서 최적이다.
• 결과는 $ \mathcal{S}_\mathcal{I}(a,b) $ 내 모든 $ z $ 에 대해 균일하게 성립하며, 수렴 속도는 일반적인 고유값 간격 $ N^{-1} $ 과 정확히 일치하므로 날카로운 결과이다.