Leonid Petrov — arXiv Feed

Joint moments of characteristic polynomials in the circular Jacobi ensemble and Painlevé equations

2026-08-25T07:37:15-04:00

In this paper, we establish a connection between joint moments of characteristic polynomials and their derivatives in the Circular Jacobi Ensemble, a generalisation of the Circular Unitary Ensemble, and solutions of nonlinear Painlevé equations. For finite $N$ , we show that these joint moments are characterised by a solution of the $\sigma$ -Painlevé V equation for all real moment exponents in their admissible range. Under an appropriate large- $N$ scaling limit, we further prove that the limiting joint moments admit a representation in terms of a solution of the $\sigma$ -Painlevé III $'$ equation for a certain range of moment exponents. As applications, we answer a question posed by Assiotis et al. in [Math. Physics. Anal. Geom. 25 (2022), no.2, Paper No. 15, 24pp, Remark 1.7] by showing that the characteristic function of a distinguished random variable is connected with the $\sigma$ -Painlevé III $'$ equation for complex parameters. Furthermore, we investigate joint moments involving higher-order derivatives of characteristic polynomials in the Circular Jacobi Ensemble. As a consequence, we extend a result of Assiotis et al. in [Comm. Pure Appl. Math. 79 (2026), no. 7, 1771-1827, Theorem 1.11] concerning the joint moments of a sequence of random variables arising from the ergodic decomposition of Hua-Pickrell measures, from real parameters to complex parameters.

The Oldest Allele in a Moran Model

2026-08-24T12:22:14-04:00

We study an infinitely many neutral alleles Moran model where the allelic frequencies are recorded in increasing order of the ages of the alleles. Using this representation, we show that in stationarity the number of individuals with the oldest allelic type evolves as a Markov process, even across times when the oldest allele goes extinct (and is immediately replaced by the formerly second oldest allele). We also find the diffusive limit of the oldest allele chain.

Disc counting statistics of the real Ginibre ensemble

2026-08-24T08:04:24-04:00

We study the disc counting statistics of the real Ginibre ensemble, whose spectrum consists of real and non-real eigenvalues, the latter occurring in complex-conjugate pairs. As the matrix dimension increases, we derive the asymptotic behaviour of the joint cumulants of the numbers of real and non-real eigenvalues contained in a centred disc in the bulk, edge, and exterior regimes of the circular law. As consequences, we establish a joint central limit theorem and obtain the conjectured asymptotics for the number variance. Furthermore, our results reveal a universality phenomenon: although the fluctuations of the real and non-real eigenvalues separately retain dependence on the symmetry class, their combined fluctuations exhibit the same limiting behaviour as in the previously established complex and symplectic Ginibre ensembles.

Strong asymptotics for Jacobi-Piñeiro orthogonal polynomials

2026-08-20T18:13:58-04:00

We investigate the asymptotic behavior of Jacobi-Piñeiro polynomials of degree $2n$ orthogonal on $[0,1]$ with respect to weights $w_j(x) = x^{\alpha_j}(1-x)^\beta$ , $j=1,2$ where $\alpha_1,\alpha_2, \beta>-1$ , and $\alpha_1-\alpha_2 \notin \mathbb{Z}$ . These polynomials are characterized by a Riemann-Hilbert problem for a $3 \times 3$ matrix-valued function. We use the Deift-Zhou steepest descent method for Riemann-Hilbert problems to obtain strong uniform asymptotics in the complex plane. The local parametrix around the origin is constructed using Meijer G-functions. We match the local parametrix around the origin with the global parametrix with a double matching, a technique that was recently introduced.

One-point fluctuations for exponential last passage percolation under upper-tail conditioning

2026-08-20T16:32:12-04:00

We study exponential directed last passage percolation conditioned on the last passage time to a specified macroscopic point being atypically large. We determine the one-point fluctuations throughout two of the three spatial regions arising under this conditioning, as well as on the boundaries between these regions, extending the work of Baik-Cordaro-Tripathi. Depending on the location of the observation point, the limiting fluctuations are governed by the GUE Tracy-Widom distribution, a Gaussian distribution, or the one-spike BBP distribution on the boundaries between these regions. Through the correspondence with tandem queues, our results also describe how an atypically late departure at one customer-station pair affects departure epochs elsewhere in the network.

Notes on Hydrodynamic Limits and Related Topics

2026-08-20T12:47:37-04:00

In these lecture notes, we discuss various `hydrodynamic LLN' and `CLT' scaling limits, among others, in types of stochastic interacting particle systems, connecting `microscopic' behaviors to continuum laws. Via `short stories', the aim is to present some of the `basics' for students and those entering the field, as a complement to books such as Kipnis-Landim 1999, Komorowski-Landim-Olla 2012, Liggett 1985, Liggett 1999. To be concrete, attention is restricted to a few `mass conservative' systems on discrete spaces, namely exclusion and zero-range processes, that have proved robust in the study of different phenomena. After preliminaries, we discuss the `entropy' and `relative entropy' methods to prove hydrodynamic limits of the bulk mass in finite volume, as well as other items such as construction of systems in infinite volume and the structure of their invariant measures, and scaling limits of local functionals, such as occupation times of sites and the motion of a tagged particle. In the last part, we also discuss equilibrium fluctuations of the bulk mass when the process starts from an invariant measure.

Positive definite, positive semidefinite and totally positive matrices over finite fields

2026-08-18T08:25:08-04:00

