The {\it{eigenvector empirical spectral distribution}} (VESD) is a useful tool in studying the limiting behavior of eigenvalues and eigenvectors of covariance matrices. In this paper, we study the convergence rate of the VESD of sample covariance matrices to the deformed Mar{\v c}enko-Pastur (MP) distribution. Consider sample covariance matrices of the form $\Sigma^{1/2} X X^* \Sigma^{1/2}$, where $X=(x_{ij})$ is an $M\times N$ random matrix whose entries are independent random variables with mean zero and variance $N^{-1}$, and $\Sigma$ is a deterministic positive-definite matrix. We prove that the Kolmogorov distance between the {\it{expected VESD}} and the deformed MP distribution is bounded by $N^{-1+\epsilon}$ for any fixed $\epsilon>0$, provided that the entries $\sqrt{N}x_{ij}$ have uniformly bounded 6th moments and $|N/M-1|\ge \tau$ for some constant $\tau>0$. This result improves the previous one obtained in [44], which gave the convergence rate $O(N^{-1/2})$ assuming $i.i.d.$ $X$ entries, bounded 10th moment, $\Sigma=I$ and $M<N$. Moreover, we also prove that under the finite $8$th moment assumption, the convergence rate of the VESD is $O(N^{-1/2+\epsilon})$ almost surely for any fixed $\epsilon>0$, which improves the previous bound $N^{-1/4+\epsilon}$ in [44].

In 1979, G. Parisi predicted a variational formula for the thermodynamic limit of the free energy in the Sherrington-Kirkpatrick model and described the role played by its minimizer. This formula was verified in the seminal work of Talagrand and later generalized to the mixed p-spin models by Panchenko. In this paper, we prove that the minimizer in Parisi's formula is unique at any temperature and external field by establishing the strict convexity of the Parisi functional.

We prove an estimate on the smallest singular value of a multiplicatively and additively deformed random rectangular matrix. Suppose $n\le N \le M \le dN$ for some constant $d\ge 1$. Let $X$ be an $M\times n$ random matrix with independent and identically distributed entries, which have zero mean, unit variance and arbitrarily high moments. Let $T$ be an $N\times M$ deterministic matrix with comparable singular values $c\le s_{N}(T) \le s_{1}(T) \le c^{-1}$ for some constant $c>0$. Let $A$ be an $N\times n$ deterministic matrix with $\|A\|=O(\sqrt{N})$. Then we prove that for any $\epsilon>0$, the smallest singular value of $TX-A$ is larger than $N^{-\epsilon}(\sqrt{N}-\sqrt{n-1})$ with high probability. If we assume further the entries of $X$ have subgaussian decay, then the smallest singular value of $TX-A$ is at least of the order $\sqrt{N}-\sqrt{n-1}$ with high probability, which is an essentially optimal estimate.

We prove that uniform random quadrangulations of the sphere with$n$faces, endowed with the usual graph distance and renormalized by$n$^{−1/4}, converge as$n$→$∞$in distribution for the Gromov–Hausdorff topology to a limiting metric space. We validate a conjecture by Le Gall, by showing that the limit is (up to a scale constant) the so-called$Brownian map$, which was introduced by Marckert–Mokkadem and Le Gall as the most natural candidate for the scaling limit of many models of random plane maps. The proof relies strongly on the concept of$geodesic stars$in the map, which are configurations made of several geodesics that only share a common endpoint and do not meet elsewhere.

No abstract uploaded!

No abstract uploaded!

The tip multifractal spectrum of a 2-dimensional curve is one way to describe the behavior of the uniformizing conformal map of the complement near the tip. We give the tip multifractal spectrum for a Schramm–Loewner evolution (SLE) curve, we prove that the spectrum is valid with probability 1, and we give applications to the scaling of harmonic measure at the tip.

