Integrable Probability arXiv Feed
au:Borodin or au:"Alexei Borodin" — filter by authorin:JAMS or in:"Ann. Probab." — filter by journalcat:math.PR — filter by arXiv categoryy:2024 or y:2020-2024 — filter by yearTOP — show top-journal papers onlyClick any author, journal, or category badge to filter. Esc to clear.
Integrable probability papers on arXiv by members of the community (broadly understood), plus additional papers discovered by semantic search across the full arXiv database.
Not all papers by tracked authors are included — only those related to integrable probability.
The first version of this feed was imported from the NSF FRG page (opens in new tab); now updated via automated arXiv queries.
Please email me if you notice a missing author or any mistakes.
Subscribe via RSS
|
Download BibTeX of all papers (updated periodically)
October 2026
Abstract
We study sample covariance matrices formed from i.i.d., centered, isotropic random vectors , whose coordinates are allowed to be dependent. We assume the random vectors satisfy a concentration bound for their quadratic forms . Under this assumption, we establish edge universality and obtain the Tracy--Widom law for the largest eigenvalue of the sample covariance matrix after appropriate centering and scaling. In particular, our result removes the unconditionality assumption in existing edge universality results for log-concave ensembles~\cite{BX}, resolving the question that originally motivated our work. The concentration assumption also allows us to extend the result beyond the log-concave class.
Abstract
Using a rational analogue of an Rogers--Ramanujan type identity discovered by Warnaar, we establish sum-like expressions for a bivariate generating function for cylindric partitions of rank three. These expressions provide the finishing touch towards conjectures of Warnaar and Kanade--Russell.
Abstract
We evaluate the characteristic polynomial of a random unitary matrix at the points on the unit circle where the polynomial takes real values. In analogy to the Riemann zeta function we call these the ``random matrix Gram points''. On the random matrix side we answer in the affirmative a question of Richard Hall, who expected that zeta on the Gram points should be larger than typical. We prove that for any fixed positive , the th moment of the characteristic polynomial averaged over the random matrix Gram points is larger than the th moment of the characteristic polynomial averaged over the whole unit circle. For large the leading-order asymptotic term is the same size, but all lower-order terms are larger. In an appendix we show how to adapt the heuristic recipe of Conrey--Farmer--Keating--Rubinstein--Snaith to conjecture moments of zeta at Gram points, and show by example that the second and fourth moments calculated by following the recipe agree with previously proven results.
Abstract
We investigate the longest increasing subsequences (LIS) of random permutations distributed under the Luce model with weights . Our first main result provides a general formula for the growth rate of the LIS. For example, we get that the LIS of Luce permutations with Sukhatme weights grows like . We also give simple conditions for the LIS to be of order , for any given . Our second main result states that, if the weights are of the form for a sufficiently regular function , then the LIS divided by converges in probability to a constant described by a variational problem. Throughout, we make extensive use of coupling techniques; in particular, to prove our second main result, we argue that Luce random permutations behave similarly to random permutations sampled from the ``Luce permuton''.
Abstract
We determine the bigraded Frobenius series of the superspace coinvariant ring for the wreath products , which is the quotient of the ring of polynomial-valued differential forms by the ideal generated by diagonal -invariants with vanishing constant term. The bigraded multiplicity of each irreducible character is a sum over standard Young -tableaux of an expression determined by the flag major index and flag descent number. This refines the Frobenius series graded by fermionic degree found by Bhattacharya--Rhoades (2026). It specializes at to the hyperoctahedral group case, and it specializes at to recover the theorem of Murai--Rhoades--Wilson (2025) on the symmetric group . We give new formulas for the bigraded multiplicities of hook characters of the superspace coinvariant ring in types and , yielding upper bounds on trigraded multiplicities of the sign character of and the determinant character of . As a corollary, we deduce a conjecture of F. Bergeron (2024) on hook components of .
Abstract
Let be fixed, and let be positive integers satisfying . Let be an random combinatorial matrix whose rows are independent and uniformly distributed over the vectors in having exactly ones. For an arbitrary deterministic invertible matrix , we establish two-sided bounds for the smallest singular value of . More precisely, for every , \[ \mathbb P\!\left( \varepsilon\le s_n(Q_nM)\cdot \|M^{-1}\|_{\mathrm{HS}} \le \varepsilon^{-2} \right) \ge 1-C\left( \varepsilon+d^{-1/2}+e^{-c\,\operatorname{sr}(M^{-1})} \right)-2e^{-cn}, \] where positive constants depend only on and denotes the stable rank. Consequently, for any sequence of deterministic invertible matrices satisfying , is typically of order . In particular, taking recovers the known scale for .
Abstract
We investigate the role of asymmetry in the stationary measures of integrable probabilistic models with boundary interactions. This paper is the first in a series of three focusing on the open asymmetric simple exclusion process (ASEP). The stationary measure of open ASEP admits a two-layer Gibbs representation in which a correction term, which we call the Bryc polynomial, encodes the dependence on the asymmetry parameter. We give a probabilistic interpretation of the Bryc polynomial in terms of a birth-death-type Markov chain. As an application, we obtain sharp total variation approximations of the stationary measure by Bernoulli product measures in the high- and low-density phases of the shock region. These approximations yield further applications, including moderate deviation estimates for the stationary measure of open ASEP, as well as asymptotic results for the related model of open ASEP with light particles.
Abstract
We study the stationary measure of the open asymmetric simple exclusion process (ASEP) on the coexistence line. Our starting point is Bryc's two-layer Gibbs representation, together with the probabilistic interpretation, established in the first paper in this series, of the associated Bryc polynomial in terms of an inhomogeneous birth--death chain. Our main result gives a sharp approximation of the stationary measure. Specifically, we identify the shock location as an explicit observable of the particle configuration and obtain optimal bounds on the total variation distance between its law and the uniform measure on the lattice. Conditional on the shock location, we approximate, in total variation distance, the joint law of the occupation variables away from the shock by Bernoulli product measures with different densities to its left and right. These approximations lead to further applications, including fluctuation limits, moderate deviation bounds for the particle density, and limits for correlation functions under the stationary measure.
Abstract
We study the stationary measure of the open asymmetric simple exclusion process (ASEP) in the fan region with general boundary rates. The stationary height fluctuations in the maximal current phase and on its boundaries within the fan region are known to have limits involving Brownian excursions or meanders. We establish discrete counterparts of these limits: on intervals away from the boundaries, the joint law of the occupation variables is approximated in total variation distance by particle configurations constructed from lazy simple random walk excursions and meanders. In the high- and low-density phases of the fan region, we establish Bernoulli product approximations with errors that decay exponentially in the distance from the relevant boundary. Our methods combine refined combinatorial estimates for Motzkin paths and the matrix product ansatz in the totally asymmetric case with a Fourier-type decomposition (based on Askey--Wilson polynomials) of the change of measure for general asymmetry. In addition, in the maximal current phase and on its boundaries within the fan region, we prove moderate deviation estimates for particle counts on arbitrary sets of sites. These estimates rely on a covariance identity for an auxiliary path measure, through which total nonnegativity (of order two) of the matrices yields negative correlations.
Abstract
In this note, we give a short probabilistic proof of the two-layer Brownian representation of the stationary measure of the open Kardar--Parisi--Zhang (KPZ) equation on a finite interval with arbitrary Neumann boundary parameters, conjectured by Barraquand and Le Doussal and recently established in full generality by Dunlap, Gu and Rosati. Our starting point is Bryc's two-layer Gibbs representation of the stationary measure of the open asymmetric simple exclusion process (ASEP) under Liggett's condition, together with our previously established probabilistic interpretation of the associated Bryc polynomials in terms of inhomogeneous birth--death Markov chains. Martingale estimates yield uniform asymptotics for these polynomials under weakly asymmetric scaling. Combined with known convergence results from open ASEP to open KPZ, these estimates recover the stationary height profile of the open KPZ equation as the first marginal of a reweighted pair of Brownian motions.
Abstract
We describe two-parametric Littlewood identities for inhomogeneous stable spin Hall-Littlewood and inhomogeneous spin -Whittaker symmetric functions. These identities generalize the Littlewood identity found by Warnaar (2006) for the ordinary Hall-Littlewood polynomials. However, instead of Rogers-Szegö coefficients used by Warnaar, coefficients in our identities are given by products of Al-Salam-Chihara polynomials. Our proof relies on integrable lattice models. A key step includes three new solutions to the reflection and boundary equations for the fully fused six-vertex model. To find these solutions we use isomorphisms between tensor products of evaluation -modules, as well as covectors stable under a -Onsager coideal subalgebra.
Abstract
We prove a conjecture proposed by Chen that, for every fixed , the sequence is log-concave, where counts permutations of whose longest increasing subsequence has length .
Abstract
We establish operator universality at the spectral edge of general Wigner matrices by showing that their suitably rescaled Lanczos tridiagonalizations, constructed from deterministic initial vectors, converge to the stochastic Airy operator with Dirichlet boundary condition. The result requires control of the entire tridiagonal structure on the scale relevant to the limiting differential operator. The Lanczos coefficients for non-Gaussian Wigner matrices are strongly dependent. We analyze their fluctuations by combining multi-resolvent local laws with quantitative control of Krylov subspaces and a comparison of the resulting Jacobi coefficients with those of the Gaussian ensembles. This approach identifies the limiting Brownian potential directly from the original matrix entries, without invoking Dyson Brownian motion in the step of Green function comparison. The resolvent limit yields joint convergence of edge eigenvalues, embedded eigenfunctions, and spectral measures. For Gaussian, Laguerre and circular -ensembles we construct sharp couplings to the stochastic Airy, Bessel and Sine operators, preserving the finite-matrix laws. For fixed , their squared Hilbert-Schmidt inverse errors have almost-sure orders and , respectively, with in the hard-edge upper bound. We also identify the positive leading order coefficient for the circular ensemble, whose second derivative jumps at .
Abstract
In this paper, we investigate the expansion of the integral form Macdonald polynomials in the elementary symmetric function basis. Motivated by the connection between chromatic symmetric functions and LLT polynomials, we conjecture the -positivity of a normalization of . We verify the conjecture when is either a hook or a two-column partition. In the hook case, we express the -coefficients as weighted sums over compositions. In the two-column case, we derive explicit formulas for the -coefficients.
Abstract
We study the edge scaling limit of the logarithms of the squared singular values of multiplicative Brownian motion on the general linear group with general initial conditions. Building on earlier work [arXiv:2605.06429], we prove that, for every initial condition in the enhanced state space, the limiting dynamics are determinantal and derive an explicit double-contour formula for their extended multitime correlation kernel. The dependence on the initial state, including information not determined by the particle coordinates, is encoded by the associated Laguerre--Pólya entire function. We further extend the construction of solutions to the corresponding infinite-dimensional system of singular stochastic differential equations to arbitrary initial configurations, including those with coinciding coordinates.
Abstract
We construct stationary measures of the multispecies ASEP on the positive half-line, with injection of highest-priority particles at the left boundary, using vertex models. Our results extend existing constructions of Martin and Aggarwal-Nicoletti-Petrov for the full-line ASEP to the half-line. For the two-species process, we compute the joint distribution of occupation colors at the first two sites. We also partially address the connection with the half-space ASEP speed process.
Abstract
We consider a classical Coulomb gas of particles in dimension , confined by a general external potential in the large- limit. In that limit the density of particles has a compact support called the droplet. We study smooth linear statistics at large . While the leading variance of is by now well understood in many settings, much less is known about higher cumulants. When and the potential is rotationally invariant, these cumulants were computed recently at leading order. We derive explicit formulas for the leading third and fourth cumulants for general smooth droplets, without assuming that either the potential or the test function is radial. The answers are expressed as intrinsic boundary integrals over the equilibrium droplet. The third cumulant depends on the cube of the jump of the normal derivative of the harmonic extension of , while the fourth cumulant also involves the exterior Dirichlet-to-Neumann operator, the mean curvature of the boundary, and normal derivatives of the background density. Our derivation uses the tilted potential , the associated deformation of the droplet, and shape derivatives of the Dirichlet energy. This deformation has a natural interpretation as the infinitesimal motion of the free boundary in the obstacle-problem formulation of the equilibrium measure. Our formulas extend previous studies by Wiegmann and Zabrodin in two dimensions for . In their setting, we provide a rigorous proof of the limiting third and fourth cumulants, supporting the proposed formulas for general Coulomb gases. We specialize the general formulas to the disk, the half-plane and the sphere, and check them against exactly solvable non-radial deformations, including translated elliptic droplets in and spheroidal deformations in .
Abstract
This work constructs a non-trivial solution to the ill-defined two-dimensional Kardar-Parisi-Zhang (KPZ) equation at criticality: a distribution-valued stochastic process on the plane, which we call the Critical 2D KPZ Process. It arises as the limit in law of centred log partition functions of two-dimensional Gaussian directed polymers at criticality, viewed as discretised KPZ solutions. It is also a measurable function of the Critical 2D Stochastic Heat Flow (SHF) via a renormalised Cole-Hopf transformation, obtained as the limit of the centred logarithm of the spatially mollified SHF as the smoothing is removed. The Critical 2D KPZ process is non-Gaussian and its covariance displays a log-log blow-up near the diagonal. Gaussian fluctuations emerge when averaging on shrinking spatial scales. The results presented here were obtained through extensive interactions with generative AI (ChatGPT 5.5, 5.6 and 6); see Appendix C for more details. The proof builds on key recent and older advances in our field, combined in novel ways with substantial AI input. We view this manuscript as a "community project" that represents more than a decade's collective effort. We have worked on making the presentation as accessible as possible, highlighting the main ideas and their relations to the literature.
Abstract
We use the partition functions of a five-vertex exactly solvable lattice model from statistical mechanics to give new combinatorial formulae for the Schubert structure constants of torus-equivariant quantum K-theory rings of Grassmannians. The proofs build on our earlier results where we identified the transfer matrices of these lattice models geometrically as multiplications by the Hirzebruch classes of the tautological subbundle and the (dual) quotient bundle. The structure constants are solutions of a system of linear equations with coefficients given by the partition function. This system may be solved explicitly by specializing the Hirzebruch parameter to an equivariant parameter, leading to explicit formulae using Cramer's rule. Furthermore, we give an alternative combinatorial construction of the partition functions in terms of sums over cylindric set-valued tableaux, which we call cylindric Grothendieck polynomials. The latter generalise McNamara's factorial skew Grothendieck and Yeliussizov's modified Grothendieck polynomials and can be thought of as an equivariant K-theoretic version of Postnikov's toric Schur function. We also provide new polynomial representatives for equivariant K-theory classes of Richardson varieties and their variants.
Abstract
We derive exact formulas for the diagonal distribution of half-space geometric last-passage percolation with a two-parameter family of stationary initial conditions throughout the phase diagram. Unlike previous results, our formulas remove the random shift at the origin. Under critical scaling, we obtain new formulas for the half-space stationary KPZ fixed-point distribution at in the high-density and maximal-current phases. In the low-density phase, we derive a new representation and show that it matches the corresponding formula from the stationary log-gamma polymer.
Abstract
The stationary increment measure of the open Kardar-Parisi-Zhang equation on a bounded interval can be represented as the independent sum of a Brownian motion and another process defined via an explicit Radon-Nikodym derivative with respect to the measure induced by a Brownian motion. In the fan region, its finite-dimensional Laplace transform was previously expressed in terms of the continuous dual Hahn process. The main contribution here consists of computing the finite-dimensional Laplace transform for a range of parameters in the shock region by an analytic continuation argument. As an application, we provide a direct proof of the convergence of the finite-dimensional distributions of the increment height function of open WASEP under the stationary measure in the shock region, by computing the convergence of finite-dimensional Laplace transform without relating to the open KPZ equation.
Abstract
We study joint moments of characteristic polynomials and their second derivatives for Haar-distributed matrices from the unitary symplectic and orthogonal groups. We establish their large- asymptotics throughout the full admissible range of the moment parameters and identify the limiting coefficients probabilistically in terms of the first inverse power sum of the Bessel point process. This extends the corresponding results of Assiotis, Gunes, Keating, and Wei [Proc. Lond. Math. Soc., 132 (2026), e70136] to the full admissible parameter range. Throughout this full range, we also express these limits in terms of a distinguished solution of the -Painlevé III' equation. Our approach reduces the problem to singular inverse-power statistics of the Jacobi unitary ensemble and combines their finite- -Painlevé V representation with a refined Riemann-Hilbert analysis yielding precise small- asymptotics and uniform derivative estimates for the finite- correction term. A key ingredient in the probabilistic identification is a kernel-based proof that the rescaled first inverse power sum of the Jacobi ensemble converges in distribution to its Bessel point process counterpart.
Abstract
In directed last-passage percolation with i.i.d. exponential vertex weights, the coalescence time of two semi-infinite geodesics started from vertices and in a fixed direction has a heavy tail with (KPZ) exponent , and hence infinite expectation. We study instead the first coalescence among triples of semi-infinite geodesics in first- and last-passage percolation models (FPP/LPP) on that arise from stationary compatible families, such as Busemann geodesics. For directed FPP and LPP with arbitrary weights, we prove that We also consider FPP and LPP models in which antidiagonal motion is allowed, and prove that for every integer , \[\mathbb{E}\bigl[\min\{\operatorname{Coal}(a,a-k(e_1-e_2)),\operatorname{Coal}(a,a+k(e_1-e_2))\}\bigr]<\infty.\] These results support the heuristic prediction that adjacent coalescence times are negatively associated.
Abstract
We characterize central limit theorems for polynomial linear statistics of continuous -particle systems on , in both the fixed-temperature regime and the high-temperature regime . In each regime, we give necessary and sufficient conditions in terms of the coefficients of the logarithmic Bessel generating function, together with an explicit combinatorial formula for the limiting covariance. At fixed temperature, this formula recovers the classical second-order free moment--cumulant relation, up to the factor ; at high temperature, it yields new second-order free -cumulants and their covariance relations. The proof combines moment extraction through Dunkl operators with a topological expansion of the logarithmic multivariate Bessel function indexed by ribbon graphs on surfaces. This expansion gives a topological interpretation of the free-cumulants: only genus zero contributes at fixed temperature, whereas arbitrarily high genera contribute at high temperature. As applications, we prove that -addition preserves the CLT in both regimes and determine its limiting covariance from the input free cumulants, and establish a suitably rescaled transition from the high-temperature covariance to its fixed-temperature counterpart as .
Abstract
In this article, we prove that the spectrum, the norm, and the largest eigenvalue of a self-adjoint polynomial in independent matrices from the GOE satisfy a large deviation principle.
Abstract
We construct an infinite-dimensional diffusion on the space of interlaced particle configurations on Z x R, a process we call the dynamical bead model. The process (self-)interacts through instantaneous reflection and preserves the bead measures of [C. Boutillier, Ann. Probab., 2009]. Through the height function correspondence, our process is equivalent to a reversible interface in (2+1) dimensions. We provide compelling evidence that at large space-time scales the interface exhibits diffusive behaviour, as we compute exactly its diffusivity matrix and show it is non-degenerate and independent of time.
Abstract
In this paper, we investigate the relation between t-embeddings and T-graphs associated with the same weighted planar bipartite graph. We prove that, given such a T-graph, the corresponding t-embeddings are classified by positive stationary measures on this T-graph and describe a reconstruction procedure of the t-embedding given such a measure. As an application, we develop an approach for constructing t-embeddings of arbitrary sufficiently non-degenerate planar bipartite graphs.
Abstract
We study the Seppäläinen-Johansson model of first-passage percolation along a region of the plane where it moves from an ordered to a disordered phase, represented by a 2nd order discontinuity in its limit shape function. The fluctuations in this transition region are of order one, and we determine that certain limiting statistics are described by Fredholm determinants and are closely related to the largest eigenvalue of the GUE.
Abstract
We consider exponential last passage percolation in half-space with a deterministic first row. This model corresponds to a version of TASEP in half-space in which particles enter at the origin and are restricted by a moving wall. We determine the possible limiting distributions of rescaled last passage times depending on the boundary strength and the choice of the wall. As key inputs to our analysis, we establish several results for point-to-point half-space last passage times, including upper tail bounds, weak convergence, and bounds on the probabilities that geodesics touch the diagonal for certain start and end points.
Abstract
A recent graph-theoretic criterion for a Hamiltonian to exhibit a free-parafermion spectrum is that its frustration graph be an indifference graph. We show that an arbitrary finite indifference graph can be realized as a frustration graph of the quantized six-vertex model at odd roots of unity, by a suitable choice of the defining wiring diagram and parameters. This provides a concrete realization of the criterion through the quantum Frobenius property of the quantized six-vertex model, recently established from its three-dimensional integrable structure.
Abstract
We study random walks with geometric or log-gamma increments, reweighted by certain exponential functionals. We establish precise asymptotics for all moments of exponential functionals, revealing a phase transition as the drifts of the walks vary. To describe the large time limit of reweighted random walks, we introduce a discrete transformation on random walk paths and establish Girsanov relations. As a corollary, we show that the stationary measures of geometric last passage percolation and the log-gamma polymer on a strip converge to their half-space counterparts as the strip width tends to infinity. In the log-gamma case, this provides an alternative proof of a recent result in arXiv:2607.17030.
Abstract
The symmetric exclusion process (SEP) is a paradigmatic model of interacting particle transport. Its current--density correlations are governed by an infinite many-body hierarchy of equations. Only recently, this hierarchy was found to collapse at large scales onto a single closed equation. Is this remarkable simplification created by coarse graining, or does it reflect a deeper property already present in the microscopic dynamics? Here we show that the second scenario is realized. At arbitrary time on the lattice, the entire hierarchy collapses onto a single linear equation which is exactly equivalent to a persistence problem for a single auxiliary random walker. Thus an interacting many-body problem reduces exactly to a genuine one-body problem before any scaling limit is taken. This reduction determines the complete finite-time correlation profile, including during nonequilibrium relaxation. More broadly, the microscopic operator uncovered here opens a route to more general correlation functions in the SEP, and may point to a general organizing principle for nonequilibrium interacting particle systems.
Abstract
Let count permutations of whose longest increasing subsequence has length . Chen conjectured that is log-concave for every positive integer . In this paper, we prove that Chen's conjecture is true when . The proof combines Gessel's determinantal formula with quantitative coefficient extraction, Tracy--Widom density estimates, and exact inequalities for the extreme lengths.
Abstract
We derive a Fredholm-Pfaffian formula for the one-time multi-point transition probabilities of half-space polynuclear growth with any nonnegative boundary nucleation rate and general deterministic upper semi-continuous initial profiles. We prove the formula by showing that it satisfies an integral Kolmogorov equation with a unique bounded solution. We also discuss the formal 1:2:3 scaling limit of the kernel and its relation to the half-space KPZ fixed point.
Abstract
We introduce a determinantal linear process on the space of -forms over a weighted graph endowed with a unitary local system, which we call the determinantal linear spanning forest (DLSF). This process is the analog of the uniform spanning tree in the presence of an arbitrary rank quenched gauge field in (or , , , ), and generalizes the determinantal random cycle-rooted spanning forest which corresponds to the case of and . We describe the support of this process by means of a polymatroid which is described in terms of the circular and bicircular matroids. We also explain how the DLSF relates to the computation of currents in an electrical network in the presence of a gauge field. The partition function of the DLSF is shown to be the determinant of the twisted Laplacian, for which we also give a matrix-tree type theorem that computes its characteristic polynomial. We initiate the study of the twisted Laplacian determinant as a functional on the space of connections, and give a first order expansion of it in a neighborhood of the trivial connection, for which we provide a probabilistic interpretation. We extend the construction to infinite graphs, exhibiting the analogs of the classical wired and free distributions on spanning forests. We give a sufficient geometric condition on the local spectral gap for the existence of a Poincaré inequality as well as the uniqueness of the infinite volume measure, and the exponential decay of correlations. We show that on any recurrent graph, the wired and free DLSF coincide. We construct examples where they do not, for example, unlike in the classical case, on for .
Abstract
In this paper we construct a one-parameter family of periodic Airy line ensembles for as weak scaling limits of Borodin's periodic Schur processes. We establish several properties of , including a determinantal point process structure with an explicit correlation kernel. We also show that, after a parabolic shift and normalization, satisfies, on its domain , the same Brownian Gibbs property enjoyed by the usual Airy line ensemble. A main technical contribution of the paper is a general tightness criterion for determinantal line ensembles, which extends the argument used by Johansson in [Comm. Math. Phys. 242 (2003), 277-329] to prove a functional limit theorem for the discrete polynuclear growth model.
Abstract
Hermite and Laguerre -ensembles are important and well studied models in random matrix theory, with the special cases corresponding to classical matrix ensembles. Although sharp moderate-deviation estimates have been established for the largest eigenvalues, substantially less is known for other edge eigenvalues. Even for the limiting distributions , the leading constants in the tail asymptotics are not known for general and . In this paper, for each fixed , we prove matching upper and lower moderate deviation bounds for the right tail of the -th largest eigenvalue, with the optimal exponential constant , for general . We also establish sharp lower tail estimates for the second largest eigenvalue of the Gaussian and Laguerre orthogonal ensembles, with exponential constant , same as for the largest eigenvalue. For , we obtain optimal upper and lower tail asymptotics for general and . We also obtain new stochastic domination results for the edge eigenvalues. Our proofs combine variational arguments from tridiagonal matrix models, last passage percolation, stochastic comparisons of sums of edge eigenvalues, and superposition-decimation identities.
Abstract
We study the geometry of semi-infinite geodesics in the Airy line ensemble. We establish sharp estimates for disjointness probabilities involving geodesics with nearby directions and endpoints in the top line of the Airy line ensemble. Moreover, we obtain a power law upper bound on the coalescence depth of two semi-infinite geodesics with same asymptotic direction and establish exponential concentration estimates for geodesic jump times more generally. We also establish pathwise Hölder -continuity of Busemann profiles associated to semi-infinite geodesics with common asymptotic directions (in terms of the asymptotic direction). We characterize exceptional behaviour of semi-infinite geodesics through their Busemann profiles and establish geometric properties of the associated sets of non-coalescence and non-uniqueness. Finally, we determine the Hausdorff dimensions of collections of semi-infinite geodesics: collections with a fixed top-line endpoint and arbitrary positive asymptotic direction almost surely have Hausdorff dimension , while collections with a fixed positive direction and arbitrary endpoint almost surely have Hausdorff dimension .
Abstract
We prove the Berezinskii-Kosterlitz-Thouless (BKT) phase transition in the large-distance behaviour of the two-point correlation function of the point process for all , confirming the first term of an asymptotic expansion predicted by Forrester.
Abstract
Permutons can be thought of as "continuum limits" of permutations of diverging size, which encode information about the number of occurrences of any given finite subpattern. We construct a natural composition operation on permutons, analogous to composition of permutations, and explore its properties via a correspondence between permutons, Markov kernels, and Markov operators. We show that no extension of permutation composition to permutons can be (jointly) continuous with respect to the usual "weak" topology on permutons, but we construct a "strong" topology with respect to which composition is jointly continuous. Moreover, we characterize strong convergence of sequences of permutations of diverging size: a sequence of permutations converges strongly if and only if it converges weakly and its limit is the transpose of the graph of a measure-preserving function. The key insight behind these results is to apply functional analysis to certain operators associated to permutons; with this shift in perspective, our results have elegant and relatively "soft" proofs. As an application, we describe the permuton limit of a product of finitely many independent uniformly sampled separable permutations and give a formula for the expectation of this limit. We also show that in this case, as the number of compositions goes to infinity, the aforementioned expectation tends to the uniform permuton.
Abstract
The strong Szegő theorem and the Fisher-Hartwig formula describe the asymptotic behavior of determinants obtained by restricting translationally invariant operators to large finite regions. In particular, Szegő's theorem states that, beyond the leading extensive contribution to such a determinant, there is a universal finite correction, where are the Fourier coefficients of the logarithm of the symbol. Although this formula appears throughout analysis and mathematical physics, the origin of its quadratic form and of the weight is often obscured in its standard analytic derivations. We develop an elementary picture in which the strong Szegő term arises from boundary leakage: paths contributing to the matrix exponential leave the finite interval and subsequently return. Counting the shortest such excursions immediately produces the factor . We then show that an isolated boundary contributes exactly one half of the strong Szegő correction, while the two boundaries of a long interval decouple. For finite-range matrices this gives a rigorous and intuitive elementary derivation of the asymptotics, with extensions to sufficiently rapidly decaying interactions. A complementary free-fermion argument provides a physical interpretation of the same cancellation mechanism. Numerical examples illustrate the asymtotics in the well controlled regime, the Fisher Hartwig logarithmic correction. The boundary formulation does not itself require translation invariance and therefore suggests natural extensions to disordered systems, higher dimensions, and symbols with internal matrix structure. To illustrate this we numerically go beyond the formal applicability region of the Szegő formula, considering slow decay and disordered cases, and showing in some cases a robustness of the approach even for slowly decaying but disordered potentials.
Abstract
While models of undirected and directed geometry in random i.i.d. disorder with rapidly decaying tails are expected to behave similarly and exhibit features of the Kardar--Parisi--Zhang universality class, the scenario when the background noise is fractal is expected to be significantly different. A canonical example of the latter is Liouville Quantum Gravity, where the planar Gaussian free field forms the fractal background. A toy model capturing some of the essential features is given by Mandelbrot's fractal percolation, a random Cantor set. In the simplest setting, dyadic cubes of different scales are retained independently with probability (a parameter of the model), and the intersection of all retained cubes forms the set of open sites. Connectivity properties of such sets have been intensely studied. Across [CCD88, Cha95], a surprising result was proven. Namely, while there is a phase transition for connectivity of the fractal percolation set (like in usual bond percolation), directed percolation never occurs. Results of a similar spirit have been obtained in geometric measure theory as well. As a step towards understanding last passage percolation (LPP) driven by the Gaussian free field (initiated in [GGN24]), we study the problem of quantifying LPP in fractal percolation. Developing a multi-scale framework, we show an almost polynomial gap from linear growth of the passage time. The obstruction to high LPP values is the Lipschitz nature of directed paths. Along these lines, [Cha96] showed that any connected subset of fractal percolation must have Hausdorff dimension strictly larger than one. We show that one can construct directed paths passing through open sites in fractal percolation provided they move exponentially fast in the spatial direction---a model we introduce and term as unrectifiable directed percolation.
Abstract
We prove that the sine beta point process is tail trivial for every beta >0. The key ingredient in the proof is a certain uniform entropy bound, that makes use of the DLR equations of Dereudre et al. (2021) and the bounded two-point correlation density obtained by Qu and Valko (2025) and Assiotis and Najnudel (2026), that allows us to show that tail events are translation invariant. Tail triviality then follows either from the uniqueness theorem of Assiotis (2026) or from translation mixing, which we derive from the correlation estimates of Dumaz and Malvy (2026). As consequences, we obtain the following applications: (a) extremality of sine beta in the full convex class of fixed kernel Gibbs measures without stationarity and finite-energy conditions; (b) L^2-ergodicity of the associated infinite-particle reversible diffusion for every beta>0; and (c) the almost-sure collision/noncollision dichotomy with threshold beta=1 for the diffusion.
Abstract
Exponential directed last passage percolation is a random maximization model connected to tandem queues and interacting particle systems. We study the line-to-point model of the percolation, in which the starting point is allowed to vary inside a line boundary, conditioned on an atypically large last-passage time to a fixed endpoint. Extending recent results for the point-to-point model where the starting point is fixed, we develop heuristic predictions for the conditional limit shape and prove corresponding fluctuation results in certain regimes through steepest-descent analysis of exact multi-point formulas. Our calculations suggest that the conditional fluctuations near diagonal are described by Brownian bridges coupled through an additional shared Gaussian component.
Abstract
We study a one-dimensional multispecies TASEP model with each particle type having an associated speed parameter and faster particles overtaking slower ones with the rate being the positive part of their speed difference. After expressing the sorting time of a finite system given a reverse-speed initial condition as the maximum of sorting times for all different fast-slow cuts, we prove that forbidding fast-fast and slow-slow overtakes does not alter the law of the sorting time. This allows for an LPP formulation and an analytic expression of the sorting time.
Abstract
We consider the adjacency matrix of a uniformly random simple -regular graph on vertices. For every fixed and every fixed bulk energy, we prove that the rescaled eigenvalue point process converges to the process with intensity . We also establish universality of consecutive gaps at deterministic bulk labels and of fixed-energy correlation measures. The eigenvector input consists of polynomial estimates for mixed fourth moments on a uniformly sampled star. These estimates prevent loss of variance in the qualitative joint Gaussian limits of Backhausz--Szegedy, yielding independent variance-one Gaussian waves for fixed tuples at nearby energies. Exact switching identities then give a stationary marked flow and the microscopic loop hierarchy, which identifies the point-process limit. The appendices prove the stronger correlation and gap statements using counting-moment bounds and conditional log-gas laws for tagged limits.
Abstract
In this paper, we show that the height function of the six-vertex model in simply connected domains with piecewise constant boundary conditions converges to the Gaussian Free Field (GFF) with Dirichlet boundary conditions. This implies the conformal invariance of the scaling limit of the model. The argument relies on the full plane GFF convergence combined with Reflexion Positivity and qualitative estimates. The proof works for a wide range of spectral parameters (among which close to , and between and ) making this paper the first instance of a conformal invariance result for a continuum of models. The paper has consequences for other models, including the FK percolation with cluster-weight and the Ashkin-Teller model, which we will pursue in future work.
Abstract
We study the full counting statistics of the integrated current in a one-dimensional hard-rod gas in equilibrium, employing two recently proposed fluctuating hydrodynamic frameworks: (i) Ballistic Macroscopic Fluctuation Theory (BMFT) and (ii) diffusive fluctuating hydrodynamic theory based on correlated Brownian motion approximation of hard-rod quasiparticles (CBFHD). We investigate how the predictions of these two approaches compare with each other and with the underlying microscopic dynamics. We show within BMFT that the current distribution admits a large-deviation form under ballistic scaling, and we derive an explicit expression for the large-deviation function, valid for an arbitrary inhomogeneous initial state factorized in point-particle coordinates. The corresponding current cumulants are also obtained. We subsequently analyze the current statistics within the CBFHD framework, which explicitly incorporates stochastic noise and dissipative terms, and demonstrate that the current distribution likewise obeys a large-deviation principle in this approach. We find that, while BMFT predictions for the first few cumulants and the rate function are in excellent agreement with direct numerical simulations of the microscopic hard-rod dynamics, CBFHD prediction is not. We further verify the predictions of the CBFHD independently using a collection of correlated Brownian motions, which provides an approximate description of the dynamics of hard-rod quasiparticles. Our findings indicate that CBFHD does not accurately capture large-scale current fluctuations in the hard-rod gas, whereas BMFT yields predictions in quantitative agreement with the microscopic dynamics, at least in the context of full counting statistics.
September 2026
Abstract
Let be a continuous-time Markov process taking values in a countable state space . We consider copies of started from distinct positions, and write for the resulting process of configurations in . We compare two forms of interaction between the particles: the first is conditioning the particles to never collide (noncolliding dynamics), and the second is simply switching off any jumps that would cause two particles to occupy the same position (exclusion dynamics). We show under mild conditions on the tail probabilities of the collision time that the dynamics are related by an explicit exponential martingale change of measure of one another. This generalises a recent result in the integrable probability literature relating TASEP on the ring to asymmetric walkers conditioned to never collide.
Abstract
In the adjacent transposition shuffle on a cycle, distinct cards are placed at the vertices of a cycle, and adjacent cards swap positions according to independent Poisson clocks of rate 1. We prove that this Markov chain exhibits total variation cutoff at time , with a window of order at most .
Abstract
For any , consider the process defined by , where is the directed geodesic from to in the directed landscape. Let . We show that the laws of and are mutually absolutely continuous if and only if . This proves Conjecture 14.5 in Dauvergne, Ortmann, and Virág, 2022 in the affirmative.
Abstract
In this paper, we establish Edgeworth expansions of extreme eigenvalue distributions for two types random unitary ensembles at the spectral edges and reveal certain universal structures for the correction terms. More precisely, we show the correlation kernels for the Gaussian-type unitary ensembles and Laguerre-type unitary ensembles admit full expansions at the soft edge and the hard edge, respectively. The coefficients of the correction terms are given by finite sums of Airy function or the Bessel functions of the first kind and their derivatives with polynomial coefficients. By lifting the kernel expansion to the associated Fredholm determinant, we further obtain full expansions for the largest eigenvalue distribu- tion of the Gaussian-type unitary ensembles and the smallest eigenvalue distribution of the Laguerre-type unitary ensembles. In addition, the first correction term therein involves the derivatives of the leading term.
Abstract
Homogeneous coordinate rings of Grassmannians are among the most fundamental instances of cluster algebras arising in ``nature.'' Despite their importance, however, the cluster variables of these cluster algebras are not well-understood, except in special cases. For the Grassmannian of 3-planes in n-space, Fomin and Pylyavskyy conjectured [FP16] that the cluster and frozen variables are exactly the indecomposable, non-elliptic web invariants which are arborizable. They proved that in particular, if an SL3 tensor diagram is a planar tree, its web invariant is a cluster or frozen variable. In this article, we give a new proof of their result for Gr(3, n) planar tree webs, and we generalize the result to Gr(4, n): if an SL4 tensor diagram for Gr(4, n) is a planar tree web, its web invariant is a cluster or frozen variable. Our proof uses some of the cluster quasihomomorphisms from [EZLP+23] and [EZPT+25] (namely, upper promotion and spurion promotion), as well as a new 6-leg pattern map.
Abstract
We study spin-current fluctuations of the folded XXZ chain, which describes the strong-anisotropy limit of the XXZ chain. We obtain an exact expression for the moment-generating function of the time-integrated spin current at infinite temperature. A long-time asymptotic analysis establishes that the rescaled time-integrated current converges to a non-Gaussian distribution described by the M-Wright function. The derivation relies on a nonlocal map to the one-dimensional t-0 model and its exact spin-charge separation. By extending this approach, we also show that the dynamical correlation function of the spins summed over two adjacent sites exhibits a diffusive Gaussian profile at long times.
Abstract
The balance constant is a poset parameter measuring how well random linear extensions can be split according to their relative order on two elements. We give an -improvement for the balance constant over the bound by Brightwell, Felsner and Trotter (1995), for some small . This is the first general result towards the -- conjecture in over 30 years.
Abstract
The dimer model and random lattice permutations are two fundamental objects at the interface of probability, combinatorics, and mathematical physics. We study these models on finite periodic boxes in within a common framework. For the dimer model, edges connecting arbitrary vertices carry a weight which depends on their relative displacement and dimer configurations are weighted through their occupied edges. Superimposing two independent perfect matchings gives the double-dimer model, whose configurations are collections of disjoint loops. Permutations, instead, are weighted through the spatial displacement of their jumps. For broad classes of weights of finite or infinite range we prove long-range order and the occurrence of macroscopic loops. This extends nearest-neighbour results to arbitrary-range edges and jumps. In particular, long-range weights yield long-range order and macroscopic loops already in dimensions . In dimension two, this behaviour is qualitatively different from that of the nearest-neighbour model. We complement these results with sharp absence criteria.
Abstract
The Anderson tight-binding model is a fundamental model of quantum transport and localization in disordered media. Completing a problem left open by Bourgain and Kenig, this paper proves Anderson localization near the bottom of the spectrum for the lattice Anderson model with Bernoulli potential, in any dimension . The proof uses the multiscale framework of Fröhlich-Spencer and Bourgain-Kenig, and the main new ingredient is a probabilistic discrete unique continuation principle (PDUC) for the discrete Schrödinger equation. This PDUC is established via a bootstrap argument and a key probabilistic lemma proved by adaptively revealing the random potential.
Abstract
We study the quantized six-vertex (Q6V) model on arbitrary admissible diagrams at , where is a primitive th root of unity with odd. For the commuting layer transfer matrices with mixed boundary conditions, we establish the quantum Frobenius property \[ \mathbb{T}(y)\mathbb{T}(\varepsilon y)\cdots \mathbb{T}(\varepsilon^{N-1}y) = \mathscr{P}_N(y^N)\mathbb{I}, \] where is a scalar Laurent polynomial obtained directly from the monomial expansion of . Our proof exploits the underlying three-dimensional integrability, in particular the tetrahedron equations and their consequences. The quantum Frobenius relation essentially determines the joint eigenvalues of the associated Hamiltonians up to multiplicities and shows that they take a free-parafermion form, thereby extending such structures from previously known one-dimensional examples to a broad class of genuinely two-dimensional quantized six-vertex models. The construction also admits independent local gauge choices and recovers the long-range-interacting model, which includes the original free-parafermion model as a special case. We also apply the same framework to the relativistic quantum Toda chain at odd roots of unity and obtain the free-parafermion spectra for the relevant Q6V transfer matrices with fixed boundary states.
Abstract
We find a new coupling between the Airy sheet and the Airy line ensemble , which we call the down-the-middle isometry. In this coupling, the Airy sheet is encoded by last passage values that start at for some and end on the top line . Using this coupling, we give a new construction of the (extended) Airy sheet as a scaling limit of Brownian last passage percolation and study coalescence of geodesics in the directed landscape. This coupling extends to discrete polymers and last passage percolation models, as well as to the KPZ line ensemble and KPZ sheet.
Abstract
The KPZ upper-tail field introduced by Liu and Zhang arises as the scaling limit of the narrow-wedge KPZ fixed point in a microscopic neighborhood of a conditioned large height. In this paper, we give a direct probabilistic identification of the upper-tail field. We identify the grounded positive-time sector with the KPZ fixed point started from the two-sided Brownian \(V\)-profile \(B_{\mathrm{ts}}(2x)-2|x|\), while the negative-time sector is represented through an explicit exponential tilt of the same evolution. The argument uses the directed-landscape variational structure and the local Brownian comparison under upper-tail conditioning. As consequences of these representations, we recover probabilistically the positive- and negative-time large-scale limits of the KPZ upper-tail field.
Abstract
We establish the large- asymptotics of Hankel determinants for unitary random matrix ensembles possessing a higher-order edge singularity. We focus on ensembles where the equilibrium measure is supported on a single interval and the limiting eigenvalue density vanishes to order for . Notably, we explicitly evaluate the constant term in the asymptotic expansion, which involves a regularized integral of the Hamiltonian associated with the Painlevé I () hierarchy. As a by-product, we also prove the universality of the eigenvalue correlation kernel near this singular edge and derive a limiting kernel expressed through functions related to a special solution of the equation. Our method relies on the Deift-Zhou nonlinear steepest descent analysis for the Riemann-Hilbert problem of orthogonal polynomials.
Abstract
We study infinite-dimensional Dyson Brownian motions obtained as limits of finite systems without rescaling the actual stochastic dynamics. For , we construct determinantal processes on an extended space of initial data and prove convergence of their finite-dimensional distributions under essentially optimal conditions. This extends the seminal results of Katori and Tanemura. The additional parameters record information at infinity and enter through an associated Laguerre-Pólya entire function. Moreover, for explicit classes of configurations, we establish convergence on path space and the Markov property. We prove rescaled long-time convergence, in finite-dimensional distributions, to the stationary extended process from arbitrary symmetric initial configurations with power-law counting exponent . This extends the integer lattice relaxation result of Katori and Tanemura which was the only such result for explicit deterministic initial conditions. For , we prove convergence of finite particle systems from regular initial data to the unique strong solution of an infinite-dimensional stochastic differential equation in a certain rigid-path-regularity class. This extends seminal works of Tsai and Osada. We finally derive a stochastic partial differential equation of Burgers-type for the Stieltjes transform of the dynamics.
Abstract
This paper establishes a mesoscopic central limit theorem for linear statistics of -ensembles or log-gas, as the dimension , in the temperature regime . For simplicity, we assume that the potential is one-cut regular and analytic. In this regime, the size of the fluctuations depends on and the mesoscopic scale. We show that there is a transition at a critical between a Random Matrix regime, where the limiting variance is given by the -norm and a Poisson regime where the limiting variance is given by the -norm. We also describe the critical regime. The proof of the CLT relies on optimal local laws at intermediate temperatures and Stein's method for -ensembles. In particular, in this regime, it is necessary to construct new correction terms to the classical equilibrium measure to obtain a suitable re-centring of linear statistics and describe their fluctuations. We also obtain a free energy expansion.
Abstract
Koornwinder moments generalize the Askey--Wilson moments arising in the asymmetric simple exclusion process. Rains conjectured that, under the specialization , the minimal numerators of Koornwinder moments have nonnegative integer coefficients. While the one-row case of this conjecture was previously proved by Corteel, Mandelshtam, and Williams using rhombic staircase tableaux, its dual counterpart, the one-column case, has remained open. In this paper, we derive a closed, manifestly positive combinatorial formula for the normalized numerators of the coefficients of the rescaled Askey--Wilson polynomials. This proves Rains' conjecture for one-column partitions, thereby providing the exact dual counterpart to the previous result.
Abstract
We define a natural orientation on a spanning-tree-decorated planar map whereby, roughly speaking, each directed edge in the map is oriented to match the direction of the contour exploration of the spanning tree. We study directed distances (lengths of shortest directed paths) with respect to this orientation. We construct the Busemann function which measures directed distances to along a natural interface in the uniform infinite spanning-tree-decorated map. We show that this Busemann function, re-scaled appropriately, converges in law to a -stable Lévy process. We also show that in a uniform spanning-tree-decorated map with edges, directed distances are typically of order . Using a strong coupling argument, we deduce analogous statements for directed distances in other random planar maps in the -Liouville quantum gravity (LQG) universality class, including uniform meandric systems and mated-CRT maps for . These results give the scaling dimension for a hypothetical directed version of the -LQG metric. Our proof strategy is inspired by work of Borga and Gwynne (2025) on directed distances in bipolar-oriented triangulations.
Abstract
Eigenvector overlaps quantify nonorthogonality and govern the response and dynamics of non-Hermitian systems. We extend the universality of non-Hermitian random matrices to the statistics of these overlaps. We obtain analytical expressions for eigenvector overlaps in the spectral bulk and near the origin, covering ten symmetry classes in the limit of large matrix size. Using fermionic replica nonlinear models, we relate these overlaps to Hermitian level statistics, symmetry class by symmetry class and topological sector by topological sector. Numerical calculations in various physical models support the universality of the normalized overlaps in the regimes studied. Our work establishes a duality between Hermitian level statistics and non-Hermitian eigenvector overlaps.
Abstract
In this paper we construct a two parameter family of generalized Gaussian fields which interpolates between the Gaussian Free Field, conditioned Gaussian Free Fields, and Gaussian Free Fields with independent noise. Depending on the values of the parameters, we characterize these generalized Gaussian fields as partial conditionings of a Gaussian Free Field with perturbations of its first two Fourier modes. We show that this family arises as the fluctuations of the height function associated with the corner process of a two parameter variant of Wigner random matrices.
Abstract
We consider a Coulomb gas on a Jordan curve in an external potential at inverse temperature and obtain an asymptotic expansion of the free energy up to and a central limit theorem for linear statistics. We focus on the one-cut regime, where the density of the weighted equilibrium measure of in is strictly positive on . The constant term in the (normalized) expansion consists of two parts: the Fredholm determinant of a generalized Grunsky operator and the Dirichlet energy of the logarithm of the density of the weighted equilibrium measure of . The coefficient of the latter vanishes for . The variance of the fluctuations of the linear statistics only depends on the Dirichlet energy of the test function and is therefore independent of . Essential in our approach is that the generalized Grunsky operator and the accompanying equilibrium parametrization allow us to transport the particles on the curve in the external potential to a reference object in a way that preserves the equilibrium measure. In our setting, the unit circle is the natural reference object.
Abstract
We prove that the height-level path ensemble of a uniform alternating sign matrix converges to the Airy line ensemble near every interior point of the north-west arctic arc. Under the standard bijection, this is equivalently a limit theorem for the osculating paths of the domain-wall six-vertex model at the ice point. The proof has two main inputs. The integrable input is a new Pfaffian formula for rectangular frozen-corner probabilities, from which we obtain one-point convergence of the top path to the GUE Tracy--Widom distribution. The probabilistic input is an exact resampling law and accompanying bridge estimates for a strictly ordered midpoint encoding of the paths; these yield tightness and the Brownian Gibbs property for every subsequential limit. Aggarwal and Huang's strong characterization then identifies each such limit as the Airy line ensemble.
Abstract
We establish the existence and uniqueness of solutions to a -dimensional fractional stochastic Burgers-type equation featuring a Caputo fractional time derivative, a fractional Laplacian, and space-time white noise forcing modified by a Riemann-Liouville fractional integral. From a physical perspective, equations of this type describe non-Markovian dynamics with long-range temporal and spatial dependencies. The main difficulty is that the solution to the linearized equation has low spatial Hölder regularity, rendering the singular nonlinearity classically ill-defined. The proof relies on two main ingredients: first, combining the Da Prato-Debussche decomposition with rough path theory to rigorously interpret the nonlinear product; and second, deriving Schauder-type estimates for fractional heat kernels given in terms of slowly decaying Mittag-Leffler functions. The nonlocality of the fractional time derivative requires a strong assumption on the initial condition.
Abstract
In this note we announce limit shape theorems for the height function of the six-vertex model in the antiferroelectric (for arbitrary ) and disordered (in a restricted range ) regimes under domain-wall boundary conditions.
Abstract
We prove two families of distinguished cases of the heat flow conjecture for random matrices of Hall-Ho: elliptic Gaussian sources with planar targets, and deterministic Hermitian sources with targets on the boundary of the covariance disk. In particular, the empirical zero measure of the heat-evolved characteristic polynomial of a complex Ginibre matrix converges almost surely to the semicircle law. The strategy of the proof is to transfer estimates for the corresponding well-understood Gaussian matrix ensembles, through exact heat identities, into control of the zeros of the heat-evolved characteristic polynomials.
Abstract
We explain how to derive the Robinson-Schensted-Knuth (RSK) correspondence, a fundamental bijection in algebraic combinatorics, from the Yang-Baxter equation. The Yang-Baxter equation arose in the study of quantum many-body systems and later became a cornerstone of the theory of solvable lattice models, particularly vertex models. In a vertex model, arrows occupy the edges of a grid, and each vertex carries a Boltzmann weight determined by the arrows on the four edges meeting at it. The weight of a configuration is the product of these local weights, and a partition function is the sum of the weights of all configurations with prescribed boundary conditions. For a vertex model whose partition functions are the Schur polynomials, the two sides of each instance of the Yang-Baxter equation admit exactly one weight-preserving matching of their summands. Carried across a grid, this forced matching is the classical RSK correspondence in the form of Fomin's growth diagrams. In natural coordinates the local matching becomes the combinatorial three-dimensional R, a set-theoretic solution of the Zamolodchikov tetrahedron equation. Its periodic closure returns the combinatorial R-matrices of one-row crystals. The forced matching is special to the Schur weights. For the Hall-Littlewood and q-Whittaker deformations and their spin versions, at generic parameter values no deterministic matching works for all boundary data. Reading each instance of the Yang-Baxter equation probabilistically, we replace the matching by a coupling of the two sides - a bijectivization, or probabilistic bijection - and obtain Markov operators that transport probability measures attached to vertex models. Iterated over the grid, these operators produce randomized RSK-type dynamics and interacting particle systems, including q-PushTASEP and the stochastic six-vertex model.
Abstract
In this note we announce the limiting shape for the height function of a uniformly distributed alternating sign matrix (equivalently, a six-vertex model with equal weights and domain-wall boundary conditions)
Abstract
We study the minimum and maximum expected scores over all adaptive strategies in sequential card guessing with yes/no feedback. A deck contains labels, each appearing times, and is shuffled uniformly. The cards are guessed one at a time, and after each guess the player is told only whether it was correct. Always guessing the same label gives exactly points. How far can the expected score depart from if all previous feedback is used either to seek or to avoid correct guesses? Diaconis, Graham, and Spiro conjectured that the minimum expectation is as with sufficiently large in terms of , and that, for each fixed , the limit superior of the maximum expectation as is at most . Writing these extrema as and , we prove \[ 0\le m-P^-_{m,n}\le3m^{3/4} \qquad(n\ge m\ge1), \] and the finite-parameter bound \[ \frac{P^+_{m,n}}m\le\sum_{k=1}^{n}\frac1{k!}<e-1 \] for all positive integers . The first bound gives a first-order asymptotic uniform in ; the universal constant in the second cannot be reduced. Thus both conjectures follow. The lower bound combines a switching comparison for restricted permutations with martingale second-moment estimates. For the upper bound, we represent the repeated cards as a random interleaving of independent decks of distinct cards and prove that the optimal expectation is additive when the component-deck labels are revealed.
Abstract
For a large class of discrete- and continuous-time models satisfying a linear flow property, we prove convergence to the critical stochastic heat flow under diffusive scaling of space and time. Convergence is in the sense of finite-dimensional distributions in the space of measure-valued stochastic flows. The class of discrete models we consider includes random walks in space-time random environments with finite-range jumps, and directed polymer models with finite-range spatial correlations. The class of continuum models we consider includes several distinct types of linear-multiplicative stochastic PDEs whose driving noise is Gaussian with a smooth and compactly supported covariance kernel.
Abstract
This paper is devoted to introduce new probabilistic tree-based models complementing already known Matsumoto-Yor and Hamza-Vallois tree models with two tree beta-type models.
Abstract
We prove the Colomo-Pronko conjecture for alternating sign matrices with a prescribed square of zeros at a corner, for all matrix sizes and freezing parameters. A known multiple-integral formula for the frozen-corner count yields determinant representations built from fixed polynomial kernels. We relate these kernels to the conjectured determinant through an inverse identity for the commutator of a signed Pascal matrix with reversal. In odd dimension, the comparison uses the one-dimensional nullspace and projection along it to eliminate the central coordinate. Combined with the asymptotic analysis of Colomo and Pronko, our result removes the conjectural assumption from their GUE Tracy-Widom fluctuation theorem for the intersection of the frozen boundary with the main diagonal in uniformly random alternating sign matrices. The finite-dimensional algebraic core of the proof has been formalized in Lean 4.
Abstract
We resolve a conjecture on the law of large numbers for coherent measures on Young diagrams associated with Macdonald-nonnegative specializations of the algebra of symmetric functions for all . We prove that row and column lengths divided by the diagram size converge to the corresponding specialization parameters. The proof combines a monomial estimate for skew Macdonald polynomials with exponential tilting and binomial decompositions obtained from the branching rule for the polynomials.
Abstract
The two-dimensional one-component plasma forms a droplet in the case that the one-body potential is proportional to the number of charges. Our interest is in studying the large form of the corresponding partition function, with the plasma confined to the surface of the cylinder. The one-body potential is (up to technical restrictions) allowed to be arbitrary in the direction of the axis of the cylinder, whereas the coupling is restricted to the exactly solvable value . It is demonstrated that a change of variables can be made to an annular droplet plasma model, up to a certain Jacobian factor. The latter can be interpreted as the characteristic function of a certain linear statistic, which generalises the original aim of our study to also analysing the large form of the characteristic function for the annular droplet plasma model. While this is well known in the case of soft wall boundary conditions, our analysis (based on Laplace's method applied to certain integrals, and the Euler-Maclaurin summation formula) covers the case of one or two hard walls at the boundary, or inside, of the droplet for which the soft edge formulas are shown to no longer hold. Use of this to study the large expansion of the cylinder partition function requires knowledge of the same expansion in the annular case. This is not available in the case of two hard walls inside the droplet, which we treat instead via a direct analysis. In distinction to the case of the annular partition function, the large expansion of the logarithm of the cylinder partition function is found to always have no term proportional to , independent of the presence of hard walls.
Abstract
We prove a generalized circular rearrangement inequality: among all circular arrangements of a finite collection of positive numbers, the greedy arrangement simultaneously maximizes the sums of products of consecutive entries for every . This resolves a 2023 conjecture of Holmes-Holroyd-Ramírez and extends the classical circular rearrangement inequality for products of adjacent entries. As an application, we show that the greedy ordering minimizes the speed of a cyclic birth-and-death chain, resolving another conjecture of the same authors.
Abstract
In this paper, we study finite dimensional representations of shifted quantum affine algebras of type A. We give an explicit description of the tensor product of simple evaluation modules of the quantum loop algebra and a one-dimensional representation of the shifted quantum affine algebra under the separation condition. As a consequence, we give the q-characters of some finite dimensional simple modules of the shifted quantum affine algebra.
Abstract
We consider a pinned harmonic chain perturbed by random velocity flips, which conserve the total energy. Previous results have established the diffusive hydrodynamic limit for the energy profile at the level of expectations. In this article, under suitable fourth-moment bounds on the initial data, we prove for the first time a law of large numbers for the empirical energy distribution. Our approach is based on the Wigner distribution and shows that its random fluctuations vanish in the macroscopic limit. Consequently, the empirical energy profile converges in probability to the deterministic solution of the corresponding heat equation.
Abstract
We give the theoretical foundations of random matrix theory through the definitions of a random matrix, a random probability measure and the corresponding empirical spectral distribution. The technical tool we use is the Stieltjes transform method through which we prove optimal convergence of the empirical spectral distribution of random sample covariance matrices to the deterministic Marchenko-Pastur distribution. We also give new results about the rigidity of the eigenvalues of this random sample covariance matrix and the rate of their convergence. We then define the Dyson equation method to prove new local laws about a random matrix model that interpolates between the Marchenko-Pastur distribution, the elliptical law and the circular law. Through our work these local laws can be considered universal.
Abstract
We consider the cokernel corners process for an i.i.d. matrix with entries in a finite local ring. When the distribution of the entries is uniform, this process is a Markov chain, and hence the ergodic theorem for Markov chains can be applied. This implies, in particular, that for uniformly distributed p-adic random matrices, the cokernels of the corners are distributed according to the Cohen-Lenstra measure, almost surely. The purpose of this note is to show that the conclusion of the ergodic theorem also holds for i.i.d matrices, provided that the distribution of the entries is not concentrated on the translate of a subring, or the translate of an ideal. This will follow from the bounds proved in a previous paper of the author.
Abstract
We study the distribution of the number of double deficiencies (DD) in permutations of length n avoiding one or two patterns of length 3. Using structural decompositions of these avoidance classes--together with a lattice-path decomposition in the 321-avoiding case--we derive functional equations and convolution-type recurrences that efficiently compute the corresponding double-deficiency generating functions in all but one single-pattern case. In the 321-avoiding permutations, the resulting generating function is algebraic; we derive exact formulas for the mean and variance and prove that the distribution is close in total variation to Bin(n-2,1/4), with an explicit convergence rate. We also identify a DD-preserving symmetry that yields DD-Wilf equivalences, reducing the number of two-pattern cases that need to be considered separately. For the resulting two-pattern classes, we obtain explicit recurrences, including C-finite relations.
Abstract
We study the large time behavior of a periodic ASEP with fixed period started from the step initial condition. We establish the limit shape theorem and the asymptotic distribution of a single second class particle. Our single-species results extend the previous results of Lam and Ayyer-Linusson about periodic TASEP to the ASEP case. These results can also be interpreted as describing the limit behavior of a random walk on affine Hecke algebra.
Abstract
For some discrete parameters, multivariate (Dunkl-)Bessel processes on Weyl chambers are projections of Brownian motions on Euclidean spaces ; the most prominent examples are Dyson Brownian motions. More generally, the projections of Brownian motions on with drift are Bessel processes with drifts on , where again the associated transition densities can be described in terms of multivariate Bessel functions. Moreover, Dunkl processes with drift can be defined in an analogous way. In this paper we prove central limit theorems for these Bessel and Dunkl processes with drift for regular drift vectors, and arbitrary root systems and multiplicities. These results improve recent results of the author.
Abstract
We consider the exclusion process with energy, a model with mass and energy conservation, and we put it in a boundary-driven setting. We are interested in its hydrodynamic and hydrostatic limits, together with the associated dynamical and stationary fluctuations. We provide a new proof for the derivation of dynamical fluctuations that does not require explicit estimates on the two-point correlation functions. While the hydrodynamic limit is given by two uncoupled heat equations, the dynamical fluctuations are described by a system of coupled generalized Ornstein-Uhlenbeck processes.
Abstract
As is already known, there is a correspondence between circle patterns and the dimer model, which is invariant under Euclidean transformations. We establish a new correspondence that is invariant under the group of Möbius transformations PO(3,1), which is the natural symmetry group of circle patterns. We show that Miquel dynamics preserve the partition function, and that convexity of the t-embedding guarantees positivity of the face weights. We also show that the correspondence has a somewhat hidden symmetry group PO(3,3), which relates to the Lorentz lift of the t-embedding. Finally, we consider isoradial graphs, Doyle spirals, the once-punctured torus and isothermic circle patterns.
Abstract
For a symmetric weight on a finite discrete lattice and its Christoffel transforms , we prove that, at each Christoffel step, the conjugated projection associated with the transformed Christoffel--Darboux kernel differs from the original orthogonal projection by a rank-one operator on the subspace of functions vanishing at the origin. This provides a general mechanism for transferring local asymptotic results from an orthogonal polynomial ensemble to its Christoffel-transformed counterpart as the lattice size tends to infinity. As a main application, we study local fluctuations of random Young diagrams arising from skew Howe duality. The corresponding particle ensemble is obtained from the Krawtchouk orthogonal polynomial ensemble on a quadratic lattice by a Christoffel transform. We identify four asymptotic regimes of local fluctuations in the limit with . Besides the universal bulk fluctuations governed by the discrete sine kernel and universal Airy fluctuations at the right edge of the limit shape, we obtain the discrete Hermite kernel in the critical regime , and the discrete hard-wall sine kernel at the left corner.
Abstract
We prove that, for all , the Dobrushin-Lanford-Ruelle (DLR) equations for the unit-intensity point process have a unique stationary solution under a particular finite electric energy condition. Uniqueness fails if this condition is dropped. This answers a question of Dereudre-Hardy-Leblé-Maïda and gives a canonical statistical physics characterisation of . The main technical ingredient is a relative entropy estimate which allows us to compare the conditional distribution in a finite box of a DLR solution to a corresponding circular ensemble.
Abstract
We investigate a special ansats that allows for an iterative solution of the constant Yang-Baxter equation. Testing this ansatz, we construct four sequences of the constant R-matrices. In each sequence the R-matrices act on the tensor squares of vector spaces of linearly growing dimensions. Each R-matrix also depends on a single complex parameter. By analyzing the spectra of the R-matrices, we conclude that the first and the third series are associated with the symmetric tensor representations of the quantum groups U_q(sl(2)) and U_q(sl(2|1)), respectively. The other two series appear to be related to representations of quantum groups in the case where q is a root of unity. The elements of the second series are associated with nilpotent representations of U_q(sl(2)). We also check the correspondence between these R-matrices and link invariants. Applying the Reshetikhin-Turaev procedure to the first, second, and fourth sequences yields, respectively, colored Jones polynomials, ADO invariants, and Alexander polynomials. To our surprise, a similar calculation for the third series yields a trivial result, identically equal to 1.
Abstract
We consider the adjacency matrix of a random directed -regular graph on vertices. For fixed , we prove that the empirical eigenvalue density converges in probability to the oriented Kesten--McKay law as . The key technical input is the small-ball probability estimate for the smallest singular value. The proof combines a fixed-rank transposition argument with finite-field anticoncentration for shifted inverse compressions. We also prove a polynomial hard-edge estimate, which allows us to deduce the global law from the vanishing small-ball probability.
Abstract
We establish Gaussian free field fluctuations for uniformly random lozenge tilings of simply connected polygonal domains with sides whose directions cycle through the three lattice directions. More precisely, assuming that the liquid region is connected and that the boundary data do not force the height at any interior point, we prove that the fluctuations of centered height function converge to the Gaussian free field in the liquid region, confirming a prediction of Kenyon and Okounkov from 2007. We introduce a tiling action function that encodes the geometry of the limit shape through its critical points. The action function has a complex conjugate pair of critical points in the liquid region, repeated real critical points on the arctic boundary, and distinct real critical points in the frozen region. Using this tiling action function, we construct an approximation to the inverse Kasteleyn matrix in terms of explicit single-contour and double-contour integrals and prove that the approximation is uniform throughout the polygonal domain. The convergence to the Gaussian free field then follows from standard kernel computations.
Abstract
For every , we establish multi-time soft-edge moment convergence for the Jack--Plancherel process, a discrete -analogue of Dyson Brownian motion. The limiting moments admit an absolutely convergent expansion in terms of nonnegative Brownian bridges decorated by pairs of equal and opposite jumps. At equal times, we identify this limit with the joint Laplace statistics of the Airy process for . The same Brownian expansion yields an asymptotic -topological expansion of the marginal -conjecture and, at , an explicit nonnegative formula for the Witten--Kontsevich intersection numbers.
Abstract
We construct isomonodromy transformations for linear -difference systems of the form , where has diagonal leading coefficient. These transformations shift eigenvalues of the leading coefficient together with roots of . They lift compatibility to the right eigenpairs of , yielding a discrete local tau function. The resulting deformation equations preserve the Birkhoff connection matrix, and reduce in their limit to the isomonodromic deformation of a meromorphic connection on with an irregular singularity of Poincaré rank one at .
Abstract
We prove uniform endpoint estimates, killed local limits, and conditioned-bridge invariance principles for triangular arrays of centered, aperiodic, integer-valued random walks with uni- form exponential tails. The endpoint estimates cover all nonnegative heights, while the local limits treat positive diffusive and zero endpoints. The increment law may vary with the scaling parameter.
Abstract
We establish an identity relating the spatial derivative of a two-point free-energy correlation to two fundamental geometric quantities: the annealed probability that two independently sampled polymers, in the same environment, meet before reaching the prescribed terminal level, and the exit-point location of a single polymer. Our result is inspired in part by recent integration-by-parts work of Gu and Quastel for the KPZ equation. A crucial step in their continuous setting relies on an ingenious application of Itô's formula, which has no direct analogue in our semi-discrete setting. Instead, we exploit intrinsic symmetries of the polymer model together with the memoryless property. We also provide two applications: (i) we recover the horizontal Burke property by proving that the anchored stationary horizontal free-energy profile is a two-sided Brownian motion; (ii) in the zero-temperature limit, we obtain the corresponding identity for Brownian last-passage percolation.
Abstract
We consider beta-Dyson Brownian motion, with , started from a deterministic configuration with uniformly bounded support. Let be the semicircular free-convolution flow issued from the initial empirical measure, and set . For every fixed , with probability at least , every particle remains within an epsilon-neighborhood of for all . The result holds from time zero, requires no regularity assumption at the initial spectral edges, and applies to multi-cut supports with macroscopic interior gaps. A key ingredient is a deterministic local resolvent exclusion principle: an comparison of Stieltjes transforms on a complex disc above a real point separated from the reference support excludes eigenvalues from the corresponding real interval. This gives a model-independent mechanism for converting local resolvent estimates into spectral confinement.
Abstract
We study the two-dimensional directed polymer model in a Gaussian environment which is independent in time and spatially correlated, with covariances either summable or with a critical decay, satisfying as for some . We determine the precise high-temperature asymptotics of the free energy, confirming a conjecture of Lacoin (Ann. Probab. 2011), later refined by Cosco, Cottini and Donadini (2025). We also establish a phase transition for the diffusively rescaled partition functions: below some critical point they converge to the Lebesgue measure, while above it they converge to zero. A key feature of our approach is that both results are obtained using only second-moment estimates and are based on a simplified change-of-measure argument in the supercritical regime, that may prove useful for other disordered models.
August 2026
Abstract
We prove that the stationary density fluctuation field of the one-dimensional weakly asymmetric simple inclusion process converges to the energy solution of the stochastic Burgers equation. The key technical ingredient is a second-order Boltzmann-Gibbs principle tailored to the microscopic current. As a consequence, we also obtain Ornstein-Uhlenbeck fluctuations for the symmetric inclusion process in a weakly condensing regime, a result of independent interest.
Abstract
The Cardy-Smirnov function of a four-marked planar domain encodes the conjectural scaling limit of the crossing probability for critical percolation. Cardy predicted a closed formula for this limit on the rectangle in the cases of both site and bond percolation, and Smirnov proved the formula together with conformal invariance for critical site percolation on the triangular lattice. Kleban and Zagier later showed that on the rectangle the function admits a modular interpretation: it is determined by a modular functional equation together with a mild analytic ansatz. We carry out the analogue on the parallelogram. We derive a closed conformal-map formula for the Cardy-Smirnov function as an incomplete beta integral, prove a closed modular formula realizing it as an integral of the modular form , and establish a uniqueness theorem showing that a single functional equation, together with a -expansion ansatz and a nondegeneracy condition, determines the function uniquely.
Abstract
We consider a class of discrete multi-time determinantal point processes whose correlation kernels admit a double-contour integral form, containing Schur processes as simple examples. We show that the form of the correlation kernel is preserved under conditioning of the point process to a restricted domain, however, with an integrand which becomes more complicated and which is characterized by a Riemann-Hilbert problem. We apply our general result to domino tilings of reduced Aztec diamonds and to lozenge tilings of hexagons with holes. Our results show a striking analogy with the Its-Izergin-Korepin-Slavnov method for integrable kernels.
Abstract
Exact cumulants of entanglement entropies of random state ensembles have traditionally been studied within the random matrix framework. In this work, we propose an alternative approach based on the intrinsic connection to integrable systems. The central idea is to embed entropic quantities into tau functions satisfying Toda-type lattice equations, which in turn yield linear difference equations for their averages. Directly solving the difference equations recovers exact entropy formulas in the literature. The integrable systems approach bypasses the case-by-case, ensemble-dependent derivations required by random matrix methods. The approach also suggests a possible route towards unified and more efficient higher-order cumulant calculations by exploring integrable hierarchies.
Abstract
We study the continuum scaling limit of the Langevin dynamics for three-dimensional U(1) lattice Yang--Mills theory. The model is defined on the discrete 3D torus with a general class of plaquette actions that are suitably normalized, including Wilson, Manton, and Villain actions. Under the weak-coupling scaling and in the DeTurck gauge, we prove that, locally in time and in probability, the rescaled logarithmic field converges to the solution of the one-form stochastic heat equation. In particular, the limiting dynamics are universal and do not depend on the higher-order details of the plaquette action.
Abstract
Let be the dihedral group of order . Consider a continuous-time random walk on driven by arbitrary symmetric rates whose support generates . For , we say the pair is rate-monotonic if for each fixed time , the -distance between the random walk's distribution at time and the uniform distribution is monotonically decreasing as a function of the rates. Lyons and White proved that and are rate-monotonic. Somewhat counterintuitively, they found several pairs with that are not rate-monotonic, and they asked whether any such pairs exist with or . We resolve their question, proving that is rate-monotonic for all positive integers and . In fact, we prove a generalization of this result to a broader family of groups that includes generalized dihedral groups, dicyclic groups, and generalized quaternion groups. In the other direction, we prove that for every real that is not an even integer, there exists a positive integer such that is not rate-monotonic. The results of this paper were formally verified in Lean by AxiomProver assuming standard literature.
Abstract
We investigate random recursive simplicial complexes growing by adding, at each step, a vertex together with a simplex formed by joining the new vertex with a randomly chosen existing simplex. We also add all faces of the new simplex to ensure that the resulting object remains a simplicial complex. If the choice of an existing simplex is uniform among simplices of dimension , the number of simplices of any admissible dimension is an asymptotically self-averaging random variable. This feature allows us to determine the asymptotic growth law of the average of when the number of vertices diverges. We also probe the degree distribution, examine the probabilities of various extreme outcomes, and analyze the characteristics of the first vertex.
Abstract
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 , we show that these joint moments are characterised by a solution of the -Painlevé V equation for all real moment exponents in their admissible range. Under an appropriate large- scaling limit, we further prove that the limiting joint moments admit a representation in terms of a solution of the -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 -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.
Abstract
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.
Abstract
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.
Abstract
We investigate the asymptotic behavior of Jacobi-Piñeiro polynomials of degree orthogonal on with respect to weights , where , and . These polynomials are characterized by a Riemann-Hilbert problem for a 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.
Abstract
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.
Abstract
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.
Abstract
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.
Abstract
We perform a large- asymptotic analysis of the Fredholm determinant of an integrable integral operator with generalized sine kernel, where 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.
Abstract
We study a general class of gradient interface models with Hamiltonian , , assuming essentially that the potential is even, on for some , and . 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 . The construction of such a representation is one of the key new contributions of this work.
Abstract
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.
Abstract
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.
Abstract
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).
Abstract
We study Toeplitz determinants for whose Fourier coefficients satisfy . This regime extends beyond 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 . We show that this quadratic term accounts for the possible growth with , 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.
Abstract
While the solution to the 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 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 , up to a vanishing error, all the mass of the point-to-plane SHF in any domain is concentrated on balls of radius , each containing 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.
Abstract
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.
Abstract
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.
Abstract
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.
Abstract
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.
Abstract
We prove, in the case , the extension to the -symmetric world of the original Macdonald positivity conjecture. This is achieved by giving a combinatorial interpretation of the Kostka coefficients in terms of standard fillings of the diagram associated to the -partition . This interpretation generalizes the one in the usual Macdonald case, which is given by a major index statistic on standard tableaux.
Abstract
We study rooted Gibbs-DLR measures in the directed polymer model on 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 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.
Abstract
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.
Abstract
We study the nested permutation classes . We relate the corresponding generating functions , to -difference equations. Solving these equations allows to show that the generation functions are not D-finite, and that , .
Abstract
In this paper, we investigate the monic orthogonal polynomials and the Hankel determinants 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}\}\), 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 satisfies the -th member of the discrete Painlevé I hierarchy. Besides, we derive the second-order differential equation satisfied by , 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 . Based on the discrete Painlevé I hierarchy and the above differential identities, we obtain new partial differential equations satisfied by and .Using the discrete Painlevé I hierarchy and the asymptotic theory of linear difference equations, we derive the full asymptotic expansions of the recurrence coefficient , the nontrivial leading coefficient , and the Hankel determinant as , for general and . Notably, while the logarithmic term 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 .
Abstract
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 , 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 we obtain closed-form expressions for the specialization.
Abstract
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.
Abstract
We prove universality results for zeros of random -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 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 . 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 -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.
Abstract
We provide limit theory for the trace of the squared sample correlation matrix , constructed from observations of a -dimensional random vector with iid components. If the entries have finite fourth moment and and grow proportionally, it is known that 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 and any growth rate of the dimension, we prove that the universal CLT remains valid for . For , we identify a critical dimension growth at which the fluctuations of become non-Gaussian. Moreover, if the dimension grows faster and we establish a non-universal CLT with norming sequences depending on the value of . Our findings are illustrated in a simulation study.
Abstract
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.
Abstract
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 and disorder cumulants with distinct power-law dependences on the density , 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.
Abstract
We study partitions which attain the largest Littlewood-Richardson coefficient. More precisely, we prove that the largest is attained at partitions such that and is a disjoint union of squares. We conjecture that for , all largest must satisfy this property. We confirm this conjecture numerically, for .
Abstract
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 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--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 case.
Abstract
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 point process. Our construction and scaling limits extend in a natural way to non-classical .
Abstract
For ordinary permutations on letters, the distribution of the number of fixed points of random permutations is well known to approach the Poisson distribution in total variation distance as 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 letters whose expected number of fixed points converges to a constant , we must have and the distribution of number of fixed points converges to the distribution as . If , then the rate of convergence in total variation distance can be as slow as .
July 2026
Abstract
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.
Abstract
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 -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 KPZ benchmarks , . Our results identify three-dimensional polariton crystals as a controllable quantum fluid route to higher-dimensional nonequilibrium universality.
Abstract
We introduce and study the perturbed -corners process, a deformation of the classical -corners process. The latter is a probability measure on interlacing arrays of real numbers that, for the classical values , describes the joint distribution of eigenvalues of principal submatrices of a Gaussian random matrix with the corresponding GOE/GUE/GSE symmetry. For classical , the perturbed process arises from adding a deterministic diagonal matrix to a GOE/GUE/GSE matrix. The eigenvalues of the resulting random matrix depend symmetrically on , but this symmetry does not extend to the whole corners process. We extend the construction to all via multivariate Bessel functions, and analyze the crystallization () 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 's grow linearly in . 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 followed by iterated ordinary derivatives.
Abstract
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.
Abstract
We give explicit formulas for the smallest non-trivial irreducible representation of quantum affine superalgebras in types D, (in the all-fermionic parity), F, (in the distinguished parity), and G, (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 -matrices.
Abstract
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.
Abstract
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.
Abstract
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 and . 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 , we construct type 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 specialization.
Abstract
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 , 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.
Abstract
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 -Mehler formula as special cases. Some new -series identities related to the multivariable Al-Salam-Carlitz polynomials are also presented.
Abstract
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 . 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 , that computes the steady-state joint occupation distribution for any set of vertices; since the steady-state distribution is very elusive, this provides a potentially useful tool to study its -wise correlations for small values of . 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.
Abstract
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.
Abstract
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 -symmetric Dirichlet form: the collision set is polar if and only if . This provides an infinite-particle counterpart of the classical collision threshold for the finite-particle Dyson Brownian motion with inverse temperature .
Abstract
KPZ-type extremal fluctuations have recently been proved for several models of random walks in space-time random environments (RWRE) in 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 -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.
Abstract
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 , further efforts are necessary to better understand these connections in dimensions . In this paper, we study this problem numerically for and , focusing on a discrete model with Dirichlet distributed transition probabilities. This model is a generalization of an integrable model in , 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 in time , is compatible with KPZ growth in dimension and . In spatial dimension , we confirm the existence of a phase transition as the angle 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 and that it is related to the extreme diffusion coefficient introduced recently.
Abstract
The KP equation is a prototypical -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.
Abstract
We study two classes of oscillatory Jacobi-type weight functions , . Other restrictions are imposed on and 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.
Abstract
Determinantal point processes (DPPs) are widely used as probabilistic models for diverse random subsets, but their approximation error under model misspecification has not been fully characterized. We study the population-level approximation of a strictly positive target distribution p* by DPPs under the forward Kullback-Leibler divergence. Using information-geometric analysis and the standard quality-diversity decomposition of an L-ensemble kernel, in which the diagonal quality component Q encodes item-specific weights and the diversity component D controls repulsive interactions among items, we show that the quality component can be chosen uniquely to match all first-order inclusion probabilities of p*. The DPP approximation problem therefore reduces to the optimization of the diversity component. This reduction yields a global optimality result for attractively dependent targets: under conditions including weak positive association, the independent product distribution with the same first-order marginals, corresponding to D=I, is an optimal DPP approximation. For more general target distributions, positively correlated pairs yield lower bounds on the approximation error. On the repulsive side, a matching of disjoint negatively correlated pairs yields an upper bound on the approximation error, or equivalently a guaranteed improvement over the independent approximation. We further study local optimality around D=I and, more generally, around block-diagonal diversity matrices by analyzing perturbations between their blocks.
Abstract
Disordered systems such as spin glasses and polymers characteristically exhibit random energy landscapes with many macroscopically separated energetic valleys corresponding to near-ground states. This high complexity renders these systems extremely sensitive to perturbations of external parameters. For instance, the support of associated Gibbs measures may change macroscopically under such perturbations, a phenomenon known as chaos in the literature. In experiments, chaotic phenomena are typically studied via temperature perturbations. In this article, we initiate the rigorous study of temperature-chaotic properties of the continuum directed random polymer (CDRP), a canonical model in the KPZ universality class. The CDRP is driven by white noise and is parametrized by inverse temperature , and is known [Wu '26, Das-Zhu '24] to converge in the zero-temperature limit to the directed landscape constructed in [Dauvergne-Ortmann-Virág '22], the putative universal scaling limit of models in the KPZ universality class. The main result of this article considers the CDRP free energies coupled through the same white noise at a pair of inverse temperatures , and shows that they decouple in the limit , converging to a pair of independent directed landscapes. This is the first such "energetic de-correlation across temperatures" result. Our key estimate measures the "pivotality" or "influence" of spatially thin strips in models of last passage percolation. As a byproduct, the proof strategy also allows to show that the directed landscape is a two-dimensional black noise (in the sense of [Tsirelson-Vershik '98]), previously conjectured by Virág. This provides the third known example of a two-dimensional black noise after critical planar percolation [Schramm-Smirnov '11] and the Brownian web [Ellis-Feldheim '16].
Abstract
We construct the periodic directed landscape, which is the conjectured scaling limit for periodic models in the Kardar-Parisi-Zhang universality class. We establish the convergence of periodic exponential last passage percolation to the periodic directed landscape. Moreover, we confirm conjectures by Baik, Liao and Liu on the local structure of the periodic KPZ fixed point, and establish the convergence of periodic ASEPs to periodic KPZ fixed points, coupled according to the same periodic directed landscape. Our main tool to construct the periodic directed landscape (and prove convergence to it) is a technique for gluing full-space directed landscapes (and their prelimits), which is of independent interest.
Abstract
We study stationary measures of the log-gamma polymer on a finite diagonal strip and in a half-space. We establish the phase diagram for the stationary measure on the strip. To this end we develop a representation of {the Laplace transform of} this stationary measure by independent random variables. We also present an analytic approach that extends Barraquand's contour integral representation of the Laplace transform. We prove that, as the strip width tends to infinity, the stationary measure of the log-gamma polymer on the strip converges to the stationary measure of the half-space log-gamma polymer. Finally, we derive a contour integral formula for the Laplace transform of the stationary measure of the half-space log-gamma polymer.
Abstract
We give a Graham-positive expansion for the product of two double Schubert polynomials indexed by two inverse Grassmannian permutations. Surprisingly, the nonzero structure constants are double Schubert polynomials in two disjoint sets of equivariant variables. We also give a positive expansion for the product of two single Schubert polynomials indexed by a -avoiding permutation (e.g., a Grassmannian permutation) and an inverse Grassmannian permutation. Unexpectedly, the nonzero structure constants are Edelman--Greene coefficients.
Abstract
A random sequence sampled from a known continuous distribution is observed with the objective to choose an item with the overall rank one. A rejected item cannot be recalled and is immediately erased from the memory. Under this memory constraint, the choice problem is not amenable to recursive methods of optimal stopping and becomes a global optimisation task. We focus on a heavy-traffic form of the problem with infinitely many choice opportunities, which we state in terms of a planar Poisson process (PPP). Symmetries of the PPP are used to derive basic structural properties of the optimal stopping rule, including the balance at the boundary equation, and two key integral identities. Throughout, we make throrough comparison to the classic full-information counterpart of the problem, revisiting both discrete- and continuous-time models. The optimal value, stopping rule and other characteristics of the problem are determined analytically and approximated numerically with high precision.
Abstract
We study the class of functions counting accepting paths of non-deterministic log-space Turing machines and construct methods to prove containment in . We prove that a large number of classical combinatorial and number theoretic functions belong to this class: classical functions from enumerative combinatorics (multinomial coefficients, Catalan numbers, linear extensions of trees, Stirling numbers, etc), algebraic combinatorics (number of standard Young tableaux, etc), discrete geometry, number theoretic functions, representation theoretic multiplicities in a large class of cases. We show that -plethysm coefficients of bounded length outer partition can be counted by log-space polytime verifiers. We pose numerous questions and conjectures on containment and its generalizations, that suggest venues for conditionally disproving -completeness. While studying which combinatorial functions are in provides a formal way of (dis)proving the existence of combinatorial interpretations, the lower class serves as an analogue for functions computable in polynomial time.
Abstract
We study the vector Dyson equation with parameter in the complex upper half-plane , where and is a nonnegative matrix, not necessarily symmetric. This equation has a unique vector solution , for which we establish a complete measure decomposition and prove regularity. We then develop a graph-theoretic approach to the singularity and stability problem for non-symmetric matrices . The graph structure of identifies the possible degeneracies of the stability operator as approaches the real axis. In particular, for non-backtracking matrices, we prove square-root growth at regular edges, cubic-root growth at regular cusps, and complete stability estimates. We also obtain the corresponding estimates for symmetric matrices in the periodic setting.
Abstract
This article associates to each positive real number a dimer model on a snake graph with activity parameter whose structure is determined by the continued fraction expansion of . When is rational, the model is finite and gives rise to a probability measure on perfect matchings. For irrational , the model is infinite, and is defined as a limit over rational approximations to ; the main technical result of the paper shows that this limit is well defined and independent of the choice of rational approximants. denotes the odds that a -random perfect matching includes a distinguished edge. When is rational, , where is the algebraic -deformation introduced by Morier-Genoud and Ovsienko. This agreement, together with evidence from the irrational case, suggests a close connection between the probabilistic and algebraic constructions.
Abstract
We study energy transport in a finite one-dimensional unpinned harmonic chain with stochastic nearest-neighbor momentum exchanges and Langevin heat baths at its endpoints. Such systems are known to exhibit superdiffusive transport driven by long-wavelength acoustic modes, leading to fractional macroscopic behavior. While fractional heat equations have been rigorously derived for infinite chains, the corresponding boundary conditions for finite systems in contact with heat baths remain unclear due to the nonlocality of the fractional Laplacian. Under the superdiffusive time scaling , where is the system size, we prove that the averaged microscopic energy profile converges, as , to a temperature field solving a fractional heat equation on , with the generator given by a Neumann fractional Laplacian and additional nonlocal boundary terms induced by the heat baths. Our results provide a rigorous derivation of macroscopic boundary conditions for superdiffusive heat transport in open chains and introduce new boundary conditions for fractional Laplacians, that are motivated by a physical model.
Abstract
We give an algebraic branching program of polynomial size which computes Cayley determinant of right quantum matrices. This is a rare example of an efficient computation of a noncommutative determinant, and the first such example for quantum groups. We extend the results to the -Cayley determinant of -right quantum matrices, as well as to their multiparameter generalization. The proofs are entirely combinatorial, as we relate Cayley, Moore and Valiant determinants using bijections/involutions on words. We then employ the celebrated determinant construction of Mahajan and Vinay (SODA'97), to obtain the results.
Abstract
Non-stationary time evolution of interacting particle systems is in general a rather difficult topic however, exceptional examples are known where hidden processes within the model make the description manageable. These hidden processes often take the form of a finite number of interacting random walks and in some cases, rather than being hidden, were very explicitly revealed as second class particles associated with the model. The examples of asymmetric exclusion and exponential bricklayers model are known in this context, where such distributional structure was demonstrated in infinite volume. Here we find boundary mechanisms that save this remarkable structure in finite volume of the model with semi-open boundaries that let ordinary particles, but not the shocks, through. This finding also allows to characterise nontrivial two-species, non-reversible stationary distributions subject to our special boundary rates.
Abstract
The spiked mixture model (SMM) has been introduced as a probabilistic model that generalizes the single-spike (Wishart) model to a mixture model form. With applications ranging from imaging mass spectrometry in the life sciences to hyperspectral imaging in computer vision, it is crucial to understand under which circumstances its signals can be recovered from noisy measurements. The highly multiplexed nature of these measurement types furthermore necessitates such analysis to hold in highdimensional settings. In this paper, we prove that the extreme eigenvalues of the covariance matrix from the SMM exhibit a phase transition in high-dimensional regimes. We show that this phase transition, and thus signal recovery by extreme eigenvalues, depends on several interacting factors: the correlation between spikes (i.e., how similar in content underlying signals are), the energy parameters (i.e., the absolute strength of each underlying signal), and the mixture probabilities (i.e., how likely it is to encounter each underlying signal). This work provides sharp information-theoretic bounds on the parameters needed to detect one or more spikes from extreme eigenvalues of the SMM covariance matrix, and these guarantees could potentially impact any application of the SMM. Understanding this interplay could serve as a tool for driving experimental design in analytical chemistry and life sciences.
Abstract
We compute the dimension of the center of the Temperley--Lieb algebra over a field of characteristic zero for every nonzero value of the parameter . The proof uses the cellular filtration by cup number, together with known facts about the representation theory of the Temperley--Lieb algebra, especially the structure of its standard modules and their radicals. Dilation and compression maps compare the induced graded pieces of the center at levels and , giving an upper bound of one for each such piece. A deformation argument gives the matching lower bound, and hence . We also prove that every central element is fixed by the canonical anti-automorphism and by the natural diagram-reflection automorphism. Finally, we give a congruence criterion for the trivial-radical case and record a Gram-matrix computation for leading terms.
Abstract
The type D ASEP is an asymmetric two--species interacting particle system on , in which two separately conserved species hop, bind into a composite ``bound pair'', and split. The model, along with its reversible measures and orthogonal polynomial duality, was constructed using the representation theory of . The reversible measures and orthogonal polynomial duality are each a product of two copies of the single-species ASEP reversible measures and orthogonal polynomial duality. In this paper, we study the long-time asymptotics of the type D ASEP. In the fixed-- regime, using an exact current--decoupling identity, we prove that the asymptotic hydrodynamic limit and Tracy--Widom fluctuations decouple, as predicted from the duality. In the weak--asymmetry (Edwards--Wilkinson) regime, when , we prove that the two density fluctuation fields \underline{decouple}: each converges to a linear stochastic heat equation, with no cross--coupling in either the drift or the noise, the limiting noises having vanishing cross--correlation. More surprisingly, we then prove that the two limiting normal random variables are \underline{correlated} with a seemingly new correlation function. The correlation is exactly equal to , with the positive parts of the normal random variables having correlations expressed by the Bessel--Struve function. This paper, with the exception of the abstract and introduction, was written entirely by Claude Opus 4.8 and Fable 5. The proofs were then formalized in Lean, using Aristotle by Harmonic AI. The human author of this paper verified the proofs manually.
Abstract
We study the global fluctuations of the height function associated with the Airy line ensemble. Using its determinantal structure and a steepest-descent analysis of the extended Airy kernel, we prove that, after a suitable rescaling, the centered height function converges to an explicit pullback of the Gaussian free field. The convergence holds in the sense of joint moments of linear statistics against compactly supported continuous test functions.
Abstract
We prove Krattenthaler's conjecture from 2001 about the -phenomenon for lozenge tilings of semiregular hexagons. In a first step we reduce the problem to the case of regular hexagons. In a second step we further reduce the question to a special case already covered in the literature. This is achieved by vast application of the celebrated Zeilberger Algorithm and the Holonomic Ansatz.
Abstract
We provide a definitive analysis of the number of correct guesses in the complete-feedback card guessing game after an asymmetric single-shelf shuffle with parameter . We explicitly describe the optimal strategy that maximizes the expected number of correct guesses. We study the number of correct guesses, under an optimal strategy, using methods from analytic combinatorics. In addition, we find the explicit distribution for the number of correct guesses, and thus find the mean and the variance for the number of correct guesses. We show that the distribution of the number of correct guesses is log-concave. We also study the limiting behaviour of the number of correct guesses as the number of cards goes to infinity. In particular, we prove a (local) central limit theorem and a large deviation principle with an explicit rate function. We also prove phase transitions for the number of correct guesses near and . Prior to this work, only the optimal strategy and the expected number of correct guesses under the optimal strategy were known for the special case of .
Abstract
We study the weakly coupled elliptic Anderson model with spatial white noise on the four-dimensional torus, which provides a basic example of a critical SPDE requiring renormalization at arbitrarily high orders. With coupling where is sufficiently small, we prove that the Green's function of the corresponding random Schrödinger operator, suitably centered and rescaled, converges to a centered Gaussian random field with explicit covariance. The main difficulty is that, for such critical models, one must expand up to order , while the perturbative expansion contains factorially many pairings and a growing number of renormalization terms. To overcome this, we construct a truncated renormalized parametrix and prove sharp high-order bounds for its remainder. A central ingredient is a multiscale analysis based on a new version of Hepp trees, combined with new estimates for summations over permutations. These estimates reveal a precise balance between logarithmic losses from scale summation and factorial gains from the structure of primitive pairings. The methods developed here are intended as a first step toward a general theory for critical SPDEs with weak couplings.
Abstract
We study a -dimensional directed polymer in a spatially correlated random environment generated by power-law tail variables: , where the variables are i.i.d. and have a regularly varying right tail with exponent . The spatial covariance of the environment has long-range decay with Hurst parameter . We identify the limiting fluctuations of the log-partition function in the intermediate disorder regime and show that the critical tail exponent is . When , the model has the same scaling limits as the corresponding Gaussian spatially correlated polymer: if , the centered log-partition function converges to the logarithm of the solution of the stochastic heat equation driven by fractional spatial noise; if , its normalized fluctuation converges to a centered Gaussian law. In the regime , at the scale , the log-partition function still satisfies Gaussian fluctuation. The main ingredient is a truncation comparison argument adapted to long-range moving-average environments, together with an invariance principle for polynomial chaos. Due to the non-locality of the environments, we perform a far-near field analysis, as well as multiscale analysis, to prove that the truncated version does not change the log-partition function at the corresponding scales.
Abstract
Recently, there has been significant interest in applying the method of moments developed by Wood and others to study distributions of finite abelian groups that arise in number theory and combinatorics. When the moments do not grow too fast, they determine a unique distribution. We construct large families of distributions that have the same moments. These families include several distributions that arise naturally in the study of sandpile groups of families of random graphs. Wood determined the distribution of Sylow -subgroups of sandpile groups of Erdős--Rényi random graphs. This was extended by Mészáros to sandpile groups of random -regular graphs, who observed an interesting special case when is even and . We study Sylow -subgroups of sandpile groups of random bipartite graphs and similarly find a special case for . Although this distribution differs from that of Mészáros, we show that they have the same moments and fit into our broader construction. To compute the moments of the distributions we study, we apply combinatorial tools from the theory of Hall--Littlewood functions.
Abstract
We show that the lower edge of the appropriately scaled size Laguerre beta-ensemble with parameter converges to the process as when and . This completes the picture of the possible edge scaling limits of the Laguerre beta-ensemble with a fixed . When our proof establishes operator level convergence of the inverse of the scaled Dumitriu-Edelman tridiagonal matrix to the inverse of the stochastic Airy operator. Our methods allow us to prove similar operator level limits for the known soft edge scaling limits of the Laguerre and Gaussian beta-ensembles. For we give a different argument that relies on coupling and a result of Dumaz-Li-Valko for the transition between the hard and soft edge limits of the Laguerre beta-ensemble.
Abstract
We prove that the universal local point processes of random matrix theory are characterized by their loop equation hierarchies. More precisely, for every rational , the point process is the unique solution of the bulk loop equation hierarchy, and the point process is the unique solution of the edge loop equation hierarchy. These uniqueness results provide a direct route to universality: it suffices to verify the corresponding approximate loop equations for the ensemble. In many models, these equations follow from local laws and integration by parts.
Abstract
We study non-exchangeable interacting diffusions with pairwise interaction strengths encoded by a sequence of matrices. Under suitable structural and denseness conditions on these matrices, we prove a universal Central Limit Theorem for the global fluctuation field. As the number of particles becomes large, it converges in distribution to the unique solution of a stochastic partial differential equation (SPDE), the same Gaussian limit as in the exchangeable mean field case. The result applies, for instance, to scaled adjacency matrices of -regular graphs when . A spatial interaction model shows that the denseness threshold is sharp. The proof proceeds with an analysis in negative Sobolev spaces, building on sharp quantitative propagation of chaos results together with functional inequalities.
Abstract
Given any , let be a random matrix over , with all entries independent and -balanced (lying in each residue class mod with probability at most ). Wood proved that as the distribution of approaches Cohen and Lenstra's conjectured distribution of class groups. Given such that , we prove that the distribution of still approaches the Cohen--Lenstra distribution even if we weaken the hypothesis by allowing up to entries per column and up to entries per row of to not be -balanced. We also weaken the independence condition by allowing certain types of dependence between the entries of each column. In addition, we prove that, for any , the cokernels of random band matrices of width with -balanced entries in the band and arbitrary entries outside of it will also approach the Cohen--Lenstra distribution, which answers a question of Kang--Lee--Yu.
Abstract
We study invariant measures for soliton systems described by Mealy automata. Motivated by recently introduced soliton models associated with 2-letter, 3-state Mealy automata, we formulate the time evolution induced by Mealy automata on bi-infinite configuration spaces. We provide sufficient conditions for the invariance of Bernoulli product measures and derive a criterion for the invariance of two-sided space-homogeneous Markov distributions. We then apply these general results to three soliton models, which can be interpreted as variants of the box-ball system (BBS). For two of these models, BBS-S(2) and BBS-V(2), we prove that Bernoulli product measures are invariant. For the remaining model, BBS-C(2), we establish a more general result: the invariance of two-sided space-homogeneous Markov distributions, which include Bernoulli product measures as a special case. Furthermore, for all three models, we compute the phase shift associated with the interaction of two solitons, as well as the velocity of an isolated soliton. Although the latter has already been studied previously, both quantities constitute fundamental characteristics for understanding the generalized hydrodynamics of these systems. These results provide a foundation for the study of invariant measures, generalized Gibbs ensembles, and generalized hydrodynamic behavior in Mealy-automaton soliton systems.
Abstract
Let be an real-valued random matrix with independent, mean-zero, variance-one entries whose fourth moments are uniformly at most . Suppose that there exists such that the entries of satisfy We prove that there are constants , depending only on and , such that for every fixed invertible matrix and every , In the Gaussian case, we also show that the above estimate is sharp in the sense that
Abstract
We give a new proof of the celebrated theorem of Knutson and Tao that the spectra of triples of Hermitian matrices exactly correspond to positions of boundary rays of honeycombs. Most importantly, our proof gives new insights into why honeycombs are related to Hermitian matrices in the first place. Our proof is axiomatic: We distill four essential properties shared by honeycombs and spectra of Hermitian triples, and show that any two objects sharing these four properties must be equivalent. In this way, we argue that honeycombs are `model organisms' for Hermitian triples: they are families of objects satisfying the same defining properties, but in more obvious ways.
Abstract
Determinantal point processes have recently emerged as a kernel-based alternative to standard independent sampling for constructing efficient minibatches, coresets, and other compact representations of large-scale datasets. In particular, sampling mechanisms based on DPPs are believed to demonstrate better approximation properties compared to classical i.i.d. samplers, even at the scale of the exponent. One of the key strengths of DPP based samplers is that they can be deployed over very general spaces, in contrast to more classical sampling methods beyond i.i.d. which tend to work in very well-structured settings, principally Euclidean spaces. In this work, we establish explicit rate guarantees for determinantal sampling in spaces that extend far beyond known Euclidean setups, focusing on spectral kernels obtained from eigenspaces of naturally associated Laplacian and other Markov diffusion operators. This includes, in particular, Riemannian manifolds and weighted networks. In determinantal sampling from compact Riemannian manifolds, we establish sampling rates that automatically pick up the intrinsic dimensionality of the underlying manifold. In the setting of networks, we investigate DPP-based samplers on the celebrated k-nearest neighbour graphs, as well as weighted random geometric graphs, and demonstrate a similar improved dependence on the intrinsic dimensionality of the data. Overall, our approach achieves guarantees of that match known rates on Euclidean spaces of comparable dimension. In terms of techniques, we connect to the celebrated Weyl's Law for manifold spectra, and leverage tools from the theory of Markov diffusions and Dirichlet forms as well as certain ingredients from the theory of pseudodifferential operators, which could be of independent interest in this area.
Abstract
We prove equality conditions for the Ahlswede--Daykin (AD) inequality and the Fortuin--Kasteleyn--Ginibre (FKG) inequality. We then present a number of applications and special cases of these equality conditions. These include Björner's and Fishburn's inequalities for linear extensions of finite posets, the Lam--Postnikov--Pylyavskyy (LPP) and the Okounkov inequalities for Schur positivity of products of Schur functions. We conclude with equality conditions for the Ahlswede--Daykin--Schur (ADS) inequality recently introduced in Chan--Chen--Pak--Soskin (2026), which is an AD type extension of the LPP inequality.
Abstract
We prove a mesoscopic central limit theorem for linear eigenvalue statistics of correlated Hermitian random matrices. The class considered here includes Wigner and Wigner-type matrices, as well as models whose entry correlations decay polynomially in the distance between index pairs. The proof combines a multivariate cumulant expansion with multi-resolvent local laws and a detailed analysis of the resulting variance kernel on the operator-level.
Abstract
We study a bi-infinite system of interacting Brownian particles on the real line with singular asymmetric interactions mediated by the collision local times. Particles perform Brownian motions, and when neighboring particles collide, the associated local time is split in proportions and . We first develop well-posedness theory for the particle system, proving pathwise uniqueness and strong existence under natural growth assumptions on the initial configuration and local times. We also identify a family of stationary distributions for the infinite-dimensional process of gaps between successive particles: for every , the product measure with i.i.d. Exp gaps is invariant. Our main result concerns the equilibrium fluctuations of the associated particle-count (height) function in a weakly asymmetric regime. Taking so that , initializing the interparticle gaps with i.i.d. Exp random variables, and applying a microscopic Hopf-Cole transform to the diffusively rescaled count function, we prove convergence, as , to the multiplicative stochastic heat equation (SHE) with Brownian exponential initial data. Equivalently, the logarithm of the limit is the Hopf-Cole solution of the KPZ equation with two-sided Brownian initial data. The proof combines localization to finite particle subsystems via chains of collisions, Brownian last-passage percolation estimates, a key local time cancellation property, and a martingale problem for a scale-adapted mollification of the Hopf-Cole field, whose space-time regularity is tuned to match that of the limiting SHE. The resulting fluctuation theorem places these Brownian particle systems with asymmetric singular collision dynamics within the KPZ universality class.
Abstract
The goal of the minimal model program for the Hilbert scheme of points on a surface aims is to describe the (stable) base loci of all divisors, their associated birational models, and the maps between them. We answer all of these questions for the Hilbert scheme of points on the blowup of the affine plane at the origin. The birational models are Brill-Noether loci in a larger Hilbert scheme, and nested variants thereof, and the wall-crossing maps are described as explicit projections. We also establish several new facts about the homogeneous coordinate ring of this Hilbert scheme, including finding the minimal set of line bundles whose sections generate the ring.
Abstract
We prove that the number of lozenge tilings of a certain triangular region is given by the formula \[T_n=\prod_{\substack{1\leq a<b\leq 3n+2\\(a,b)\not=(n+1,2n+2)}}\left|1+ζ^a+ζ^b\right|^{1/3},\] where . This answers a question of Ciucu and Krattenthaler, both by finding the exact formula and by explaining why has many prime factors. The proof reduces the lozenge tiling enumeration problem to evaluating the determinant of the bipartite adjacency matrix of the dual graph of , and then evaluates this determinant by diagonalising .
Abstract
We present an embedded (in ) degeneration of Lusztig varieties (which generalize type Hessenberg varieties) to certain unions of Richardson varieties, giving a simultaneous reproof (and extension) of results of Anderson--Tymoczko, Harada--Horiguchi--Masuda--Park, and Kim. Although torus-equivariant, the degeneration is not Gröbner. In the case that the Lusztig variety is the permutahedral toric variety, this degeneration provides a subdivision of the permutahedron into Bruhat interval polytopes, and we prove a more general result showing equivariant degenerations of projective toric varieties produce subdivisions of the moment polytope (as was shown in the Gröbner case by Sturmfels). A Gröbner degeneration would result in a {\em regular} subdivision, and despite our degeneration not being Gröbner we show in types that our subdivisions of the permutahedron are indeed regular.
Abstract
Bergeron, Garsia, Haiman and Tesler conjectured in 1999 that, for all partitions , the polynomial has nonnegative integer coefficients, where is the Bergeron--Garsia nabla operator, which acts diagonally on the modified Macdonald basis, and is the monomial symmetric function. In this article, we prove this conjecture, and more generally that for all . We establish a recursion showing that has an expansion with coefficients in in the symmetric functions , where denotes the operator introduced by Haglund, Morse and Zabrocki. Combining this expansion with the compositional shuffle theorems of Carlsson--Mellit and Mellit, and with the Schur positivity of LLT polynomials, completes the proof. The same method, using the -positivity of column LLT polynomials after the substitution , also gives an -positive analogue.
June 2026
Abstract
In previous work \cite{GW}, we developed a theory of modulated \(2j-k\) bi-orthogonal polynomial systems \(\{P_n(z;r),Q_n(z;r)\}\) and \(j-2k\) bi-orthogonal polynomial systems \(\{R_n(z;r),S_n(z;r)\}\), which generalize the classical \(j-k\) Toeplitz systems. In the present paper, we further develop this theory in several directions. We derive simplified and unified recurrence relations for both families of polynomials, prove a more transparent Christoffel--Darboux formula, and give Riemann--Hilbert characterizations of the \(2j-k\) and \(j-2k\) systems.
Abstract
We establish the explicit leading order asymptotics, with a quantitative error bound, of tail probabilities of collision times for a class of integrable space-time inhomogeneous Markov chains, in discrete and continuous time. The corresponding process conditioned not to intersect arises in interacting particle systems with local push-block interactions thereby confirming a recent prediction. The generic discrete nature of the spatial inhomogeneities rules out powerful coupling-with-Brownian-motion techniques, so our proof strategy proceeds instead via a novel steepest-descent analysis combined with a Karlin--McGregor semigroup expansion in terms of dominant-index contributions.
Abstract
We study the macroscopic fluctuation theory (MFT) of a multispecies generalization of the symmetric simple exclusion process (mSSEP) on the infinite line, in which particles of species -- including, possibly, a vacancy species -- exchange positions at unit rate. Working with the full, redundant set of coarse-grained densities $\bfrho=\{\rho_0,\dots,\rho_N\}$ keeps the relabelling symmetry of the model manifest throughout. We first extend the argument of Derrida and Gerschenfeld to the multispecies setting, showing that the cumulant generating function of the multispecies current between two regions of an arbitrary graph depends on the boundary densities $ \bfrho_L,\bfrho_R$ and on the fugacities $\bflambda$ only through a single scalar variable . We then formulate the MFT saddle-point equations for the mSSEP on the infinite line and show that they define an integrable system: they are of Landau--Lifshitz type, and a Zakharov--Takhtajan gauge transformation recasts them in AKNS form. Solving the resulting linear scattering problem by the inverse scattering method, we recover the cumulant generating function for the multispecies current, as well as the initial and final density profiles conditioned on a prescribed current fluctuation. In particular, we show that coincides with the function obtained by Derrida and Gerschenfeld for the single-species SSEP, now derived for an arbitrary number of species directly from the integrable structure of the multispecies MFT equations.
Abstract
We prove that the differential of the log of the partition function for the -matrix model with quartic interactions converges in a certain double-scaling regime to the differential of the -function for the string equation. This confirms the convergence of the critical Ising model on random surfaces to the topological minimal model, which was stated in the works of Douglas and Shenker, Brézin and Kazakov, and Gross and Migdal. Our analysis is based on a steepest-descent analysis of a Riemann-Hilbert problem associated to a family of biorthogonal polynomials. New features in the matching problem in the construction of local parametrices appear.
Abstract
Let be a whole plane Gaussian free field, and let be a bounded domain in two dimensions. We study the asymptotics as of the Liouville quantum gravity (LQG) heat trace, defined as the integral over of the on-diagonal LQG heat kernel. Our main result is to show that the second term in the spectral expansion as of the expected heat trace is governed by a nontrivial exponent, given by the KPZ (Knizhnik--Polyakov--Zamolodchikov) relation. A similar but stronger (almost sure) result applies to the related notion of heat content. Along the way we obtain various results on the short-term behaviour of the heat kernel, notably solving a conjecture of \cite{BW} concerning its annealed asymptotics, and showing the finiteness of all moments of the properly rescaled heat kernel.
Abstract
For every , we prove global existence for the two-dimensional generalized parabolic Anderson model on the whole space with nonlinearity , driven by an enhanced noise . Here the noise has polynomially weighted spatial Besov--Hölder regularity , and is the corresponding renormalized second-order object. If is globally Lipschitz, the solution is unique. The proof combines a weight-compatible annular high-low decomposition with a paracontrolled transport representation. The final remainder is estimated simultaneously in weighted norm and in a higher order weighted parabolic Hölder norm, with strictly different polynomial weights. This weight gap absorbs the polynomial losses arising from the enhanced noise and the transport coefficient. Several key refinements of the earlier works allow the maximum-principle and Schauder estimates to yield a global a priori bound for a larger range of . Uniqueness is proved in a time-dependent exponentially weighted topology.