Motivated by the equivalent definitions of positive definite (resp. positive semidefinite) matrices over real and complex fields, we give four (resp. five) inequivalent definitions for these matrices over finite fields. Our starting point is the recent definition due to Cooper--Hanna--Whitlatch (RMJ. Math., 2024) of positive elements in finite fields. We also use this definition to study totally positive matrices over finite fields. For all of these cases, we give explicit enumeration formulae or give bounds. Most of our formulas are new, but we summarize results from the existing literature for completeness. For positive semidefinite matrices of type 5 and totally positive matrices, we give structural formulas using the rationality of the Weil zeta function, i.e. Dwork's theorem, and conjecture a quasipolynomial-type formula.

Asymptotic analysis of a Fredholm determinant occurring in the description of the dynamical correlation functions of the Lieb--Liniger Bose gas

2026-08-17T11:41:16-04:00

We perform a large- $x$ asymptotic analysis of the Fredholm determinant of an integrable integral operator with generalized sine kernel, where $x$ controls the strength of the oscillations along the integration contour of the operator and will play the role of the distance variable in applications to the correlation functions of integrable quantum systems. Our generalized sine kernel involves a number of functional parameters that will allow us to adapt it to the analysis of the dynamical correlation functions of the Lieb--Liniger model at finite temperature and for all positive values of the coupling constant. It will also allow us to consider a class of equilibrium correlation functions that are governed by generalized Gibbs ensembles. Our work is based on the analysis of a matrix Riemann--Hilbert problem that is canonically connected with our integrable integral operator.

Gradient Gibbs measures with non-convex potentials and the universality class of the Gaussian Free Field

2026-08-14T13:45:42-04:00

We study a general class of gradient interface models with Hamiltonian $H=\beta \sum V(\nabla\varphi)$ , $\beta>0$ , assuming essentially that the potential $V$ is even, $V'(s)\ge \alpha s$ on $[0,\infty)$ for some $\alpha>0$ , and $-M\leq V''\le C$ . We establish a Helffer-Sjöstrand representation for these models, and use it to prove that their scaling limits are Gaussian Free Fields (GFFs), and that their covariances decay at the same rate as the GFF. This extends results for strictly convex potentials to a large class of non-convex potentials and to arbitrary temperatures. Additionally, we prove Brascamp-Lieb and dimension-free Poincaré inequalities for the models. We obtain these results by representing the interface as a mixture of gradient interface models with strictly convex potentials, extending an idea by Biskup-Spohn who had considered mixtures of Gaussians at moderate inverse temperature $\beta=1$ . The construction of such a representation is one of the key new contributions of this work.

Form factors of field exponentials in the 1+1 dimensional Sine-Gordon model

2026-08-14T11:04:39-04:00

This work provides explicit expressions, in all sectors of the Fock space: soliton, antisoliton and breather, for the fundamental building blocks -the form factors- of the generalised integral kernels which represent the exponential of the field operator in the 1+1 dimensional Sine-Gordon quantum field theory. The obtained form factors reproduce the various results obtained earlier in the literature using, in particular, the free field approach.

Hitting-time mixing for the star transposition shuffle

2026-08-13T15:41:38-04:00

We prove a hitting-time analogue of cutoff for the star transposition shuffle on the symmetric group S_n. Let tau be the first time at which every non-top card has been selected. We show that the shuffle is asymptotically mixed at time tau: more precisely, the total variation distance between the law of Y_tau and the uniform distribution on S_n is at most exp(-(log n)^(1/2+o(1))). Our proof compares the star transposition shuffle with the random transposition shuffle using simultaneous diagonalization of the two transition kernels, and then adapts the hitting-time strategy of Jain and Sawhney. This introduces a technique that can be applied to card shuffles that are not necessarily conjugacy invariant.

Fractional moments of the Stochastic Heat Flow and 2D Directed Polymers

2026-08-13T11:22:59-04:00

We estimate the fractional moments of the normalized mass assigned by the Critical 2D Stochastic Heat Flow to small balls. Our results also cover the discrete case corresponding to the 2D directed polymer model and provide estimates that are uniform in all parameters. One key takeaway of our results is that the vanishing of the fractional moments is completely governed by the divergence of the second moment. We use a quite robust method, by refining the change of measure argument and introducing a novel coarse-graining procedure, reducing the proof to essentially second moment estimates (in fact, we also provide sharp second moment estimates for directed polymers, of independent interest).

On Toeplitz determinants with slow Fourier decay

2026-08-13T08:46:25-04:00

We study Toeplitz determinants $\det T_n(e^f)$ for $f$ whose Fourier coefficients satisfy $f_k=O(|k|^{-1})$ . This regime extends beyond $H^{1/2}$ and includes symbols with Fisher-Hartwig singularities. We develop an operator-theoretic approach based on the Baker-Campbell-Hausdorff formula that separates the quadratic term \[ \sum_{k=1}^{\infty}\min(k,n)f_kf_{-k} \] from the higher-order terms in the expansion of $\log\det T_n(e^{tf})$ . We show that this quadratic term accounts for the possible growth with $n$ , while every fixed higher-order coefficient remains bounded. For symbols with bounded positive and negative Fourier parts, our estimates yield two-sided bounds for the determinant after removal of the quadratic contribution. For a broader admissible class, including Fisher-Hartwig-type symbols, we obtain uniform higher-order coefficient bounds and a central limit theorem for the associated CUE linear statistics. We also obtain bounds on mixed exponential moments for CUE-derived random fields beyond the characteristic polynomial.

Logarithmic intermittency of the critical 2D SHF

2026-08-12T13:10:03-04:00