This paper analyzes the random fluctuations obtained by a heterogeneous multi-scale first-order finite element method applied to solve elliptic equations with a random potential. Several multiscale numerical algorithms have been shown to correctly capture the homogenized limit of solutions of elliptic equations with coefficients modeled as stationary and ergodic random fields. Because theoretical results are available in the continuum setting for such equations, we consider here the case of a secondorder elliptic equations with random potential in two dimensions of space. We show that the random fluctuations of such solutions are correctly estimated by the heterogeneous multi-scale algorithm when appropriate fine-scale problems are solved on subsets that cover the whole computational domain. However, when the fine-scale problems are solved over patches that do not cover the entire domain, the random fluctuations may or may not be estimated accurately. In the case of random potentials with short-range interactions, the variance of the random fluctuations is amplified as the inverse of the fraction of the medium covered by the patches. In the case of random potentials with long-range interactions, however, such an amplification does not occur and random fluctuations are correctly captured independent of the (macroscopic) size of the patches. These results are consistent with those obtained in [9] for more general equations in the one-dimensional setting and provide indications on the loss in accuracy that results from using coarser, and hence computationally less intensive, algorithms.

It is well known that the spectral measure of eigenvalues of a rescaled square non-Hermitian random matrix with independent entries satises the circular law. We consider the product $T X$, where $T$ is a deterministic $N\times M$ matrix and $X$ is a random $M\times  N$ matrix with independent entries having zero mean and variance $(N \wedge M) ^{-1}$. We prove a general local circular law for the empirical spectral distribution (ESD) of $TX$ at any point $z$ away from the unit circle under the assumptions that$ N \sim M$, and the matrix entries $X_{ij} have suciently high moments. More precisely, if $z$ satisfes $| |z|
- 1| \ge \tau$  for arbitrarily small $\tau > 0, the ESD of $TX$ converges to $\tilde{X}_D (z) d A(z)$, where $\tilde{X}_D$ is a rotation-invariant function determined by the singular valuau of $T$ and $dA$ denotes the Lebesgue measure on $\mathbb{C}$. The local circular law is valid around $z$ up to scale $(N \wedge M) ^{-1/4+\epsilon}$ for any  $\epsilon >0$. Moreover, if $|z| >1$ or the matrix entries of $X$ have vanishing third moments, the local circular law is valid around $z$ up to scale $(N \wedge M) ^{-1/2+\epsilon}$ for any  $\epsilon >0$.

Let $K$ be a self-similar set satisfying the open set condition. Following Kaimanovich’s elegant idea, it has been proved that on the symbolic space $X$ of $K$ a natural augmented tree structure $E$ exists; it is hyperbolic, and the hyperbolic boundary $\partial_H X$ with the Gromov metric is H\"older equivalent to $K$. In this paper we consider certain reversible random walks with return ratio $0 < \lambda < 1$ on $(X,E)$. We show that the Martin boundary ${\mathcal M}$ can be identified with $\partial_H X$ and $K$. With this setup and a device of Silverstein, we obtain precise estimates of the Martin kernel and the Na\"im kernel in terms of the Gromov product. Moreover, the Naïm kernel turns out to be a jump kernel satisfying the estimate $\Theta (\xi,\eta) \asymp |\xi-\eta|^{-(\alpha+\beta)}$, where $\alpha$ is the Hausdorff dimension of $K$ and $\beta$ depends on $\lambda$. For suitable $\beta$, the kernel defines a regular non-local Dirichlet form on $K$. This extends the results of Kigami concerning random walks on certain trees with Cantor-type sets as boundaries.

A fundamental and very well studied region of the Erdős–Rényi process is the phase transition at$m$∼$n$/2 edges in which a giant component suddenly appears. We examine the process beginning with an initial graph. We further examine the Bohman–Frieze process in which edges between isolated vertices are more likely. While the positions of the phase transitions vary, the three processes belong, roughly speaking, to the same universality class. In particular, the growth of the giant component in the barely supercritical region is linear in all cases.

We prove the Wigner-Dyson-Mehta conjecture at fixed energy in the bulk of the spectrum for generalized symmetric and Hermitian Wigner matrices. Previous results concerning the universality of random matrices either require an averaging in the energy parameter or they hold only for Hermitian matrices if the energy parameter is fixed. We develop a homogenization theory of the Dyson Brownian motion and show that microscopic universality follows from mesoscopic statistics

