No abstract uploaded!

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].

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}$.

The goal of this paper is to prove rigorous results for the behavior of genealogies in a one-dimensional long range biased voter model introduced by Hallatschek and Nelson [Theor. Pop. Biol. 73 (2008) 158–170]. The first step, which is easily accomplished using results of Mueller and Tribe [Probab. Theory Related Fields 102 (1995) 519–545], is to show that when space and time are rescaled correctly, our biased voter model converges to a Wright–Fisher SPDE. A simple extension of a result of Durrett and Restrepo [Ann. Appl. Probab. 18 (2008) 334–358] then shows that the dual branching coalescing random walk converges to a branching Brownian motion in which particles coalesce after an exponentially distributed amount of intersection local time. Brunet et al. [Phys. Rev. E (3) 76 (2007) 041104, 20] have conjectured that genealogies in models of this type are described by the Bolthausen–Sznitman coalescent, see [Proc. Natl. Acad. Sci. USA 110 (2013) 437–442]. However, in the model we study there are no simultaneous coalescences. Our third and most significant result concerns “tracer dynamics” in which some of the initial particles in the biased voter model are labeled. We show that the joint distribution of the labeled and unlabeled particles converges to the solution of a system of stochastic partial differential equations. A new duality equation that generalizes the one Shiga [In Stochastic Processes in Physics and Engineering (Bielefeld, 1986) (1988) 345–355 Reidel] developed for the Wright–Fisher SPDE is the key to the proof of that result.

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 introduce an interacting particle system in which two families of reflected diffusions interact in a singular manner near a deterministic interface I. This 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. A related interacting random walk model with discrete state spaces has recently been introduced and studied in Chen and Fan (2014). In this paper, we establish the functional law of large numbers for this new system, thereby extending the hydrodynamic limit in Chen and Fan (2014) to reflected diffusions in domains with mixed-type boundary conditions, which include absorption (harvest of electric charges). We employ a new and direct approach that avoids going through the delicate BBGKY hierarchy.

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].

We consider ASEP on a bounded interval and on a half‐line with sources and sinks. On the full line, Bertini and Giacomin in 1997 proved convergence under weakly asymmetric scaling of the height function to the solution of the KPZ equation. We prove here that under similar weakly asymmetric scaling of the sources and sinks as well, the bounded interval ASEP height function converges to the KPZ equation on the unit interval with Neumann boundary conditions on both sides (different parameter for each side), and likewise for the half‐line ASEP to KPZ on a half‐line. This result can be interpreted as showing that the KPZ equation arises at the triple critical point (maximal current / high density / low density) of the open ASEP.

We consider the KPZ equation in one space dimension driven by a sta- tionary centred space–time random field, which is sufficiently integrable and mixing, but not necessarily Gaussian. We show that, in the weakly asymmetric regime, the solution to this equation considered at a suitable large scale and in a suitable reference frame converges to the Hopf–Cole solution to the KPZ equation driven by space–time Gaussian white noise. While the limit- ing process depends only on the integrated variance of the driving field, the diverging constants appearing in the definition of the reference frame also depend on higher order moments.

We introduce the dynamical sine-Gordon equation in two space dimensions with parameter β, which is the natural dynamic associated to the usual quantum sine- Gordon model. It is shown that when β2 ∈ (0, 16π ) the Wick renormalised equation 3 is well-posed. In the regime β2 ∈ (0,4π), the Da Prato–Debussche method [J Funct Anal 196(1):180–210, 2002; Ann Probab 31(4):1900–1916, 2003] applies, while for β2 ∈ [4π, 16π ), the solution theory is provided via the theory of regularity structures 3 [Hairer, Invent Math 198(2):269–504, 2014]. We also show that this model arises naturally from a class of 2 + 1-dimensional equilibrium interface fluctuation models with periodic nonlinearities. The main mathematical difficulty arises in the construction of the model for the associated regularity structure where the role of the noise is played by a non-Gaussian random distribution similar to the complex multiplicative Gaussian chaos recently analysed in Lacoin et al.

No abstract uploaded!

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.

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 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].

We derive the alternating arm exponents of the critical Ising model. We obtain six different patterns of alternating boundary arm exponents which correspond to the boundary conditions (⊖⊕), (⊖ free) and (free free), and the alternating interior arm exponents.

We compute the almost-sure Hausdorff dimension of the double points of chordal SLE$_\kappa$ for $\kappa > 4$, confirming a prediction of Duplantier--Saleur (1989) for the contours of the FK model. We also compute the dimension of the cut points of chordal SLE$_\kappa$ for $\kappa > 4$ as well as analogous dimensions for the radial and whole-plane SLE$_\kappa(\rho)$ processes for $\kappa > 0$. We derive these facts as consequences of a more general result in which we compute the dimension of the intersection of two flow lines of the formal vector field $e^{ih/\chi}$, where $h$ is a Gaussian free field and $\chi > 0$, of different angles with each other and with the domain boundary.

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.

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.

Via a Dirichlet form extension theorem and making full use of two-sided heat kernel estimates, we establish quenched invariance principles for ran- dom walks in random environments with a boundary. In particular, we prove that the random walk on a supercritical percolation cluster or among random conductances bounded uniformly from below in a half-space, quarter-space, etc., converges when rescaled diffusively to a reflecting Brownian motion, which has been one of the important open problems in this area. We establish a similar result for the random conductance model in a box, which allows us to improve existing asymptotic estimates for the relevant mixing time. Furthermore, in the uniformly elliptic case, we present quenched invariance prin- ciples for domains with more general boundaries.

No abstract uploaded!

No abstract uploaded!

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.

We develop mathematical models describing the evolution of stochastic age-structured populations. After reviewing existing approaches, we formulate a complete kinetic framework for age-structured interacting populations undergoing birth, death and fission processes in spatially dependent environments. We define the full probability density for the population-size age chart and find results under specific conditions. Connections with more classical models are also explicitly derived. In particular, we show that factorial moments for non-interacting processes are described by a natural generalization of the McKendrick-von Foerster equation, which describes mean-field deterministic behavior. Our approach utilizes mixed-type, multidimensional probability distributions similar to those employed in the study of gas kinetics and with terms that satisfy BBGKY-like equation hierarchies.

In this paper, we consider particle systems with interaction and Brownian motion. We prove that when the initial data is from the sampling of Chorin's method, i.e., the initial vertices are on lattice points hi∈Rd with mass ρ0(hi)hd, where ρ0 is some initial density function, then the regularized empirical measure of the interacting particle system converges in probability to the corresponding mean-field partial differential equation with initial density ρ0, under the Sobolev norm of L∞(L2)∩L2(H1). Our result is true for all those systems when the interacting function is bounded, Lipschitz continuous and satisfies certain regular condition. And if we further regularize the interacting particle system, it also holds for some of the most important systems of which the interacting functions are not. For systems with repulsive Coulomb interaction, this convergence holds globally on any interval [0,t]. And for systems with attractive Newton force as interacting function, we have convergence within the largest existence time of the regular solution of the corresponding Keller-Segel equation.

Abstract In this paper we develop a new martingale method to show the convergence of the regularized empirical measure of many particle systems in probability under a Sobolev norm to the corresponding mean field PDE. Our method works well for the simple case of Fokker Planck equation and we can estimate a lower bound of the rate of convergence. This method can be generalized to more complicated systems with interactions.