While the solution to the $1+1$ dimensional stochastic heat equation with multiplicative noise is closely related to the exponential of a Brownian motion, the two-dimensional picture exhibits an additional weak-to-strong disorder transition. In [CSZ '23], the critical two-dimensional stochastic heat flow (SHF) was constructed as the scaling limit of the partition function of $2+1$ dimensional directed polymers under the logarithmic intermediate-disorder scaling at criticality. The SHF is a random measure and, like many naturally occurring random measures, it is expected to exhibit rich intermittency. [CSZ '25] established that it is almost surely singular with respect to the Lebesgue measure. More recently, [GT '26] showed that the logarithm of the SHF averaged over small balls is asymptotically Gaussian, with both its mean and variance diverging as the ball radius tends to zero. In this paper we prove a sharp result quantifying the singularity of the support of the SHF as well as its intermittency. In particular, we show that, almost surely, for all small $\varepsilon>0$ , up to a vanishing error, all the mass of the point-to-plane SHF in any domain is concentrated on ${1}/{\big(\varepsilon^2\log^{1/2+o(1)}(1/\varepsilon)\big)}$ balls of radius $\varepsilon$ , each containing $\varepsilon^2{\log^{1/2+o(1)}(1/\varepsilon)}$ mass, thus precisely establishing its logarithmic fractal behavior. A key ingredient in the proof is a refined large-deviations theory, which allows access to conditional distributions, by taking advantage of the Gaussian-like behavior of the SHF at quasi-critical scales. A further useful observation that features prominently is that conditioning a Brownian motion on its endpoint being unusually large essentially induces a shift in the mean of its increments, and consequently, at small enough scales, their distributions do not alter significantly.

Normalized skew Schur polynomials are Lorentzian

2026-08-12T13:06:13-04:00

We prove the conjecture of Huh, Matherne, Mészáros, and St.~Dizier that the normalization of every skew Schur polynomial in finitely many variables is Lorentzian. We first realize every nonzero skew Schur polynomial in finitely many variables as a specialization of a Schubert polynomial and prove that it is dually Lorentzian. The dual Jacobi--Trudi identity then identifies its normalization with the finite dual of a skew Schur polynomial obtained by rectangular complementation. As a consequence, skew Kostka numbers satisfy log-concavity inequalities along the root directions.

Continuous-time directed polymers to the Critical SHF via moments

2026-08-11T23:23:01-04:00

We consider a class of two-dimensional continuous-time directed random polymers in a correlated Brownian environment. The partition functions of these polymers form a discrete approximation of the Critical Stochastic Heat Flow (SHF), constructed by Caravenna, Sun and Zygouras (2023). Using a recent moment based axiomatic characterization of the SHF by Tsai (2024), we prove that the point-to-point partition functions, viewed as random measures, converge to the critical SHF. The proof proceeds by first establishing convergence of all moments through the resolvent method introduced by Gu, Quastel and Tsai (2021) and based on Rajeev (1999) and Dimock and Rajeev (2004). This is followed by establishing tightness and verifying the remaining axioms in the characterization.

On the most likely geodesic in last passage percolation

2026-08-11T15:34:58-04:00

We consider the problem of identifying the paths most likely to occur as geodesics in last passage percolation. Heuristics suggest that these modal paths should be the extreme corner paths, going straight between the corners of the cube of accessible vertices. We identify three mechanisms which favour such corner paths, and show that they are more likely to appear than all but a vanishingly small portion of paths. We make more specific comparisons in exponential last passage percolation, where moderate deviation estimates may be used to show that corner paths are nearly modal in a precise sense. Finally, we show a form of monotonicity of geodesic probabilities in a special case and conjecture that this holds in general.

Matrix Szegő Function and Matrix Orthogonal Polynomials for Multiple Cuts

2026-08-11T11:16:45-04:00

Following recent developments on matrix-valued orthogonal polynomials (MVOPs), we construct the matrix Szegő factorisation of the matrix weight function when it is supported on multiple intervals. We find that on the gaps the matrix function has unitary multiplicative jumps. In the next part, performing the Deift-Zhou steepest descent analysis we study the associated global parametrix which comes from the Riemann-Hilbert problem of the MVOPs. A solution is constructed visualising the rows of the parametrix as sections of vector bundles on a hyper-elliptic Riemann surface. An alternate view point to the global parametrix is also presented using a vanishing lemma. As an application we obtain strong asymptotics of MVOPs for multiple cuts which extends the works of Deaño, Kuijlaars, and Román.

A proof of the m-Symmetric Macdonald positivity at t=1

2026-08-10T16:47:45-04:00

We prove, in the case $t=1$ , the extension to the $m$ -symmetric world of the original Macdonald positivity conjecture. This is achieved by giving a combinatorial interpretation of the Kostka coefficients $K_{\Omega \Lambda}(q,1)$ in terms of standard fillings of the diagram associated to the $m$ -partition $\Omega$ . This interpretation generalizes the one in the usual Macdonald case, which is given by a major index statistic on standard tableaux.

Rooted Gibbs-DLR Measures in Planar Directed Polymers

2026-08-10T13:33:42-04:00

We study rooted Gibbs-DLR measures in the directed polymer model on $\mathbb{Z}^2$ with an ergodic disorder distribution which satisfies an additional mild hypothesis. We prove that the set of extremal rooted Gibbs-DLR measures is closed and totally ordered, characterize extremality in terms of path coalescence, and show that each fully supported extremal rooted Gibbs measure canonically generates a globally consistent and coalescing family of extremal rooted Gibbs measures indexed by all lattice sites. These families are, moreover, totally ordered. Building on this structure, we prove strong existence and strong uniqueness of the associated Busemann process, together with an $L^1$ continuity theorem for the shift-covariant Busemann cocycles which is joint in the inverse temperature, the tilt parameter, and the random environment. This yields, as a corollary, in-probability continuity of the generated extremal Gibbs measures corresponding to directions of differentiability under bounded i.i.d. perturbations of the weights. In positive temperature, it shows in-probability convergence of the generated extremal Gibbs measures. At zero temperature, it also yields quenched subsequential large deviation principles for the corresponding positive-temperature rooted Gibbs-DLR measures on path space. The rate functions are determined by a zero-temperature Busemann cocycle and vanish precisely on the infinite geodesics generated by that cocycle.

Limiting eigen-structure of spiked sample covariance matrices under missing observations

2026-08-10T01:27:53-04:00

High-dimensional Principal Component Analysis (PCA) has become an essential tool in modern data analysis, offering dimensionality reduction and feature extraction. However, the presence of missing data introduces significant challenges, distorting the performance of PCA and complicating statistical inference. In this paper, we study the asymptotic behavior of PCA under a spiked population model with missing observations, leveraging recent advances in random matrix theory. We demonstrate that while the spiked sample eigenvalues exhibit asymptotic normality, the limiting parameters differ substantially from those in the complete data case, reflecting the non?trivial influence of the missing data mechanism. As an application of our results, we propose a test to evaluate the independent structure of a spiked population.

q-Difference Equations for two classes of Pattern Avoiding Permutations

2026-08-09T23:27:07-04:00

We study the nested permutation classes $Av(4123,4231,4312)\subset Av(4123,4312)$ . We relate the corresponding generating functions $D^{(2)}(z)$ , $D^{(3)}(z)$ to $q$ -difference equations. Solving these equations allows to show that the generation functions are not D-finite, and that $[z^n]D^{(3)}(z)\sim 0.0189672325071988\ldots\, (4.46840899160393\ldots)^n$ , $[z^n]D^{(2)}(z)\sim0.0849833632890843\ldots\, (3+2\sqrt2)^n n^{-3/2}$ .

Generalized Freud weight, discrete Painlevé I hierarchy and full asymptotics of Hankel determinants

2026-08-08T23:52:54-04:00

In this paper, we investigate the monic orthogonal polynomials $P_{n}(x;T_{m};\lambda)$ and the Hankel determinants $D_{n}(T_{m}; \lambda)$ associated with the generalized Freud weight \[w(x;T_{m};λ) = |x|^{2λ+1}\exp\biggl(-\sum_{k=1}^m t_k x^{2k}\biggr),\quad m \in \mathbb{Z}^+,\; t_{k} \in \mathbb{R} , x\in\mathbb{R}\setminus\{0\},\] where $T_{m}=\{t_{1},t_{2},\cdots, t_{m}\}$, $t_{m}>0$ and $λ>-1$.By employing ladder operators and compatibility conditions, we find that all members of the discrete Painlevé I hierarchy have a unified structure and the recurrence coefficient $β_n$ of $P_{n}(x;T_{m};\lambda)$ satisfies the $m$ -th member of the discrete Painlevé I hierarchy. Besides, we derive the second-order differential equation satisfied by $P_{n}(x;T_{m};\lambda)$ , the partial derivatives of the recurrence coefficients $β_n$ with respect to parameters $t_1, t_2, \dots, t_{m-1}$ and the corresponding differential identities for $D_{n}(T_{m}; \lambda)$ . Based on the discrete Painlevé I hierarchy and the above differential identities, we obtain new partial differential equations satisfied by $\ln\beta_n$ and $\ln D_{n}(T_{m}; \lambda)$ .Using the discrete Painlevé I hierarchy and the asymptotic theory of linear difference equations, we derive the full asymptotic expansions of the recurrence coefficient $\beta_n$ , the nontrivial leading coefficient $\mathrm{p}(n; T_m; \lambda)$ , and the Hankel determinant $D_n(T_m; \lambda)$ as $n\to\infty$ , for general $T_m$ and $\lambda>-1$ . Notably, while the logarithmic term $\ln n$ appears in the leading-order contributions, it is absent from the remainder terms in these expansions.We illustrate our results under the specific decic Freud weight $w(x;t_1,t_2;\lambda)=|x|^{2\lambda+1} \exp\bigl(-x^{10}-t_2x^4-t_1x^2\bigr)$ .

Universal Hilbert series coefficients of the superspace coinvariant ring

2026-08-08T11:28:08-04:00

The coefficients that determine the Hilbert series of the superspace coinvariant ring are indexed by hook-shaped partitions. We give a manifestly positive combinatorial interpretation of these coefficients, together with several generating functions for them. Specializing this Hilbert series at $u=-q^2$ , we show that its coefficients are differences of binomial coefficients. Consequently, this proves a conjecture of Sagan--Swanson (2024) that these coefficients are palindromic up to sign. More generally, for every $m \geq 1$ we obtain closed-form expressions for the $u = -q^m$ specialization.

Negative association of Busemann functions in exponential last-passage percolation

2026-08-07T09:52:03-04:00

One hallmark of exactly solvable KPZ random growth models is product-form invariant measures. In the setting of exponential last-passage percolation (LPP), this corresponds to the independence of Busemann increments along any down-right path. However, this independence breaks down when multiple asymptotic directions are considered simultaneously, owing to the fact that jointly invariant measures are not jointly product-form. This paper shows that the failure of independence is one-sided: Busemann increments across arbitrary directions are negatively associated. As an application, we derive an exponential concentration inequality for sums of Busemann increments on the diffusive scale, even when the increments are not independent. While our argument relies on a Burke property that is special to exponential weights, all other proof ingredients$\unicode{x2014}$including hidden LPP monotonicities and braid relations for queueing maps$\unicode{x2014}$hold for arbitrary weights.

Universality for Random p-adic Polynomials and Random Matrices Via the Resultant Distribution Method

2026-08-06T16:38:06-04:00

We prove universality results for zeros of random $p$ -adic polynomials with independent, sufficiently non-concentrated coefficients. We show that when the degree goes to infinity, the limiting joint root statistics in arbitrary finite extensions of $\mathbb{Q}_p$ are universal and coincide with those of the Haar coefficient model, whose limiting root statistics were determined by Caruso (arXiv:2110.03942). In particular, our results extend the Haar-model formulas to a broad class of coefficient distributions and to general joint root statistics over finite extensions of $\mathbb{Q}_p$ . Our approach is based on a new method, which we call the resultant distribution method. Instead of studying the roots of a random polynomial directly, we analyze the distributions of the resultant valuations against suitable fixed test polynomials. We prove that these resultant distributions determine the limiting law of the polynomial and, together with suitable degree estimates, yield convergence of general root statistics. As a further application, we apply the same method to characteristic polynomials of random $p$ -adic matrices with independent entries. We prove that the distinguished factors of the characteristic polynomials converge to universal limiting distributions independent of the entry distribution, thereby showing that such a random matrix model is remarkably robust.

Non-Gaussian fluctuations for traces of squared sample correlation matrices in high dimensions

2026-08-06T06:49:45-04:00

We provide limit theory for the trace of the squared sample correlation matrix $\mathbf R$ , constructed from $n$ observations of a $p$ -dimensional random vector with iid components. If the entries have finite fourth moment and $p$ and $n$ grow proportionally, it is known that $\operatorname{tr}({\mathbf R}^2)$ satisfies a central limit theorem (CLT) and the centering and scaling sequences are universal in the sense that they do not depend on the entry distribution. Under symmetry and regular variation assumption with index $\alpha$ and any growth rate of the dimension, we prove that the universal CLT remains valid for $\alpha>3$ . For $\alpha<3$ , we identify a critical dimension growth at which the fluctuations of $\operatorname{tr}({\mathbf R}^2)$ become non-Gaussian. Moreover, if the dimension $p$ grows faster and $\alpha \le 3$ we establish a non-universal CLT with norming sequences depending on the value of $\alpha$ . Our findings are illustrated in a simulation study.

Asymptotics of Polynomials Orthogonal on an Interval with Varying Weights

2026-08-05T09:41:42-04:00

We study strong asymptotics of orthogonal polynomials on an interval with varying weights, extending the framework of Totik's theorem on weighted orthogonal polynomials. Our results generalize this theorem in several directions. In particular, we allow the supports of the associated equilibrium measures to converge to a proper subinterval of the original interval of orthogonality or even collapse to a point. These extensions are motivated by applications to multiple orthogonality where such varying measures arise naturally.

A Spectral Route to Directed-Polymer Glasses

2026-08-04T10:26:28-04:00

A finite density of mutually avoiding directed polymers in a quenched random medium is a minimal model of glassy line matter. The dilute theory, solved by replica Bethe ansatz, predicts an interaction free energy proportional to $\rho^2$ and disorder cumulants with distinct power-law dependences on the density $\rho$ , but direct numerical tests have been hindered by the combinatorially large many-polymer transfer matrix. We recast the problem as filling logarithmic eigenvalues of a single-polymer transfer-matrix product, obtaining the quenched free energy, its cumulants, and a disorder-induced linear spectral edge consistent with the replica prediction.

On the largest Littlewood--Richardson coefficient

2026-08-04T03:56:15-04:00

We study partitions which attain the largest Littlewood-Richardson coefficient. More precisely, we prove that the largest $c^\lambda_{\mu \nu}$ is attained at partitions such that $\mu \subseteq \nu$ and $\nu/\mu$ is a disjoint union of squares. We conjecture that for $n \ge 16$ , all largest $c^\lambda_{\mu \nu}$ must satisfy this property. We confirm this conjecture numerically, for $16 \le n \le 45$ .

Formulas for Koornwinder polynomials

2026-08-03T15:05:42-04:00

This paper provides formulas for Koornwinder polynomials in analogy with the creation formula, the alcove walk formula and the non-attacking fillings formula for the type $GL_n$ Macdonald polynomials. We state the creation formula in terms of the divided-difference operators used in Schubert calculus, and we use a box-greedy reduced word to reformulate the alcove walk formula in terms of uncompressed set-valued tableaux. Then two types of compression, ``around-the-end compression'' and ``across-the- $0$ -gap compression'', are used to derive a formula for Koornwinder polynomials in terms of compressed set-valued tableaux. Throughout we work in the full generality of relative Koornwinder polynomials, which are the analogues of the permuted basement Macdonald polynomials used in the type $GL_n$ case.

Multilevel limits of spiked random matrix minors

2026-08-02T18:51:16-04:00

We consider the largest eigenvalues of minors of the classical Gaussian random matrices with a finite rank spike in the critical BBP regime. We show that they converge to a multi-level interlacing particle system whose lowest level is the Airy $_\beta$ point process. Our construction and scaling limits extend in a natural way to non-classical $\beta>0$ .

Poisson approximations of the number of fixed points in random multiset permutations

2026-08-01T11:16:14-04:00

For ordinary permutations on $n$ letters, the distribution of the number of fixed points of random permutations is well known to approach the Poisson $(1)$ distribution in total variation distance as $n\to\infty$ super-exponentially quickly. We use Stein's method to get related results for the number of fixed points of random permutations of multisets. Given a sequence of multisets on $n$ letters whose expected number of fixed points converges to a constant $c$ , we must have $c\geq 1$ and the distribution of number of fixed points converges to the ${\rm Poisson}(c)$ distribution as $n\to\infty$ . If $c>1$ , then the rate of convergence in total variation distance can be as slow as $n^{-1/2}$ .

Interacting Quantum Symmetric Exclusion Process

2026-07-30T10:15:35-04:00

We introduce and solve the Interacting Quantum Symmetric Exclusion Process (IQSEP), a family of models describing the stochastic quantum hopping of charged particles along the edges of a lattice, with hopping amplitudes that depend on the occupations of neighbouring sites. In the absence of interactions, they reduce to the standard quantum simple symmetric exclusion process, exhibiting coherent diffusive transport. For interactions of order one, they capture incoherent diffusive transport and its fluctuations, characterized by density-dependent diffusivity and mobility, making contact with the macroscopic fluctuation theory. By rescaling the interaction strength appropriately with the lattice mesh, we define a mesoscopic scaling regime that retains a finite coherence length in the continuous thermodynamic limit. This regime interpolates between coherent behavior at small length scales and incoherent behavior at large scales. The resulting scaling theory accounts for fluctuations of quantum coherences in interacting diffusive systems, going beyond the scope of standard fluctuating hydrodynamics.

Three-Dimensional Kardar--Parisi--Zhang Scaling in Polariton Condensates

2026-07-30T08:13:36-04:00

Kardar--Parisi--Zhang (KPZ) universality provides an example of macroscopic scaling generated by microscopic violation of detailed balance. While one- and two-dimensional realizations have been explored in driven condensates and growing interfaces, demonstrating KPZ scaling in three spatial dimensions remains a major challenge. Here we propose a three-dimensional exciton-polariton crystal as a platform for observation of 3D KPZ universality. Starting from a stochastic driven-dissipative Gross-Pitaevskii equation for a condensate formed in a three-dimensional photonic-crystal lower-polariton band, we eliminate the massive density and reservoir modes and obtain an effective $3+1$ -dimensional KPZ equation for the condensate phase. Numerical simulations of both the KPZ equation and the full driven-dissipative polariton model show an intermediate-asymptotic regime in which the first-order coherence obeys $-\ln |\gone(0,\Delta t)|\propto |\Delta t|^{2\beta}$ and $-\ln |\gone(\Delta r,0)|\propto |\Delta r|^{2\chi}$, with exponents consistent with the $3+1$ KPZ benchmarks $\beta= 0.1845$ , $\chi= 0.3135$ . Our results identify three-dimensional polariton crystals as a controllable quantum fluid route to higher-dimensional nonequilibrium universality.

Perturbed Beta Corners Process

2026-07-30T00:00:00-04:00

We introduce and study the perturbed $\beta$ -corners process, a deformation of the classical $\beta$ -corners process. The latter is a probability measure on interlacing arrays of real numbers that, for the classical values $\beta=1,2,4$ , describes the joint distribution of eigenvalues of principal submatrices of a Gaussian random matrix with the corresponding GOE/GUE/GSE symmetry. For classical $\beta$ , the perturbed process arises from adding a deterministic diagonal matrix $A=diag(a_1,...,a_N)$ to a GOE/GUE/GSE matrix. The eigenvalues of the resulting random matrix depend symmetrically on $a_1,...,a_N$ , but this symmetry does not extend to the whole corners process. We extend the construction to all $\beta>0$ via multivariate Bessel functions, and analyze the crystallization ( $\beta \to\infty$ ) of the resulting interlacing array, proving a law of large numbers and a central limit theorem in two regimes. For fixed perturbation, the eigenvalue repulsion dominates, and the array crystallizes on the same lattice of polynomial-derivative roots, with the same discrete Gaussian free field fluctuations, as in the unperturbed case treated by Gorin and Marcus (arXiv:1706.07393). New phenomena appear in the second regime, where the $a_i$ 's grow linearly in $\beta$ . Then the external source competes with the repulsion at leading order. Here the array freezes on a deformed lattice characterized by a coupled system of optimality equations. The fluctuations are governed by the same discrete Gaussian free field, now attached to the deformed lattice. The deformed lattice equations decouple in the special case of a single spike in the last coordinate. Then the deformed lattice is obtained explicitly by applying one shifted derivative $D_c f = f' + cf$ followed by iterated ordinary derivatives.

Free versions of the strong Szegő limit theorem

2026-07-28T13:01:35-04:00

The Strong Szegő Limit Theorem is a theorem about the asymptotics of the determinants of large Toeplitz matrices. It can be reformulated as a probabilistic statement about eigenvalue statistics of random unitary matrices. We prove a multivariate generalization of the theorem in this latter form, replacing a single unitary with a system of independent random unitaries. It turns out that the standard proofs of the classical theorem do not generalize to this setting, and instead we must import tools from random matrix theory, free probability, and noncommutative function theory. In addition, we obtain results about the averaged determinants of random unitary pencils, as the size of the unitaries tends to infinity; and an auxiliary result giving a formula for the spectral radius of a matricial sum of free Haar unitaries.

Intertwiners of representations of exceptional type quantum affine superalgebras

2026-07-27T08:55:36-04:00

We give explicit formulas for the smallest non-trivial irreducible representation $V$ of quantum affine superalgebras in types D $_{2\vert 1;\alpha}$ , $\dim V=18$ (in the all-fermionic parity), F $_{3\vert 1}$ , $\dim V=41$ (in the distinguished parity), and G $_{2\vert 1}$ , $\dim V=32$ (in the distinguished parity), both in the Drinfeld-Jimbo and in the new Drinfeld realizations. We use this information to obtain an explicit expression for the corresponding $R$ -matrices.

On non-monotonicity of logarithmic energy for random matrices

2026-07-27T04:53:40-04:00

We construct a finite-energy Wigner counterexample and a concrete one-parameter family of Gaussian-regularised Bernoulli entry laws yielding counterexamples to the conjectural dimensional monotonicity of the quadratically penalised logarithmic energy for mean empirical spectral distributions by Chafaï, Dadoun, and Youssef.

Jeu de taquin forests and the inverse infinite RSK correspondence

2026-07-27T01:33:37-04:00

The Plancherel-random infinite Young tableau arises from applying the infinite Robinson-Schensted-Knuth (RSK) correspondence to a sequence of i.i.d. random variables distributed uniformly on the unit interval [0,1]. Building on the isomorphism of dynamical systems established in prior work, we provide an explicit geometric characterization of the inverse map. Our approach makes use of a previously unexplored structure, the jeu de taquin forest on the infinite tableau, where edges connect boxes according to local comparison rules. We prove that each tree in this forest almost surely extends toward infinity with a well-defined asymptotic direction, establishing a canonical bijection between trees and values in the i.i.d. input sequence. The ordering is recovered through tree lifetimes under iterated jeu de taquin transformations.

Type C multiline queues and the open-boundary TASEP

2026-07-27T00:00:00-04:00

The totally asymmetric simple exclusion process (TASEP) with open boundaries is a finite Markov chain describing particles hopping between adjacent sites on a one-dimensional lattice with left and right boundary transitions governed by parameters $\alpha$ and $\beta$ . The multispecies TASEP is a higher-rank generalization in which particles have different species. Multiline queues were introduced by Ferrari and Martin (2007) to compute the stationary distribution of the multispecies TASEP on a circle. It has remained an open problem to find a combinatorial formula for the stationary distribution of the multispecies open-boundary TASEP. Using Kirillov--Reshetikhin crystals of type $C$ , we construct type $C$ multiline queues and a corresponding Ferrari--Martin pairing algorithm that projects them to TASEP configurations. This yields a combinatorial formula for the stationary distribution of the multispecies open-boundary TASEP for the $\alpha=\beta=1$ specialization.

Non-equilibrium fluctuations of a two-species exclusion process with slow boundary

2026-07-26T09:10:23-04:00

We study the non-equilibrium density fluctuations of a one-dimensional two-species symmetric simple exclusion process in contact with critical slow boundary reservoirs. The boundary rates are of order $1/n$ , which leads to Robin boundary conditions at the macroscopic level. We prove that the coupled fluctuation field associated with the empirical densities of two species converges to a generalized Ornstein-Uhlenbeck process. The limiting process is characterized by a linear martingale problem whose coefficients reflect the combined effects of bulk diffusion, species conversion, and the boundary reservoirs.

Bilinear generating functions of the multivariable Al-Salam-Carlitz polynomials and applications

2026-07-25T04:28:59-04:00

In this paper, by the method of comparing coefficients, we establish a new generating function of multivariable Al-Salam-Carlitz polynomials which contains both Rogers' and Bowman's symmetirc expansion formulas and the classical $q$ -Mehler formula as special cases. Some new $q$ -series identities related to the multivariable Al-Salam-Carlitz polynomials are also presented.

Mixing times and spectra of non-equilibrium symmetric exclusion processes on general graphs

2026-07-24T22:08:16-04:00

The symmetric exclusion process (SEP) is a classical model of interacting particles on a graph, in which multiple particles execute random walks subject to the constraint that no two particles may occupy the same vertex. The version of the process in which the number of particles is fixed, is by now very well understood, through the study of reversible Markov chains. In this paper, we study the non-equilibrium version of the SEP, in which some vertices of the graph interact with external heat baths, held at different temperatures, that mediate transport through the graph. The resulting steady-state distribution is still the stationary distribution of a natural Markov chain, but the chain is no longer reversible. This means that even the steady-state distribution is very hard to describe, and indeed is the subject of a rich literature, which has focused almost exclusively on graphs which are finite subsets of the lattice $Z^d$ . Very little is known about other important quantities, notably the mixing time. We derive various quantitative results about the non-equilibrium SEP for arbitrary graphs with arbitrary heat baths. Our main result is a bound on the mixing time in terms of the worst-case hitting time of a single particle to any of the heat baths; this bound is tight up to two logarithmic factors in the size of the graph. We then use this result to derive a simple algorithm, running in time roughly $n^{O(k)}$ , that computes the steady-state joint occupation distribution for any set of $k$ vertices; since the steady-state distribution is very elusive, this provides a potentially useful tool to study its $k$ -wise correlations for small values of $k$ . Finally, we prove that the spectrum of the Markov chain associated with the non-reversible dynamics for the SEP is independent of the heat bath temperatures; thus in particular the spectrum is always real.

Effective field theories of nonlinear fluctuating hydrodynamics in one dimension

2026-07-24T13:57:48-04:00

Describing emergent macroscopic phenomena in low spatial dimensions is known to be notoriously challenging, primarily due to strong interactions that render perturbative approaches inapplicable. On the other hand, low-dimensional systems host a wealth of unorthodox phenomena. A prominent example is the emergence of superdiffusive transport in one-dimensional interacting systems, traditionally studied in the framework of nonlinear fluctuating hydrodynamics. After identifying and discussing internal inconsistencies in the previous formulations, in this work we develop a general and systematic approach for constructing effective field theories of one-dimensional hydrodynamic systems in the form of coupled stochastic Langevin-type equations compatible with the physical requirements of local equilibrium states such as the thermodynamic Maxwell relation and fluctuation-dissipation symmetry. We implemented a general numerical integration scheme and exemplified our construction on a simple model of two interacting hydrodynamic modes with a non-Gaussian stationary equilibrium measure.

Collision and non-collision for diffusions on configuration space

2026-07-24T10:12:33-04:00

We develop criteria for collision and non-collision of reversible infinitely many interacting diffusion processes in the real line. The approach is potential-theoretic and is based on capacity estimates for symmetric Dirichlet forms on the configuration space. Our main results are model-independent in the sense that no determinantal or Pfaffian structure, prescribed interaction potential, or explicit labelled stochastic differential equation is required. The non-collision criterion involves only the second and third correlation functions of the reversible measure, whereas the collision criterion is based on a local lower bound for a finite-volume conditional density of the reversible measure. As an application, we identify the sharp collision threshold for the diffusion associated with the $\mathsf{Sine}_\beta$ -symmetric Dirichlet form: the collision set is polar if and only if $\beta \ge 1$ . This provides an infinite-particle counterpart of the classical collision threshold for the finite-particle Dyson Brownian motion with inverse temperature $\beta$ .

The critical KPZ scale for the Averaging Process

2026-07-23T11:50:08-04:00

KPZ-type extremal fluctuations have recently been proved for several models of random walks in space-time random environments (RWRE) in $1+1$ dimensions. A general moment criterion predicts the spatial scale at which this behavior should occur, but does not by itself guarantee non-trivial fluctuations at that scale. In this paper, we show that the averaging process provides an instance in which this criterion is not sharp: because of a degeneracy in the update mechanism, the actual KPZ limit appears only beyond the predicted scale. The relevant critical contribution arises from the interplay of two distinct fluctuation mechanisms, a phenomenon that appears to be rather special within the RWRE setting. Our proof builds on Dobrushin-type local-time limit theorems for zero-sum additive functionals, refined estimates for tilted $k$ -point motions, and the recent moment-based axiomatic characterization of Cole--Hopf solutions to the one-dimensional KPZ equation. We also identify the behavior on the two sides of the critical scale, thereby sharply separating the subcritical, critical, and supercritical regimes of the model.

Random walks in Dirichlet random environment in dimension d+1

2026-07-22T11:27:25-04:00

The atypical behaviour of random walks in time-dependent random environment was recently related to Kardar-Parisi-Zhang (KPZ) growth. While this is now well-understood in spatial dimension $d=1$ , further efforts are necessary to better understand these connections in dimensions $d>1$ . In this paper, we study this problem numerically for $d=1, 2$ and $3$ , focusing on a discrete model with Dirichlet distributed transition probabilities. This model is a generalization of an integrable model in $d=1$ , and it has the advantage of admitting an explicit, product-form, stationary measure. We verify that the growth of the variance of the logarithm of point-to-point probabilities, namely from the origin to position $x$ in time $t$ , is compatible with KPZ growth in dimension $d=1$ and $d=2$ . In spatial dimension $d=3$ , we confirm the existence of a phase transition as the angle $\vert x\vert /t$ increases and we obtain a lower bound based on an exact second moment calculation. We find that in the weak disorder phase the point-to-point probability acquires a heavy tailed distribution, and that in the strong disorder phase the cumulants of its logarithm grow with time. Further, we show that for this model, we can compute exactly the sample to sample variance of the thermal average $\overline{ \langle x \rangle^2}$ and that it is related to the extreme diffusion coefficient introduced recently.

Combinatorial geometry of the 2D Toda lattice and Davey Stewartson equation

2026-07-22T09:10:20-04:00

The KP equation is a prototypical $(2+1)$ -dimensional integrable PDE. Its soliton solutions are famously parametrized by the Sato Grassmannian. In seminal work, Kodama and Williams made the surprising discovery that the combinatorics of soliton solutions are intimately related to the combinatorics of the totally positive Grassmannian as pioneered by Postnikov. They introduced novel algorithmic methods inspired by polyhedral structures arising from tropical geometry. Soliton solutions to the 2D Toda lattice and the Davey--Stewartson equation, two closely related integrable systems with soliton solutions, are also classified by the Sato Grassmannian. Kodama suggested that the methods of his work with Williams could generalize to these two integrable equations. In this work, we show that this is indeed the case. We derive algorithms to produce contour plots from elements in the totally nonnegative Grassmannian in both cases. In the asymptotic setting, we recover and refine previous work of Biondini and Wang; as well as Biondini, Kireyev and Maruno.

Recurrence Coefficients of the Orthogonal Polynomials for Oscillatory Jacobi-type Weight Functions

2026-07-22T02:21:21-04:00

We study two classes of oscillatory Jacobi-type weight functions $x^c(1-x^2)^{\lambda-1/2}\exp(i\zeta x)$ , $x\in[-1,1], \lambda>-1/2, c\in\{0,1\}$ . Other restrictions are imposed on $\lambda$ and $\zeta$ to guarantee the existence of the associated orthogonal polynomials. By using the ladder operators established in the recent literature for monic orthogonal polynomials associated with Jacobi-type weight functions and three compatibility conditions, we derive two coupled difference equations satisfied by the three-term recurrence coefficients. Compared with the existing results, these equations are structurally simpler and of lower order. Once the initial values are determined, the recurrence coefficients can be computed at any stage through the difference equations. The obtained expressions enable us to conjecture the symbolic forms for the recurrence coefficients.