Consider an age-dependent, single-species branching process defined by a progeny number distribution, and a lifetime distribution associated with each independent particle. In this paper, we focus on the associated inverse problem where one wishes to formally solve for the progeny number distribution or the lifetime distribution that defines the Bellman-Harris branching process. We derive results for the existence and uniqueness (the identifiability) of these two distributions given one of two types of information: the extinction time probability of the entire process (extinction time distribution), or the distribution of the total number of particles at one fixed time. We demonstrate that perfect knowledge of the distribution of extinction times allows us to formally determine either the progeny number distribution or the lifetime distribution. Furthermore, we show that these constructions are unique. We then consider “data” consisting of a perfectly known total number distribution given at one specific time. For a process with known progeny number distribution and exponentially distributed lifetimes, we show that the rate parameter is identifiable. For general lifetime distributions, we also show that the progeny distribution is globally unique. Our results are presented through four theorems, each describing the constructions in the four distinct cases.

In this paper, we prove a necessary and sufficient condition for the edge universality of sample covariance matrices with general population. We consider sample covariance matrices of the form $\mathcal Q = TX(TX)^{*}$, where $X$ is an $M_2\times N$ random matrix with $X_{ij}=N^{-1/2}q_{ij}$ such that $q_{ij}$ are $i.i.d.$ random variables with zero mean and unit variance, and $T$ is an $M_1 \times M_2$ deterministic matrix such that $T^* T$ is diagonal. We study the asymptotic behavior of the largest eigenvalues of $\mathcal Q$ when $M:=\min\{M_1,M_2\}$ and $N$ tend to infinity with $\lim_{N \to \infty} {N}/{M}=d \in (0, \infty)$. We prove that the Tracy-Widom law holds for the largest eigenvalue of $\mathcal Q$ if and only if $\lim_{s \rightarrow \infty}s^4 \mathbb{P}(\vert q_{ij} \vert \geq s)=0$ under mild assumptions of $T$. The necessity and sufficiency of this condition for the edge universality was first proved for Wigner matrices by Lee and Yin [34].

In this paper, we consider the universality of the local eigenvalue statistics of random matrices. Our main result shows that these statistics are determined by the first four moments of the distribution of the entries. As a consequence, we derive the universality of eigenvalue gap distribution and$k$-point correlation, and many other statistics (under some mild assumptions) for both Wigner Hermitian matrices and Wigner real symmetric matrices.

No abstract uploaded!

No abstract uploaded!

We consider the edge statistics of large dimensional deformed rectangular matrices of the form $Y_t=Y+\sqrt{t}X,$ where $Y$ is a $p \times n$ deterministic signal matrix whose rank is comparable to $n$, $X$ is a $p\times n$ random noise matrix with centered i.i.d. entries with variance $n^{-1}$, and $t>0$ gives the noise level. This model is referred to as the interference-plus-noise matrix in the study of massive multiple-input multiple-output (MIMO) system, which belongs to the category of the so-called signal-plus-noise model. For the case $t=1$, the spectral statistics of this model have been studied to a certain extent in the literature \cite{DOZIER20071099,DOZIER2007678,VLM2012}. In this paper, we study the singular value and singular vector statistics of $Y_t$ around the right-most edge of the singular value spectrum in the harder regime $n^{-2/3}\ll t \ll 1$. This regime is harder than the $t=1$ case, because on one hand, the edge behavior of the empirical spectral distribution (ESD) of $YY^\top$ has a strong effect on the edge statistics of $Y_tY_t^\top$ since $t\ll 1$ is ``small", while on the other hand, the edge statistics of $Y_t$ is also not merely a perturbation of those of $Y$ since $t\gg n^{-2/3}$ is ``large". Under certain regularity assumptions on $Y,$ we prove the edge universality, eigenvalues rigidity and eigenvector delocalization for the matrices $Y_tY_t^\top$ and $Y_t^\top Y_t$. These results can be used to estimate and infer the massive MIMO system. To prove the main results, we analyze the edge behavior of the asymptotic ESD for $Y_tY_t^\top$, and establish some sharp local laws on the resolvent of $Y_tY_t^\top$. These results can be of independent interest, and used as useful inputs for many other problems regarding the spectral statistics of $Y_t$.

