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)
September 2026
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.
Abstract
The generalized Yang-Baxter equations are multi-site versions of the standard Yang-Baxter equation. When spectral parameters are included, such equations are expected to lead to integrable Hamiltonians with local interactions involving multiple degrees of freedom. In this work we characterize both the continuous and discrete symmetries of these equations required to establish an equivalence class of solutions. We find that the set of such symmetries depend on the number of sites on which the equation is supported. In several cases there are more symmetries than the standard Yang-Baxter equation, thus placing heavy constraints on the number of inequivalent solutions and the associated integrable models.
Abstract
We prove an off-diagonal expansion of the kernel of the Toeplitz operator whose symbol is the indicator function of a compact domain with smooth boundary in a complete symplectic manifold of bounded geometry. Using our approach, we extend two results to the non-compact setting: the first concerns the asymptotics of the trace of polynomials in this operator, and the second establishes a Weyl law for this Toeplitz operator.
Abstract
This paper is motivated by a recent result of Pitman and Yakubovich stating that a partially exchangeable, stationary (PES) random partition of is exchangeable. This echoes an earlier theorem of Kallenberg on the equivalence of spreadability (contractability) and exchangeability for infinite partitions. We revise the hierarchy of symmetries with a focus on partitions of finite sets , and ask about the extent to which these relaxed symmetry types differ from exchangeability, framing the question in terms of the geometry of the polytope of distributions defined by the symmetry constraints. We show that single-orbit exchangeable partitions remain extreme among spreadable and PES distributions, and that is the first case where non-exchangeable extreme partitions occur, causing the associated polytopes to deviate from a simplex structure.
Abstract
We derive a determinantal expression for the inhomogeneous asymmetric five vertex model in a rectangular geometry with arbitrary boundary conditions at the bottom and top. Standard non-intersecting lattice path, or free fermion, approaches are not applicable and the determinantal form thus is not immediate.
Abstract
We study a piecewise constant function with a finite number of discontinuities in any interval. We assume that the associated walk satisfying , pinned by , has finite length excursions over past minima. This is the continuous generalization of an initial ball configuration in the discrete Box Ball System introduced by Takahashi and Satsuma, where solitons of integer sizes are identified. We extend the slot decomposition developed by Ferrari, Nguyen, Rolla and Wang in the discrete setting to the continuous case. Each soliton of is represented by a point in two dimensional space, one coordinate for position and the other for the soliton height, mapping to a point configuration. We consider a distribution on walks given by a product measure on the decomposition of the path into excursions over past minima. Excursions are distributed as products of their solitons weights, which are determined by the soliton heights. We show that when the weight function is in the slot decomposition of is a Poisson process. This extends to the continuous case an approach of Ferrari and Gabrielli. As an example, we compute the intensity measure of the Poisson process associated to the asymmetric telegraph process introduced by Kac. In a forthcoming paper we discuss the dynamic properties.
Abstract
A fundamental identity for the number of vacillating tableaux was originally obtained through the representation theory of partition algebras. We extend this identity to arbitrary differential posets and show that it, together with several analogous identities, follows directly from the structural properties of differential posets. Specializing to Young's lattice and its Cartesian powers, we further obtain new bijective proofs via a simple deletion--insertion process.
Abstract
Motivated by the surge of interest in statistics of non-Hermitian random matrices as a framework for description of universal characteristics of dissipative chaotic quantum many-body systems, we address the problem of characterizing the parametric correlations of spectral densities. Considering parameter-dependent ensemble of complex Ginibre matrices we derive an explicit, closed-form expression for the parametric number covariance in the systems of symmetry class for eigenvalues in a circular domain containing on average a finite number of eigenvalues in the spectral bulk. This behavior is expected to be universal, as further supported by numerical evidence for the real Ginibre ensemble, non-Hermitian Bernoulli Wigner matrices and bi-unitarily invariant ensembles. We also discuss a relation between parametric correlations of spectral densities and the distribution of the so-called eigenvector non-orthogonality factor, which attracted considerable interest in recent years.
Abstract
In this work we consider an Aztec diamond model split into two unequal regions which are asymptotically fixed in size. Each region is weighted with a distinct two-periodic weighting. We refer to this model as the t-split two-periodic Aztec diamond, to signify its difference from the previous work title Split Two-Periodic Aztec Diamond, where the model was split into two equal regions. We derive an integral expression for the correlation kernel of the model and give a partial description of the scaling limit behavior, along with a conjecture for the remainder. We refer to the larger and smaller sides of the model as the dominant and non-dominant sides, and to the location of the weight change as the interface. The dominant side exhibits a limit shape that depends only on its own weighting and is identical to that of the two-periodic Aztec diamond, while the non-dominant side appears to have a novel limit shape that depends on both weightings and the location of the interface. Lastly, we consider the complete limit shape in the case where the dominant side two-periodic parameter goes to 0.
Abstract
We present a new classification of the class of all the hypergeometric orthogonal polynomial sequences. The classification uses properties of the coefficients and in the three-term recurrence relation satisfied by the orthogonal polynomial sequences. Such coefficients are rational functions of and are determined by a set of six parameters. The partial fraction decomposition of is a linear combination of five linearly independent functions of with coefficients , for , that are polynomials in our six parameters. For each value of a parameter , associated with the eigenvalues, we classify the according to the set of coefficients that are nonzero. There are certain particular values of that must be considered separately. Our classification yields a collection of 53 disjoint classes and it is different from the Askey scheme in several aspects. It is similar to the classifications proposed recently by Koornwinder.
Abstract
In an asymmetric shelf shuffle, a deck of cards is dealt sequentially from the bottom and assigned one of the shelves uniformly at random. The card is placed at the top of the assigned shelf with probability , and at the bottom of the assigned shelf with probability . Analysis of the shelf shuffle has gained much attention recently, and the case was first treated by Diaconis--Fulman--Holmes [Ann. Appl. Prob. 23 (2013), no. 4, 1692--1720]. In this paper, we extend the analysis of the shelf shuffle to general . In particular, we study the distribution of cycles, cycle lengths, number of descents, number of valleys, number of inversions, and the RSK shape of a permutation obtained from an asymmetric shelf shuffle. Our results extend the analysis of Diaconis--Fulman--Holmes to arbitrary . Furthermore, our analysis of the distribution of descents and inversions is new even for .
Abstract
A mechanical shuffler consists of shelves. A deck of cards, arranged in increasing order, is dealt from the bottom sequentially. Each card is assigned a shelf uniformly at random and placed on the top (bottom) of the existing pile with probability () independently. We refer to this as asymmetric shelf-shuffle. We find the law of the permutation induced by the asymmetric shelf-shuffle and show that the pair consisting of the number of descents and the number of valleys is a sufficient statistic. This generalizes a result of Diaconis, Fulman, and Holmes (Ann. Appl. Prob., 2013) corresponding to the case . For , Chen and Ottolini (ECP, 2025) established the cutoff in the total variation distance near . We establish the cutoff for the asymmetric shelf shuffle. Let be the uniform measure on the set of all permutations of . For a fixed and , we show that \[\TV\left(ν_{n, \lfloor cn^{3/2}\rfloor }^{(p)}, ν_n\right)=1-2Φ\left(-\frac{|2p-1|}{4\sqrt{3}c}\right)+O_{c, p}(n^{-1/2})\;.\] We also establish the cutoff in the separation distance near and in the relative entropy near . In both cases, we also obtain the cutoff profile explicitly.
Abstract
The Fomin local rules for Schensted column insertion can be seen as a two-lane box-ball system, in which a carrier moves particles forward or laterally. Running such two-lane dynamics in parallel on a periodic lattice gives rise to a two-dimensional generalization of the box-ball system, which we call the \emph{skew column RSK dynamics}. Equivalently, this is a deterministic dynamics on pairs of skew semistandard Young tableaux . We prove that this dynamics exhibits solitonic behavior and construct an explicit bijection that linearizes the time evolution. The resulting coordinates consist of two horizontally weak tableaux recording the asymptotic soliton data, integer riggings , and a weakly decreasing sequence of integers . A key feature of the construction is an explicit projection from the skew column RSK dynamics to the classical box-ball system; under this projection, the rigging is precisely the Kerov--Kirillov--Reshetikhin rigging of the associated box-ball configuration. Our proof uses two commuting affine crystal structures on pairs of skew tableaux and a novel connectivity theorem for distinguished subgraphs of tensor products of Kirillov--Reshetikhin crystals. We also derive Greene-type formulas for the soliton lengths in terms of last-passage percolation on the associated cylindrical environment. Finally, by taking generating functions in the linearizing coordinates, we obtain bijective proofs of Cauchy and Kawanaka--Littlewood-type identities for transformed Hall--Littlewood polynomials.
Abstract
In this paper, we establish a general comparison theorem for the order statistics of the edge eigenvectors for generalized Wigner matrices. Consequently, we derive the Gumbel law for the maximal edge eigenvector component and prove the universality of the Gaussian fluctuations of the order statistics in an intermediate regime close to the maximum. In addition, our comparison result also implies a quantitative first order estimate for moderately small order statistics.
Abstract
This work investigates the stability of outliers from the circular law, via the convergence of their associated diagonal overlaps between eigenvectors - also known as the squared eigenvalue condition numbers. We consider and compare two paradigmatic cases, namely: 1) the Complex Ginibre Ensemble conditioned on the existence of an outlier, and 2) the outlier induced by a rank-one Hermitian perturbation of a Complex Ginibre matrix. In both cases, we prove almost sure convergence towards a specific constant that only depends on the radius of the outlier and its status - either conditioned or induced. These results can be generalized to other complex integrable ensembles with the same techniques, and complement our understanding of eigenvalue stability in non-Hermitian ensembles.
Abstract
In this paper, we consider the Cauchy problem for the defocusing modified Korteweg-de Vries (mKdV) equation with a step-like initial data. Based on the Riemann-Hilbert problem associated with the mKdV equation, we derive the long-time asymptotic expansion of the solution to the defocusing mKdV equation in two transition regions using the nonlinear steepest descent method. It comes out that the leading term in the expansion is shown to match the corresponding background constants. The subleading term, however, decays at the order , and its coefficient is derived from the associated Painlevé \uppercase\expandafter{\romannumeral34\relax} model.
Abstract
We prove that the Opdam-Cherednik kernel, also known as the nonsymmetric Opdam hypergeometric function, can be written as the Laplace transform of a positive measure supported on the convex hull of the Weyl group orbit of its argument. As a consequence, the trigonometric Dunkl intertwining operator is positivity preserving. The main ingredient in the proof is a new formula for the Opdam-Cherednik kernel as a degeneration of nonsymmetric Macdonald polynomials. As a further application, we prove majorization inequalities for Macdonald polynomials and Heckman-Opdam hypergeometric functions associated with arbitrary root systems.
Abstract
Many-body integrable systems can be understood as a gas of quasiparticles. They propagate ballistically and drive large-scale transport. However, with the exception of the hard rods system, no tools have been available to numerically track such quasiparticles. Focusing on the Toda fluid, whose integrability relies on the availability of a Lax pair, we present a numerical scheme to track quasiparticle trajectories as determined by the time-dependent eigenvectors of the Lax matrix. Simulating the Toda fluid in thermal equilibrium, this tracking scheme is used to numerical confirm Brownian motion of a quasiparticle. Simulated is also the motion of a tagged particle. Our numerical results for the diffusion constant matches with a novel TBA prediction. We believe our numerical scheme can be extended to other classical many-particle models possessing a Lax matrix.
Abstract
The asymmetric simple exclusion process (ASEP) is an important particle model with deep connections to orthogonal polynomials. Motivated by this connection, Corteel and Williams introduced the Koornwinder moments at , which generalize the moments of Askey--Wilson polynomials. They showed that the partition function of the two-species ASEP is equal to for a one-row partition . In this paper, we investigate a conjecture of Rains on the positivity of the minimal numerator of the Koornwinder moment . We derive the first explicit formula for this moment, thereby obtaining a precise formulation of the conjecture by determining the minimal denominator of . We also propose a generalization of the conjecture for the more general Koornwinder moments indexed by two partitions at special parameter values. We prove the generalized Rains' conjecture in two special cases: and . For , we construct a lattice path model and obtain a combinatorial formula for in terms of non-intersecting lattice paths. For , we establish an explicit product formula for and give a combinatorial interpretation using lecture hall tableaux.
Abstract
We prove a local bijection between alternating sign matrices (ASMs) and descending plane partitions (DPPs) in the refined fibres of Behrend--Di Francesco--Zinn-Justin. The map is constructed by moving a boundary line through the ASM and DPP pictures. At each step only a bounded neighbourhood of that boundary is changed. In growth-diagram language, the boundary states are ordinary ASM row frontiers and shifted DPP path frontiers, and the elementary moves are reversible Fomin-type square rules. The proof checks that overlapping local moves are confluent, that every square is invertible, and that the completed boundary again satisfies the ordinary ASM and DPP coordinate inequalities. The only terminal choice is resolved by a radius-two fan rule on the visible height contour. The same local framework gives a statistic on refined ASMs whose generating function is the DPP degree enumerator. The statistic has a Striker--Fulmek quadratic core \[ Q_{\rm SF}(A)= \sum_{1\le r<s\le n}\ \sum_{1\le b<a\le n} (n-r+1)A_{r,a}A_{s,b}, \] plus a bounded fan-boundary correction and, outside the normal boundary gauge, an explicit compact-potential endpoint term. In the normal gauge this is \[ W_{\rm ASM}^{\rm SF}(A)=Q_{\rm SF}(A)+B_{\rm fan}(A). \] Thus the paper gives both a local ASM--DPP bijection and a statistic-preserving refined \(q\)-enumeration. The supplementary files reproduce selected local checks and the displayed arithmetic examples.
Abstract
A generalised inverse first passage problem asks whether, given a probability measure on , one can find a boundary such that the stopping time:\[τ:=\inf\left\{t:Λ\int_0^t ω(W_s-b(s))ds \geq U\right\}\] has distribution , where , and is a monotonic decreasing function. We construct a branching-selection particle system whose hydrodynamic limit is governed by a free boundary problem and connect this to the generalised inverse first passage problem. In the -particle system, particles move as independent Brownian motions, branch at a prescribed rate, and are removed at a rate proportional to their location relative to a position which is a function of the empirical distribution. We identify the limit of as the solution of the inverse first passage problem.
Abstract
A symmetric function is called Schur-positive if it admits an expansion in the Schur basis with nonnegative coefficients. In this paper, we study the Schur-positivity of symmetric functions naturally associated with set partitions, with respect to a descent set function that considers i as descent, if i and i+1 share a block in the partition. The Schur expansion involves hook-shaped Young diagrams, and the corresponding coefficients are given by Touchard-Riordan polynomials, which enumerate matchings by their number of crossings.
Abstract
Balanced words on a finite alphabet are those words in which every letter of the alphabet occurs the same number of times. The notion of descents and major index extends in a natural way to words. It is known that the bivariate generating polynomials for descents and major index over balanced words on the alphabet with occurrences each is palindromic, but a bijective proof has been missing even for balanced binary words. We give an explicit bijection proving this result. For permutations (which are also balanced), our bijection is different from the complementation map. We also show that for balanced binary words, this bijection simultaneously flips the ascent and comajor index as well.
Abstract
In this paper, we introduce and prove a set of rank recurrence relations for -Whittaker and Macdonald operators. We also derive an explicit expression for the -th power of the -Whittaker operators in terms of the -deformed binomial probability distribution, and we express the rank recurrence relations for the Macdonald operators in terms of the Cauchy determinant.
Abstract
We introduce a censored version of the stochastic six-vertex model. We show that for parameters , this model started from the initial condition is stochastically dominated at any time by the blocking measure. This is a partial analog of the censoring inequality for monotone spin systems. In particular, this result allows us to control the behavior of second-class particles. The proof uses parabolic Kazhdan--Lusztig -polynomials, whose appearance is explained using a connection between the stochastic six-vertex model and the Iwahori--Hecke algebras of symmetric groups. Furthermore, we find an intertwining relation for this process using normalized parabolic Kazhdan--Lusztig -polynomials as an intertwining kernel.
Abstract
We prove the Harish-Chandra isomorphism theorem for supersymmetric spaces, describing the polynomial algebra of eigenvalues of invariant differential operators. The polynomials obtained satisfy novel invariance conditions, which remain somewhat mysterious. We also prove the Harish-Chandra isomorphism for ghost distributions, which satisfy a `square root' of the invariance conditions coming from invariant differential operators. All proofs are algebraic, and rely on a rank-one reduction argument and the Chevalley restriction theorem.
Abstract
We study the large-time asymptotic behavior of a new KdV soliton gas. We first introduce a pure-soliton Riemann--Hilbert(RH) problem with \(2N\) poles and two different types of residue conditions. We show that, as \(N\to\infty\), this discrete problem converges to primitive-potential RH problem introduced by Dyachenko, Zakharov, and Zakharov, and the jump matrix of this soliton gas RH problem has two nonzero reflection coefficients. To analyze the large-time behavior, we apply the Deift--Zhou nonlinear steepest descent method together with an appropriate \(g\)-function mechanism. Through a sequence of transformations, the original RH problem is reduced to explicitly solvable model problems on an associated hyperelliptic Riemann surface. This allows us to derive an explicit leading-order asymptotic formula for the solution in terms of Jacobi elliptic function. The result provides a rigorous asymptotic description of a new KdV soliton gas and extends the available analysis beyond the previously studied case \(r_2\equiv 0\).
Abstract
In this paper, we introduce and study a family of diffusion processes associated with the Pitman-Yor process and the two-parameter Poisson-Dirichlet distribution. The diffusion coefficients indexed by a non-negative parameter are smaller than the corresponding one-parameter models in terms of quadratic forms or bilinear forms when is positive. The well known Petrov's diffusion corresponds to among the unlabelled diffusions. If is the same as the stable parameter in the Pitman-Yor process, we obtain both labelled and unlabelled reversible diffusion processes with the Pitman-Yor process and the two-parameter Poisson-Dirichlet distribution as the corresponding reversible measures. We construct these processes analytically through Dirichlet forms. In comparison with existing models in the literature, our models possess two fundamental new features. Firstly, our labelled model is characterized by an explicit generator, which is the first among all models studied so far. This makes it possible to establish the crucial integration-by-parts formula. Secondly, a novel foundational structure of the two-parameter distributions is the existence of a diversity index (a multiple of local time) for the positive stable parameter. By slowing down the diffusion in our model, the essential role of the diversity index is revealed in the evolution of the population. Additionally, we also obtain properties including ergodicity, path behaviour, and finite dimensional approximations.
Abstract
We prove that the Gromov--Witten invariants of the elliptic orbifold lines , , and satisfy a certain system of Hirota Quadratic (or Bilinear) Equations. Our result is the analogue of the so-called Toda conjecture in the Gromov-Witten theory of or more precisely its non-extended version. A new feature in our constructions is a certain bilinear operator whose principal symbol can be expressed in terms of elliptic theta functions.
Abstract
We study the KPZ equation on a circle with an additive spatial perturbation , where is a spacetime white noise and is a smooth spatial function. When , it is well-known that the unique invariant measure is the Brownian bridge. In the presence of the perturbation, we show that the equation admits a unique invariant measure that is absolutely continuous with respect to the Brownian bridge. We further prove the measure has a finite relative entropy with respect to the law of the bridge and that, for any , the corresponding Radon-Nikodym derivative belongs to , provided that is sufficiently small. The proof uses the discretization and mollification scheme of \cite{FQ}, together with an application of the log-Sobolev and spectral gap inequalities for the underlying Gaussian measure.
Abstract
We define and study the rectangular finite free heat flow, a dynamical system on polynomials that plays the role of the heat equation in the setting of rectangular finite free probability. We show several equivalent characterizations of the evolution (including PDE and gradient flow formulations), establish basic properties of the dynamics, and determine the asymptotic distributions of the polynomial roots in the long-time and high-degree limits. We also discuss connections with Calogero-Moser systems and Dunkl processes, and we show that the rectangular finite free heat flow describes the mean curvature expansion of a family of compact Lie group orbits.
Abstract
We generalize the Ahlswede--Daykin inequality (1978) to a Schur positive \emph{ADS inequality}, which also contains the Lam--Postnikov--Pylyavskyy inequality (2007) as a special case. We then present a number of further generalizations and applications. Notably, we resolve Mihalcea's conjecture on log-supermodularity of stable Grothendieck polynomials.
Abstract
In this manuscript, we explore the application of neural networks to predict the natural parameter of Schramm-Loewner Evolution (SLE) theory. SLE is a family of random fractal curves that has significant implications in Statistical Mechanics and Conformal Field Theory. This parameter plays an important role in the theory as there are models of Planar Statistical Physics that are proven to have SLE as scaling limits as well as others that are conjectured to have this limit for various choices of the parameter . In addition, there are three different statistical behaviors of the SLE curves as the parameter changes in Leveraging the powerful pattern recognition capabilities of neural networks, this study aims to develop a predictive model that can estimate the parameter with good accuracy.
Abstract
Permuted-basement Macdonald polynomials are nonsymmetric generalizations of symmetric Macdonald polynomials that form a basis for the polynomial ring for each fixed . There are combinatorial formulas for them as generating functions over composition-shaped non-attacking fillings. In this extended abstract, we bijectively prove identities for the relationship between , , , and . These identities correspond to two combinatorial operations on non-attacking fillings: (1) swapping adjacent entries in the basement, generalizing a result of Alexandersson (2019), and (2) swapping adjacent parts in the shape, which yields a straightening rule for expanding in the polynomials .
Abstract
In this paper, we prove universality of the distribution of the cokernels of a random -adic matrix with inhomogeneously balanced columns. More precisely, let be an integer and be a random matrix over whose -th column is -balanced. We prove that if as for every , then the cokernels of converge in distribution, as , to the same limiting law as the cokernels of Haar-random matrices over . This extends a universality theorem of Nguyen and Wood to random -adic matrices with inhomogeneously balanced columns.
Abstract
We introduce a novel exclusion process with two conservation laws, mass and energy, designed to mimic the essential features of continuous systems like interacting oscillators within the framework of interacting particle systems. This distinguishes our model from conventional multi-species processes where only particle numbers are conserved. As a basis for our fluctuation analysis, we first show that applying nonlinear fluctuating hydrodynamics (NFH) to this model reveals a wide variety of universality classes depending on the parameter choices. The main objective of this work is to study the stationary fluctuations of these conserved quantities. For a suitable choice of parameters, we rigorously show that the fluctuation fields converge to uncoupled stochastic Burgers equations (SBE) in the scaling limit. The proof relies on the second-order Boltzmann-Gibbs principle that we establish for this model, along with the spectral gap estimate and the equivalence of ensembles. Of independent interest is our general proof of the diagonalizability of the Jacobian matrix for the macroscopic current with distinct real eigenvalues. While this property is often taken as given in the physics literature, we establish it rigorously for multi-component systems even when the eigenvectors cannot be explicitly computed, offering a firm mathematical foundation for a broad class of models.
Abstract
Fix a multivariate polynomial in non-commuting variables of arbitrary degree, and consider independent complex Ginibre matrices . We prove that the empirical spectral distribution of converges as tends to infinity to the so-called Brown measure of evaluated at free circular variables. For polynomials of degree at most 2, the convergence was proven by Cook, Guionnet, and Husson \cite{cook2022spectrum}, and we prove that the convergence in fact holds for polynomials of any degree. The main step in the proof is a least singular value lower bound for for almost all complex shifts , and we prove this via a least singular value lower bound for a wide class of tensorized Ginibre matrices of finite type with a deterministic shift, which is of independent interest. We further show that the Brown measure convergence holds beyond Gaussians: the same convergence holds when the entry law has mean 0, variance 1, bounded density on and finite moments of all orders.
May 2026
Abstract
We introduce the \emph{tandem exclusion process} (TEP), a one-dimensional stochastic lattice model motivated by tandem running in ants. Particles evolve through two cooperative local transitions, at rate (leader advancement) and at rate (follower recovery). We prove that the stationary measure on the dynamically active sector is the Gibbs measure , where and counts neighboring occupied pairs, and derive exact closed-form expressions for the stationary current and spatial correlations using transfer-matrix methods. The current is asymmetric under particle--hole exchange , with its maximum occurring at densities strictly larger than , in contrast to the symmetric current of the totally asymmetric simple exclusion process (TASEP). For , cooperative dynamics enhances the current above the TASEP value and generates strong spatial clustering; in the limit , the current approaches , corresponding to nearly unconstrained collective transport. These results suggest that tandem coordination alone can substantially enhance collective transport efficiency at moderate and high densities, even without pheromone-mediated long-range communication.
Abstract
Gasper gave a two-step extension of the Askey-Roy formula for a beta-integral type defined on the unit circle. This involved increasing the number of parameters and adding a balancing condition. We derive a multidimensional generalization of the final formula in Gasper's extension. By considering a two-step degeneration involving parameters our generalization becomes Tarasov-Varchenko's formula -- itself a multidimensional generalization of the original Askey-Roy formula. The proof is done by Aomoto's method. This generalized integration by parts strategy makes use of a class of functions known as the fundamental invariants of type .
Abstract
We consider the directed polymer model on , in an i.i.d.\ random environment , focusing on the critical dimension . Our main contribution is to give a sharp lower bound on the free energy in the high-temperature regime. Our proof uses a percolation argument inspired by Lacoin (2010), for which we introduce a key property of bounded ``-energy'': this property quantifies the regularity of the polymer measures at diffusive scales and we show that it propagates along open paths. Writing for the quenched free energy, and setting and , our lower bound combined with Theorem 2.8 of Berger, Caravenna, and Turchi (2025) gives
Abstract
We consider matrices with independent, identically distributed entries, and prove that the sequence of measures converge to the Gaussian Multiplicative Chaos in the full subcritical regime as . Our result holds for both symmetry classes and in particular is new even for real Ginibre matrices, and is the first such convergence for any non-invariant ensemble of random matrices. We also establish the asymptotics for the -point function of at any collection of mesoscopically separated points . Our methods are analytic and probabilistic in nature, relying in part on the dynamical approach based on Dyson Brownian motion.
Abstract
We study combinatorial structures arising from finite-time transition probabilities of the Totally Asymmetric Simple Exclusion Process with open boundary conditions. While much of the existing combinatorial theory regarding the TASEP concerns the steady-state distribution, we focus instead on the transient dynamics. We first show that the enumeration of transition sequences between two configurations of the open TASEP is equivalent to the enumeration of standard Young tableaux of a family of non-classical shapes which have been of recent interest in the combinatorial literature. This extends to the open-boundary setting the correspondence between the TASEP with periodic boundaries and cylindric tableaux. We then introduce a family of tableau-like objects associated with Young diagrams in which repetitions of cells are allowed, subject to the partial order induced by the diagram. For each diagram, we collect the numbers of these objects into an exponential generating function. We prove that the entries of the homogeneous open TASEP transition matrix can be expressed as signed sums of such generating functions over suitable families of diagrams. This gives a combinatorial and order-theoretic interpretation of finite-time transition probabilities for the open TASEP, analogous to the combinatorial mappings known for steady-state probabilities.
Abstract
We prove asymptotic equivalents for finite-level representations of symmetric groups, that is, for Young diagrams having all but finitely many boxes on their first row. We deduce that random walks on symmetric groups generated by conjugacy classes with a macroscopic number of fixed points have a Poissonian cutoff profile. We also prove that the random involution walk exhibits cutoff and find its cutoff profile. Finally, we obtain numerics for the random transposition walk on a deck of 52 cards, giving concrete estimates on the question that originally motivated Diaconis and Shahshahani.
Abstract
While Macroscopic Fluctuation Theory (MFT) has been highly successful in analyzing non-equilibrium steady states, its application to non-steady-state processes remains limited. In this study, we apply MFT to the relaxation process of one-dimensional boundary-driven diffusive systems coupled to particle reservoirs at both ends. We exactly derive the current variance for systems with a constant diffusion coefficient and arbitrary mobility, as well as the cumulant generating function for the current in Reflective Brownian Motion (RBM). Our results demonstrate that non-steady current fluctuations during the approach to a steady state can be quantitatively described within the MFT framework.
Abstract
This paper establishes new upper bounds for the sum of the largest eigenvalues of symmetric matrices. When applied to the adjacency matrix of a graph, our results improve upon a related bound due to Mohar {\bf [On the sum of k largest eigenvalues of graphs and symmetric matrices, J. Combin. Theory Ser. B 99 (2009) 306--313]}. Furthermore, in the case of the Laplacian matrix, we prove that the well-known Brouwer's conjecture {\bf [Spectra of Graphs, Springer, New York, 2012]} holds for small values of for almost all graphs, thereby taking a significant step toward its complete resolution.
Abstract
We prove that the continuous-time, single-flip Glauber dynamics for lozenge tilings of the size- hexagon mix in time . This was predicted to hold on fairly general domains of diameter (on the basis of the ``Lifshitz law'' heuristic) but had previously only been established in domains such that the associated limit shape has no frozen facets. To access the hexagon, we introduce a multi-scale comparison argument between the height function of the random tiling evolving under the Glauber dynamics and the limit shape of a volume-tilted tiling (whose tilting parameter varies suitably in time).
Abstract
The symmetric simple exclusion process (SSEP) is a paradigmatic model of classical non-equilibrium dynamics. Exact results for large deviations of particle current in the SSEP have been obtained in various settings using integrability-based methods. In this Article, we discuss how these results are modified in the presence of localized slow bonds. We consider three conventional geometries: (a) a finite one-dimensional lattice weakly coupled to unequal reservoirs at its boundaries, (b) a semi-infinite one-dimensional lattice weakly coupled to a boundary reservoir, and (c) an infinite one-dimensional lattice with localized slow bonds near the origin. For each case, we present exact expressions for the large deviation function of current and validate them through rare-event simulations based on the cloning algorithm. In connection with our results, we present an elementary derivation of the exact large deviation function for the current in the semi-infinite SSEP, complementing recent results obtained through more elaborate techniques.
Abstract
We give a complete classification of the eternal solutions for the KPZ fixed point. Each of these is a (possibly infinite) patching together of the known eternal solutions, called Busemann functions. The resulting evolution of the KPZ fixed point exhibits a shock at each of the boundaries between the different Busemann functions. Moving forward in time, the shocks coalesce, while moving backwards in time, additional shocks can form. We describe several geometric properties of this tree of shocks.
Abstract
We prove that the maximal dimension of an irreducible representation of the symmetric group satisfies for some constant . This answers a question raised by Vershik--Kerov in 1985.
Abstract
We prove a Fredholm determinantal identity for the tilted Toeplitz minor generalizing the Borodin-Okounkov-Geronimo-Case (BOGC) identity to oblique splittings of the Hardy space. The tilts enter only through an oblique projection that multiplies the trace-class kernel inside the Fredholm determinant; the BOGC operator constructed from is unchanged. Baik-Liao-Liu (arXiv:2603.01964) and Liu-Tripathi (arXiv:2604.24747) have recently shown that the same tilted Toeplitz minor admits a contour Fredholm-determinantal representation, in connection with the periodic Totally Asymmetric Simple Exclusion Process (TASEP). In the periodic TASEP application of Baik-Liao-Liu, the formula plays an important role in identifying the periodic KPZ fixed point with general initial data. Our formula is a companion to their Fredholm determinant and readily reduces to the original BOGC identity. The one-sided tilted Toeplitz minor (that is, when all ) admits a bialternant form recovering Schur and Grothendieck polynomials as special cases. A Cauchy-Binet expansion realizes as a restricted sum over partitions of products of Jacobi-Trudi type determinants, generalizing Gessel's theorem. In the pure-shift setting this specializes to a skew Schur expansion. Finally, for finite Laurent exponential symbols, we record explicit resolvent-block flow identities and formulate the associated finite-dimensional closure problem. We also illustrate a possible asymptotic application leading to finite-rank perturbations of the Airy kernel.
Abstract
This manuscript studies a special case of the Hurwitz enumeration problem: for branched covers from genus g compact Riemann surface to the Riemann sphere, with three branch points, and require the branching data at one of the branch points to be of the partition (m 1^{n-m}), obtain a formula of Hurwitz number. The Hurwitz enumeration problem can be transformed into enumeration of a class of unrooted hypermaps. We first provide a enumeration formula for rooted hypermaps, thereby obtaining the weighted Hurwitz numbers. Next give the quantitative relationship between the enumeration of unrooted hypermaps and that of rooted hypermaps in the general setting. Finally, combining these two results, we obtain a formula for the unrooted hypermaps or the Hurwitz numbers in the special setting.
Abstract
Pipedreams and bumpless pipedreams are two combinatorial models that compute double Grothendieck polynomials. While studying matrix Schubert varieties, Pechenik, Speyer, and Weigandt defined a permutation statistic that captures the leading monomial of the top-degree components of a Grothendieck polynomial. Combinatorially, their result implies that there exists a unique pipedream (or bumpless pipedream) with row weight and column weight . A construction of such a pipedream was subsequently given by Chou and Yu. In this paper, we resolve the bumpless pipedream version of this problem by providing an explicit algorithm.
Abstract
The genus--g Fuss--Catalan (FC) number counts the number of ways to obtain a genus-g surface by identifying the edges of a pn--gon via p-valent hyperedges. For p=2 these are the genus--g Catalan numbers which are generated as the trace correlations in the Gaussian matrix model (GUE). Here we construct a simple two-matrix model which generates the higher-genus Fuss--Catalan numbers for any p as the coefficients of its 1/N-expansion. We obtain exact sum rules and an explicit formula for the higher-genus Fuss--Catalan numbers which generalises the Harer--Zagier formula to p>2. We discuss the relation of the higher-genus FC numbers to the intersection numbers and the Euler characteristic of the moduli space of spin-p curves.
Abstract
We study geodesics in the Brochette first-passage percolation model, where edges on the same axis-parallel line share a common random passage time, inducing long-range dependence. We focus on the maximal transversal deviation H n of geodesics from the origin to ne 1 . We prove existence of geodesics under mild assumptions and establish the order of magnitude of H n depending on the behavior of the passage-time distribution near its infimum. These results yield explicit wandering exponents in this dependent setting.
Abstract
We prove a complete characterization of the extremal invariant measures for half-space geometric last-passage percolation with an arbitrary boundary parameter. This is the first result of its kind for a model in the KPZ universality class that has boundary effects and an unbounded domain. A description of a class of invariant measures was previously given in a work of Barraquand and Corwin, where it was conjectured that these should comprise all extremal invariant measures. To complete the classification, we prove a one force--one solution principle: when started in the distant past from an arbitrary initial condition with a given asymptotic slope at , the recentered solution at time converges to a process which is distributed as the associated invariant measure with the specified slope. This limiting process is called the Busemann process, the first of its kind constructed for a half-space model. The Busemann process across all slopes is distributed as the joint invariant measure for geometric half-space LPP, recently constructed by Dauvergne and Zhang. There, it was conjectured that the constructed family of jointly invariant measures comprises all extremal jointly invariant measures; our analysis also confirms this conjecture. When the model has a strong (attractive) boundary, the collection of slopes for the invariant measures has a discontinuity, which does not arise in the full-space case. To handle this difficulty, we combine the control of the directions of semi-infinite geodesics with techniques from the theory of half-space Gibbsian line ensembles. Along the way, we classify the set of directions of semi-infinite geodesics for half-space geometric LPP, confirming a recent conjecture of Dauvergne and Zhang.