BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Percy Deift appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
This is a survey of Harold Widom's work in random matrices. We start with his pioneering papers on the sine-kernel determinant, continue with his and Craig Tracy's groundbreaking results concerning the distribution functions of random matrix theory, touch on the remarkable universality of the Tracy-Widom distributions in mathematics and physics, and close with Tracy and Widom's remarkable work on the asymmetric simple exclusion process.
We study the determinant $\det(I-γK_s), 0<γ<1$, of the integrable Fredholm operator $K_s$ acting on the interval $(-1,1)$ with kernel $K_s(λ, μ)= \frac{\sin s(λ- μ)}{π(λ-μ)}$. This determinant arises in the analysis of a log-gas of interacting particles in the bulk-scaling limit, at inverse temperature $β=2$, in the presence of an external potential $v=-\frac{1}{2}\ln(1-γ)$ supported on an interval of length $\frac{2s}π$. We evaluate, in particular, the double scaling limit of $\det(I-γK_s)$ as $s\rightarrow\infty$ and $γ\uparrow 1$, in the region $0\leqκ=\frac{v}{s}=-\frac{1}{2s}\ln(1-γ)\leq 1-δ$, for any fixed $0<δ<1$. This problem was first considered by Dyson in \cite{Dy1}.
In this paper we continue our analysis \cite{BDIK} of the determinant $\det(I-γK_s),γ\in(0,1)$ where $K_s$ is the trace class operator acting in $L^2(-1,1)$ with kernel $K_s(λ,μ)=\frac{\sin s(λ-μ)}{π(λ-μ)}$. In \cite{BDIK} various key asymptotic results were stated and utilized, but without proof: Here we provide the proofs (see Theorem 1.2 and Proposition 1.3 below).
We consider the Laguerre Unitary Ensemble (aka, Wishart Ensemble) of sample covariance matrices $A = XX^*$, where $X$ is an $N \times n$ matrix with iid standard complex normal entries. Under the scaling $n = N + \lfloor \sqrt{ 4 c N} \rfloor$, $c > 0$ and $N \rightarrow \infty$, we show that the rescaled fluctuations of the smallest eigenvalue, largest eigenvalue and condition number of the matrices $A$ are all given by the Tracy--Widom distribution ($β= 2$). This scaling is motivated by the study of the solution of the equation $Ax=b$ using the conjugate gradient algorithm, in the case that $A$ and $b$ are random: For such a scaling the fluctuations of the halting time for the algorithm are empirically seen to be universal.
The authors present evidence for universality in numerical computations with random data. Given a (possibly stochastic) numerical algorithm with random input data, the time (or number of iterations) to convergence (within a given tolerance) is a random variable, called the halting time. Two-component universality is observed for the fluctuations of the halting time, i.e., the histogram for the halting times, centered by the sample average and scaled by the sample variance, collapses to a universal curve, independent of the input data distribution, as the dimension increases. Thus, up to two components, the sample average and the sample variance, the statistics for the halting time are universally prescribed. The case studies include six standard numerical algorithms, as well as a model of neural computation and decision making. A link to relevant software is provided in for the reader who would like to do computations of his'r own.
We present the results of an empirical study of the performance of the QR algorithm (with and without shifts) and the Toda algorithm on random symmetric matrices. The random matrices are chosen from six ensembles, four of which lie in the Wigner class. For all three algorithms, we observe a form of universality for the deflation time statistics for random matrices within the Wigner class. For these ensembles, the empirical distribution of a normalized deflation time is found to collapse onto a curve that depends only on the algorithm, but not on the matrix size or deflation tolerance provided the matrix size is large enough (see Figure 4, Figure 7 and Figure 10). For the QR algorithm with the Wilkinson shift, the observed universality is even stronger and includes certain non-Wigner ensembles. Our experiments also provide a quantitative statistical picture of the accelerated convergence with shifts.
Sigma-Delta modulation is a popular method for analog-to-digital conversion of bandlimited signals that employs coarse quantization coupled with oversampling. The standard mathematical model for the error analysis of the method measures the performance of a given scheme by the rate at which the associated reconstruction error decays as a function of the oversampling ratio $λ$. It was recently shown that exponential accuracy of the form $O(2^{-rλ})$ can be achieved by appropriate one-bit Sigma-Delta modulation schemes. By general information-entropy arguments $r$ must be less than 1. The current best known value for $r$ is approximately 0.088. The schemes that were designed to achieve this accuracy employ the "greedy" quantization rule coupled with feedback filters that fall into a class we call "minimally supported". In this paper, we study the minimization problem that corresponds to optimizing the error decay rate for this class of feedback filters. We solve a relaxed version of this problem exactly and provide explicit asymptotics of the solutions. From these relaxed solutions, we find asymptotically optimal solutions of the original problem, which improve the best known exponential error decay rate to $r \approx 0.102$. Our method draws from the theory of orthogonal polynomials; in particular, it relates the optimal filters to the zero sets of Chebyshev polynomials of the second kind.
We consider the one-dimensional focusing nonlinear Schrödinger equation (NLS) with a delta potential and even initial data. The problem is equivalent to the solution of the initial/boundary problem for NLS on a half-line with Robin boundary conditions at the origin. We follow the method of Bikbaev and Tarasov which utilizes a Bäcklund transformation to extend the solution on the half-line to a solution of the NLS equation on the whole line. We study the asymptotic stability of the stationary 1-soliton solution of the equation under perturbation by applying the nonlinear steepest-descent method for Riemann-Hilbert problems introduced by Deift and Zhou. Our work strengthens, and extends, earlier work on the problem by Holmer and Zworski.
We present new and streamlined proofs of various formulae for products and ratios of characteristic polynomials of random Hermitian matrices that have appeared recently in the literature.
In this announcement we present a general and new approach to analyzing the asymptotics of oscillatory Riemann-Hilbert problems. Such problems arise, in particular, in evaluating the long-time behavior of nonlinear wave equations solvable by the inverse scattering method. We will restrict ourselves here exclusively to the modified Korteweg de Vries (MKdV) equation, $$y_t-6y^2y_x+y_{xxx}=0,\qquad -\infty<x<\infty,\ t\ge0, y(x,t=0)=y_0(x),$$ but it will be clear immediately to the reader with some experience in the field, that the method extends naturally and easily to the general class of wave equations solvable by the inverse scattering method, such as the KdV, nonlinear Schrödinger (NLS), and Boussinesq equations, etc., and also to ``integrable'' ordinary differential equations such as the Painlevé transcendents.