The variance of the number of particles in a set is an important quantity in understanding the statistics of non-interacting fermionic systems in low dimensions. An exact map of their ground state in a harmonic trap in one and two dimensions to the classical Gaussian unitary and complex Ginibre ensemble, respectively, allows to determine the counting statistics at finite and infinite system size. We will establish two new results in this setup. First, we uncover an interpolating central limit theorem between known results in one and two dimensions, for linear statistics of the elliptic Ginibre ensemble. We find an entire range of interpolating weak non-Hermiticity limits, given by a two-parameter family for the mesoscopic scaling regime. Second, we considerably generalize the proportionality between the number variance and the entanglement entropy between Fermions in a set A and its complement in two dimensions. Previously known only for rotationally invariant sets and external potentials, we prove a holographic principle for general non-rotationally invariant sets and random normal matrices. It states that both number variance and entanglement entropy are proportional to the circumference of A.

We study the individual digits for the absolute value of the characteristic polynomial for the Circular beta-Ensemble. We show that, in the large N limit, the leading digits obey Benford's Law and furthermore that the digits of position increasing with the size of the ensemble become uniformly distributed. The key to the proofs is a bound on the rate of convergence in total variation norm in the CLT for the logarithm of the absolute value of the characteristic polynomial.

We compute the genus 0 free energy for the 2-matrix model with quartic interactions, which acts as a generating function for the Ising model’s partition function on a random, 4-regular, planar graph. This is consistent with the predictions of Kazakov and Boulatov on this model, as well as subsequent confirmation of this formula using combinatorial methods. We also provide a new parametric formula for the free energy and give a characterization of the phase space. Our analysis is based on a steepest descent Riemann–Hilbert analysis of the associated biorthogonal polynomials and the corresponding isomonodromic τ-function. A key ingredient in the analysis is a parametrization of the spectral curve. This analysis lays the groundwork for the subsequent study of the multicritical point, which we will study in a forthcoming work.

We study the limiting behavior of random lozenge tilings of the hexagon with a q-Racah weight as the size of the hexagon grows large. Based on the asymptotic behavior of the recurrence coefficients of the q-Racah polynomials, we give a new proof for the fact that the height function for a random tiling concentrates near a deterministic limit shape and that the global fluctuations are described by the Gaussian free field. These results were recently proved using (dynamic) loop equation techniques. In this paper, we extend the recurrence coefficient approach that was developed for (dynamic) orthogonal polynomial ensembles to the setting of q-orthogonal polynomials. An interesting feature is that the complex structure is easily found from the limiting behavior of the (explicitly known) recurrence coefficients. A particular motivation for studying this model is that the variational characterization of the limiting height function has an inhomogeneous term. The study of the regularity properties of the minimizer for general variational problems with such inhomogeneous terms is a challenging open problem. In a general setup, we show that the variational problem gives rise to a natural complex structure associated with the same Beltrami equation as in the homogeneous situation. We also derive a relation between the complex structure and the complex slope. In the case of the q-Racah weighting of lozenge tilings of the hexagon, our representation of the limit shape and their fluctuations in terms of the recurrence coefficients allows us to verify this relation explicitly.

We study random domino tilings of the Aztec diamond with a biased 2 × 2 periodic weight function and associate a linear flow on an elliptic curve to this model. Our main result is a double integral formula for the correlation kernel, in which the integrand is expressed in terms of this flow. For special choices of parameters the flow is periodic, and this allows us to perform a saddle point analysis for the correlation kernel. In these cases we compute the local correlations in the smooth disordered (or gaseous) region. The special example in which the flow has period six is worked out in more detail, and we show that in that case the boundary of the rough disordered region is an algebraic curve of degree eight.

This paper is motivated by computing correlations for domino tilings of the Aztec diamond. It is inspired by two of the three distinct methods that have recently been used in the simplest case of a doubly periodic weighting, that is, the two-periodic Aztec diamond. One of the methods, powered by the domino shuffle, involves inverting the Kasteleyn matrix giving correlations through the local statistics formula. Another of the methods, driven by a Wiener–Hopf factorization for two-by-two matrix-valued functions, involves the Eynard–Mehta Theorem. For arbitrary weights, the Wiener–Hopf factorization can be replaced by an LU- and UL-decomposition, based on a matrix refactorization, for the product of the transition matrices. This paper shows that, for arbitrary weightings of the Aztec diamond, the evolution of the face weights under the domino shuffle and the matrix refactorization is the same. In particular, these dynamics can be used to find the inverse of the LGV matrix in the Eynard–Mehta Theorem.

