Abstract We investigate a system of N Brownian particles with the Coulomb interaction in any dimension \(d\ge 2\), and we assume that the initial data are independent and identically distributed with a common density \(\rho _0\) satisfying \(\int _{\mathbb {R}^{d}}\rho _0\ln \rho _0\,\hbox {d}x<\infty \) and \(\rho _0\in L^{\frac{2d}{d+2}} (\mathbb {R}^{d}) \cap L^1(\mathbb {R}^{d}, (1+|x|^2)\,\hbox {d}x)\). We prove that there exists a unique global strong solution for this interacting partsicle system and there is no collision among particles almost surely. For \(d=2\), we rigorously prove the propagation of chaos for this particle system globally in time without any cutoff in the following sense. When \(N\rightarrow \infty \), the empirical measure of the particle system converges in law to a probability measure and this measure possesses a density which is the unique weak solution to the mean-field Poisson–Nernst–Planck equation of single component.

We study the random field Ising model on Z^2 where the external field is given by i.i.d. Gaussian variables with mean zero and positive variance. We show that the effect of boundary conditions on the magnetization in a finite box decays exponentially in the distance to the boundary.

We consider Green's functions $G(z):=(H-z)^{-1}$ of Hermitian random band matrices $H$ on the $d$-dimensional lattice $(\Z/L\Z)^d$. The entries $h_{xy}=\overline h_{yx}$ of $H$ are independent centered complex Gaussian random variables with variances $s_{xy}=\mathbb E|h_{xy}|^2$, which satisfy a banded profile so that $s_{xy}$ is negligible if $|x-y|$ exceeds the band width $W$. For any fixed $n\in \N$, we construct an expansion of the $T$-variable, $T_{xy}=|m|^2 \sum_{\alpha}s_{x\alpha}|G_{\alpha y}|^2$, with an error $\OO(W^{-nd/2})$, and use it to prove a local law on the Green's function. This $T$-expansion was the main tool to prove the delocalization and quantum diffusion of random band matrices for dimensions $d\ge 8$ in part I \cite{PartI_high} of this series.

We analyze the random fluctuations of several multiscale algorithms, such as the multiscale finite element method (MsFEM) and the finite element heterogeneous multiscale method (HMM), that have been developed to solve partial differential equations with highly heterogeneous coefficients. Such multiscale algorithms are often shown to correctly capture the homogenization limit when the highly oscillatory random medium is stationary and ergodic. This paper is concerned with the random fluctuations of the solution about the deterministic homogenization limit. We consider the simplified setting of the one-dimensional elliptic equation, where the theory of random fluctuations is well understood. We develop a fluctuation theory for the multiscale algorithms in the presence of random environments with short-range and long-range correlations. For a given mesh size h, we show that the fluctuations converge in distribution in the space of continuous paths to Gaussian processes as the correlation length ε→0. We next derive the limit of such Gaussian processes as h→0 and compare this limit with the distribution of the random fluctuations of the continuous model. When such limits agree, we conclude that the multiscale algorithm captures the random fluctuations accurately and passes the corrector test. This property serves as an interesting benchmark to assess the behavior of the multiscale algorithm in practical situations where the assumptions necessary for the theory of homogenization are not met. What we find is that the computationally more expensive methods MsFEM, and HMM with a choice of parameter δ=h, correctly capture the random fluctuations both for short-range and long-range oscillations in the medium. The less expensive method HMM with δ<h correctly captures the fluctuations for long-range oscillations and strongly amplifies their size in media with short-range oscillations. We present a modified scheme with an intermediate computational cost that captures the random fluctuations in all cases.

We present a generalization of the method of the local relaxation flow to establish the universality of local spectral statistics of a broad class of large random matrices. We show that the local distribution of the eigenvalues coincides with the local statistics of the corresponding Gaussian ensemble provided the distribution of the individual matrix element is smooth and the eigenvalues ${x_j}_{j=1}^N$ are close to their classical location ${\gamma_j}_{j=1}^N$ determined by the limiting density of eigenvalues. Under the scaling where the typical distance between neighboring eigenvalues is of order 1/N, the necessary apriori estimate on the location of eigenvalues requires only to know that $\E |x_j - \gamma_j |^2 \le N^{-1-\e}$ on average. This information can be obtained by well established methods for various matrix ensembles. We demonstrate the method by proving local spectral universality for Wishart matrices.