Integrable Probability arXiv Feed
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
March 2026
Abstract
We give an almost purely combinatorial expression for Wilson loop expectations of the Yang-Mills holonomy process with values in the unitary group on a compact oriented surface, possibly with boundary and arbitrary boundary conditions. Our main result computes the non-normalized expectation of products of traces of holonomies along an arbitrary family of immersed curves with transverse self-intersections and no triple points. It is expressed as a sum over assignments of highest weights of the unitary group to the connected components of the complement of the curves. Each term is a product of a Gaussian exponential factor, dimensions of unitary representations, and local contributions at the intersection points given by the sine or cosine of an angle determined by the surrounding highest weights. As an application, we obtain a new and short proof of the Makeenko-Migdal equations on arbitrary compact surfaces.
Abstract
The Central Limit Theorem (CLT) establishes that sufficiently large sequences of independent and identically distributed random variables converge in probability to a normal distribution. This makes the CLT a fundamental building block of statistical reasoning and, by extension, in reasoning about computing systems that are based on statistical inference such as probabilistic programing languages, programs with optimisation, and machine learning components. However, there is no general theory of CLT-like results currently, which forces practitioners to redo proofs without having a good handle on the essential ingredients of CLT-type results. In this paper, we introduce dilated seminorm-enriched category theory as a unifying framework for central limits, and we establish an abstract central limit theorem within that framework. We illustrate how a strengthened version of the classical CLT and the law of large numbers can be obtained as instances of our framework. Moreover, we derive from our framework a novel central limit theorem for symplectic manifolds, the CLT for observables, which finds applications in statistical mechanics.
Abstract
In 2000, Kadell gave an orthogonality conjecture for a symmetric function generalization of the Zeilberger--Bressoud -Dyson constant term identity. The non-zero part of Kadell's conjecture is a constant term identity indexed by a weak composition . This conjecture was first proved by Károlyi, Lascoux and Warnaar in 2015. They further formulated a closed-form expression for the above constant term when all parts of the composition are distinct. In 2021, Zhou obtained a recursion for this constant term for an arbitrary composition . In this paper, by categorizing the variables into two parts, we generalize Zhou's result.
Abstract
We introduce and analyze a broad class of continuous directed polymers in driven by Gaussian environments that are white in time and spatially correlated, under Dalang's condition. Using an Itô-renormalized stochastic-heat-equation representation, we establish structural properties of the partition function, including positivity, stationarity, scaling, homogeneity, and a Chapman--Kolmogorov relation. On finite time intervals, we prove Brownian-type pathwise behavior, namely Hölder continuity and identification of the quadratic variation. We then obtain a sharp measure-theoretic dichotomy: the quenched polymer measure is singular with respect to Wiener measure if and only if (equivalently, the noise is non-trace-class), and it is equivalent otherwise. Finally, in dimension , we prove diffusive behavior at large times in the high-temperature regime. This extends the Alberts--Khanin--Quastel framework from the white-noise setting to higher-dimensional Gaussian environments with general spatial covariance.
Abstract
In this paper, we consider the near-critical dimer model in the setup of isoradial superpositions with Temperleyan boundary conditions. We show that the centered height function converges as the mesh size tends to zero to a limiting field which agrees with the (electromagnetically tilted) sine-Gordon model, whose derivative correlations are described by Grassmann variables (or equivalently determinants involving a massive Dirac operator). This answers a longstanding question in the field. A crucial part of the work is to develop a notion of discrete massive holomorphic functions and the tools to study such functions, in particular finding an exact discrete form of the massive Cauchy--Riemann equations, which is satisfied by the inverse Kasteleyn matrix. In comparison with previous studies, a key novelty of this part of our work is that the mass is not only allowed to be non-constant but can be complex-valued.
Abstract
We study the isotropic six-vertex model on with spectral parameter , that is, with weights and . We show that the associated height function converges, in the scaling limit, to a properly scaled full-plane Gaussian free field. The result extends to anisotropic weights upon using a suitable embedding of the lattice.
Abstract
We consider the Sherrington--Kirkpatrick spin glass model with zero external field and at inverse temperature . Let be the corresponding log-partition function. Under the assumption that is bounded away from , we prove that Var As a consequence, we obtain Var for any fixed constant . We also prove a Gaussian central limit theorem for the centered and scaled .
Abstract
We prove the equivalent of the Baik, Ben Arous, Péché (2004) phenomenon for a novel, doubly sparse model where both the Wigner noise matrix and signal vector(s) are sparse. Specifically, we consider a deformed sub-Gaussian sparse Wigner ensemble with a fixed number of sub-Gaussian spike vectors of the same-order sparsity added. We show that spike vectors with signals greater than one are correlated with the top eigenvectors of the deformed ensemble and that each spike vector of signal greater than one induces an outlier eigenvalue. Notably, our results hold in the supercritical sparsity regime for the Wigner matrix () and for any sparse spike vector with an unbounded number of entries (). No further relationship between the sparsities of the noise matrix () and spike vectors () is necessary. This generalizes the work of Benaych-Georges and Nadakuditi (2010) and Péché (2005).
Abstract
This paper establishes a comparison theorem for the maximum eigenvalue of a sum of independent random symmetric matrices. The theorem states that the maximum eigenvalue of the matrix sum is dominated by the maximum eigenvalue of a Gaussian random matrix that inherits its statistics from the sum, and it strengthens previous results of this type. Corollaries address the minimum eigenvalue and the spectral norm. The comparison methodology is powerful because of the vast arsenal of tools for treating Gaussian random matrices. As applications, the paper improves on existing eigenvalue bounds for random matrices arising in spectral graph theory, quantum information theory, high-dimensional statistics, and numerical linear algebra. In particular, these techniques deliver the first complete proof that a sparse random dimension reduction map has the injectivity properties conjectured by Nelson & Nguyen in 2013.
Abstract
We study the stochastic heat flow with constant initial data and analyze its spatial average on the scale of . We prove that the logarithm of the averaged process satisfies a pointwise central limit theorem: After being centered by and scaled down by , it converges in distribution to a standard Gaussian.
Abstract
We consider the periodic totally asymmetric simple exclusion process with a general initial condition that properly approximates a periodic upper-semicontinuous function. We find the large time limit of the rescaled space-time multipoint distribution of the height function in the relaxation time scale. The limiting functions form a consistent family of finite-dimensional distributions; thus, they define the periodic KPZ fixed point with a general upper-semicontinuous initial condition. The main technical novelty of the paper is a hitting expectation representation of the energy function and the characteristic function in the finite-time multipoint distribution formula obtained in arXiv:1912.10143. The representation of the characteristic function is partly inspired by the work of arXiv:1701.00018, arXiv:2509.03246, while the representation of the energy function is based on an extensive guess-and-check exploration.
February 2026
Abstract
We consider stochastic lattice gases with stationary product weights and a polynomial perturbation vanishing with the system size that leads to condensation. If the density of particles exceeds a critical value the system phase separates into a bulk with homogeneous distribution of particles and a condensed phase. Depending on parameter values, the latter consists of a single macroscopic cluster or a diverging number of independent clusters on a smaller scale. We establish the condensation transition via the equivalence of ensembles and the main novelty is a derivation of the cluster size distribution using size-biased sampling, generalizing previous work on zero-range and inclusion processes. Simulations of zero-range processes illustrate our theoretical results on the condensate scale and size distribution.
Abstract
We study the rough-smooth boundary in the two-periodic Aztec diamond, a random domino tiling model exhibiting three types of macroscopic regions. We show that the height function at this boundary converges to an independent sum of an Airy surface and an i.i.d. noise field with fluctuations governed by the full-plane smooth phase. Going further, we prove convergence of a family of Temperleyan backbone paths to the Airy line ensemble. This gives the first convergence result for a family of undirected paths converging to the Airy line ensemble, as well as Airy convergence at a noisy boundary.
Abstract
Fluctuations of integrated currents have attracted considerable interest over the past decades in the context of statistical mechanics. Recently, anomalous current fluctuations, characterized by the M-Wright function, were obtained exactly in a classical automaton [. Krajnik et al., Phys. Rev. Lett. 128, 160601 (2022)], and previous studies have shown that the anomalous behavior can arise in a variety of classical systems. Despite the rapidly growing interest in such anomalous behaviors, which capture a universal aspect of one-dimensional many-body transport, the exact derivation of the M-Wright function in quantum many-body systems has remained elusive. In this Letter, we present the first exact microscopic derivation of the M-Wright function in quantum many-body dynamics by analyzing the integrated spin current in a one-dimensional Fermi-Hubbard model with infinitely strong repulsive interactions. Our results lay the groundwork for exploring anomalous integrated currents in a broad class of quantum many-body systems.
Abstract
We consider a probabilistic cellular automaton (PCA) of evaporation-deposition on the one-dimensional lattice having sites with periodic boundary conditions, in which each site, during each epoch, can be in one of two states: and . Fix a positive integer . There are two types of transitions at each discrete time, which are as follows: (i) the first site in every contiguous block of s becomes a with probability , and (ii) the first site in every contiguous block of s followed immediately by a also becomes a with probability . As in a PCA, all of these transitions occur simultaneously. We show that the resulting discrete-time Markov chain is ergodic, and we give an explicit formula for its limiting distribution, the partition function and the density. We also propose necessary and sufficient conditions for this Markov chain to be reversible. For , we provide a fully analytical expression for the free energy of this model.
Abstract
Coalescing particles on a line merge when they meet. When every site is initially occupied, only finitely many particles survive at any positive time, and their positions form a Pfaffian point process: all correlation functions are determined by pairwise quantities arranged in antisymmetric matrices. Previous proofs of this structure relied on analytic methods specific to time-homogeneous dynamics. We identify the checkerboard duality as the structural reason for the Pfaffian: on a discrete planar graph, binary random choices create two complementary non-crossing forests, one tracing ancestral lineages backward and the other carrying the coalescing particles forward as domain boundaries. This duality converts the absence of particles in an interval into coalescence of ancestral lineages at its endpoints. A cancellative labeling then converts coalescence into annihilation, for which a companion paper provides a Pfaffian formula. The resulting empty-interval formula holds for any discrete graph with the checkerboard structure and arbitrary inhomogeneous edge probabilities, requiring no symmetry or specific distributions. This covers settings beyond existing methods, including totally asymmetric dynamics and position-dependent transition rules, and yields an explicit Pfaffian point process in each setting. For Brownian motion, the formula recovers the known Pfaffian point process and the empty-interval probabilities previously derived by PDE methods.
Abstract
The numerical range of a non-normal matrix plays a central role as a descriptor of non-normal effects beyond spectral information. We study a class of fundamental non-Hermitian random matrix ensembles that interpolate between the Hermitian and non-Hermitian regimes. Our analysis focuses on the elliptic Ginibre ensemble and its chiral counterpart, as well as on non-Hermitian Wishart matrices. For each of these models, we explicitly characterise the geometry of the numerical range in the large-system limit. In particular, we show that for the elliptic Ginibre ensemble and its chiral version, the limiting numerical range is an ellipse, whereas for the non-Hermitian Wishart ensemble it is described by a non-elliptic envelope.
Abstract
We construct a new class of infinite-dimensional diffusions taking values in a generalized Kingman simplex. Our model describes the temporal evolution of the relative frequencies of infinitely-many types which are "labeled" by an arbitrary finite number of marks or colors, but "unlabeled" within each mark. We start with a finite-dimensional construction which extends to Wright-Fisher diffusions a self-similarity property known for Dirichlet distributions, and corresponds to a multiple skew-product representation of the Wright-Fisher diffusion relative to the marks in the population. After ranking decreasingly the frequencies within each mark, we identify the limit in distribution of the resulting diffusion when the number of types for each mark goes to infinity, and describe its infinitesimal operator. The limiting process reduces to a diffusion in the Thoma simplex in the special case of only two marks, whereas the infinitely-many-neutral-alleles model is recovered when all frequencies have the same mark. The stationary measure of the limit diffusion is shown to be the recently introduced multiple Poisson-Dirichlet distribution, which extends Kingman's Poisson-Dirichlet distribution and is the de Finetti representing measure for a family of random partitions whose elements are marked.
Abstract
In this note, we prove convergence of the half-space exponential last passage percolation (LPP) model, away from the boundary, to the directed landscape. Our approach couples the half-space and full-space LPP models and constructs two barrier events based on the monotonicity of last passage paths. Combining this coupling with moderate deviation estimates for both models and the known convergence of full-space LPP to the directed landscape, we establish the desired convergence.
Abstract
The quantum symmetric simple exclusion process (QSSEP) is a recent extension of the symmetric simple exclusion process, designed to model quantum coherent fluctuating effects in noisy diffusive systems. It models stochastic nearest-neighbor fermionic hopping on a lattice, possibly driven out-of-equilibrium by boundary processes. We present a direct formulation in the continuum, and establish how this formulation captures the scaling limit of the discrete version. In the continuum, QSSEP emerges as a non-commutative process, driven by free increments, conditioned on the algebra of functions on the ambiant space to encode spatial correlations. We actually develop a more general framework dealing with conditioned orbits with free increments which may find applications beyond the present context. We view this construction as a preliminary step toward formulating a quantum extension of the macroscopic fluctuation theory.
Abstract
In this paper, we consider five models of heavy-tailed queues involving Mittag-Leffler distributions that generalize the classical queues. These models are suitable modifications of previously defined models in such a way that the classical queue can be recovered by a suitable selection of parameters. We provide the distribution of inter-arrival and service times of both the original and modified queueing models. We then study the scaling limits of all the proposed models and we argue that the behaviour of the limiting processes can be used to characterise the traffic regime of the queues.
Abstract
We study the nonlinear multiplicative stochastic heat equation on Dirichlet spaces with white in time noise under weak disorder. We show that positive invariant random fields with uniformly bounded second moments are in one-to-one correspondence with bounded positive harmonic functions. The proof combines a remote past pullback construction with a uniqueness argument that applies a contraction inspired by chaos expansion. As a consequence, the space of invariant measures inherits geometric structure from the Martin Boundary. We further establish a small-noise Gaussian fluctuation result within each harmonic sector and show that the long-time behavior of solutions is completely determined by the Martin boundary data of the initial condition. These results reveal a direct connection between stochastic PDE dynamics and boundary theory in potential analysis.
Abstract
We obtain some results on the asymptotic behavior and zero distribution of the so-called geometric polynomials. The asymptotics is given both on compact subsets of and on compact subsets of the interval . The zeros of these polynomials are simple and lie in ; moreover, the zeros of consecutive polynomials interlace. Its zero distribution is a measure whose density is similar to Cauchy weight. Some orthogonality properties of these polynomials are also proved.
Abstract
Let be a fixed sequence of real numbers. At each stage, pick integers uniformly at random without replacement and then for each replace by . It is easy to observe that all the co-ordinates converge to . In this article, we extend the result of \cite{chatterjee2019note} by establishing order of decay of the expected distance. Furthermore, we establish the mixing time to be in between and .
Abstract
In two widely circulated manuscripts from the 1980s, I. G. Macdonald introduced certain multivariate hypergeometric series and and their -analogs and . These series are given by explicit expansions in Jack and Macdonald polynomials, and they generalize the hypergeometric functions of one and two matrix arguments from statistics. In a recent joint paper with Siddhartha Sahi, we constructed differential operators that characterize the Jack series thereby answering a question of Macdonald. In this paper we construct analogous -difference operators that characterize the Macdonald series . More precisely, we construct three -difference operators , , . The equation characterizes , while the equations and each characterize . These characterizations are subject to certain symmetry, boundary and stability condition. In the special case of , our operator was previously constructed by Kaneko in 1996.
Abstract
When particles on a line collide, they may annihilate-both are destroyed. Computing exact annihilation probabilities has been difficult because collisions reduce the particle count, while determinantal methods require a fixed count throughout. The ghost particle method, introduced in a companion paper for coalescence, keeps destroyed particles walking as invisible ghosts that restore the missing dimension. We apply this method to annihilation: when two particles annihilate, both trajectories continue as invisible walkers, yielding an exact determinantal formula that specifies the number of annihilations, where survivors end up, and where ghosts end up. For complete annihilation (no survivors), the formula simplifies to a Pfaffian-an algebraic relative of the determinant built from pairwise quantities-connecting to Pfaffian point process theory. The annihilation formula also yields results about coalescence: pairwise coalescence can be reinterpreted as complete annihilation, producing a Pfaffian coalescence formula. These formulas are exact for any finite initial configuration and apply to discrete lattice paths, birth-death chains, and continuous diffusions including Brownian motion.
Abstract
We consider a card guessing game with complete feedback. An ordered deck of cards labeled up to is shelf-shuffled exactly one time. One after the other a single card is drawn from the shuffled deck. The guesser makes has guess and the card is shown until no cards remain. We provide a distributional analysis of the number of correct guesses under the optimal strategy. We re-obtain the previously derived expectation and add a complete description of the distribution. We also obtain a central limit theorem for the number of cards tending to infinity. Furthermore, we discuss an unbalanced, biased shelf shuffle and show how to derive the extend our analysis, also adding the complete position matrix. Finally, a refined analysis of the number of correct guesses is carried out, distinguishing between pure luck guesses and certified correct guesses.
Abstract
In this paper, we derive new combinatorial formulas for symmetric Macdonald polynomials and integral Macdonald polynomials , in terms of several new statistics and the major index for a partition . Compared to previous formulas, these new ones contain the fewest terms. Moreover, three existing formulas for symmetric Macdonald polynomials established by Corteel--Mandelshtam--Williams (2022), Corteel--Haglund--Mandelshtam--Mason--Williams (2022) and Mandelshtam (2025) are recovered. Our proof relies on a new statistic on super fillings, employing the superization formulas of Haglund--Haimain--Loehr (2005) and Ayyer--Mandelshtam--Martin (2023), together with our recent approach to modified Macdonald polynomials.
Abstract
We present an exact integral formula for the multi-particle propagator of the one-dimensional Fermi--Hubbard model on an infinite lattice. The proof is based on the nested Bethe ansatz without relying on the string hypothesis. Our formula enables an explicit integral representation of the time evolution of arbitrary finite-particle wave functions and thereby provides a foundation for the exact analysis of nonequilibrium dynamics in the Hubbard model. It can further be applied to related open quantum models.
Abstract
The numerical range of a non-normal matrix plays a central role as a descriptor of non-normal effects beyond spectral information. We study a class of fundamental non-Hermitian random matrix ensembles that interpolate between the Hermitian and non-Hermitian regimes. Our analysis focuses on the elliptic Ginibre ensemble and its chiral counterpart, as well as on non-Hermitian Wishart matrices. For each of these models, we explicitly characterise the geometry of the numerical range in the large-system limit. In particular, we show that for the elliptic Ginibre ensemble and its chiral version, the limiting numerical range is an ellipse, whereas for the non-Hermitian Wishart ensemble it is described by a non-elliptic envelope.
Abstract
We construct a new class of infinite-dimensional diffusions taking values in a generalized Kingman simplex. Our model describes the temporal evolution of the relative frequencies of infinitely-many types which are "labeled" by an arbitrary finite number of marks or colors, but "unlabeled" within each mark. We start with a finite-dimensional construction which extends to Wright-Fisher diffusions a self-similarity property known for Dirichlet distributions, and corresponds to a multiple skew-product representation of the Wright-Fisher diffusion relative to the marks in the population. After ranking decreasingly the frequencies within each mark, we identify the limit in distribution of the resulting diffusion when the number of types for each mark goes to infinity, and describe its infinitesimal operator. The limiting process reduces to a diffusion in the Thoma simplex in the special case of only two marks, whereas the infinitely-many-neutral-alleles model is recovered when all frequencies have the same mark. The stationary measure of the limit diffusion is shown to be the recently introduced multiple Poisson-Dirichlet distribution, which extends Kingman's Poisson-Dirichlet distribution and is the de Finetti representing measure for a family of random partitions whose elements are marked.
Abstract
We study the local statistics of orthogonal polynomial ensembles near a hard edge, subject to a multiplicative deformation of the measure. Probabilistically, this deformation corresponds to a position-dependent conditional thinning of the particles. We prove that, under critical hard edge scaling and for a large class of potentials and deformation symbols, the correlation kernel of the conditional ensemble converges to a universal limit, which we identify as the conditional thinned Bessel point process. We derive an explicit expression for this limiting kernel in terms of the solution to a nonlocal integrable system depending on a parameter. For a special choice of the parameter, this system was recently identified in the study of multiplicative statistics of the Bessel point process. Our results establish that this system governs the full correlation structure of the conditional Bessel point process, extending the classical connection between the standard Bessel kernel and the Painlev\'e V equation.
Abstract
In this note, we prove convergence of the half-space exponential last passage percolation (LPP) model, away from the boundary, to the directed landscape. Our approach couples the half-space and full-space LPP models and constructs two barrier events based on the monotonicity of last passage paths. Combining this coupling with moderate deviation estimates for both models and the known convergence of full-space LPP to the directed landscape, we establish the desired convergence.
Abstract
We study an -species -PushTASEP, an integrable long-range stochastic process, on a one-dimensional periodic lattice with inhomogeneities and arbitrary capacity at each lattice site. The Markov matrix is identified with an alternating sum of commuting transfer matrices over all fundamental representations of . Stationary probabilities are expressed in a matrix product form involving a fusion of quantized corner transfer matrices for the strange five-vertex model introduced by Okado, Scrimshaw, and the second author. The resulting partition function, which serves as the normalization factor of the stationary probabilities, is obtained from the case by a finite plethystic substitution of length .
Abstract
Systems of instantaneously annihilating or coalescing Brownian motions on the line are considered. The extreme points of the set of entrance laws for this process are shown to be Pfaffian point processes at all times and their kernels are identified.
Abstract
The quantum symmetric simple exclusion process (QSSEP) is a recent extension of the symmetric simple exclusion process, designed to model quantum coherent fluctuating effects in noisy diffusive systems. It models stochastic nearest-neighbor fermionic hopping on a lattice, possibly driven out-of-equilibrium by boundary processes. We present a direct formulation in the continuum, and establish how this formulation captures the scaling limit of the discrete version. In the continuum, QSSEP emerges as a non-commutative process, driven by free increments, conditioned on the algebra of functions on the ambiant space to encode spatial correlations. We actually develop a more general framework dealing with conditioned orbits with free increments which may find applications beyond the present context. We view this construction as a preliminary step toward formulating a quantum extension of the macroscopic fluctuation theory.
Abstract
In this paper, we consider five models of heavy-tailed queues involving Mittag-Leffler distributions that generalize the classical queues. These models are suitable modifications of previously defined models in such a way that the classical queue can be recovered by a suitable selection of parameters. We provide the distribution of inter-arrival and service times of both the original and modified queueing models. We then study the scaling limits of all the proposed models and we argue that the behaviour of the limiting processes can be used to characterise the traffic regime of the queues.
Abstract
We study the nonlinear multiplicative stochastic heat equation on Dirichlet spaces with white in time noise under weak disorder. We show that positive invariant random fields with uniformly bounded second moments are in one-to-one correspondence with bounded positive harmonic functions. The proof combines a remote past pullback construction with a uniqueness argument that applies a contraction inspired by chaos expansion. As a consequence, the space of invariant measures inherits geometric structure from the Martin Boundary. We further establish a small-noise Gaussian fluctuation result within each harmonic sector and show that the long-time behavior of solutions is completely determined by the Martin boundary data of the initial condition. These results reveal a direct connection between stochastic PDE dynamics and boundary theory in potential analysis.
Abstract
We obtain some results on the asymptotic behavior and zero distribution of the so-called geometric polynomials. The asymptotics is given both on compact subsets of and on compact subsets of the interval . The zeros of these polynomials are simple and lie in ; moreover, the zeros of consecutive polynomials interlace. Its zero distribution is a measure whose density is similar to Cauchy weight. Some orthogonality properties of these polynomials are also proved.
Abstract
Let be a fixed sequence of real numbers. At each stage, pick integers uniformly at random without replacement and then for each replace by . It is easy to observe that all the co-ordinates converge to . In this article, we extend the result of \cite{chatterjee2019note} by establishing order of decay of the expected distance. Furthermore, we establish the mixing time to be in between and .
Abstract
We describe a simple Monte Carlo method for estimating by tossing a coin. Although the underlying Catalan-number series identities appear implicitly in the probability theory literature, the probabilistic interpretation of presented here seems to be new.
Abstract
Previous work of Ayyer, Martin, and Williams gave a probabilistic interpretation of the Macdonald polynomials at in terms of a Markov chain called the multispecies -Push TASEP, a Markov chain involving particles of types hopping around a ring. In particular, they showed that for each composition obtained by permuting the parts of , the stationary probability of being in state is proportional to the ASEP polynomial , and the normalizing constant (or partition function) is . There is an inhomogeneous generalization of Macdonald polynomials due to Knop and Sahi called interpolation Macdonald polynomials , as well as an inhomogeneous generalization of ASEP polynomials called interpolation ASEP polynomials that we introduced in previous work. In this article we introduce a new Markov chain called the interpolation -Push TASEP, and show that its steady state probabilities and partition function are given by the interpolation ASEP polynomials and the interpolation Macdonald polynomial, evaluated at . This generalizes the previous result of Ayyer, Martin, and Williams.
Abstract
In two widely circulated manuscripts from the 1980s, I. G. Macdonald introduced certain multivariate hypergeometric series and and their -analogs and . These series are given by explicit expansions in Jack and Macdonald polynomials, and they generalize the hypergeometric functions of one and two matrix arguments from statistics. In a recent joint paper with Siddhartha Sahi, we constructed differential operators that characterize the Jack series thereby answering a question of Macdonald. In this paper we construct analogous -difference operators that characterize the Macdonald series . More precisely, we construct three -difference operators , , . The equation characterizes , while the equations and each characterize . These characterizations are subject to certain symmetry, boundary and stability condition. In the special case of , our operator was previously constructed by Kaneko in 1996.
Abstract
When particles on a line collide, they may annihilate-both are destroyed. Computing exact annihilation probabilities has been difficult because collisions reduce the particle count, while determinantal methods require a fixed count throughout. The ghost particle method, introduced in a companion paper for coalescence, keeps destroyed particles walking as invisible ghosts that restore the missing dimension. We apply this method to annihilation: when two particles annihilate, both trajectories continue as invisible walkers, yielding an exact determinantal formula that specifies the number of annihilations, where survivors end up, and where ghosts end up. For complete annihilation (no survivors), the formula simplifies to a Pfaffian-an algebraic relative of the determinant built from pairwise quantities-connecting to Pfaffian point process theory. The annihilation formula also yields results about coalescence: pairwise coalescence can be reinterpreted as complete annihilation, producing a Pfaffian coalescence formula. These formulas are exact for any finite initial configuration and apply to discrete lattice paths, birth-death chains, and continuous diffusions including Brownian motion.
Abstract
We consider a card guessing game with complete feedback. An ordered deck of cards labeled up to is shelf-shuffled exactly one time. One after the other a single card is drawn from the shuffled deck. The guesser makes has guess and the card is shown until no cards remain. We provide a distributional analysis of the number of correct guesses under the optimal strategy. We re-obtain the previously derived expectation and add a complete description of the distribution. We also obtain a central limit theorem for the number of cards tending to infinity. Furthermore, we discuss an unbalanced, biased shelf shuffle and show how to derive the extend our analysis, also adding the complete position matrix. Finally, a refined analysis of the number of correct guesses is carried out, distinguishing between pure luck guesses and certified correct guesses.
Abstract
In this paper, we derive new combinatorial formulas for symmetric Macdonald polynomials and integral Macdonald polynomials , in terms of several new statistics and the major index for a partition . Compared to previous formulas, these new ones contain the fewest terms. Moreover, three existing formulas for symmetric Macdonald polynomials established by Corteel--Mandelshtam--Williams (2022), Corteel--Haglund--Mandelshtam--Mason--Williams (2022) and Mandelshtam (2025) are recovered. Our proof relies on a new statistic on super fillings, employing the superization formulas of Haglund--Haimain--Loehr (2005) and Ayyer--Mandelshtam--Martin (2023), together with our recent approach to modified Macdonald polynomials.
Abstract
The misanthrope process is an interacting particle system where particles move between neighbouring sites with hop rates depending only on the number of particles at the departure and arrival sites. Motivated by a discretised version of the Hammersley--Aldous--Diaconis process, we introduce a partially asymmetric long range misanthrope process (PALRMP) on a finite one-dimensional lattice with periodic boundary conditions where particles can move between sites that are not necessarily neighbours, as long as there are no particles in between the departure and arrival sites. In this model, each site has an inhomogeneous rate parameter associated to it, and the hop rate of a particle moving from site to site depends upon the parameter associated to the target site , the direction the particle moves, and the number of particles at sites and . We also consider the homogeneous PALRMP, where all the 's are 1. We find necessary and sufficient conditions on the hop rates under which the stationary distribution is of factorised form for both the PALRMP and the homogeneous PALRMP, as well as the extreme variants, namely the ones where the particle motion is totally asymmetric (TALRMP) and symmetric (SLRMP). As an illustrative example, we study in detail the discrete Hammersley--Aldous--Diaconis process.
Abstract
The directed landscape is a random directed metric on the plane that arises as the scaling limit of metric models in the KPZ universality class. For a pair of points p, q, the disjointness gap G(p; q) measures the shortfall when we optimize length over pairs of disjoint paths from p to q versus optimizing over all pairs of paths. Any spatial marginal of G is simply the gap between the top two lines in an Airy line ensemble. In this paper, we show that when the start and end time are fixed, the disjointness gap fully encodes the set of exceptional geodesic networks. The correspondence uses simple features of the disjointness gap, e.g. zeroes, local minima. We give a similar correspondence relating semi-infinite geodesic networks to a Busemann gap function. The proofs are deterministic given a list of soft properties related to the coalescent geometry of the directed landscape.
Abstract
We prove a conjecture of Fyodorov and Keating on the supercritical moments of the partition function of the CE field or equivalently the supercritical moments of moments of the characteristic polynomial of the CE ensemble for general and general real moment exponents. Moreover, we give the first expression for all correlation functions of the point process for all . The main object behind both results is the Hua-Pickrell stochastic zeta function introduced by Li and Valk\'{o}.
Abstract
For , let be the KPZ line ensemble with parameter , satisfying the homogeneous -Brownian Gibbs property with . We prove quantitative concentration estimates for the th line which yield the asymptotics as . A key step in the proof is a general integration by parts formula for -Brownian Gibbs line ensembles which yields the identity for any .
Abstract
We investigate the long-time behavior of the Airy wanderer line ensembles, an infinite-parameter family of Brownian Gibbsian line ensembles arising as edge-scaling limits of inhomogeneous models in the Kardar--Parisi--Zhang universality class. These ensembles are governed by sequences of nonnegative parameters that encode the asymptotic slopes of the curves at positive and negative infinity. Our main results characterize the fluctuations around this leading-order behavior and establish functional limit theorems for the ensembles near both ends of the spatial axis. We show that, at a macroscopic level, an Airy wanderer line ensemble organizes into groups of finitely many curves sharing a common asymptotic slope. After appropriate centering and scaling, each such group converges to a Dyson Brownian motion whose dimension equals the size of the group. In the case where only finitely many slope parameters are positive, we further prove a curve separation phenomenon: the upper curves follow deterministic parabolic trajectories, while the remaining lower curves remain globally flat and converge to the classical Airy line ensemble.
Abstract
We derive finite-particle rates for the regularized Stein variational gradient descent (R-SVGD) algorithm introduced by He et al. (2024) that corrects the constant-order bias of the SVGD by applying a resolvent-type preconditioner to the kernelized Wasserstein gradient. For the resulting interacting -particle system, we establish explicit non-asymptotic bounds for time-averaged (annealed) empirical measures, illustrating convergence in the \emph{true} (non-kernelized) Fisher information and, under a condition on the target, corresponding convergence for a large class of smooth kernels. Our analysis covers both continuous- and discrete-time dynamics and yields principled tuning rules for the regularization parameter, step size, and averaging horizon that quantify the trade-off between approximating the Wasserstein gradient flow and controlling finite-particle estimation error.
Abstract
Colored interlacing triangles, introduced by Aggarwal-Borodin-Wheeler (2024), provide the combinatorial framework for the Central Limit Theorem for probability measures arising from the Lascoux-Leclerc-Thibon (LLT) polynomials. Colored interlacing triangles depend on two key parameters: the number of colors and the depth of the triangle . Recent work of Gaetz-Gao (2025) connects these objects to Schubert calculus and resolves the enumeration for and arbitrary depth . However, the enumerative behavior for general has remained open. In this paper, we analyze the complementary regime: fixed depth and arbitrary number of colors . We prove that in this setting, colored interlacing triangles are in bijection with Dumont derangements, identifying their enumeration with the Genocchi medians. This connects the probabilistic model to a rich hierarchy of classical combinatorial objects. Furthermore, we introduce a -deformation of this enumeration arising naturally from the LLT transition energy. This yields new -analogs of the Genocchi medians. Finally, we present computational results and sampling algorithms for colored interlacing triangles with higher or , which suggests the limits of combinatorial tractability in the parameter space.
Abstract
We present an exact method for calculating the large deviation function describing rare fluctuations in the number of particles for product-kernel aggregation. Starting from the master equation, we derive an exact integral representation for the probability of observing particles at time starting from monomers for any finite . From this, we obtain an exact expression for the exponential moment for integer . Employing a replica conjecture -- numerically validated by finite- scaling -- we extend this result to real . The convex envelope of the large deviation function, obtained via a Legendre-Fenchel transform of the exponential moment, shows singular behavior. The singular structure allows us to construct the full phase diagram of product-kernel aggregation, which contains a tricritical point, separating continuous and discontinuous transitions. We also compute the asymptotic form of the LDF for small .
Abstract
We study equilibrium measures for Riesz gases in dimension with pairwise interaction kernel , subject to radially symmetric external fields. We characterise broad classes of confining potentials for which the equilibrium measure is supported on the unit ball and admits an explicit density. Our main contribution is a converse construction: starting from a prescribed radially symmetric equilibrium density given as a power series in the squared radius, we determine the associated external potential and establish the corresponding Euler-Lagrange variational conditions. A key ingredient in the proof is an identity between two hypergeometric functions evaluated at unit argument, which is of independent interest. As applications, we identify the external potentials corresponding to equilibrium densities proportional to , , and show that these potentials can be expressed in terms of Gauss hypergeometric functions , reducing to polynomials for special values of . We also determine the equilibrium measure associated with purely power-type external potentials, often referred to as Freud or Mittag--Leffler potentials in the context of log gases, for which the equilibrium density admits an explicit representation. Furthermore, we apply our framework to a Coulomb gas in dimension confined by a harmonic potential to the half-space. We derive a necessary condition under which the equilibrium measure is fully supported on the boundary hyperplane of dimension , with the induced density corresponding to that of a Riesz gas with exponent .
January 2026
Abstract
The purpose of this note is to explain the structure, general strategy, and main ideas of the proof in the work of Huang, McKenzie, and Yau (2024) on the Ramanujan property and edge universality of random regular graphs. The core of the argument is the derivation of self-consistent equations and a microscopic version of the loop equations for random -regular graphs. We first recall the local law for random -regular graphs, and then illustrate the main ideas behind the derivation of the self-consistent equations and the first loop equation.
Abstract
We study the totally asymmetric simple exclusion process (TASEP) on with step initial condition, in which all particles have distinct types. Our main object of interest is the type of the rightmost particle -- the leader -- at large time . We prove a central limit theorem for this random variable. Somewhat unexpectedly, the problem is closely connected to certain observables of voter and coalescing processes on ; we therefore derive their asymptotics as well. We also analyze the large-time behavior of a few other related observables, including certain multi-particle ones.
Abstract
The Tsetlin library is a random shuffling process on permutations of letters, where each letter can be interpreted as a book; book is brought to the front of the bookshelf with an assigned probability . We define a -deformation of the Tsetlin library by replacing the symmetric group action on permutations by the action of the type Iwahori-Hecke algebra. We compute the stationary distribution and spectrum of this Markov chain by relating it to a Markov chain on complete flags over the finite field vector space and applying techniques from semigroup theory. We also generalize the -Tsetlin library to words (with repeated letters), and compute its stationary distribution and spectrum.
Abstract
The fluctuations of lozenge tilings of hexagons with one or several cuts (nonconvexities) along opposite sides are governed by the (discrete-continuous) tacnode kernel , upon letting the hexagon become very large (or in other terms, keeping the hexagon fixed, with the tiles becoming very small). This is a point process with a finite number of (continuous) points along a discrete set of parallel lines within a specific region (see \cite{AJvM1,AJvM2}). Letting , one finds a liquid phase inscribed in the polygon, whose boundary (arctic curve) has a cusp near each cut, with two solid phases descending into the cusp (split-cusp). Duse-Johansson-Metcalfe \cite{DJM} show that in this situation the tile-fluctuations should obey the cusp-Airy statistics. It would have seem natural to expect to see the same cusp-Airy kernel in the neighborhood of the cut, for the limit () of the tacnode kernel . As it turns out, another statistics appears: the {\em transversal cusp-Airy} statistics, which was a puzzling fact to all of us. This statistics is derived and fully explained in this paper.
Abstract
In this paper, we prove that biorthogonal ensembles on the real line with a specific derivative structure admit an explicit correlation kernel of double contour integral form. We will demonstrate that this expression is a valuable starting point for asymptotic analysis and that our class of biorthogonal ensembles admits a large variety of limit kernels, by proving that two new classes of limit kernels can occur. The first type is a deformation of the hard edge Bessel kernel which arises in polynomial ensembles describing the eigenvalues of the sum of two random matrices, while the second type arises for Muttalib-Borodin type deformations of polynomial ensembles.
Abstract
We carry out the asymptotic analysis of repulsive ensembles of N particles which are discrete analogues of continuous 1d log-gases or beta-ensembles of random matrix theory. The ensembles that we study have several groups of particles which can have different intensities of repulsion. They appear naturally in models of random domino and lozenge tilings, random partitions, supersymmetric gauge theory, asymptotic representation theory, discrete orthogonal polynomial ensembles, etc. We allow filling fractions to be either fixed, or free, or to vary while respecting affine constraints. We are interested in the macroscopic behavior of the distribution of particles, captured by linear statistics, partition functions, and their finite-size corrections as N is large. We prove the law of large numbers and large deviations for the empirical measure around the equilibrium measure. To reach finite-size correction we assume off-criticality. For fixed filling fractions, we prove an asymptotic expansion for the partition function and for the cumulants of linear statistics, in particular establishing a central limit theorem. For varying filling fractions, we prove that the central limit theorem is perturbed by an additional discrete Gaussian component oscillating with N. We apply our general results to the study of uniformly random lozenge tilings on a large class of domains -- not necessarily planar, simply-connected, nor orientable. When the analogues of filling fractions are fixed and this domain is orientable, we show that the Gaussian fluctuations on the vertical extend to the whole liquid region and are governed there by the Gaussian free field, as predicted by the Kenyon-Okounkov conjecture. We also establish a modification of the Kenyon-Okounkov conjecture in the non-orientable case. Complementarily, we prove discrete Gaussian fluctuations for filling fractions, when they are not fixed.
Abstract
We introduce a family of generalized Rokhsar-Kivelson (RK) Hamiltonians, which are reverse-engineered to have an arbitrary edge-weighted superposition of dimer coverings as their exact ground state at the RK point. We then focus on a quantum dimer model on the triangular lattice, with doubly-periodic edge weights. For simplicity we consider a periodic model in which all weights are set to one except for a tunable horizontal edge weight labeled . We analytically show that the model exhibits a continuous quantum phase transition at , changing from a topological quantum spin liquid () to a columnar ordered state (). The dimer-dimer correlator decays exponentially on both sides of the transition with the correlation length and as a power-law at criticality. The vison correlator exhibits an exponential decay in the spin liquid phase, but becomes a constant in the ordered phase. We explain the constant vison correlator in terms of loops statistics of the double-dimer model. Using finite-size scaling of the vison correlator, we extract critical exponents consistent with the 2D Ising universality class.
Abstract
In this note, we briefly introduce the background and motivation of the collaborative work [arXiv:2508.20797], and provide an outline of the main results. The latter relates to matrix and higher order scalar differential equations satisfied by certain Hankel and Toeplitz determinants involving I-Bessel functions, or equivalently certain unitary matrix integrals, and moreover puts this property in a broader context. We also investigate large gaps between zeros of the derivatives of the Hardy -function, assuming the validity of a certain joint moments conjecture in random matrix theory.
Abstract
We present evidence for a conjectural relationship between absorption times for discrete Whittaker processes and maximal heights of non-intersecting Brownian bridges.
Abstract
We develop the basic theory of eigenvalues of -adic random matrices, analogous to the classical theory for random matrices over and . Such eigenvalue statistics were proposed as a model for the zeroes of -adic -functions by Ellenberg-Jain-Venkatesh, who computed the limiting distribution of the number of eigenvalues in a unit disc. We compute the full joint distribution of the eigenvalues of an matrix with Haar distribution, obtaining Coulomb gas type formulas as in the archimedean case, with Vandermonde terms leading to eigenvalue repulsion. From these Coulomb gas density functions we derive asymptotics of eigenvalue statistics as . These include exact computations, such as a closed form for the limiting pair correlation of eigenvalues in , and similar results in quadratic extensions. Such formulas yield concrete numerical predictions on zeroes of -adic -functions. For eigenvalues in arbitrary extensions of we also give precise estimates on their pair-repulsion and expected number of eigenvalues in each extension. Finally, we compute the asymptotic probability that all eigenvalues lie in . Our proofs combine results from several distinct areas: -adic orbital integrals, roots of random -adic polynomials, the Sawin-Wood moment method for random modules, and Markov chains associated with measures on integer partitions.
Abstract
We consider random integer partitions that follow the Poissonized Plancherel measure of parameter . Using RiemannHilbert techniques, we establish the asymptotics of the multiplicative averages for fixed in the regime and . We compute the large- expansion of expressing the rate function and the subsequent divergent and oscillatory contributions explicitly in terms of elliptic theta functions. The associated equilibrium measure presents, in general, nontrivial saturated regions and it undergoes two third-order phase transitions of different nature which we describe. Applications of our results include an explicit characterization of tail probabilities of the height function of the -deformed polynuclear growth model and of the edge of the positive-temperature discrete Bessel process and asymptotics of radially symmetric solutions to the 2D Toda equation with step-like initial data.
Abstract
Half-space models in the Kardar-Parisi-Zhang (KPZ) universality class exhibit rich boundary phenomena that alter the asymptotic behavior familiar from their full-space counterparts. A distinguishing feature of these systems is the presence of a boundary parameter that governs a transition between subcritical, critical, and supercritical regimes, characterized by different scaling exponents and fluctuation statistics. In this paper we construct the pinned half-space Airy line ensemble on -- a natural half-space analogue of the Airy line ensemble -- expected to arise as the universal scaling limit of supercritical half-space KPZ models. The ensemble is obtained as the weak limit of the critical half-space Airy line ensembles introduced in arXiv:2505.01798 as the boundary parameter tends to infinity. We show that has a Pfaffian point process structure with an explicit correlation kernel and that, after a parabolic shift, it satisfies a one-sided Brownian Gibbs property described by pairwise pinned Brownian motions. Far from the origin, converges to the standard Airy line ensemble, while at the origin its distribution coincides with that of the ordered eigenvalues (with doubled multiplicity) of the stochastic Airy operator with .
Abstract
It is well known that the -Whittaker polynomials, which are specializations of the Macdonald polynomials , expand positively as the sum of Schur polynomials. Macdonald polynomials have a quasisymmetric refinement: the quasisymmetric Macdonald polynomials , and a nonsymmetric refinement: the ASEP polynomials . We study the specializations of both these families of polynomials and show analogous properties: the quasisymmetric Macdonald polynomials expand positively as a sum of quasisymmetric Schur functions, , and the ASEP polynomials expand positively as a sum of Demazure atoms, . As a corollary of the latter, we prove more generally that any permuted basement Macdonald polynomial has a positive expansion in the Demazure atoms at . We give a description of the structure coefficients of and in both cases in terms of the charge statistic on a restricted set of semistandard tableaux.
Abstract
We study pattern densities in binary sequences, finding optimal limit sequences with fixed pattern densities.
Abstract
Studies of density matrices for random quantum states lead naturally to the fixed trace Laguerre ensemble in random matrix theory. Previous studies have uncovered explicit rational function formulas for moments of purity statistic (trace of the squared density matrix), and also a third order linear differential equation satisfied by the eigenvalue density. We further probe the origin of these results from the viewpoint of integrability, which is taken here to mean wider classes of recursions and differential equations, and give extensions. Prominent in our study are first order linear matrix differential equations. One application given is to the derivation of the third order scalar equation for the density. Another is to obtain the explicit rational function formula for the variance of the purity statistic in the generalised fixed trace Laguerre ensemble. In the original case (), the purity cumulants are expressed in terms of the large argument expansion of a particular -Painlev\'e IV transcendent. In a different but related direction, the exact computation of the two-point correlation for the fixed determinant circular unitary ensemble SU is given the Appendix.
December 2025
Abstract
We study fixed point biased involutions that avoid a pattern. For every pattern of length three we obtain limit theorems for the asymptotic distribution of the (appropriately centered and scaled) number of fixed points of a random fixed point biased involution avoiding that pattern. When the pattern being avoided is either , , or , we find a phase transition depending on the strength of the bias. We also obtain a limit theorem for distribution of fixed points when the pattern is for any and partial results when the pattern is .
Abstract
Let be a sequence of probabilities on the nonnegative integers, and be a sequence of independent random variables with law . For denote and , and assume . For , define the tilted probability , and let be a sequence of independent variables with law , and denote , with . Choose and denote . The Gibbs Conditioning Principle (GCP) holds if converges weakly to the law of , as . We prove the GCP for log-concave 's, meaning , subject to a technical condition that prevents condensation. The canonical measures are the distributions of the first variables, conditioned on their sum being . Efron's theorem states that for log-concave 's, the canonical measures are stochastically ordered with respect to . This, in turn, leads to the ordering of the conditioned tilted measures in terms of . This ordering is a fundamental component of our proof.
Abstract
We study the sine-Gordon measure defined on each homotopy class. The energy space decomposes into infinitely many such classes indexed by the topological degree . Even though the sine-Gordon action admits no minimizer in homotopy classes with , we prove that the Gibbs measure on each class nevertheless concentrates and exhibits Ornstein-Uhlenbeck fluctuations near the multi-soliton manifold in the joint low-temperature and infinite-volume limit. Moreover, we show that soliton collisions are unlikely events, so that typical states consist of solitons separated at an appropriate scale. Finally, we identify the joint distribution of the multi-soliton centers as the ordered statistics of independent uniform random variables, so that each soliton's location follows a Beta distribution.
Abstract
We consider the five-vertex model on a rectangular domain of the square lattice, with the so-called `scalar-product' boundary conditions. We address the evaluation of the free-energy density of the model in the scaling limit, that is when the number of sites is sent to infinity and the mesh of the lattice to zero, while keeping the size of the domain constant. To this aim, we reformulate the partition function of the model in terms of a discrete log-gas, and study its behaviour in the thermodynamic limit. We reproduce previous results, obtained by using a differential equation approach. Moreover, we provide the explicit form of the resolvent in all possible regimes. This work is preliminary to further studies of limit shape phenomena in the model.
Abstract
In this paper we introduce a new class of integrable 3D lattice models, possessing continuous families of commuting layer-to-layer transfer matrices. Algebraically, this commutativity is based on a very special construction of local Boltzmann weights in terms of quantum dilogarithms satisfying the inversion and pentagon identities. We give three examples of such quantum dilogarithms, leading to integrable 3D lattice models. The partition function per site in these models can be exactly calculated in the limit of an infinite lattice by using the functional relations, symmetry and factorization properties of the transfer matrix. The results of such calculations for 3D models associated with the Faddeev modular quantum dilogarithm are briefly presented.
Abstract
Stochastic growth models in the Kardar-Parisi-Zhang (KPZ) universality class exhibit remarkable fluctuation phenomena. While a variety of powerful methods have led to a detailed understanding of their typical fluctuations or large deviations, much less is known about behavior on intermediate, or moderate deviation, scales. Addressing this problem requires refined asymptotic control of the integrable structures underlying KPZ models. Motivated by this perspective, we study multiplicative statistics of discrete orthogonal polynomial ensembles (dOPEs) in different scaling regimes, with a particular focus on applications to tail probabilities of the height function in the stochastic six-vertex model. For a large class of dOPEs, we obtain robust singular asymptotic estimates for multiplicative statistics critically scaled near a saturated-to-band transition. These asymptotics exhibit universal crossover behavior, interpolating between Airy, Painlev\'e II, and Bessel-type regimes. Our proofs employ the Riemann-Hilbert Problem (RHP) approach to obtain asymptotics for the correlation kernel of a deformed version of the dOPE across the critical scaling windows. These asymptotics are then used on a double integral formula relating this kernel to partition function ratios, which may be of independent interest. At the technical level, the RHP analysis employs a novel parameter-dependent local parametrix, which requires a separate asymptotic analysis of its own. Using these results, together with a known identity relating a Laplace-type transform of the stochastic six-vertex model height function to a multiplicative statistic of the Meixner point process, we derive moderate deviation estimates for the height function in both the upper and lower tail regimes, with sharp exponents and constants.
Abstract
We construct a new solution to the three-dimensional reflection equation, a boundary analogue of the tetrahedron equation. The -operator is the one obtained by Sun, Terashima, Yagi, and the authors in 2024, involving four quantum dilogarithms with arguments in the -Weyl algebra. The new -operator similarly involves ten such quantum dilogarithms. Our approach is based on the quantum cluster algebra associated with the symmetric butterfly quiver on the wiring diagram of type C.
Abstract
The Airy wanderer line ensembles are infinite-parameter generalizations of the classical Airy line ensemble that arise naturally as scaling limits of inhomogeneous (spiked) models in the Kardar-Parisi-Zhang universality class. In this paper, we establish several structural properties of these ensembles. Our results show their laws depend continuously on the parameters, which encode the asymptotic slopes of the ensemble's curves near positive and negative infinity. We further prove that these ensembles admit multiple monotone couplings with respect to their parameters. Finally, we show that the Airy wanderer line ensembles are extreme points in the space of all Brownian Gibbsian line ensembles on the real line.
Abstract
We obtain precise asymptotics for the weighted number of domino tilings of an L-shaped subset of the Aztec diamond, obtained by removing an approximate rectangle in a corner of the Aztec diamond. By tuning the size of the removed corner, we observe different types of asymptotics. For a small removed corner, the number of tilings is close to that of the full Aztec diamond. Enlarging the removed corner to a critical size, a phase transition described in terms of the Tracy-Widom distribution occurs. Further increasing the size of the removed region, we observe a sharp decrease of the number of tilings, until it is finally approximated by the number of tilings of two smaller disjoint Aztec diamonds. We obtain uniform asymptotics for the number of domino tilings which fully describe these transitions.
Abstract
In classical percolation theory, the van den Berg-Kesten (BK) inequality is a fundamental tool that shows that disjoint events induce negative conditionings on each other. The inequality also holds in the context of last passage percolation (LPP), which is the zero temperature limit of polymer models and an important subclass in the Kardar-Parisi-Zhang (KPZ) universality class. Recently, an analog of the BK inequality was discovered in the context of zero temperature line ensembles and the scaling limit of LPP, where it was used to study upper tail probabilities of the weight and the scaling limit of geodesics under such upper tail conditionings. However, while it has become apparent that such an inequality in the positive temperature setting would have a number of applications, it seems likely that a direct generalization of the zero temperature inequality would not hold. In this work we prove a version of the BK inequality for the KPZ line ensemble and the continuum directed random polymer. We do so by working with the log gamma polymer, making use of its integrability and the geometric RSK correspondence. Our inequality serves as a key input in analyzing the KPZ line ensemble and proving sharp upper tail estimates of the KPZ equation in arXiv:2208.08922, and proving convergence of the continuum directed random polymer to Brownian bridge under the upper tail event in arXiv:2311.12009. The crucial role of integrability in the validity of such an inequality is highlighted via a counter-example for a non-integrable model.
Abstract
We prove a large deviations principle for the largest eigenvalue of Gaussian Kronecker matrices, namely matrices defined as the sum of tensors of independent Gaussian matrices in the regime where the dimension of the Gaussian matrices goes to infinity.
Abstract
We study the asymmetric simple exclusion process (ASEP) on a segment and are interested in its total variation distance to equilibrium when started from an initial configuration . We provide a general result which gives the cutoff window and profile whenever a KPZ-type limit theorem is available for an extension of to . We apply this result to obtain the cutoff window and profile of ASEP on the segment with flat, half-flat and step initial data. Our arguments are entirely probabilistic and make no use of Hecke algebras.
Abstract
We prove local convergence of the -PNG model with zero boundary to the stationary -PNG model, confirming a recent conjecture of Drillick and Lin (2024). The stationary -PNG model is the one with both left and bottom boundaries of Poisson nucleations with rate parameters and , respectively, for some . In the proof, we consider the trajectories of certain second class particles via a basic monotone coupling of three -PNG processes, and adapt microscopic concavity ideas used in particle models (e.g., Bal\'azs and Sepp\"al\"ainen (2009)), as well as blocking measure bounds like in Ferrari, Kipnis and Saada (1991).
Abstract
We consider a random interface model on the discrete torus with sites, obtained from the classical corner flip dynamics but with a weak global perturbation, namely an asymmetry of order of the direction of growth that switches direction based on the sign of the total area under the interface. The slopes of this model can be viewed as a non-simple exclusion process at half filling with globally dependent rates. We show that, for , the hydrodynamic equation of the empirical density is given by a time concatenation of the viscous Burgers equation and the heat equation. Moreover, for prime and , we establish convergence in law of the equilibrium fluctuations to an infinite-dimensional Ornstein-Uhlenbeck process.
Abstract
These are a set of lecture notes for a mini-course I gave at The University of Warwick from October 30th to November 1st, 2024. Recordings of the lectures are available on Oleg Zaboronski's webpage at https://warwick.ac.uk/fac/sci/maths/people/staff/oleg_zaboronski/jeffrey_kuan_visit/ . The main body of the notes covers the content of the lectures, and provides an introduction to Drinfel'd-Jimbo quantum groups and the Yang-Baxter equation, with a probabilist as the target audience. The appendix contains several topics, requested by colleagues during my visit to the United Kingdom, which all depend on the main set of notes. The notes begin by defining what it means for the asymmetric simple exclusion process (ASEP) to be integrable, in the sense of satisfying the Yang-Baxter equation. It then provides the algebraic background for the Yang-Baxter equation, by defining Drinfel'd-Jimbo groups as a quasi-triangular Hopf algebra. The algebraic background motivates generalizations of ASEP to stochastic vertex models and "fused" models. Each section corresponds to approximately an hour of lecture time. The appendix covers the F.R.T. construction, Hecke algebras, the matrix product ansatz, and orthogonal polynomial vertex weights. The topics in the appendix can be read independently of each other. Accessibility Statement: This PDF meets the technical standards of WCAG2.1AA, which complies with Ohio Administrative Policy IT-09 , Texas Administrative Code 206.70 and Title II of the Americans with Disabilities Act (effective April 24, 2026) . A webpage version of this PDF, typeset in MathML, is also available at https://go.osu.edu/QuantumKuan . To block web crawlers, the webpage is password protected. The password is TaySwift13. As an additional benefit, the webpage will have space for public comments and a list of updated errata, without the need to update the arXiv version.
Abstract
We introduce a new family of symmetric polynomials arising from exactly solvable lattice models associated with the quantised loop algebra . The polynomials unify -Whittaker polynomials, inhomogeneous -Whittaker polynomials, Grothendieck polynomials and their duals. Using Yang--Baxter equation, we derive Cauchy identities and combinatorial formulas for the transition coefficients.
Abstract
In this paper, we find necessary and sufficient conditions for the Law of Large Numbers of averaged empirical measures of -particle ensembles, in terms of the asymptotics of their Bessel generating functions, in the fixed temperature regime. This settles an open problem posed by Benaych-Georges, Cuenca and Gorin. For one direction, we use the moment method through Dunkl operators, and for the other we employ a special case of the formula of Chapuy--Dolega for the generating function of infinite constellations. As applications, we prove that the LLN for -sums and -corners of random matrices are given by the free convolution and free projection, respectively, regardless of the value of inverse temperature parameter . We also prove the LLN for a time-slice of the -Dyson Brownian motion.
Abstract
The ergodic theory of the open KPZ equation has seen significant progress in recent years, with explicit invariant measures described in a series of works by Corwin--Knizel, Barraquand--Le Doussal, and Bryc--Kuznetsov--Wang--Weso{\l}owski. In this paper, we provide a stochastic analytic proof of the formula for the invariant measures. Our approach starts from the Gaussian invariant measure for the case of homogeneous boundary conditions. We approximate the inhomogeneous problem by a homogeneous one with a singular boundary potential. Using tools including change of measure, time reversal for Markov processes, and It\^o's formula, we then reduce the problem to analyzing the KPZ nonlinearity in a thin boundary layer. Finally, using the theory of regularity structures, we establish a central limit theorem for the time-integrated nonlinearity near the boundary, which completes the proof of the invariance. Although it is known that different boundary parameters give rise to distinct physical regimes for the invariant measures, our method is robust and does not rely on any particular choice of boundary parameters.
Abstract
We introduce a multi-parameter family of random edge weights on the Aztec diamond graph, given by certain Gamma variables, and prove several results about the corresponding random dimer measures. Firstly, we show there is no phase transition at the level of the free energy. This provides rigorous backing for the physics predictions of Zeng-Leath-Hwa and later works that dimer models with random weights are in the glassy `super-rough' phase at all temperatures with no phase transition. Secondly, we show that the random dimer covers themselves enjoy exact distributional equalities of certain marginals with path locations in new `hybrid' integrable polymers. These reduce to the stationary log-Gamma, strict-weak, and Beta polymer in random environment in certain cases, allowing transfer of known results from integrable polymers to dimers with random weights. As an example application, we prove that the turning points at the boundaries of the Aztec diamond exhibit fluctuations of order , in contrast to the fluctuations for deterministic weights. Underlying all these is a key integrability property of the weights: they are the unique family for which independence is preserved under the shuffling algorithm.
Abstract
We present the transition probability for the asymmetric simple exclusion process on the half-space for general initial conditions and particle insertion at the boundary. In the limit of total asymmetry, where particles only jump to the right, we show that the transition probability reduces to a single Pfaffian for a large class of initial conditions. We further derive a Fredholm Pfaffian formula for the multipoint distribution of the totally asymmetric model conditioned on the number of particles in the system using a connection with a Pfaffian point process supported on Gelfand-Tsetlin patterns.
Abstract
We study the free boundary -Whittaker and Hall--Littlewood processes, two probability measures on sequences of partitions. We prove that a certain observable of the free boundary -Whittaker process exhibits a symmetry after a random shift, generalizing a previous result of Imamura, Mucciconi, and Sasamoto, and an extension of that result due to the first author. Our proof is completely different, and as part of our proof, we find contour integral formulas for the free boundary -Whittaker process. We also show a matching between certain observables in the free boundary Hall--Littlewood process and a quasi-open six vertex model, and explain how work of Finn and Vanicat gives an evaluation of a bounded sum over skew Hall--Littlewood functions as a rectangular Koornwinder polynomial.
November 2025
Abstract
We derive a form of master loop equations for the lattice Yang-Mills-Higgs theory with structure group , or . Compared to the pure Yang-Mills setting, several new operations arise. In fact, to obtain a closed recursion we must broaden the class of observables to include open Wilson lines. Our approach is based on the conditional Langevin dynamic and yields a concise proof via It\^o's formula.
Abstract
We conjecture affine or Hall-Littlewood analogues of the dual Jacobi-Trudi formulas for orthogonal and symplectic Schur functions indexed by rectangular partitions of maximal height. These conjectures are then used to derive -analogues of many known Rogers-Ramanujan identities for the characters of standard modules of affine Lie algebras. This includes -analogues of the classical Rogers-Ramanujan identities, (some of) the Andrews-Gordon identities and the , and GOW identities. We also prove an affine analogue of the dual Jacobi-Trudi formula for Schur functions indexed by rectangular partitions of arbitrary height.
Abstract
A solvable model of directed polymer with matrix-valued disorder is introduced in arXiv:2203.14868. The disorder is made of inverse-Wishart random matrices, so that the model nicely generalizes the well-studied log-gamma polymer, recovered when . Much of the features of the log-gamma polymer seem to have analogues for higher , although the integrability needs to be better understood. In this paper, we introduce stationary inverse-Wishart polymer models on a quadrant or a strip of . In each setting, we identify stationary measures, characterized explicitly in terms of random walks with inverse-Wishart increments in special cases, or more complicated two-layer Gibbs measures for generic choices of boundary parameters. We also make conjectures about asymptotics of the free energy, and explain important differences between matrix-valued polymer models and their scalar counterpart, due to non-commutativity.
Abstract
Actions in the Airy line ensemble represent distances from an infinitely far object. We characterize the Airy sheet by S(x,.)=T^x(.,1), where T^x is the unique action in the Airy line ensemble satisfying a growth condition depending on x. This provides a new simple framework for establishing convergence to the Airy sheet. We present simple conceptual proofs of such results in the case of Brownian last passage, the O'Connell-Yor semidiscrete polymer, the log-gamma polymer and the KPZ equation.
Abstract
We introduce a multi-colour multi-urn generalisation of the Bernoulli-Laplace urn model, consisting of urns, colours, and balls, with balls of each colour and balls in each urn. At each step, one ball is drawn uniformly at random from each urn, and the chosen balls are redistributed among the urns based on a permutation drawn from a distribution on the symmetric group . We study the mixing time of this Markov chain for fixed , , and , as . We show that there is cutoff whenever the chain on corresponding to the evolution of a single ball is irreducible, and that the same holds for a labeled version of the model. As an application, we also obtain partial results on cutoff for a card shuffling version of the model in which the cards are labeled and their ordering within each stack matters.
Abstract
We review our joint work on the scaling limits of disordered systems, linking the notion of disorder relevance/irrelevance to that of sub/super-criticality of singular SPDEs. This line of research culminated in the construction of the Critical 2D Stochastic Heat Flow (SHF), a universal process which provides a non-trivial solution to the Stochastic Heat Equation in dimension 2, a critical singular SPDE that lies beyond the reach of existing solution theories. The SHF also offers a rare example of a non-Gaussian scaling limit for a disordered system at its phase transition point in the critical dimension.
Abstract
We prove convergence of solutions of Dirichlet problems and Green's functions on Tutte's harmonic embeddings to those of the linearized Monge-Amp\`ere equation . The potential appears as the limit of piecewise linear potentials associated with the embeddings and the only assumption that we use is the uniform convexity of . Even if is quadratic, this setup significantly generalizes known results for discrete harmonic functions on orthodiagonal tilings. Motivated by potential applications to the analysis of 2d lattice models on irregular graphs, we also study the situation in which the limits are harmonic in a different complex structure.
Abstract
A Poisson line process is a random set of straight lines contained in the plane, as the image of the map , for each point of a Poisson process in the space-velocity plane. By associating a step with each line of the process, a random surface called multitime walk field is obtained. The diffusive rescaling of the surface converges to the multitime Brownian motion, a classical Gaussian field also called L\'evy-Chentsov field. A cut of the multitime fields with a perpendicular plane, reveals a one dimensional continuous time random walk and a Brownian motion, respectively. A hard rod is an interval contained in that travels ballistically until it collides with another hard rod, at which point they interchange positions. By associating each line with the ballistic displacement of a hard rod and associating surface steps with hard rod jumps, we obtain the hydrodynamic limits of the hard rods in the Euler and diffusive scalings. The main tools are law of large numbers and central limit theorems for Poisson processes. When rod sizes are zero we have an ideal gas dynamics. We describe the relation between ideal gas and hard-rod invariant measures.
October 2025
Abstract
We study longest and shortest directed paths in the following natural model of directed random planar maps: the uniform infinite bipolar-oriented triangulation (UIBOT), which is the local limit of uniform bipolar-oriented triangulations around a typical edge. We construct the Busemann function which measures directed distance to along a natural interface in the UIBOT. We show that in the case of longest (resp.\ shortest) directed paths, this Busemann function converges in the scaling limit to a -stable L\'evy process (resp.\ a -stable L\'evy process). We also prove up-to-constants bounds for directed distances in finite bipolar-oriented triangulations sampled from a Boltzmann distribution, and for size- cells in the UIBOT. These bounds imply that in a typical subset of the UIBOT with edges, longest directed path lengths are of order and shortest directed path lengths are of order . These results give the scaling dimensions for discretizations of the (hypothetical) -directed Liouville quantum gravity metrics. The main external input in our proof is the bijection of Kenyon-Miller-Sheffield-Wilson (2015). We do not use any continuum theory. We expect that our techniques can also be applied to prove similar results for directed distances in other random planar map models and for longest increasing subsequences in pattern-avoiding permutations.
Abstract
This is the second paper in a series studying the global asymptotics of discrete -particle systems with inverse temperature parameter in the high temperature regime. In the first paper, we established necessary and sufficient conditions for the Law of Large Numbers at high temperature in terms of Jack generating functions. In this paper, we derive a functional equation for the moment generating function of the limiting measure, which enables its analysis using analytic tools. We apply this functional equation to compute the densities of the high temperature limits of the pure Jack measures. As a special case, we obtain the high temperature limit of the large fixed-time distribution of the discrete-space -Dyson Brownian motion of Gorin-Shkolnikov. Two special cases of our densities are the high temperature limits of discrete versions of the GE, computed by Allez-Bouchaud-Guionnet in [Phys. Rev. Lett. 109 (2012), 094102; arXiv:1205.3598], and LE, computed by Allez-Bouchaud-Majumdar-Vivo in [J. Phys. A, vol. 46, no. 1 (2013), 015001; arXiv:1209.6171]. Moreover, we prove the following crystallization phenomenon of the particles in the high temperature limit: the limiting measures are uniformly supported on disjoint intervals with unit gaps and their locations correspond to the zeros of explicit special functions with all roots located in the real line. We also show that these zeros correspond to the spectra of certain unbounded Jacobi operators.
Abstract
We study a problem with three equivalent formulations: describing Gibbs measures for five-vertex model in quadrant; classifying coherent systems on a p-deformation of the Gelfand-Tsetlin graph related to Grothendieck polynomials; finding the Martin boundary for discrete time TASEP with p-geometric jumps. We find a wide family of the Gibbs measures, parameterized by certain analytic functions. A subset of our measures have probabilistic interpretation as interacting particle systems with fixed particles speeds. In contrast to previous related boundary problems, we find that admissible speeds are not arbitrary, but must be larger than . For this subset we further establish Law of Large Numbers and Central Limit Theorem, connecting the fluctuations to families of independent GUE eigenvalues. As a consequence, the measures from the subset are extreme points of the Martin boundary. It remains open whether our list of measures is exhaustive.
Abstract
In arbitrary spatial dimension , we study a generalized model of random walks in a time-varying random environment (RWRE) defined by a stochastic flow of kernels. We consider the quenched probability distribution of the random walker under a scaling where the time is of order and the spatial window is of size . This spatial window may not necessarily be centered close to the origin. We show that as there are Gaussian fluctuations up to a certain specific spatial centering radius in the tail of the quenched probability distribution, which we call the critical scale. This critical scale depends on the spatial dimension of the underlying random walk, specifically when , when , and when . In the particular case of centering the fluctuation window at the origin, our results recover and generalize some known fluctuation results for related models. However, farther from the origin, the previous literature is more sparse. The noise coefficient in the limiting Gaussian field is nontrivial and depends on the invariant measure of the two-point motion of the underlying RWRE model. We furthermore reconcile some of these coefficient formulas with previous works. As part of the proof, we introduce a general class of Markov chains with short-range interactions that admit nice estimates and limit formulas. One of the key technical results for such Markov chains is that in , one can propagate test functions backwards in time to obtain precise limiting moment formulas.
Abstract
We consider a one-dimensional gas of hard rods, one of the simplest examples of an interacting integrable model. It is well known that the hydrodynamics of such integrable models can be understood by viewing the system as a gas of quasiparticles. Here, we explore the dynamics of individual quasiparticles for a variety of initial conditions of the background gas. The mean, variance, and two-time correlations are computed exactly and lead to a picture of quasiparticles as drifting Brownian particles. For the case of a homogeneous background, we show that the motion of two tagged quasiparticles is strongly correlated, and they move like a rigid rod at late times. Apart from a microscopic derivation based on the mapping to point particles, we provide an alternate derivation which emphasizes that quasiparticle fluctuations are related to initial phase-space fluctuations, which are carried over in time by Euler scale dynamics. For the homogeneous state, we use the Brownian motion picture to develop a Dean-Kawasaki-type fluctuating hydrodynamic theory, formally having the same structure as that derived recently by Ferrari and Olla. We discuss differences with existing proposals on the hydrodynamics of hard rods and some puzzles.
Abstract
Although for the most part classical, the topic of electrostatics finds to this day new applications. In this review we highlight several theoretical results on electrostatics, chosen to both illustrate general principles, and for their application in statistical mechanics and random matrix settings. The theoretical results include electrostatic potentials and energies associated with balls and hyperellipsoids in general dimension, the use of conformal mappings in two-dimensions, and the balayage measure. A number of explicit examples of their use in predicting the leading asymptotic form of certain configuration integrals and particle density in particular statistical mechanical systems are given, as well as with regards to questions relating to fluctuation formulas and (conditioned) gap probabilities.
Abstract
We establish an alternative, ``perpendicular" collection of generating functions for the coefficients of Gaussian polynomials, . We provide a general characterization of these perpendicular generating functions. For small values of , unimodality of the coefficients of Gaussian polynomials is easily proved from these generating functions. Additionally, we uncover new and surprising identities for the differences of Gaussian polynomial coefficients, including a very unexpected infinite family of congruences for coefficients of .
Abstract
We present some review material relating to the topic of optimal asymptotic expansions of correlation functions and associated observables for ensembles in random matrix theory. We also give an introduction to a related line of study that we are presently undertaking.
Abstract
We consider the asymmetric version of the Popkov-Sch\"{u}tz two-lane lattice gas with general jump rates, subject to the stationary measure being of product form. This still leaves five free parameters. At density 1/2 the eigenvalues of the flux Jacobian are degenerate. We compute the second order expansion of the average fluxes at density 1/2 and thereby identify the universality classes.
Abstract
A Littlewood identity is an identity equating a sum of Schur functions with an infinite product. A bounded Littlewood identity is one where the sum is taken over the partitions with a bounded number of rows or columns. The price to pay is that the infinite product has to be replaced by a determinant. The focus of this article is on refinements of such bounded Littlewood identities where one also prescribes the number of odd-length rows or columns of the partitions. Goulden [{\it Discrete Math.} {\bf99} (1992), 69--77] had given such a refinement in which the number of columns is bounded and the number of odd-length rows is prescribed. We provide refinements where the number of columns is bounded and the number of odd-length columns is prescribed. Furthermore, we present new formulations of such bounded Littlewood identities involving skewing operators. As corollaries we obtain non-standard formulas for numbers of standard Young tableaux with restricted shapes as above. In the last part of the article we discuss combinatorial interpretations of such identities in terms of up-down tableaux. As corollaries, we obtain identities between numbers of standard Young tableaux and numbers of (marked) vacillating tableaux.
Abstract
In a widely circulated manuscript from the 1980s, now available on the arXiv, I.~G.~Macdonald introduced certain multivariable hypergeometric series and in one and two sets of variables and . These two series are defined by explicit expansions in terms of Jack polynomials , and for they specialize to the hypergeometric series of matrix arguments studied by Herz (1955) and Constantine (1963) that admit analogous expansions in terms of zonal polynomials. In this paper we determine explicit partial differential equations that characterize , thereby answering a question posed by Macdonald. More precisely, for each we construct three differential operators , , , and we show that and are the unique series solutions of the equations and , respectively, subject to certain symmetry and boundary conditions. We also prove that the equation characterizes , but only after one restricts the domain of to the set of series satisfying an additional stability condition with respect to . Special cases of the operators and have been constructed previously in the literature, but only for a small number of pairs , namely for and in the zonal case by Muirhead (1970), Constantine--Muirhead (1972), and Fujikoshi (1975); and for and in the general Jack case by Macdonald (1980s), Yan (1992), Kaneko (1993), and Baker--Forrester (1997). However the operator seems to be new even for these special cases.
Abstract
Brownian Boost is a one-parameter family of stochastic differential games played on the real line in which players spend at rates of their choosing in an ongoing effort to influence the drift of a randomly diffusing point particle~. One or other player is rewarded, at time infinity, according to whether~ tends to plus or minus infinity. Each player's net receipt is the final reward (only for the victor) minus the player's total spend. We characterise and explicitly compute the time-homogeneous Markov-perfect Nash equilibria of Brownian Boost, finding the derivatives of the players' expected payoffs to solve a pair of coupled first-order non-linear ODE. Brownian Boost is a high-noise limit of a two-dimensional family of player-funded tug-of-war games, one of which was studied in~\cite{LostPennies}. We analyse the discrete games, finding them, and Brownian Boost, to exemplify key features studied in the economics literature of tug-of-war initiated by~\cite{HarrisVickers87}: a battlefield region where players spend heavily; stakes that decay rapidly but asymmetrically in distance to the battlefield; and an effect of discouragement that makes equilibria fragile under asymmetric perturbation of incentive. Tug-of-war has a parallel mathematical literature derived from~\cite{PSSW09}, which solved the scaled fair-coin game in a Euclidean domain via the infinity Laplacian PDE. By offering an analytic solution to Brownian Boost, a game that models strategic interaction and resource allocation, we seek to build a bridge between the two tug-of-war literatures.
Abstract
We study a quantum version of the -dimer model from statistical mechanics, based on the formalism from quantum topology developed by Reshetikhin and Turaev (the latter which, in particular, can be used to construct the Jones polynomial of a knot in ). We apply this machinery to construct an isotopy invariant polynomial for knotted bipartite ribbon graphs in , giving, in the planar setting, a quantum -dimer partition function. As one application, we compute the expected number of loops in the (classical) double dimer model for planar bipartite graphs.
Abstract
We study line ensembles arising naturally in symmetrized/half-space geometric last passage percolation (LPP) on the square. The weights of the model are geometrically distributed with parameter off the diagonal and on the diagonal, where and . In the supercritical regime , we show that the ensembles undergo a phase transition: the top curve separates from the rest and converges to a Brownian motion under fluctuations and spatial scaling, while the remaining curves converge to the Airy line ensemble under fluctuations and spatial scaling. Our analysis relies on a distributional identity between half-space LPP and the Pfaffian Schur process. The latter exhibits two key structures: (1) a Pfaffian point process, which we use to establish finite-dimensional convergence of the ensembles, and (2) a Gibbsian line ensemble, which we use to extend convergence uniformly over compact sets.
Abstract
In this note we provide a short proof of the distributional equality between last passage percolation with geometric weights along a general down-right path and Schur processes. We do this in both the full-space and half-space settings, and for general parameters. The main inputs for our arguments are generalizations of the Robinson-Schensted-Knuth correspondence and Greene's theorem due to Krattenthaler, which are based on Fomin's growth diagrams.
Abstract
In 1996, Knop and Sahi introduced a remarkable family of inhomogeneous symmetric polynomials, defined via vanishing conditions, whose top homogeneous parts are exactly the Macdonald polynomials. Like the Macdonald polynomials, these interpolation Macdonald polynomials are closely connected to the Hecke algebra, and admit nonsymmetric versions, which generalize the nonsymmetric Macdonald polynomials. In this paper we give a combinatorial formula for interpolation Macdonald polynomials in terms of signed multiline queues; this formula generalizes the combinatorial formula for Macdonald polynomials in terms of multiline queues given by Corteel-Mandelshtam-Williams.
September 2025
Abstract
We study partition functions with domain-wall like boundary conditions for path models issued from colored vertex models. These models display an arctic phenomenon, as attested by numerical simulations. We show that the colored vertex model is equivalent to a certain single-color ``colorblind" vertex model. In a special case of the weights for the colorblind touching paths, we derive the arctic curve using a bijective sliding map to non-intersecting paths, for which arctic curves were previously derived using the tangent method. The resulting arctic curves are only piecewise analytic, as in the known non-free fermion cases of Six vertex model with domain-wall boundaries and its relatives. We also prove a shear phenomenon, that some portions of the arctic curve are sheared versions of the analytic continuation of other portions, as already observed in the uniformly weighted Six and Twenty vertex models.
Abstract
We consider a semilinear stochastic heat equation in spatial dimension at least , forced by a noise that is white in time with a covariance kernel that decays like as . We show that in an appropriate diffusive scaling limit with a logarithmic attenuation of the noise, the pointwise statistics of the solution can be approximated by the solution to a forward-backward stochastic differential equation (FBSDE). The scaling and structure of the problem is similar to that of the two-dimensional stochastic heat equation forced by an approximation of space-time white noise considered by the first author and Gu (Ann. Probab., 2022). However the resulting FBSDE is different due to the long-range correlations of the noise.
Abstract
We say a sequence of symmetric functions is Schur log-concave if is Schur positive for all . We conjecture that a very general class of sequences of Schur functions satisfies this property, and show it for sequences of Schur functions indexed by partitions with growing first part and column. Our findings are related to work of Lam, Postnikov and Pylyavskyy on Schur positivity, and of Butler, Sagan, and the second author on -log-concavity.
Abstract
In this paper, we study run-and-tumble particles moving on two copies of the discrete torus (referred to as layers), where the switching rate between layers depends on a mean-field interaction among the particles. We derive the hydrodynamic limit of this model, as well as the large deviations from the hydrodynamic limit. Our main tool is the introduction of a weakly perturbed version of the system, whose hydrodynamic equations precisely characterize the trajectories associated with large deviations.
Abstract
We introduce the Stochastic Box-Ball System (SBBS), a probabilistic cellular automaton that generalizes the classic Takahashi-Satsuma Box-Ball System. In SBBS, particles are transported by a carrier with a fixed capacity that may fail to pick up any given particle with a fixed probability . This model interpolates between two known integrable systems: the Box-Ball System (as ) and the PushTASEP (as ). We show that the long-term behavior of SBBS is governed by isolated particles and the occasional emergence of short solitons, which can form longer solitons but are more likely to fall apart. More precisely, we first show that all particles are isolated except for a -fraction of times in any given steps, and solitons keep forming for this fraction of times. We then show that under diffusive scaling, both SBBS (for any carrier capacity) and PushTASEP converge weakly to semimartingale reflecting Brownian Motions (SRBMs) on the Weyl chamber with explicit covariance and reflection matrices, which are consistent with the microscale relations between these systems. The reflection matrix for SBBS is determined by how 2-solitons behave and exhibit ``solitonic bias'' visible in the diffusive scale. Our proof relies on a new, extended SRBM invariance principle that we develop in this work. This principle can handle processes with complex boundary behavior that can be written as "overdetermined" Skorokhod decompositions, which is crucial for analyzing the complex solitonic interaction in SBBS. We believe this tool may be of independent interest.
Abstract
Majorization inequalities have a long history, going back to Maclaurin and Newton. They were recently studied for several families of symmetric functions, including by Cuttler--Greene--Skandera (2011), Sra (2016), Khare--Tao (2021), McSwiggen--Novak (2022), and Chen--Sahi (2024+) among others. Here we extend the inequalities by these authors to Jack and Macdonald polynomials, and obtain conjectural characterizations of majorization and of weak majorization of the underlying partitions. We prove these characterizations for two variables. In fact, we upgrade -- and prove in the above cases -- the characterization of majorization, to containment of Jack and Macdonald differences lying in the Muirhead semiring.
Abstract
We show that the spatial increments of the KPZ fixed point starting from arbitrary initial data, exhibit strong quantitative comparison against rate two Brownian motion on compacts. The above estimates are uniform for uniformly bounded continuous, compactly supported initial data and countably many narrow wedges with supports contained in a fixed compact set.
Abstract
We show that the Radon-Nikodym derivative of the law of the spatial increments (with endpoints away from the origin) of inhomogeneous Brownian last passage percolation (LPP) with non-decreasing initial data against the Wiener measure on compacts is in ; and for any fixed , the norm is at most of the order for some -dependent constant . Furthermore, when the initial data is homogeneous, we establish optimal growth on norms () of the Radon-Nikodym derivative of the Brownian LPP (i.e. top line of an -level Dyson Brownian motion) away from the origin, as the number of curves tends to infinity, for all sufficiently large. As an application of our framework, we show that the Radon-Nikodym derivative of certain toy models for the KPZ fixed point lies in , inspired by its variational characterisation in terms of the directed landscape.
Abstract
In this paper, we study sequences of perfect t-embeddings of a uniformly weighted family of graphs we call generalized tower graphs. We show that the embeddings of these graphs satisfy certain technical assumptions, in particular, the rigidity assumption of Berggren-Nicoletti-Russkikh. As a result, we confirm the convergence of the gradients of the height function fluctuations of these graphs to those of the Gaussian free field.
Abstract
We will describe solvable lattice models whose partition functions depend on two sets of variables, and that have different connections with the representation theory of where is a nonarchimedean local field. If the boundary conditions are chosen in one way, they are essentially the Motivic Chern classes that were used very effectively by Aluffi, Mihalcea, Sch\"urmann and Su (AMSS) to study such problems. In particular, using this specialization we can obtain deformations of the Kazhdan-Lusztig R-polynomials that were used by Bump, Nakasuji and Naruse to study matrix coefficients of intertwining operators (introduced by Casselman). Thus we are able see that the recursion formula for the is a reflection of the Yang-Baxter equation. On the other hand, with more general boundary conditions, specializing the parameters we recover colored lattice models that were previously used by Brubaker, Buciumas, Bump and Gustafsson to represent Iwahori Whittaker functions on . Thus we term the resulting two-variable-set family of functions as ``double Whittaker polynomials.''
Abstract
Understanding how a system loses memory of its initial state is a central problem in probability and statistics. In this manuscript, we introduce the notion of abrupt decorrelation, which explicitly characterises a sharp and sudden loss of correlation over time. We study this phenomenon within a class of linear stochastic differential equations (LSDEs), where explicit descriptions are available under various statistical distances. Our main focus is on the multivariate Ornstein-Uhlenbeck process, while in the one-dimensional case we extend the analysis to LSDEs with time-dependent drifts. The results highlight strong parallels with the cut-off phenomenon in Markov processes and contribute to a broader understanding of decorrelation in stochastic systems.
Abstract
We prove an optimal order delocalization estimate for the eigenvectors of general non-Hermitian matrices : with very high probability, for any right or left eigenvector of . This improves upon the previous tightest bound of Rudelson and Vershynin [arXiv:1306.2887] of , and holds under weaker assumptions on the tail of the matrix elements. In addition to the coordinate basis, our bound holds for the norm in any deterministic orthonormal basis. Our result is proven via a dynamical method, by studying the flow of the resolvent of the Hermitization of and proving local laws on short scales.
Abstract
We prove a homogenization result for the difference of two coupled Dyson Brownian motions started from generalized Wigner matrix initial data. We prove an optimal order, high probability estimate that is valid throughout the spectrum, including up to the spectral edges. Prior homogenization results concerned only the bulk of the spectrum. We apply our estimate to address the question of quantifying edge universality. Here, we have two results. We show that the Kolmogorov-Smirnov distance of the distribution of the gap between the largest two eigenvalues of a generalized Wigner matrix (with smooth entry distribution) and its GOE/GUE counterpart is . On the other hand, we show that, for the distribution of the largest eigenvalue, there are Wigner matrices so that the analogous Kolmogorov-Smirnov distance is bounded below by .
Abstract
We derive an exact formula for the probability that a Brownian path on an annulus does not disconnect the two boundary components of the annulus. The leading asymptotic behavior of this probability is governed by the disconnection exponent obtained by Lawler-Schramm-Werner (2001) using the connection to Schramm-Loewner evolution (SLE). The derivation of our formula is based on this connection and the coupling with Liouville quantum gravity (LQG). As byproducts of our proof, we obtain a precise relation between Brownian motion on a disk stopped upon hitting the boundary and the SLE loop measure on the disk; we also obtain a detailed description of the LQG surfaces cut by the outer boundary of stopped Brownian motion on a -LQG disk.
Abstract
An Alternating Sign Matrix (ASM) is a square matrix with entries in , and such that: in each row and columns, nonzero entries alternate in sign; for any given row or column, entries sum up to 1. We define the frozen-square enumeration as the enumeration of ASMs under the refinement of having, located in a corner, an square of entries that are all zeroes. We state a conjectural formula for such enumeration, in terms of a Fredholm type determinant of some matrix whose entries are given explicitely. We provide numerical support in favour of our conjecture. We also illustrate the relevance of the conjectured formula in connection with the limit shape observed in large ASMs, its fluctuations, and the Tracy--Widom distribution.
Abstract
We investigate a novel variant of the exclusion process in which particles perform asymmetric nearest-neighbor jumps across a bond \((k, k+1)\) only if the preceding site \((k-1)\) is unoccupied. This next-nearest-neighbor constraint significantly enriches the system's dynamics, giving rise to long-range correlations and a mixed-order transition controlled by the asymmetry parameter. We focus on the critical case of half filling, where the system splits into two ergodic components, each associated with an invariant reversible measure. The combinatorial structure of this equilibrium distribution is intimately connected to the \(q\)-Catalan numbers, enabling us to derive rigorously the asymptotic behavior of key thermodynamic quantities in the strongly asymmetric regime and to conjecture their behavior in the weakly asymmetric limit. Even though the system is one-dimensional and has short-range interactions, an equilibrium phase transition occurs between a clustered phase -- characterized by slow dynamics, long-range correlations with thermodynamic additivity, and spontaneous breaking of translational symmetry -- and a fluid phase where the correlations are short-range and which is thermodynamically additive. This equilibrium phase transition features characteristics of a first-order transition, such as a discontinuous order parameter as well as characteristics of a second-order transition, namely a divergent susceptibility at the transition point. We also briefly discuss density higher than one half where ergodicity is broken.
Abstract
In this article, we continue the development of the Riemann-Hilbert formalism for studying the asymptotics of Toeplitz+Hankel determinants with non-identical symbols, which we initiated in \cite{GI}. In \cite{GI}, we showed that the Riemann-Hilbert problem we formulated admits the Deift-Zhou nonlinear steepest descent analysis, but with a special restriction on the winding numbers of the associated symbols. In particular, the most natural case, namely zero winding numbers, is not allowed. A principal goal of this paper is to develop a framework that extends the asymptotic analysis of Toeplitz+Hankel determinants to a broader range of winding-number configurations. As an application, we consider the case in which the winding numbers of the Szeg\H{o}-type Toeplitz and Hankel symbols are zero and one, respectively, and compute the asymptotics of the norms of the corresponding system of orthogonal polynomials.
Abstract
We consider Jack measures on partitions with homogeneous defining specializations. For each of the six distinct classes of measures obtained this way we prove a global law of large numbers with an explicit limiting particle density. We also demonstrate for one of these classes how to obtain a global central limit theorem, global and edge large deviation principles, and edge universality using the results of the paper. Our argument is based on explicitly evaluating Jack symmetric functions at homogeneous specializations, relating the Jack measures to the discrete -ensembles from (Publications math{\' e}matiques de l'IH{\' E}S 125, 1-78, 2017) and using the discrete loop equations to substantially reduce computations.
Abstract
We investigate the asymptotic properties of permutations drawn from the Luce model, a natural probabilistic framework in which permutations are generated sequentially by sampling without replacement, with selection probabilities proportional to prescribed positive weights. These permutations arise in applications such as ranking models, the Tsetlin library, and related Markov processes. Under minimal assumptions on the weights, we establish a permuton limit theorem describing the global behavior of Luce-distributed permutations and derive an explicit density of the limiting permuton. We further compute limiting pattern densities and analyze the differences between exact Luce permutations and their permuton approximations. We also study the local convergence of these permutations, proving a quenched Benjamini--Schramm limit and a central limit theorem for consecutive pattern occurrences. Finally, we prove a central limit theorem for the number of inversions.
Abstract
The KPZ fixed point is a universal limiting space-time random field for the Kardar-Parisi-Zhang universality class. While the joint law of the KPZ fixed point at a fixed time has been studied extensively, the multipoint distributions of the KPZ fixed point in the general space-time plane are much less well understood. More explicitly, formulas were only available for the narrow wedge initial condition arXiv:1906.01053, arXiv:1907.09876 and the flat initial condition arXiv:1907.09876 for the multipoint distributions, and the half-Brownian and Brownian initial conditions arXiv:2010.07357v1, arXiv:2504.19975 for the two-point distributions. In this paper, we obtain the first formula for the space-time joint distributions of the KPZ fixed point with general initial conditions of compact support. The formula is obtained through taking KPZ scaling limit of the multipoint distribution formulas for the totally asymmetric simple exclusion process (TASEP). A key ingredient is a probabilistic representation, inspired by arXiv:1701.00018, of the kernel encoding the initial condition for TASEP, which was first defined through an implicit characterization in arXiv:1907.09876. Moreover, we also verify that the equal time version of our formula matches the path integral formula in arXiv:1701.00018 for the KPZ fixed point when the initial condition is of compact support.
Abstract
We consider the limiting fluctuations of the geodesic in the directed landscape, conditioning on its length going to infinity. It was shown in \cite{Liu22b,Ganguly-Hegde-Zhang23} that when the directed landscape becomes large, the geodesic from to lies in a strip of size and behaves like a Brownian bridge if we zoom in the strip by a factor of . Moreover, the length along the geodesic with respect to the directed landscape fluctuates of order and its limiting one-point distribution is Gaussian \cite{Liu22b}. In this paper, we further zoom in a smaller neighborhood of the endpoints when or , and show that there is a critical scaling window for the time, geodesic location, and geodesic length, respectively. Within this scaling window, we find a nontrivial limit of the one-point joint distribution of the geodesic location and length as . This limiting distribution, if we tune the time parameter to infinity, converges to the joint distribution of two independent Gaussian random variables, which is consistent with the results in \cite{Liu22b}. We also find a surprising connection between this limiting distribution and the one-point distribution of the upper tail field of the KPZ fixed point recently obtained in \cite{Liu-Zhang25}.
Abstract
The main result of this paper is the factorization of the holonomy -matrix for quantum at a root of unity into a product of four quantum holonomy dilogarithms. This factorization extends previously known results in this direction. We collect many existing results needed for the factorization. We use the holonomy -matrices to define representations of a groupoid of braids with flat connections, which also define invariants of knots with such connections.
August 2025
Abstract
In this paper, we consider characterisations of the class of unitary matrix integrals in terms of a first-order matrix linear differential equation for a vector function of size , and in terms of a scalar linear differential equation of degree . It will be shown that the latter follows from the former. The matrix linear differential equation provides an efficient way to compute the power series expansion of the matrix integrals, which with and are of relevance to the enumeration of longest increasing subsequences for random permutations, and to the question of the moments of the first and second derivative of the Riemann zeta function on the critical line, respectively. This procedure is compared against that following from known characterisations involving the -Painlev&\'e III second-order nonlinear differential equation. We show too that the natural generalisation of the unitary group integral permits characterisation by the same classes of linear differential equations.
Abstract
A hyperbinary partition of the nonnegative integer n is a partition where every part is a power of 2 and every part appears at most twice. We give three applications of the length generating function for such partitions, denoted by h_q(n). Morier-Genoud and Ovsienko defined the q-analogue of a rational number [r/s]_q in various ways, most of which depend directly or indirectly on the continued fraction expansion of r/s. As our first application we show that [r/s]_q = q h_q(n-1)/h_q(n) where r/s occurs as the nth entry in the Calkin-Wilf enumeration of the non-negative rationals. Next we consider fence posets which are those which can be obtained from a sequence of chains by alternately pasting together maxima and minima. For every n we show there is a fence poset F(n) whose lattice of order ideals is isomorphic to the poset of hyperbinary partitions of n ordered by refinement. For our last application, Morier-Genoud and Ovsienko also showed that [r/s]_q can be computed by taking products of certain matrices which are q-analogues of the standard generators for the special linear group SL(2,R). We express the entries of these products in terms of the polynomials h_q(n).
Abstract
We construct a family of invariant measures from the perspective of a shock in the KPZ fixed point. These measures are parameterized by a positive number , and are supported on functions satisfying . Each can be described as the sum of a Brownian motion and an independent Bessel- process with drift. We show that these measures appear as the limit of the measures constructed by Barraquand, Corwin, and Yang for the conjectural open KPZ fixed point on , after recentering by an appropriately defined shock location. Furthermore, we show that, with respect to the standard, deterministic recentering at , all extremal invariant measures for the KPZ fixed point are Brownian motions with drift. To do this, we first show that any extremal invariant measures must be supported on functions having fixed asymptotic slopes at . Using a one-force-one-solution principle from the work of Busani, Sepp\"al\"ainen, and the second author, this rules out all other invariant measures except those having left slope and right slope for some . To handle this case, we derive the limiting fluctuations of the shock for a special choice of initial condition. Additionally, we derive the limiting fluctuations of the shock for the case of the invariant measure from the perspective of a shock, and for the case of initial data .
Abstract
We consider the exponential last passage percolation (LPP) with thick two-sided boundary that consists of a few inhomogeneous columns and rows. Ben Arous and Corwin previously studied the limit fluctuations in this model except in a critical regime, for which they predicted that the limit distribution exists and depends only on the most dominant columns and rows. In this article, we prove their conjecture and identify the limit distribution explicitly in terms of a Fredholm determinant of a matrix kernel. This result leads in particular to an explicit variational formula for the one-point marginal of the KPZ fixed point for a new class of initial conditions. Our limit distribution is also a novel generalization of and provides a new numerically efficient representation for the Baik--Rains distribution.
Abstract
We consider the open KPZ equation on the interval with Neumann boundary conditions depending on parameters (the so-called maximal current phase). For and stationary initial conditions, we obtain matching upper and lower bounds on the variance of the height function for . Our proof combines techniques from arXiv:2111.03650, which treated the periodic KPZ equation, with Gibbsian line ensemble methods based on the probabilistic structure of the stationary measures developed in arXiv:2103.12253, arXiv:2105.15178, arXiv:2105.03946, arXiv:2306.05983, arXiv:2404.13444.
Abstract
For random integer matrices , the cokernels of the partial products naturally define a random flag of abelian -groups. We prove that as , this flag converges universally, for any nondegenerate entry distribution, to the Cohen-Lenstra type measure which weights each flag inversely proportional to the size of its automorphism group. As a corollary, we prove universality of certain formulas for the limiting conditional distribution of given in terms of Hall-Littlewood structure constants, which were previously obtained only for Haar matrices over . Our proofs combine the general technology of Sawin-Wood, matrix product moment computations following those of Nguyen-Van Peski, and the computation done previously for Haar -adic matrices by Huang.
Abstract
We establish asymptotic formulae for general joint moments of characteristic polynomials and their higher-order derivatives associated with matrices drawn randomly from the groups and in the limit as . This relates the leading-order asymptotic contribution in each case to averages over the Laguerre ensemble of random matrices. We uncover an exact connection between these joint moments and a solution of the -Painlev\'{e} V equation, valid for finite matrix size, as well as a connection between the leading-order asymptotic term and a solution of the -Painlev\'{e} III equation in the limit as . These connections enable us to derive exact formulae for joint moments for finite matrix size and for the joint moments of certain random variables arising from the Bessel point process in a recursive way. As an application, we provide a positive answer to a question proposed by Altu\u{g} et al.
Abstract
This paper explores the connection between perfect t-embeddings and the octahedron equation in the setting of the two-periodic Aztec diamond. In particular, we show that the positions of both the t-embedding and the corresponding origami map can be expressed as sums of density functions arising from solutions to the octahedron equation with appropriate flat initial conditions.
Abstract
Binomial formulas for Schur polynomials and Jack polynomials were studied by Lascoux in 1978, and Kaneko, Okounkov--Olshanski and Lassalle in the 1990s. We prove that the associated binomial coefficients are monotone and derive some symmetric function inequalities, in particular, a Schur positivity and Jack positivity result. These inequalities are similar to those studied by Newton, Muirhead, Gantmacher, Cuttler--Greene--Skandera, Sra and Khare--Tao.
Abstract
We develop the theory of the discrete periodic Pitman transform, first introduced by Corwin, Gu, and the fifth author. We prove that, for polymers in a periodic environment, single-path and multi-path partition functions are preserved under the action of this transform on the weights in the polymer model. As a corollary, we prove that the discrete periodic Pitman transform satisfies the same braid relations that are satisfied for the full-line Pitman transform shown by Biane, Bougerol, and O'Connell. Combined with a new inhomogeneous Burke property for the periodic Pitman transform, we prove a multi-path invariance result for the periodic inverse-gamma polymer under permutations of the column parameters. In the limit to the full-line case, we obtain a multi-path extension of a recent invariance result of Bates, Emrah, Martin, Sepp\"al\"ainen, and the fifth author, in both positive and zero-temperature.
Abstract
Under typical scaling, the last passage time field of the directed last passage percolation model with exponential site distributions converges to the KPZ fixed point. In this paper, we consider an atypical scenario in which the last passage time to a specific site is unusually large, and we explore how the last passage time field changes under this one-point upper large deviation event. We prove a conditional law of large numbers and compute the limiting fluctuations in certain regimes. Our proofs rely on an analysis of explicit multi-point distributions.
Abstract
We study the large-scale geometry of t-surfaces -- pairs of perfect t-embeddings and their associated origami maps -- arising from dimer models on Aztec diamonds with periodic edge weights. We prove that these t-surfaces converge to space-like maximal surfaces in the Minkowski space . We observe that the frozen and gas regions influence the geometry of the limiting surface in striking ways: all frozen regions collapse to four boundary points, regardless of the number of frozen regions, while each gas region collapses to a distinct light-like cusp in the interior of the surface. In the absence of gas regions, the limiting surface lies entirely within ; in the general case, however, this is no longer true. The limiting surface is sensitive to the detailed structure of the model: both the positions of the cusps, and the placement of the boundary vertices, depend on the precise way the edge weights are distributed on the Aztec diamond. Nevertheless, we show that the global conformal structure remains robust and coincides with the Kenyon-Okounkov conformal structure. We further conjecture that the cusp locations encode the shift in the discrete Gaussian component that appears in the global fluctuations of the dimer model.
Abstract
We describe some recently discovered connections between one-dimensional interacting particle models and Macdonald polynomials. The first such model is the multispecies asymmetric simple exclusion process (ASEP) on a ring, linked to the symmetric Macdonald polynomial through its partition function. Through this connection, a new formula was found for by generalizing multiline queues, which were introduced by Martin in 2018 to compute stationary probabilities of the ASEP. The second particle model is the multispecies totally asymmetric zero range process (TAZRP) on a ring, which was very recently found to have an analogous connection to the modified Macdonald polynomial through its partition function. This discovery coincided with a new formula for , this time in terms of tableaux with a queue inversion statistic, which also compute stationary probabilities of the TAZRP. We explain the plethystic relationship between multiline queues and queue inversion tableaux, and along the way, derive a new formula for using the queue inversion statistic. This plethystic correspondence is closely related to fusion in the setting of integrable systems.
Abstract
The real Ginibre spherical ensemble consists of random matrices of the form , where are independent standard real Gaussian matrices. The expected number of real eigenvalues is known to be of order . We consider the probability that there are real eigenvalues in various regimes. These are when is proportional to (large deviations), when is proportional to (intermediate deviations), and when is in the neighbourhood of the mean (local central limit theorem). This is done using a Coulomb gas formalism in the large deviations case, and by determining the leading asymptotic form of the generating function for the probabilities in the case of intermediate deviations (the local central limit regime was known from earlier work). Moreover a matching of the left tail asymptotics of the intermediate deviation regime with that of the right tail of the large deviation regime is exhibited, as is a matching of the right tail intermediate deviation regime with the leading order form of the probabilities in the local central limit regime. We also give the leading asymptotic form of , i.e. the probability of no real eigenvalues.
Abstract
Consider the random bipartite Erd\H{o}s-R\'{e}nyi graph , where each edge with one vertex in and the other vertex in is connected with probability , and for a constant aspect ratio . It is well known that the empirical spectral measure of its centered and normalized adjacency matrix converges to the Mar\v{c}enko-Pastur (MP) distribution. However, largest and smallest singular values may not converge to the right and left edges, respectively, especially when . Notably, it was proved by Dumitriu and Zhu (2024) that there are almost surely no singular value outside the compact support of the MP law when . In this paper, we consider the critical sparsity regime where for some constant . We quantitatively characterize the emergence of outlier singular values as follows. For explicit and as functions of , we prove that when , there is no outlier outside the bulk; when , outliers are present only outside the right edge of the MP law; and when , outliers are present on both sides, all with high probability. Moreover, the locations of those outliers are precisely characterized by a function depending on the largest and smallest degree vertices of the random graph. We estimate the number of outliers as well. Our results follow the path forged by Alt, Ducatez and Knowles (2021), and can be extended to sparse random rectangular matrices with bounded entries.
July 2025
Abstract
We investigate the asymptotic behavior of the q-Racah probability measure on lozenge tilings of a hexagon whose side lengths scale linearly with a large parameter , while the parameters and remain fixed. This regime differs fundamentally from the traditional case , in which random tilings are locally governed by two-dimensional translation-invariant ergodic Gibbs measures. In the fixed-q regime we uncover a new macroscopic phase, the waterfall (previously only observed experimentally), where the two-dimensional Gibbs structure collapses into a one-dimensional random stepped interface that we call a barcode. We prove a law of large numbers and exponential concentration, showing that the random tilings converge to a deterministic waterfall profile. We further conjecture an explicit correlation kernel of the one-dimensional barcode process arising in the limit. Remarkably, the limit is invariant under shifts by but not by , exhibiting an emergent period-two structure absent from the original weights. Our conjectures are supported by extensive numerical evidence and perfect sampling simulations. The kernel is built from a family of functions orthogonal in both spaces and , that may be of independent interest. Our proofs adapt the spectral projection method of Borodin-Gorin-Rains (arXiv:0905.0679) to the regime with fixed~q. The resulting asymptotic analysis is substantially more involved, and leads to non-self-adjoint operators. We overcome these challenges in the exponential concentration result by a separate argument based on sharp bounds for the ratios of probabilities under the q-Racah orthogonal polynomial ensemble.
Abstract
The parabolic Airy process is the Airy process minus a parabola, initially defined by its finite-dimensional distributions, which are given by a Fredholm determinant formula with the extended Airy kernel. This process is also the one-time spatial marginal of the KPZ fixed point with the narrow wedge initial condition. There are two formulas for the space-time multipoint distribution of the KPZ fixed point with the narrow wedge initial condition obtained by arXiv:1906.01053 and arXiv:1907.09876. Especially, the equal-time case of arXiv:1907.09876 gives a different formula of the multipoint distribution of the parabolic Airy process. In this paper, we present a direct proof that this formula matches the one with the extended Airy kernel. Some byproducts in the proof include several new formulas for the parabolic Airy process, and a generalization of the Andreief's identity.
Abstract
We prove that the Kipnis-Marchioro-Presutti (KMP) process converges to the Kardar-Parisi-Zhang (KPZ) equation, as time goes to infinity, in a properly scaled observation window shifted by . Our proof is based on identifying the KMP process with a stochastic flow of kernels describing transition probabilities in a certain model of random walk in space-time random environment. This allows to apply a recent result of arXiv:2401.06073 proving convergence of the density field of random walks in random environment to the KPZ equation in a suitably general sense.
Abstract
We introduce the notion of a Young generating function for a probability measure on integer partitions. We use this object to characterize probability distributions over integer partitions satisfying a law of large numbers and those that satisfy a central limit theorem. We further establish a multilevel central limit theorem, which enables the study of random standard Young tableaux. As applications of these results, we describe the fluctuations of height functions associated with (i) the Plancherel growth process, (ii) random standard Young tableaux of fixed shape, and (iii) probability distributions induced by extreme characters of the infinite symmetric group . In all cases, we identify the limiting fluctuations as a conditioned Gaussian Free Field.
Abstract
The paper analyzes four classical signal-plus-noise models: the factor model, spiked sample covariance matrices, the sum of a Wigner matrix and a low-rank perturbation, and canonical correlation analysis with low-rank dependencies. The objective is to construct confidence intervals for the signal strength that are uniformly valid across all regimes - strong, weak, and critical signals. We demonstrate that traditional Gaussian approximations fail in the critical regime. Instead, we introduce a universal transitional distribution that enables valid inference across the entire spectrum of signal strengths. The approach is illustrated through applications in macroeconomics and finance.
Abstract
Last passage percolation (LPP) in an lower triangular domain has nice connections with various generalizations of Schur measures. LPP along an anti-diagonal, from to , gives a distribution of a highest column of a random composition with respect to a Demazure measure (a non-symmetric analog of a Schur measure). LPP along a main diagonal, from to , is distributed as a marginal of a Pfaffian Schur process. In the first case we show that the asymptotics for the constant specialization is governed by the GOE Tracy-Widom distribution, in the second case - by the GSE Tracy-Widom distribution. In the latter case we were also able to study the truncated lower triangular case, obtaining an interesting generalization of the GSE Tracy-Widom distribution.
Abstract
We establish that the family of polymer measures associated with the Stochastic Heat Flow (SHF), indexed by , has a conditional Gaussian Multiplicative Chaos (GMC) structure. Namely, taking the random measure as the reference measure, we construct the path-space GMC with noise strength and prove that the resulting random measure is equal in law to . As two applications, we prove that the polymer measure and SHF tested against general nonnegative functions are almost surely strictly positive and that the SHF converges to as .
Abstract
We study the dimer model on the square grid, with quenched random edge weights. Randomness is chosen to have a layered structure, similar to that of the celebrated McCoy-Wu disordered Ising model. Disorder has a highly non-trivial effect and it produces an essential singularity of the free energy, with decay of dimer-dimer correlations, at a point of the ``liquid'' (or ``massless'') phase where the homogeneous dimer model has instead a real analytic free energy and correlations decaying like . Moreover, at a point where the homogeneous model has a transition between a massive (gaseous) and massless (liquid) phase, the critical exponent 3/2 (Pokrovsky-Talapov law), characteristic of the transition between the two regimes, is modified by disorder into an exponent that ranges continuously between 3/2 and infinity.
Abstract
We consider the half-space geometric Last Passage Percolation model starting with stationary measures. We obtain exact formulas for LPP value along the diagonal across the entire phase diagram. We also obtain the limits of these distributions under critical scaling which should yield the one-point distribution of the half-space KPZ fixed point starting from stationary initial conditions.
Abstract
We show that the Busemann process indexed by tilts in the super-differential of the limit shape exists and is unique in the strong sense in the i.i.d.\ planar corner growth model. This means that every probability space that supports the field of i.i.d.~weights supports a copy of the process and any two realizations of the process are equal almost surely.
Abstract
We study the asymptotic behavior of random domino tilings of the Aztec diamond of size in a random environment, where the environment is a one-periodic sequence of i.i.d. random weights attached to domino positions (i.e., to the edges of the underlying portion of the square grid). We consider two cases: either the variance of the weights decreases at a critical scale , or the distribution of the weights is fixed. In the former case, the unrescaled fluctuations of the domino height function are governed by the sum of a Gaussian Free Field and an independent Brownian motion. In the latter case, we establish fluctuations on the much larger scale , given by the Brownian motion alone. To access asymptotic fluctuations in random environment, we employ the method of Schur generating functions. Moreover, we substantially extend the known Law of Large Numbers and Central Limit Theorems for particle systems via Schur generating functions in order to apply them to our setting. These results might be of independent interest.
Abstract
We study lattice spin systems and analyze the evolution of Gaussian concentration bounds (GCB) under the action of probabilistic cellular automata (PCA), which serve as discrete-time analogues of Markovian spin-flip dynamics. We establish the conservation of GCB and, in the high-noise regime, demonstrate that GCB holds for the unique stationary measure. Additionally, we prove the equivalence of GCB for the space-time measure and its spatial marginals in the case of contractive probabilistic cellular automata. Furthermore, we explore the relationship between (non)-uniqueness and GCB in the context of space-time Gibbs measures for PCA and illustrate these results with examples.
Abstract
The QSSEP, short for quantum symmetric simple exclusion process, is a paradigm model for stochastic quantum dynamics. Averaging over the noise, the quantum dynamics reduce to the well-studied SSEP (symmetric simple exclusion process). These notes provide an introduction to quantum exclusion processes, focusing on the example of QSSEP and its connection to free probability, with an emphasis on mathematical aspects.
Abstract
We consider an -soliton solution of the focusing nonlinear Schr\"{o}dinger equations. We give conditions for the synchronous collision of these solitons. When the solitons velocities are well separated and the solitons have equal amplitude, we show that the local wave profile at the collision point scales as the function. We show that this behaviour persists when the amplitudes of the solitons are i.i.d. sub-exponential random variables. Namely the central collision peak exhibits universality: its spatial profile converges to the function, independently of the distribution. We derive Central Limit Theorems for the fluctuations of the profile in the near-field regime (near the collision point) and in the far-regime.
Abstract
Motivated by the dissipative abelian sandpile model, we analyze the trajectories of a one-dimensional random walk in a landscape of soft traps. These traps, placed at increasing distances from each other, correspond to dissipative sites in the associated dissipative abelian sandpile model. We identify a critical growth rate of the sizes of intervals between successive traps where there is a transition between finiteness and non-finiteness of the expected survival time of the random walk. This corresponds to a transition between non-criticality and criticality of the associated dissipative abelian sandpile model. Therefore, in this setting, we thus identify precisely how much dissipation can be added to the original abelian sandpile model in order to disrupt its criticality.
June 2025
Abstract
We study the spatial decorrelation of the solution to the KPZ equation with narrow wedge initial data. For fixed , we determine the decay rate of the spatial covariance function, showing that as . In addition, we prove that the finite-dimensional distributions of the properly rescaled spatial average of the height function converge to those of a Brownian motion.
Abstract
We establish a scaling limit result for the length of the longest increasing subsequence of a permutation of size sampled from the Brownian separable permuton of parameter , which is the universal limit of pattern-avoiding permutations. Specifically, we prove that \[\frac{\operatorname{LIS}(\sigma_n)}{n^\alpha}\;\underset{n\to\infty}{\overset{\mathrm{a.s.}}{\longrightarrow}}\; X,\] where is the unique solution in the interval to the equation \[\frac{1}{4^{\frac{1}{2\alpha}}\sqrt{\pi}}\,\frac{\Gamma\big(\tfrac{1}{2}-\tfrac{1}{2\alpha}\big)}{\Gamma\big(1-\tfrac{1}{2\alpha}\big)}=\frac{p}{p-1},\] and is a non-deterministic and a.s. positive and finite random variable, which is a measurable function of the Brownian separable permuton. Notably, the exponent is an increasing continuous function of with , and , which corresponds to the permuton limit of uniform separable permutations. We prove analogous results for the size of the largest clique of a graph sampled from the Brownian cographon of parameter .
Abstract
In this article we prove a sprinkled decoupling inequality for the stationary Hammersley's interacting particle process. Inspired by the work of Baldasso and Texeira (2018), and Hil\'ario, Kious and Texeira (2020), we apply this inequality to study two distinct problems on the top of this particle process. First, we analyze a detection problem, demonstrating that a fugitive can evade particles, provided that their jump range is sufficiently large. Second, we show that a random walk in a dynamic random environment exhibits ballistic behavior with respect to the characteristic speed of the particle system, under a weak assumption on the probability of being away of this critical speed.
Abstract
We prove that the stochastic heat flow Caravenna Sun Zygouras (2023), Tsai (2024) is a black noise in the sense of Tsirelson (2004). As a corollary, the 2d stochastic heat equation driven by a mollified spacetime white noise becomes asymptotically independent of that noise in the scaling limit.
Abstract
The large asymptotic expansion of the partition function for the normal matrix model is predicted to have special features inherited from its interpretation as a two-dimensional Coulomb gas. However for the latter, it is most natural to include a hard wall at the boundary of the droplet. We probe how this affects the asymptotic expansion in the solvable case that the potential is radially symmetric and the droplet is a disk or an annulus. We allow too for the hard wall to be strictly inside the boundary of the droplet. It is observed the term of order , has then a different rational number prefactor to that when the hard wall is at the droplet boundary. Also found are certain universal (potential independent) numerical constants given by definite integrals, both at order , and in the constant term.
Abstract
We provide a basis transformation that inverts the coordinate Bethe Ansatz. It is widely believed that the Bethe Ansatz is complete, based on numerical evidence and combinatorial arguments. We present a constructive and comprehensive approach to show that the Bethe Ansatz is complete. The transformation formula is verified numerically with excellent accuracy for non-zero anisotropy parameters and rings of odd length. Additionally, assuming the Bethe Ansatz is complete, we derive a novel exact formula for the one-point function systems through special identities for the Izergin-Korepin determinant.
Abstract
In i.i.d. exponential last-passage percolation, we describe the joint distribution of Busemann functions, over all edges and over all directions, in terms of a joint last-passage problem in a finite inhomogeneous environment. More specifically, the Busemann increments within a grid, and associated to different directions, are equal in distribution to a particular collection of last-passage increments inside a grid. The joint Busemann distribution was previously described along a horizontal line by Fan and the fourth author, using certain queuing maps. By contrast, our new description explicitly gives the joint distribution for any collection of edges (not just along a horizontal line) using only finitely many random variables. Our result thus provides an exact and accessible way to sample from the joint distribution. In the proof, we rely on one-directional marginal distributions of the inhomogeneous Busemann functions recently studied by Janjigian and the second and fourth authors. Another innovation of our argument is a joint invariance of inhomogeneous last-passage times under permutation of the inhomogeneity parameters. This result does not appear to be contained among the last-passage invariances recently described by Dauvergne. Our proof of invariance is also rather different, using the Burke property instead of the RSK correspondence, and leading to an explicit coupling of the weights before and after the permutation of the parameters.
Abstract
The dimer model is a classical statistical mechanics model which is exactly solvable in two dimensions, but about which little is known in higher dimensions. In analogy with large limits in lattice gauge theory, we study a large limit of the dimer model in any dimension . The dependence on comes from the multinomial tiling model introduced by Kenyon and Pohoata, which gives a general framework for adding a dependence on to a tiling model. We study the behavior of this model on periodic bipartite graphs in , in the scaling limit as the multiplicity and then the size of the graph go to infinity. In this iterated limit, in any dimension , we prove a variational principle and show that random configurations concentrate on a limit shape which is the unique solution to an associated system of Euler-Lagrange equations. The rate function of the variational principle is the integral of a surface tension function, which we can compute explicitly for lattices in any dimension as the Legendre dual of the free energy for the model on the torus. We give a unified methodology for computing the surface tension and Euler-Lagrange equations in any dimension . A new structure called the critical gauge also emerges in the large limit. We show that the critical gauge functions converges in the scaling limit to a limiting gauge function which is the unique solution to a dual Euler-Lagrange equation. This limiting gauge function determines the limit shape and vice versa. We further use our techniques to compute explicit limit shapes in some two and three dimensional examples, such as the Aztec diamond and "Aztec cuboid". This is one of the first stat mech models in dimensions where limit shapes can be computed explicitly.
Abstract
For each , , we show that there exists a unique -indexed line ensemble of random continuous curves with the following properties: (1) The top curve is distributed as the time- Cole--Hopf solution to the half-space KPZ equation with narrow wedge initial condition and Neumann boundary condition with parameter . (2) The line ensemble satisfies a one-sided resampling invariance property, involving softly non-intersecting Brownian motions with an attractive potential between pairs at the boundary. We call this object the half-space KPZ line ensemble. For with fixed (critical regime) and for fixed (supercritical regime), we show that the half-space KPZ line ensemble is tight under 1:2:3 KPZ scaling as . Moreover, all subsequential limits approximate a parabola and enjoy a one-sided Brownian Gibbs property, described by non-intersecting Brownian motions with pairwise interaction at the boundary. In the critical case this agrees with the half-space Airy line ensemble recently constructed by Dimitrov and Yang. In the supercritical case, we demonstrate a novel structure involving pairwise pinned Brownian motions, one of the main technical contributions of this paper.
Abstract
We study orthogonal polynomial ensembles whose weights are deformations of exponential weights, in the limit of a large number of particles. The deformation symbols we consider affect local fluctuations of the ensemble around a bulk point of the limiting spectrum. We identify the limiting kernel in terms of a solution to an integrable non-local differential equation. This novel kernel is the correlation kernel of a conditional thinned process starting from the Sine point process, and it is also related to a finite temperature deformation of the Sine kernel as recently studied by Claeys and Tarricone. We also unravel the effect of the deformation on the recurrence coefficients of the associated orthogonal polynomials, which display oscillatory behavior even in a one-cut regular situation for the limiting spectrum.
Abstract
We introduce a class of Markov processes conditioned to avoid intersection over a moving time window of length T>0, a setting we refer to as myopic non-intersection. In particular, we study a system of myopic non-intersecting Brownian motions subject to a periodic potential. Our focus lies in understanding the interplay between the confining effect of the potential and the repulsion induced by the non-intersection constraint. We show that, in the long time limit, and as both T and the strength of the potential become large, the model converges to a system of myopic non-intersecting random walks, which transitions between standard non-intersection dynamics and exclusion behavior. The main technical contribution of the paper is the introduction of an algorithm, based on a modification of the acceptance-rejection sampling scheme, that provides an explicit construction of myopically constrained systems.
May 2025
Abstract
In this paper, we study the asymptotic behaviour of plane partitions distributed according to a weighted Muttalib--Borodin ensemble. Specifically, we establish a Large Deviation principle for the discrete Muttalib--Borodin process, characterising the rate function. Furthermore, through Riemann--Hilbert analysis, we give an explicit expression for the asymptotic shape of the partition.
Abstract
We construct and analyze several new families of permutons arising from random processes involving the Demazure product on the symmetric group. First, we consider Demazure products associated to random pipe dreams, generalizing the Grothendieck permutons introduced by Morales, Panova, Petrov, and Yeliussizov by replacing staircase shapes with arbitrary order-convex shapes. Using the totally asymmetric simple exclusion process (TASEP) with geometric jumps, we prove precise scaling limit and fluctuation results for the associated height functions, showing that these models belong to the Kardar--Parisi--Zhang (KPZ) universality class. We then consider permutons obtained by applying deterministic sequences of bubble-sort operators to random initial permutations. We again provide precise descriptions of the limiting permutons. In a special case, we deduce the exact forms of the standard bubble-sort permutons, the supports of which were computed by DiFranco. A crucial tool in our analysis is a formulation, due to Chan and Pflueger, of the Demazure product as matrix multiplication in the min-plus tropical semiring. This allows us to define a Demazure product on the set of permutons. We discuss further applications of this product. For instance, we show that the number of inversions of the Demazure product of two independent uniformly random permutations of size is .
Abstract
We study the coupling of pairs of reverse plane partitions of the same shape by assigning a certain local interaction between the reverse plane partitions. We show that they are in bijection with a certain Yang-Baxter integrable colored vertex model. By utilizing the Yang-Baxter equation for this colored vertex model, we are able to compute the generating function for the interacting pairs of reverse plane partitions. We also give a bijection between the coupled pairs of reverse plane partitions with the interaction strength set to zero and a single reverse plane partition of the same shape.
Abstract
In the present context, superintegrability is a property of certain probability density functions coming from matrix models, which relates to the average over a distinguished basis of symmetric functions, typically the Jack or Macdonald polynomials. It states that the average can be computed according a certain combination of those same polynomials, now specialised by specific substitutions when expressed in terms of the power sum basis. For a particular generalisation of the Gaussian ensemble from random matrix theory, known independently from the consideration of certain integrable gauge theories, we use results developed in a theory of multivariable Al-Salam and Carlitz polynomials based on Macdonald polynomials to prove the superintegrability identity. This then is used to deduce a duality formula for these same averages, which in turn allows for a derivation of a functional equation for the spectral moments.
Abstract
Efficient matrix trace estimation is essential for scalable computation of log-determinants, matrix norms, and distributional divergences. In many large-scale applications, the matrices involved are too large to store or access in full, making even a single matrix-vector (mat-vec) product infeasible. Instead, one often has access only to small subblocks of the matrix or localized matrix-vector products on restricted index sets. Hutch++ achieves optimal convergence rate but relies on randomized SVD and assumes full mat-vec access, making it difficult to apply in these constrained settings. We propose the Block-Orthonormal Stochastic Lanczos Quadrature (BOLT), which matches Hutch++ accuracy with a simpler implementation based on orthonormal block probes and Lanczos iterations. BOLT builds on the Stochastic Lanczos Quadrature (SLQ) framework, which combines random probing with Krylov subspace methods to efficiently approximate traces of matrix functions, and performs better than Hutch++ in near flat-spectrum regimes. To address memory limitations and partial access constraints, we introduce Subblock SLQ, a variant of BOLT that operates only on small principal submatrices. As a result, this framework yields a proxy KL divergence estimator and an efficient method for computing the Wasserstein-2 distance between Gaussians - both compatible with low-memory and partial-access regimes. We provide theoretical guarantees and demonstrate strong empirical performance across a range of high-dimensional settings.
Abstract
We study the hydrodynamic limit for three gradient spin models: generalized Kipnis-Marchioro-Presutti (KMP), its discrete version and a family of harmonic models, under symmetric and nearest-neighbor interactions. These three models share some universal properties: occupation variables are unbounded, all these processes are of gradient type, their invariant measures are product with spatially homogeneous weights, and, notably, they are all attractive, meaning that the process preserves the partial order of measures along the dynamics. In view of hydrodynamics of large-scale interacting systems, dealing with processes taking values in unbounded configuration spaces is known to be a challenging problem. In the present paper, we show the hydrodynamic limit for all three models listed above in a comprehensive way, and show as a main result, that, under the diffusive time scaling, the hydrodynamic equation is given by the heat equation with model-dependent diffusion coefficient. Our novelty is showing the attractiveness for each model, which is crucial for the proof of hydrodynamics.
Abstract
The circular ensemble for and 4 corresponds to circular orthogonal, unitary and symplectic ensemble respectively as introduced by Dyson. The statistical state of the eigenvalues is then a determinantal point process () and Pfaffian point process (). The explicit functional forms of the correlation kernels then imply that the general -point correlation functions exhibit an asymptotic expansion in , which moreover can be lifted to an asymptotic in for the spacing distributions and their generating function. We use -Painlev\'e characterisations to show that the functional form of the first correction is related to the leading term via a second derivative. In the case this finding has immediate consequence in interpreting the empirical Riemann zeros spacing distribution at large height, and that of their thinning. Explicit functional forms are used to show that the spectral form factors for and 4 also admit an asymptotic expansion in . Differential relations are identified expressing the first and second correction in terms of the limiting functional form, and evidence is presented that they hold for general . For even it is proved that the two-point correlation function permits an asymptotic expansion in , and moreover that the leading correction relates to the limiting functional form via a second derivative.
Abstract
We give a complete characterization of all eternal solutions of the KPZ fixed point satisfying the asymptotic slope condition . For fixed , there is exactly one eternal solution with probability one. However, in the second and third authors' work with Sepp\"al\"ainen, it was shown that there exists a random, countably infinite set of slopes, for which there exist at least two eternal solutions. These correspond to two non-coalescing families of infinite geodesics in the same direction for the directed landscape. We denote the two eternal solutions as and . In the present paper, we show that, for the exceptional slopes, there are in fact uncountably many eternal solutions. To give the characterization, we show that these eternal solutions are in bijection with a certain set of bi-infinite competition interfaces. Each bi-infinite interface separates the plane into two connected components -- a left component and a right component. A general eternal solution with slope is equal to on the left component and equal to on the right component. For these bi-infinite interfaces in the exceptional directions, we uncover new geometric phenomena that is not present for directed landscape geodesics. Additionally, we show that this set of eternal solutions appears as the Busemann limits for sequences going to in direction .
Abstract
We study the integrability of the quantized six-vertex model with four parameters on a torus. It is a three-dimensional integrable lattice model in which a layer transfer matrix, depending on two spectral parameters associated with the homology cycles of the torus, can be defined not only on the square lattice but also on more general graphs. For a class of graphs that we call admissible, we establish the commutativity of the layer transfer matrices by introducing four types of tetrahedron equations and two types of inversion relations. Expanding in the spectral parameters yields a family of commuting quantum Hamiltonians. The quantized six-vertex model can also be reformulated in terms of (quantized) dimer models, and encompasses known integrable systems as special cases, including the free parafermion model and the relativistic Toda chain.
Abstract
We define the log-gamma sheet and the log-gamma landscape in terms of the 2-parameter and 4-parameter free energy of the log-gamma polymer model and prove that they converge to the Airy sheet and the directed landscape, which are central objects in the Kardar-Parisi-Zhang (KPZ) universality class. Our proof of convergence to the Airy sheet relies on the invariance of free energy through the geometric RSK correspondence and the monotonicity of the free energy. To upgrade the convergence to the directed landscape, tail bounds in both spatial and temporal directions are required. However, due to the lack of scaling invariance in the discrete log-gamma polymer--unlike the Brownian setting of the O'Connell-Yor model--existing on-diagonal fluctuation bounds are insufficient. We therefore develop new off-diagonal local fluctuation estimates for the log-gamma polymer.
Abstract
Santosh Kumar was an active researcher on the topic of exact results in random matrix theory and their various applications, particularly to quantum chaos and information theory. Barely entering his mid-career, he died unexpectedly on the 18th October 2024. The present article gives an account of some of his research directions and findings. As well as serving as a tribute to his work, this is done also for the purpose of providing a resource for those who may continue along related lines in the future.
Abstract
We construct a one-parameter family of infinite line ensembles on that are natural half-space analogues of the Airy line ensemble. Away from the origin these ensembles are locally described by avoiding Brownian bridges, and near the origin they are described by a sequence of avoiding reverse Brownian motions with alternating drifts, that depend on the parameter of the model. In addition, the restrictions of our ensembles to finitely many vertical lines form Pfaffian point processes with the crossover kernels obtained by Baik-Barraquand-Corwin-Suidan (Ann. Probab., 46(6), 3015-3089, 2018)
Abstract
We prove that a parabolically rescaled and suitably renormalised height function of a weakly asymmetric simple exclusion process on a circle converges to the Cole-Hopf solution of the KPZ equation. This is an analogue of the celebrated result by Bertini and Giacomin from 1997 for the exclusion process on a circle with any particle density. The main goal of this article is to analyse the interacting particle system using the framework of regularity structures without applying the Gaertner transform, a discrete version of the Cole-Hopf transform which linearises the KPZ equation. Our analysis relies on discretisation framework for regularity structures developed by Erhard and Hairer as well as estimates for iterated integrals with respect to cadlag martingales derived by Grazieschi, Matetski and Weber. The main technical challenge addressed in this work is the renormalisation procedure which requires a subtle analysis of regularity preserving discrete convolution operators.
April 2025
Abstract
We study the relationship between Yang-Baxter maps and the independence preserving (IP) property, motivated by their role in integrable systems, from the perspective of ultra-discretization. Yang-Baxter maps satisfy the set-theoretic Yang-Baxter equation, while the IP property ensures independence of transformed random variables. The relationship between these two seemingly unrelated properties has recently started to be studied by Sasada and Uozumi (2024). Ultra-discretization is a concept primarily used in the context of integrable systems and is an area of active research, serving as a method for exploring the connections between different integrable systems. However, there are few studies on how the stationary distribution for integrable systems changes through ultra-discretization. In this paper, we introduce the concept of ultra-discretization for probability distributions, and prove that the properties of being a Yang-Baxter map and having the IP property are both preserved under ultra-discretization. Applying this to quadrirational Yang-Baxter maps, we confirm that their ultra-discrete versions retain these properties, yielding new examples of piecewise linear maps having the IP property. We also explore implications of our results for stationary distributions of integrable systems and pose several open questions.
Abstract
This article studies the temporal law of the KPZ fixed point. For the stationary geometry, we find the two-time law, which extends the single time law due to Baik-Rains and Ferrari-Spohn. For the droplet geometry, we find a relatively simpler formula for the multi-time law compared to a previous formula of Johansson and the author. These formulas are derived as the scaling limit of corresponding multi-time formulas for geometric last passage percolation.
Abstract
We study a planar Coulomb gas confined to a sufficiently smooth Jordan arc in the complex plane, at inverse temperature . Let \[\bar{Z}_{n}^\beta(\gamma) = Z_{n}^\beta(\gamma)/\left(2 \textrm{cap}(\gamma)\right)^{\beta n^2/2+(1-\beta/2)n}\] be the normalized partition function. We compute the free energy as the number of particles tends to infinity, including the constant term: \[ \lim_{n\to \infty}\log \frac{\bar{Z}_{n}^\beta(\gamma)}{\bar{Z}_{n}^\beta([-1,1])} = \frac{1}{24}J^A(\gamma) + \frac{1}{8} \left( \sqrt{ \frac{\beta}{2} } - \sqrt{ \frac{2}{\beta}} \right)^2 J^{F}(\gamma). \] Here is an explicit Selberg integral, is half the Loewner energy of a certain Jordan curve associated to plus a covariance term, and is the Fekete energy, related to the zero temperature limit of the model. We also prove an asymptotic formula for the Laplace transform of linear statistics for sufficiently regular test functions. As a consequence, the centered empirical measure converges to a Gaussian field with explicit asymptotic mean, and asymptotic variance given by the Dirichlet energy of the bounded harmonic extension of the test function outside of the arc. A key tool in our analysis is the arc-Grunsky operator associated to , reminiscent to but different from the classical Grunsky operator. We derive several basic properties the arc-Grunsky operator, including an estimate analogous to the strengthened Grunsky inequality and the relation to the Dirichlet integral. In our proofs arises from the Fredholm determinant and a substantial part of the paper is devoted to relating this to the Loewner energy.
Abstract
We study Lam's reduced random walk in a hyperbolic triangle group, which we view as a random walk in the upper half-plane. We prove that this walk converges almost surely to a point on the extended real line. We devote special attention to the reduced random walk in (i.e., the triangle group). In this case, we provide an explicit formula for the cumulative distribution function of the limit. This formula is written in terms of the interrobang function, a new function that shares several of the remarkable analytic and arithmetic properties of Minkowski's question-mark function.
Abstract
We provide a combinatorial formula for the structure constants of spin Hall--Littlewood functions. This is achieved by representing these functions and the structure constants as the partition function of a lattice model and applying the underlying Yang--Baxter equation. Our combinatorial expression is in terms of generalised honeycombs; the latter were introduced by Knutson and Tao for ordinary Littlewood--Richardson coefficients and applied to the computation of Hall polynomials by Zinn--Justin.
Abstract
We consider the Kardar-Parisi-Zhang equation on the interval with Neumann type boundary conditions and boundary parameters . We show that the -th order cumulant of the height behaves as in the large time limit , and we compute the coefficients . We obtain an expression for the upper tail large deviation function of the height. We also consider the limit of large , with , , which should give the same quantities for the two parameter family KPZ fixed point on the interval. We employ two complementary methods. On the one hand we adapt to the interval the replica Bethe ansatz method pioneered by Brunet and Derrida for the periodic case. On the other hand, we perform a scaling limit using previous results available for the open ASEP. The latter method allows to express the cumulants of the KPZ equation in terms a functional equation involving an integral operator.
Abstract
We consider the lozenge tiling model on a large, finite hexagon. It is well-known that random lozenge tilings of the hexagon correspond to a two-dimensional determinantal point process via a bijection with ensembles of non-intersecting paths. The starting point of our analysis is a formula for the correlation kernel due to Duits and Kuijlaars which involves the Christoffel-Darboux kernel of a particular family of non-Hermitian orthogonal polynomials. Our main results are split into two parts: the first part concerns the family of orthogonal polynomials, and the second concerns the behavior of the boundary of the so-called arctic curve. In the first half, we identify the orthogonal polynomials as a non-standard instance of little -Jacobi polynomials and compute their large degree asymptotics in the regime. A consequence of this analysis is a proof that the zeros of the orthogonal polynomials accumulate on an arc of a circle and an asymptotic formula for the Christoffel-Darboux kernel. In the second half, we use these asymptotics to show that the boundary of the liquid region converges to the Airy process, in the sense of finite dimensional distributions, away from the boundary of the hexagon. At inflection points of the arctic curve, we show that we do not need to subtract/add a parabola to the Airy line ensemble, and this effect persists at distances which are in the tangent direction.
Abstract
We relate the combinatorics of Hall-Littlewood polynomials to that of abelian -groups with alternating or Hermitian perfect pairings. Our main result is an analogue of the classical relationship between the Hall algebra of abelian -groups (without pairings) and Hall-Littlewood polynomials. Specifically, we define a module over the classical Hall algebra with basis indexed by groups with pairings, and explicitly relate its structure constants to Hall-Littlewood polynomials at different values of the parameter . We also show certain expectation formulas with respect to Cohen-Lenstra type measures on groups with pairings. In the alternating case this gives a new and simpler proof of previous results of Delaunay-Jouhet.
Abstract
A minimal surface in a random environment (MSRE) is a -dimensional surface in -dimensional space which minimizes the sum of its elastic energy and its environment potential energy, subject to prescribed boundary values. Apart from their intrinsic interest, such surfaces are further motivated by connections with disordered spin systems and first-passage percolation models. In this work, we consider the case of strongly correlated environments, realized by the model of harmonic MSRE in a fractional Brownian environment of Hurst parameter . This includes the case of Brownian environment ( and ), which is commonly used to approximate the domain walls of the -dimensional random-field Ising model. We prove that surfaces of dimension delocalize with power-law fluctuations, and determine their precise transversal and minimal energy fluctuation exponents, as well as the stretched exponential exponents governing the tail decay of their distributions. These exponents are found to be the same in all codimensions , depending only on and . The transversal and minimal energy fluctuation exponents are specified by two scaling relations. We further show that surfaces of dimension delocalize with sub-power-law fluctuations, with their height and minimal energy fluctuations tied by a scaling relation. Lastly, we prove that surfaces of dimensions localize. These results put several predictions from the physics literature on mathematically rigorous ground.