We consider N × N Gaussian random matrices, whose average density of eigenvalues has the Wigner semicircle form over [-√2],√2]. For such matrices, using a Coulomb gas technique, we compute the large N behavior of the probability P(N,L)(N(L)) that N(L) eigenvalues lie within the box [-L,L]. This probability scales as P(N,L)(N(L) = κ(L)N) ≈ exp(-βN(2)ψ(L)(κ(L))), where β is the Dyson index of the ensemble and ψ(L)(κ(L)) is a β-independent rate function that we compute exactly. We identify three regimes as L is varied: (i) N(-1)≪L < √2 (bulk), (ii) L∼√2 on a scale of O(N(-2/3)) (edge), and (iii) L > sqrt[2] (tail). We find a dramatic nonmonotonic behavior of the number variance V(N)(L) as a function of L: after a logarithmic growth ∝ln(NL) in the bulk (when L∼O(1/N)), V(N)(L) decreases abruptly as L approaches the edge of the semicircle before it decays as a stretched exponential for L > sqrt[2]. This "dropoff" of V(N)(L) at the edge is described by a scaling function V(β) that smoothly interpolates between the bulk (i) and the tail (iii). For β = 2 we compute V(2) explicitly in terms of the Airy kernel. These analytical results, verified by numerical simulations, directly provide for β = 2 the full statistics of particle-number fluctuations at zero temperature of 1D spinless fermions in a harmonic trap.
Let P_{β}^{(V)}(N_{I}) be the probability that a N×Nβ-ensemble of random matrices with confining potential V(x) has N_{I} eigenvalues inside an interval I=[a,b] on the real line. We introduce a general formalism, based on the Coulomb gas technique and the resolvent method, to compute analytically P_{β}^{(V)}(N_{I}) for large N. We show that this probability scales for large N as P_{β}^{(V)}(N_{I})≈exp[-βN^{2}ψ^{(V)}(N_{I}/N)], where β is the Dyson index of the ensemble. The rate function ψ^{(V)}(k_{I}), independent of β, is computed in terms of single integrals that can be easily evaluated numerically. The general formalism is then applied to the classical β-Gaussian (I=[-L,L]), β-Wishart (I=[1,L]), and β-Cauchy (I=[-L,L]) ensembles. Expanding the rate function around its minimum, we find that generically the number variance var(N_{I}) exhibits a nonmonotonic behavior as a function of the size of the interval, with a maximum that can be precisely characterized. These analytical results, corroborated by numerical simulations, provide the full counting statistics of many systems where random matrix models apply. In particular, we present results for the full counting statistics of zero-temperature one-dimensional spinless fermions in a harmonic trap.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.