Universal scaling and criticality of extremes in random matrix theory

We present a random matrix realization of a two-dimensional percolation model with the occupation probability p. We find that the behavior of the model is governed by the two first extreme eigenvalues. While the second extreme eigenvalue resides on the moving edge of the semicircle bulk distribution with an additional semicircle functionality on p, the first extreme exhibits a disjoint isolated Gaussian statistics which is responsible for the emergence of a rich finite-size scaling and criticality. Our extensive numerical simulations along with analytical arguments unravel the power-law divergences due to the coalescence of the first two extreme eigenvalues in the thermodynamic limit. We develop a scaling law that provides a universal framework in terms of a set of scaling exponents uncovering the full finite-size scaling behavior of the extreme eigenvalue's fluctuation. Our study may provide a simple practical approach to capture the criticality in complex systems and their inverse problems with a possible extension to the interacting systems.