The elliptic Ginibre ensemble of complex non-Hermitian random matrices allows us to interpolate between the rotationally invariant Ginibre ensemble and the Gaussian unitary ensemble of Hermitian random matrices. It corresponds to a two-dimensional one-component Coulomb gas in a quadrupolar field at inverse temperature beta = 2. Furthermore, it represents a determinantal point process in the complex plane with the corresponding kernel of planar Hermite polynomials. Our main tool is a saddle point analysis of a single contour integral representation of this kernel. We provide a unifying approach to rigorously derive several known and new results of local and global spectral statistics, including in higher dimensions. First, we prove the global statistics in the elliptic Ginibre ensemble first derived by Forrester and Jancovici [Int. J. Mod. Phys. A 11, 941 (1996)]. The limiting kernel receives its main contribution from the boundary of the limiting elliptic droplet of support. In the Hermitian limit, there is a known correspondence between non-interacting fermions in a trap in d real dimensions R-d and the d-dimensional harmonic oscillator. We present a rigorous proof for the local d-dimensional bulk (sine) and edge (Airy) kernel first defined by Dean et al. [Europhys. Lett. 112, 60001 (2015)], complementing the recent results by Deleporte and Lambert [arXiv:2109.02121 (2021)]. Using the same relation to the d-dimensional harmonic oscillator in d complex dimensions C-d, we provide new local bulk and edge statistics at weak and strong non-Hermiticity, where the former interpolates between correlations in d real and d complex dimensions. For C-d with d = 1, this corresponds to non-interacting fermions in a rotating trap.

We study the fluctuations of linear statistics with polynomial test functions for Multiple Orthogonal Polynomial Ensembles. Multiple Orthogonal Polynomial Ensembles form an important class of determinantal point processes that include random matrix models such as the GUE with external source, complex Wishart matrices, multi-matrix models and others. Our analysis is based on the recurrence matrix for the multiple orthogonal polynomials, that is constructed out of the nearest neighbor recurrences. If the coefficients for the nearest neighbor recurrences have limits, then we show that the right-limit of this recurrence matrix is a matrix that can be viewed as representation of a Toeplitz operator with respect to a non-standard basis. This will allow us to prove Central Limit Theorems for linear statistics of Multiple Orthogonal Polynomial Ensembles. A particular novelty is the use of the Baker-Campbell-Hausdorff formula to prove that the higher cumulants of the linear statistics converge to zero. We illustrate the main results by discussing Central Limit Theorems for the Gaussian Unitary Ensembles with external source, complex Wishart matrices and specializations of Schur measure related to multiple Charlier, multiple Krawtchouk and multiple Meixner polynomials. (C) 2021 The Authors. Published by Elsevier Inc.

We analyze domino tilings of the two-periodic Aztec diamond by means of matrix valued orthogonal polynomials that we obtain from a reformulation of the Aztec diamond as a non-intersecting path model with periodic transition matrices. In a more general framework we express the correlation kernel for the underlying determinantal point process as a double contour integral that contains the reproducing kernel of matrix valued orthogonal polynomials. We use the Riemann-Hilbert problem to simplify this formula for the case of the two-periodic Aztec diamond. In the large size limit we recover the three phases of the model known as solid, liquid and gas. We describe the fine asymptotics for the gas phase and at the cusp points of the liquid-gas boundary, thereby complementing and extending results of Chhita and Johansson.

We study a one-parameter family of probability measures on lozenge tilings of large regular hexagons that interpolates between the uniform measure on all possible tilings and a particular fully frozen tiling. The description of the asymptotic behavior can be separated into two regimes: the low and the high temperature regime. Our main results are the computations of the disordered regions in both regimes and the limiting densities of the different lozenges there. For low temperatures, the disordered region consists of two disjoint ellipses. In the high temperature regime the two ellipses merge into a single simply connected region. At the transition from the low to the high temperature a tacnode appears. The key to our asymptotic study is a recent approach introduced by Duits and Kuijlaars providing a double integral representation for the correlation kernel. One of the factors in the integrand is the Christoffel-Darboux kernel associated to polynomials that satisfy non-Hermitian orthogonality relations with respect to a complex-valued weight on a contour in the complex plane. We compute the asymptotic behavior of these orthogonal polynomials and the Christoffel-Darboux kernel by means of a Riemann-Hilbert analysis. After substituting the resulting asymptotic formulas into the double integral we prove our main results by classical steepest descent arguments.