← Back to arXiv
arXivProbabilityarXiv:2609.21140

Inverse of a Sum of Random Matrices: A Dynamical Mean-Field Approach

The paper tackles a fundamental problem in random matrix theory: when you add two large random matrices together and then invert the result, what does that inverse look like, and how well can you approximate it with something simpler? Specifically, it considers a structured sum where one matrix is randomly rotated (in a uniformly random way) before being added to another, and asks how close the inverse of this sum is to a simpler "deterministic equivalent" matrix that is much easier to work with. This kind of problem comes up throughout statistics, signal processing, and machine learning, wherever large covariance matrices or regularized least-squares problems appear.

The main results come in two flavors. The first says that the difference between the true inverse and its deterministic approximation, when applied to a fixed vector, behaves statistically like a Gaussian random vector whose scale is controlled by quantities from free probability theory, specifically the R-transform and free additive convolution. These are tools that describe how the eigenvalue distributions of large random matrices combine. The second result says that even when you look at traces involving a matrix built from up to a polynomial number of rank-one pieces, the approximation error becomes negligible. Together, these results rigorously justify using the simpler deterministic formula in place of the complicated random inverse, and they quantify exactly how good that approximation is.

The proofs rely on an unusual technique called a dynamical mean-field approach, borrowed from statistical physics. Instead of attacking the inverse matrix directly with standard random matrix tools, the authors introduce an artificial time evolution whose long-time behavior converges to the quantity they care about. They then analyze fluctuations around this dynamics using objects called free Appell polynomials, which are algebraic structures adapted to free probability. This approach is notable because it sidesteps many of the combinatorial difficulties that traditional methods face, and it may be useful for studying other complex random matrix problems where direct approaches stall.

Read original →