This is the second part of a three part series of papers. In this paper, we consider a general class of $N\times N$ random band matrices $H=(H_{ij})$ whose entries are centered random variables, independent up to a symmetry constraint. We assume that the variances $\mathbb E |H_{ij}|^2$ form a band matrix with typical band width $1\ll W\ll N$. We consider the generalized resolvent of $H$ defined as $G(Z):=(H - Z)^{-1}$, where $Z$ is a deterministic diagonal matrix such that $Z_{ij}=\left(z\mathds{1}_{1\le i \le W}+\wt z\mathds{1}_{ i > W} \right) \delta_{ij}$, with two distinct spectral parameters $z\in \mathbb C_+:=\{z\in \mathbb C:\im z>0\}$ and $\wt z\in \mathbb C_+\cup \mathbb R$. In this paper, we prove a sharp bound for the local law of the generalized resolvent $G$ for $W\gg N^{3/4}$. This bound is a key input for the proof of delocalization and bulk universality of random band matrices in [2]. Our proof depends on a fluctuations averaging bound on certain averages of polynomials in the resolvent entries, which will be proved in [10].

No abstract uploaded!

We establish necessary and sufficient conditions for consistent root reconstruction in continuous-time Markov models with countable state space on bounded-height trees. Here a root state estimator is said to be consistent if the probability that it returns to the true root state converges to 1 as the number of leaves tends to infinity. We also derive quantitative bounds on the error of reconstruction. Our results answer a question of Gascuel and Steel [GS10] and have implications for ancestral sequence reconstruction in a classical evolutionary model of nucleotide insertion and deletion [TKF91].

No abstract uploaded!

A new non-conservative stochastic reaction–diffusion system in which two families of random walks in two adjacent domains interact near the interface is introduced and studied in this paper. Such a system can be used to model the transport of positive and negative charges in a solar cell or the population dynamics of two segregated species under competition. We show that in the macroscopic limit, the particle densities converge to the solution of a coupled nonlinear heat equations. For this, we first prove that propagation of chaos holds by establishing the uniqueness of a new BBGKY hierarchy. A local central limit theorem for reflected diffusions in bounded Lipschitz domains is also established as a crucial tool.

We consider a general class of symmetric or Hermitian random band matrices $H=(h_{xy})_{x,y \in \llbracket 1,N\rrbracket^d}$ in any dimension $d\ge 1$, where the entries are independent, centered random variables with variances $s_{xy}=\mathbb E|h_{xy}|^2$. We assume that $s_{xy}$ vanishes if $|x-y|$ exceeds the band width $W$, and we are interested in the mesoscopic scale with $1\ll W\ll N$. Define the {\it{generalized resolvent}} of $H$ as $G(H,Z):=(H - Z)^{-1}$, where $Z$ is a deterministic diagonal matrix with entries $Z_{xx}\in \mathbb C_+$ for all $x$. Then we establish a precise high-probability bound on certain averages of polynomials of the resolvent entries. As an application of this fluctuation averaging result, we give a self-contained proof for the delocalization of random band matrices in dimensions $d\ge 2$. More precisely, for any fixed $d\ge 2$, we prove that the bulk eigenvectors of $H$ are delocalized in certain averaged sense if $N\le W^{1+\frac{d}{2}}$. This improves the corresponding results in \cite{HeMa2018} under the assumption $N\ll W^{1+\frac{d}{d+1}}$, and in \cite{ErdKno2013,ErdKno2011} under the assumption $N\ll W^{1+\frac{d}{6}}$. For 1D random band matrices, our fluctuation averaging result was used in \cite{PartII,PartI} to prove the delocalization conjecture and bulk universality for random band matrices with $N\ll W^{4/3}$.

We investigate the energy landscape of the mixed even p-spin model with Ising spin configurations. We show that for any given energy level between zero and the maximal energy, with overwhelming probability there exist exponentially many distinct spin configurations such that their energies stay near this energy level. Furthermore, their magnetizations and overlaps are concentrated around some fixed constants. In particular, at the level of maximal energy, we prove that the Hamiltonian exhibits exponentially many orthogonal peaks. This improves the results of Chatterjee [20] and Ding-Eldan-Zhai [29], where the former established a logarithmic size of the number of the orthogonal peaks, while the latter proved a polynomial size. Our second main result obtains disorder chaos at zero temperature and at any external field. As a byproduct, this implies that the fluctuation of the maximal energy is superconcentrated when the external field vanishes and obeys a Gaussian limit law when the external field is present.