Let $\widetilde{X}_{M\times N}$ be a rectangular data matrix with independent real-valued entries $[\widetilde{x}_{ij}]$ satisfying $\mathbb {E}\widetilde{x}_{ij}=0$ and $\mathbb {E}\widetilde{x}^2_{ij}=\frac{1}{M}$, $N,M\to\infty$. These entries have a subexponential decay at the tails. We will be working in the regime $N/M=d_N,\lim_{N\to\infty}d_N\neq0,1,\infty$. In this paper we prove the edge universality of correlation matrices ${X}^{\dagger}X$, where the rectangular matrix $X$ (called the standardized matrix) is obtained by normalizing each column of the data matrix $\widetilde{X}$ by its Euclidean norm. Our main result states that asymptotically the $k$-point ($k\geq1$) correlation functions of the extreme eigenvalues (at both edges of the spectrum) of the correlation matrix ${X}^{\dagger}X$ converge to those of the Gaussian correlation matrix, that is, Tracy-Widom law, and, thus, in particular, the largest and the smallest eigenvalues of ${X}^{\dagger}X$ after appropriate centering and rescaling converge to the Tracy-Widom distribution. The asymptotic distribution of extreme eigenvalues of the Gaussian correlation matrix has been worked out only recently. As a corollary of the main result in this paper, we also obtain that the extreme eigenvalues of Gaussian correlation matrices are asymptotically distributed according to the Tracy-Widom law. The proof is based on the comparison of Green functions, but the key obstacle to be surmounted is the strong dependence of the entries of the correlation matrix. We achieve this via a novel argument which involves comparing the moments of product of the entries of the standardized data matrix to those of the raw data matrix. Our proof strategy may be extended for proving the edge universality of other random matrix ensembles with dependent entries and hence is of independent interest.

In this paper, we derive two subgradient estimates of the CR heat equation in a closed pseudohermitian 3-manifold which are served as the CR version of the Li-Yau gradient estimate. With its applications, we first get a subgradient estimate of the logarithm of the positive solution of the CR heat equation. Secondly, we have the Harnack inequality and upper bound estimate for the heat kernel. Finally, we obtain Perelman-type entropy formulae for the CR heat equation.

This paper is motivated by the problem of subdividing a prismatic mesh to a tetrahedral mesh (without inserting Steiner points) so as to not only match arbitrarily prescribed boundary conditions but also allow arbitrary topologies in the base mesh. We explore all possible combinations of these two factors, and propose a complete solution to this 3D problem by converting it to an equivalent 2D graph problem, called cutting flow problem. For each case, we not only prove the sufficient and necessary condition for the existence of solutions, but also provide linear and provable algorithms to compute a solution whenever there is one.

The aim of the present paper is to suggest that statistical physics provides the correct language to understand the practical behavior of the LLL algorithm, most of which are left unexplained to this day. To this end, we propose sandpile models that imitate LLL with compelling accuracy, and prove for these models some of the most desired statements regarding LLL. We also formulate a few conjectures that formally capture our heuristics and would serve as milestones for further development of the theory.

Robustness of numerical methods has attracted an increasing interest in the community of computational fluid dynamics. One mathematical aspect of robustness for numerical methods is the positivity-preserving property. At high Mach numbers or for flows near vacuum, solving the conservative Euler equations may generate negative density or internal energy numerically, which may lead to nonlinear instability and crash of the code. This difficulty is particularly profound for high order methods, for multi-material flows and for problems with moving meshes, such as the Lagrangian methods. In this paper, we construct both first order and uniformly high order accurate conservative Lagrangian schemes which preserve positivity of physically positive variables such as density and internal energy in the simulation of compressible multi-material flows with general equations of state (EOS). We first develop a positivity-preserving approximate Riemann solver for the Lagrangian scheme solving compressible Euler equations with both ideal and non ideal EOS. Then we design a class of high order positivity-preserving and conservative Lagrangian schemes by using the essentially non-oscillatory (ENO) reconstruction, the strong stability preserving (SSP) high order time discretizations and the positivity-preserving scaling limiter which can be proven to maintain conservation and uniformly high order accuracy and is easy to implement. One-dimensional and two-dimensional numerical tests for the positivity-preserving Lagrangian schemes are provided to demonstrate the effectiveness of these methods.