[논문 리뷰] The eigenvectors of Gaussian matrices with an external source

Romain Allez, Joël Bun외 1명|arXiv (Cornell University)|2014. 12. 22.

Random Matrices and Applications참고 문헌 15인용 수 14

한 줄 요약

이 논문은 대칭 브라운 운동을 통해 외부 결정론적 소스에 의해 흐트러진 큰 가우시안 행렬의 고유벡터 역학을 다이슨 브라운 운동을 통해 연구한다. 고유벡터 오버랩에 대한 자율적인 진화 방정식을 유도하며, 고립된 고유벡터에 대해 그 오버랩이 적절히 중심화되고 스케일링된 후 √N으로 정규분포로 수렴함을 증명한다. 이때 분산은 시간 t에 따라 명시적으로 결정되며, 고유값이 스펙트럼의 밀도 영역 내에 위치한 정도에 따라 달라진다.

ABSTRACT

We consider a diffusive matrix process $(X_t)_{t\ge 0}$ defined as $X_t:=A+H_t$ where $A$ is a given deterministic Hermitian matrix and $(H_t)_{t\ge 0}$ is a Hermitian Brownian motion. The matrix $A$ is the "external source" that one would like to estimate from the noisy observation $X_t$ at some time $t>0$. We investigate the relationship between the non-perturbed eigenvectors of the matrix $A$ and the perturbed eigenstates at some time $t$ for the three relevant scaling relations between the time $t$ and the dimension $N$ of the matrix $X_t$. We determine the asymptotic (mean-squared) projections of any given non-perturbed eigenvector $|ψ_j^0 angle$, associated to an eigenvalue $a_j$ of $A$ which may lie inside the bulk of the spectrum or be isolated (spike) from the other eigenvalues, on the orthonormal basis of the perturbed eigenvectors $|ψ_i^t angle,i eq j$. We derive a Burgers type evolution equation for the local resolvent $(z-X_t)_{ii}^{-1}$, describing the evolution of the local density of a given initial state $|ψ_j ^0 angle$. We are able to solve this equation explicitly in the large $N$ limit, for any initial matrix $A$. In the case of one isolated eigenvector $|ψ_j^0 angle$, we prove a central limit Theorem for the overlap $\langle ψ_j^0|ψ_j^t angle$. When properly centered and rescaled by a factor $\sqrt{N}$, this overlap converges in law towards a centered Gaussian distribution with an explicit variance depending on $t$. Our method is based on analyzing the eigenvector flow under the Dyson Brownian motion.

연구 동기 및 목표

• 결정론적 행렬 A의 고유벡터가 헤르미트 브라운 운동에 의해 어떻게 변화하는지 이해하기 위해
• 큰 N 극한에서 비편향 상태와 편향 상태 사이의 고유벡터 오버랩의 渐진적 행동을 분석하기 위해
• 국소 리졸베이트와 고유벡터 사영에 대한 폐쇄형 진화 방정식을 유도하기 위해
• 고립된 고유벡터와 그 편향된 대응 고유벡터 사이의 오버랩에 대한 중심극한정리 수립하기 위해

제안 방법

• 행렬 과정을 $ X_t = A + H_t $로 모델링하며, 여기서 $ H_t $는 $ 1/√N $로 스케일된 i.i.d. 가우시안 성분을 가진 헤르미트 브라운 운동이다.
• 고유값과 고유벡터에 대한 확률적 미분방정식을 사용하여 다이슨 브라운 운동 하에서 고유벡터의 흐름을 분석한다.
• 고정된 비편향 고유벡터 $ |\psi_j^0\rangle $가 편향된 고유벡터 기저 $ |\psi_i^t\rangle $에 사영되는 데 대한 자율적 진화 방정식 (2.5)를 유도한다.
• 국소 리졸베이트 $ (z - X_t)^{-1}_{ii} $를 사용하여 주어진 초기 상태의 국소 밀도의 시간 진화를 묘사하는 버거스형 방정식을 유도한다.
• 스티엘티jes 변환 기법을 적용하고 일반적인 초기 행렬 A에 대해 큰 N 극한에서 해당 방정식을 해결한다.
• 재귀적 적분 표현과 샤클레의 관계를 활용하여 오버랩 분포의 모멘트가 $ g^{(2n)}(t) = (2n-1)!! \cdot g^{(2)}(t)^n $를 만족함을 증명하며, 이는 고립된 고유벡터에 대한 중심극한정리로 이어진다.

실험 결과

연구 질문

• RQ1 시간에 따라 변화하는 헤르미트 브라운 운동이 첨가될 때 결정론적 행렬 A의 고유벡터는 어떻게 변화하는가?
• RQ2 큰 N 극한에서 비편향 고유벡터가 편향된 고유벡터 기저에 투영되는 행동은 어떻게 되는가?
• RQ3 고립된 고유벡터와 그 편향된 대응 고유벡터 사이의 오버랩은 시간이 지남에 따라 어떻게 행동하며, 알려진 분포로 수렴하는가?
• RQ4 국소 리졸베이트에 대한 폐쇄형 진화 방정식을 유도할 수 있는가? 이는 고유벡터 흐름의 역학을 기술한다.
• RQ5 고립된 고유벡터에 대해 오버랩 $ \langle \psi_j^0 | \psi_j^t \rangle $의 분포는 무엇이며, 큰 N 극한에서의 분산은 무엇인가?

주요 결과

• 고립된 고유벡터에 대해 오버랩 $ \langle \psi_j^0 | \psi_j^t \rangle $는 적절히 중심화되고 $ \sqrt{N} $으로 스케일링된 후 법칙적으로 중심 정규분포로 수렴한다.
• 오버랩의 극한 분산은 시간 t에 의해 명시적으로 결정되며, 초기 고유값이 스펙트럼의 밀도 영역 내에 위치한 정도에 따라 달라진다.
• 논문은 국소 리졸베이트 $ (z - X_t)^{-1}_{ii} $에 대한 버거스형 진화 방정식을 유도하였으며, 이는 주어진 초기 상태의 국소 밀도의 시간 진화를 지배한다.
• 일반적인 초기 행렬 A에 대해, 비편향 고유벡터의 점근적(평균제곱) 사영은 고유벡터 흐름 역학에서 유도된 자율 방정식 (2.5)를 통해 계산된다.
• 오버랩 분포의 모멘트는 $ g^{(2n)}(t) = (2n-1)!! \cdot g^{(2)}(t)^n $를 만족하며, 이는 정규분포 극한과 완전한 모멘트 구조를 제공한다.
• 이 방법은 밀도 영역 내의 고유값과 스파이크 고유값을 모두 다룰 수 있으며, 초기 고유값이 밀도 영역 내에 있는지 여부에 따라 다른 행동을 보임을 보여준다.