Integrable Probability arXiv Feed
au:Borodin or au:"Alexei Borodin" — filter by authorin:JAMS or in:"Ann. Probab." — filter by journalcat:math.PR — filter by arXiv categoryy:2024 or y:2020-2024 — filter by yearTOP — show top-journal papers onlyClick any author, journal, or category badge to filter. Esc to clear.
Integrable probability papers on arXiv by members of the community (broadly understood), plus additional papers discovered by semantic search across the full arXiv database.
Not all papers by tracked authors are included — only those related to integrable probability.
The first version of this feed was imported from the NSF FRG page (opens in new tab); now updated via automated arXiv queries.
Please email me if you notice a missing author or any mistakes.
Subscribe via RSS
|
Download BibTeX of all papers (updated periodically)
May 2026
Abstract
We investigate the saturation regime of the condensing symmetric inclusion process on the discrete one-dimensional torus in the thermodynamical limit. In this regime, the total mass concentrates on a finite number of sites, forming condensates. Our main result establishes that, under appropriate scaling, the positions of the condensates converge to a system of coalescing Brownian motions on the continuum torus. In particular, condensates perform diffusive motion until they meet, at which point they merge and their masses coagulate. This provides a rigorous derivation of a macroscopic coalescing diffusion from an underlying interacting particle system with condensation. The main technical difficulty arises from the complicated coalescence mechanism of two condensates of particles, whose trajectories are very difficult to track completely. The key idea is to control the coalescing time instead and prove that it is negligible compared to the time-scale of condensate movement. By combining this with precise estimates of movements without coalescence, we can prove its convergence to coalescing Brownian motions.
April 2026
Abstract
We investigate the behavior of discrete interface growth models belonging to the Edwards--Wilkinson (EW) and Kardar--Parisi--Zhang (KPZ) universality classes, when defined on a complete graph, a topology commonly used to probe the infinite-dimensional limit of statistical mechanical systems. Our aim is to assess to what extent discrete lattice models reproduce the behavior of their corresponding continuum equations in this highly connected setting. After assessing the trivial behavior shown by some well known cases (like random deposition with surface relaxation or the etching model), we focus on two paradigmatic models associated with the KPZ universality class, the Restricted Solid-on-Solid (RSOS) and Ballistic Deposition (BD) models, and assess non-trivial behavior through several observables including the roughness, height fluctuations, power spectra, and two-time autocorrelation functions. Still, despite similarities with continuum equations, important differences arise in the fluctuations and long-time dynamics. In both discrete models the rescaled height fluctuations display a pronounced left tail, indicating the presence of lagging nodes. While the RSOS model exhibits a saturation roughness that decreases with system size, similarly to the EW and KPZ equations, the BD model exhibits a saturation roughness that increases with system size and an additional ultrafast growth regime, placing it outside the KPZ universality class on a complete graph.
Abstract
We construct the multilevel correlation kernel for the rising GUE eigenvalue process starting from a fixed initial configuration , and show that it converges on short time scales (as quickly as ) to the extended semi-discrete sine kernel. As an application, we show fixed-energy universality of bulk local statistics of complex Hermitian Wigner matrices matching the covariance structure of GUE and with a finite moment for . This application demonstrates that it is possible to obtain universality of bulk local statistics under near-optimal moment assumptions without using a Dyson Brownian motion relaxation step, which was a key ingredient in many results on this topic.
Abstract
Random matrices now play a role in many parts of computational mathematics. To advance these applications, it is desirable to have tools that are flexible, easy to use, and powerful. Over the last 25 years, researchers have developed a remarkable family of results, called matrix concentration inequalities, that meet the criteria. This paper offers an invitation to the field of matrix concentration and its multifarious applications.
Abstract
For a permutation , let be the permutation with scaled Lehmer code. For given and integer , the stretched Schubert coefficients are defined as . Our main result is that the function is eventually quasi-polynomial. This proves Kirillov's conjecture (2004), that the generating function for the sequence is rational. For the proof, we use combinatorics of pipe dreams to show that Schubert coefficients are given as an alternating sum of the numbers of integer points in certain polytopes. These polytopes behave nicely under stretching, and we use Ehrhart theory to obtain the result. As a consequence of the proof, we also present new counterexamples to the saturation conjecture for Schubert coefficients, and give computational applications.
Abstract
We study the central limit theorem (CLT) for linear eigenvalue statistics of several types of matrix models, whose entries are having exploding moments, i.e., moments of the entries are increasing with the size of the matrix. In particular, we study elliptic, centrosymmetric, circulant, and inter-correlated block matrices. The CLTs are established using asymptotic Wick formula.
Abstract
We study random matrices whose entries are obtained by applying consistent rank correlations, such as Hoeffding's , pairwise to a high-dimensional random vector with mutually independent components. Prior work has shown that, in the proportional high-dimensional regime, the empirical spectral distributions of large Kendall's tau and Spearman's rho matrices converge weakly almost surely to the Marchenko--Pastur law. By contrast, we prove that for consistent rank correlations such as Hoeffding's , the limiting spectral distribution is given by the semicircle law. Our result thus generalizes a recent work of Dong, Han, and Yao (2025), who considered the special case of Chatterjee's rank correlation and established the first semicircle law for a large correlation matrix in the proportional regime.
Abstract
Polynomial-in-time algorithms for computing classical Hurwitz numbers were given in [4] based on the Pandharipande equation. The paritition function of double Hurwitz numbers was proved [21] to satisfy the 2-Toda hierarchy. In this paper, similar to [21] we derive Pandharipande-type equations for double Hurwitz numbers from 2-Toda hierarchy. Based on these equations and a method from [4], we study large genus as well as large degree asymptotics of double Hurwitz numbers.
Abstract
This paper investigates the asymptotics of eigenstructure of sample covariance matrix under the spiked covariance matrix model in ultra-high-dimensional settings, where the dimensionality can grow much faster than the sample size with , . We establish the first-order convergence limits of eigenvalue locations and eigenvector projections of properly scaled sample covariance matrix. Our results are extensions of \cite{bloemendal16,ding21}.
Abstract
We study a natural analogue of Ulam's problem for random rooted trees distributed according to a Plancherel-type measure. This probability measure is closely related to the classical Plancherel measure on integer partitions. For a Plancherel random tree with vertices, we investigate the asymptotic behavior of its height , defined as the maximal distance from the root to a leaf. We prove that this height grows logarithmically. More precisely, there is a one-parameter family of random trees indexed by such that converges in probability to , where is an explicit constant depending on the parameter . The case of Plancherel trees corresponds to the parameter . The proof is based on the fact that the Plancherel random trees can be viewed as Ewens fragmentation trees, for which the height exhibits a sharp threshold phenomenon. An upper bound is obtained via -mass functionals and contraction estimates, while the lower bound is derived by embedding the model into a branching random walk with logarithmic displacements governed by a Poisson--Dirichlet distribution. The constant is characterized through a variational principle associated with this branching random walk.
Abstract
The evaluation homomorphisms from the super Yangian $\Ymn$ to the universal enveloping algebra allows one to regard the covariant tensor module of as $\Ymn$ modules. We study simple quotients of the submodules generated by a tensor product of highest weight vectors inside the tensor products of covariant evaluation modules. In the case , this recover all finite-dimensional simple modules of $\Y(\mathfrak{gl}_m)$. We give a necessary and sufficient condition for such modules to be tame, which generalizes the earlier work of Nazarov and Tarasov for $\Y(\mathfrak{gl}_m)$ to the super case.
Abstract
We introduce a linear map on symmetric functions that 'divides' a partition by a positive integer , sending a Schur function indexed by a partition of to a symmetric function indexed by partitions of . We prove that the image of this map is always Schur-positive, meaning it expands in the Schur basis with nonnegative integer coefficients. These coefficients are enumerated by a new family of combinatorial objects, called -Yamanouchi tableaux, which generalize the classical ballot (Yamanouchi) tableaux appearing in the Littlewood-Richardson rule. We also obtain expansions in the power-sum basis and establish power-sum positivity. A further application establishes a connection to work of Tewodros Amdeberhan, John Shareshian, and Richard Stanley on alternating permutations and Euler numbers.
Abstract
The Riesz gas in one-dimension consists of particles interacting via a pair potential, , and for . In the infinite density limit, with the particle support the interval , we apply a functional derivative method due to Beenakker to compute the covariance of two smooth linear statistics for the Riesz gas with exponent , . This we give in terms of a sum over Fourier components of the linear statistics with respect to a Gegenbauer polynomial basis, which generalises a known form in the case involving a cosine expansion. For the power sum linear statistic, our general formula can be reduced to a product of gamma function form, and compared against recent exact results in the literature for this case.
Abstract
We derive an identity for the determinant of the sum of two matrices, and , whose entries are defined via contour integrals. Specifically, we consider and . Under suitable assumptions on the functions , we show that can be expressed as a Fredholm determinant , where is an integral kernel acting on the contour . This result generalizes a recent identity obtained in \cite{Baik-Liao-Liu26}.
Abstract
For a fixed integer , consider representations of positive integers as sums of binomial coefficients of the form . While exact minimal bounds for the number of required summands are known only in a few low-dimensional cases, general existence results have received less explicit treatment. This paper provides: explicit elementary proofs for the cases () and (), a comparison with classical polygonal number theory, an explanation of why naive counting arguments fail for general (), conditional and unconditional existence results for general (), and a discussion of quantitative bounds and computational evidence. Together these give a unified and transparent framework for understanding additive representations by binomial coefficients.
Abstract
The random transposition shuffle on repeated cards induces a Markov chain on the quotient space of arrangements with multiplicities, and is equivalent to the many-urn mean-field Bernoulli-Laplace model introduced by Scarabotti. Writing , where there are card types and each type appears times, we determine the limiting profile for the total variation distance to stationarity at times , under the assumption . Scarabotti previously established that this process exhibits cutoff at time ; our result refines this by identifying the precise asymptotic shape of convergence inside the cutoff window. We show that the limiting profile is asymptotically Gaussian, with different explicit forms in the regimes fixed and . Together with our previous work on the fixed- regime, where the limiting profile is of Poisson type, this yields the cutoff profile for the random transposition shuffle on repeated cards for the full range of parameters and . Our argument has two main steps. First, we combine Scarabotti's Fourier-analytic framework for the many-urn Bernoulli-Laplace model with the approximation method of Jain-Sawhney (arXiv:2410.23944). More precisely, we compare the original shuffling measure with an explicitly tractable auxiliary measure directly on the repeated card quotient space, rather than passing through an intermediate comparison on the full symmetric group; this step relies in particular on our new estimates for Kostka numbers. Second, we reduce the limiting-profile problem to quotient fixed-point statistics and analyze them via Hoeffding-type combinatorial central limit theorems.
Abstract
We study Mallows random permutations conditioned to avoid a given pattern of length~. When the bias parameter is of the form , we prove that these permutations converge to a non-trivial explicit deterministic permuton that depends on the pattern and on the parameter . Along the way, we provide parametrizations for -avoiding permutons, and establish a large deviation principle for uniform -avoiding permutations. As a byproduct of the proof, we also obtain asymptotic estimates of two versions of -Catalan numbers in the regime .
Abstract
In [arXiv:2409.08465], Quastel and Gu use Stein's equation and integration by parts to give a direct proof that drifted Brownian motions are stationary (modulo height shifts) for the full-line KPZ equation. In this article, we consider the open KPZ equation with boundary conditions for a general real parameter , and emulate the approach of Quastel and Gu to provide a similar proof that Brownian motion with constant drift is invariant (modulo height shifts) in this case.
Abstract
In this paper, we propose and analyze a novel one-dimensional inhomogeneous random walk model that combines spatial decay of transition probabilities with a temporal renewal structure for each excursion. In this model, the probability of moving to the right from each state creats a spatial inhomogeneity that causes a stronger pull-back toward the origin as the process moves farther away. Furthermore, it features a random environment aspect where the parameter of each transition probability is independently resampled from a uniform distribution at the beginning of each excursion. We rigorously derive the hitting probability to an upper boundary using a scale function. Furthermore, by solving linear difference equations, we provide the probability generating function of the first hitting time, the expected occupation time for each state during an excursion (discrete Green's function), and the distribution and expectation of the maximum penetration depth.
Abstract
This work is a thorough investigation of skew-orthogonal polynomials with respect to a quartic Freud weight. We provide an explicit method to evaluate skew-orthogonal polynomials of any degree as linear combinations of orthogonal polynomials. The coefficients of these combinations can be evaluated via novel recursive relations. Moreover, we observe that skew-orthogonal polynomials with even and odd degree constitute two families of quasi-orthogonal polynomials with respect to two different semi-classical Laguerre weights, and we provide the first instance of closed recursive relations involving skew-orthogonal polynomials only.
Abstract
We study the vanishing-regularization limit of entropically regularized optimal transport (EOT) for the Euclidean distance cost in dimension . We develop a comprehensive variational convergence framework that entails two main results. First, we resolve the longstanding entropic selection problem: the EOT minimizer converges to a distinguished optimal transport plan that is characterized explicitly as the solution of a constrained EOT problem on each transport ray. Denoting by the regularization parameter, this selection holds for all -approximate minimizers, with sharp failure at the scale. Second, we establish an explicit second-order expansion of the entropic transport cost. The second-order term encodes the geometry of the regularization and reveals the optimal asymptotic tradeoff between entropy and transport cost.
Abstract
Internal DLA is a discrete random growth model describing growing clusters of particles. Its limiting shape and fluctuations are well understood when the underlying graph is the -dimensional lattice or the cylinder . In the latter geometry, the average fluctuations of IDLA have been shown to converge to the GFF. In this note we generalise this result by showing that, for any vertex-transitive base graph satisfying an eigenvalue convergence condition, the average fluctuations of IDLA on the cylinder are given by a GFF. On the way, we present an improved bound on the clusters' maximal fluctuations, which is of independent interest and which implies a shape theorem for IDLA on for any vertex-transitive base graph .
Abstract
The standard Monte Carlo estimator of relies on independent samples from and has variance of order . Replacing the samples with a determinantal point process (DPP), a repulsive distribution, makes the estimator consistent, with variance rates that depend on how the DPP is adapted to and . We examine two existing DPP-based estimators: one by Bardenet & Hardy (2020) with a rate of for smooth , but relying on a fixed DPP. The other, by Ermakov & Zolotukhin (1960), is unbiased with rate of order , like Monte Carlo, but its DPP is tailored to . We revisit these estimators, generalize them to continuous settings, and provide sampling algorithms.
Abstract
Linek's 1989 problem asks whether the numbers of independent sets of trees avoid infinitely many positive integers. We show that the set of natural numbers realized as the number of independent sets of a tree has a lower growth exponent of . We further prove that the set of positive integers representable by connected planar graphs has asymptotic density one. Lastly, we establish a phase transition: the number of independent sets of graphs with fewer than edges for any is contained in a set of density zero, whereas, following Shkredov's recent breakthrough on Zaremba's conjecture in continued fraction theory, there exists a constant such that the number of independent sets of graphs with at most edges covers all positive integers.
Abstract
The spiked Wigner ensemble is a prototypical model for high-dimensional inference. We study the spectral properties of an inhomogeneous rank-one spiked Wigner model in which the variance of each entry of the noise matrix is itself a random variable. In the high-dimensional limit, we derive exact equations for the spectral edges, the outlier eigenvalue, and the distribution of the components of the outlier eigenvector. These equations determine the BBP transition line that separates the gapped phase, where the signal is detectable, from the gapless phase. In the gapped regime, the distribution of the outlier eigenvector provides a natural estimator of the spike. We solve the equations for a noise matrix whose variances are generated from a truncated power-law distribution. In this case, the BBP transition line is non-monotonic, showing that an inhomogeneous noise can enhance signal detectability.
Abstract
In this paper we study detailed fluctuation results for a class of non-equilibrium steady states. The main example is the boundary driven harmonic model \cite{frassek2022exact}. In this model, the non-equilibrium steady state (NESS) is a mixture of products of geometric distributions, of which the local parameters are in turn distributed as uniform order statistics. For such a NESS, we prove law of large numbers, central limit theorem and large deviation results for fields of a general local functions (generalizing the density field). We also obtain quantitative results on the deviation from local equilibrium.
Abstract
The totally asymmetric simple exclusion process (TASEP) is a stochastic model for the unidirectional flow of interacting particles on a 1D-lattice that is much used in systems biology and statistical physics. Its master equation describes the evolution of the probability distribution on the configuration space. The size of the master equation grows exponentially with the length of the lattice. It is known that the complexity of the system may be reduced using mean-field approximations. We provide a rigorous definition of a family of such models using moments of any order and an extension to the pair approximation for obtaining closures for the system. The dimension of these models grows linearly with the lattice size and exponentially in the order of the approximation. Moreover, we show that the states of these models still have a probabilistic interpretation and that basic structural properties of the master equation are preserved. This extends known results on the Ribosome Flow Model which can be viewed as the first order approximation for TASEP.
Abstract
Kardar-Parisi-Zhang (KPZ) scaling has been observed in discrete polariton lattices, enabled by engineered band structures that stabilize the condensate. Whether this universality extends to intrinsically continuous systems with natural noise regularization remains an open question. We propose and numerically demonstrate KPZ scaling in a continuous quasi-one-dimensional polariton condensate stabilized by optical confinement in the transversal direction. Large-scale simulations of the stochastic Gross-Pitaevskii equation, with experimentally relevant parameters, reveal temporal and spatial scaling exponents of the two-point phase correlation function betaC = 0.30(5) and alfaC =0.46(8), and Tracy-Widom one-point phase fluctuation statistics, yielding robust KPZ dynamics intrinsic to the continuous polariton fluid.
Abstract
In this paper we consider the Toda lattice at thermal equilibrium, meaning that its variables and are independent Gaussian and Gamma random variables, respectively. We show under diffusive scaling that the space-time fluctuations for the model's currents converge to an explicit Gaussian limit. As consequences, we deduce, (i) the scaling limit for the trajectory of a single particle is a Brownian motion; (ii) space-time two-point correlation functions for the model decay inversely with time, with explicit scaling distributions predicted by Spohn (Spohn, J. Phys. A 53 (2020), 265004). Our starting point is the notion that the Toda lattice can be thought of as a dense collection of many ``quasi-particles'' that interact through scattering. The core of our work is to establish that the full joint scaling limit of the fluctuations for these quasi-particles is given by a Gaussian process, called a dressed Lévy-Chentsov field.
Abstract
For stochastic Hamilton-Jacobi (SHJ) equations, instability points are the space-time locations where two eternal solutions with the same asymptotic velocity differ. Another fundamental structure in such equations is shocks, which are the space-time locations where the velocity field is discontinuous. In this work, we study the KPZ fixed point, the central object of the KPZ universality class, which can be viewed as a prototype--albeit degenerate--of an inviscid SHJ equation in one spatial dimension. We describe the geometric structure of the instability region and give a detailed and precise analysis of its interplay with the shock structures of the two eternal solutions. We show that these shock structures allow one to reconstruct the instability region. Along the way, we obtain a complete classification of all possible configurations of semi-infinite geodesics emanating from arbitrary space-time points, in the directed landscape--the random environment in which the KPZ fixed point evolves.
Abstract
Hurwitz theory provides a large variety of enumerative problems related to algebraic geometry, mathematical physics, and combinatorics. We give a general framework to approach the large genus asymptotics of Hurwitz theory using only elementary methods and apply it to several types of Hurwitz numbers: single, double, or with an arbitrary numbers of fixed ramifications; simple and / or including completed cycles type of ramification and / or finitely many blocks of weakly monotone and / or strictly monotone types of ramifications. These, to the best of our knowledge, cover most of the Hurwitz numbers studied, and include for instance correlators of the HCIZ matrix model, Grothendieck dessins d'enfant, weighted Hurwitz numbers, and Gromov-Witten invariants of the Riemann sphere. We also apply our method to b-content Hurwitz numbers. As a specialisation, we recover some previously known about the large genus asymptotics of Hurwitz theory, namely classical results by Hurwitz and recent results of Do-He-Robertson, C. Yang, and results connected to recent work of X. Li.
Abstract
We consider the product of \(k_{n}\) independent \(n\times n\) complex Ginibre matrices and denote its eigenvalues by \(Z_{1},\ldots ,Z_{n}\). Let \(α= \lim_{n\to\infty} n / k_{n}\). Using the determinantal point process method, we reduce the study of extremal eigenvalues to the evaluation of determinants of certain \(n\times n\) matrices. In the modulus case, rotational invariance makes the relevant matrix diagonal, which yields a product representation in terms of Gamma tail probabilities. In the real-part case, the matrix is no longer diagonal; we handle this by a polar-coordinate reduction that introduces an independent uniform angle and leads to explicit formulas involving Gamma variables and trigonometric integrals. After appropriate rescaling, the spectral radius \(\max_{1\leq j\leq n}|Z_{j}|\) converges weakly to a nontrivial distribution \(Φ_α\) when \(α\in (0, +\infty)\), to the Gumbel distribution when \(α= +\infty\), and to the standard normal distribution when \(α= 0\). The family \(\{Φ_α\}_{α>0}\) extends continuously to the boundary regimes: \(Φ_α\) converges weakly to the standard normal law as \(α\to 0^{+}\) and to the Gumbel law as \(α\to +\infty\). Thus the three limiting regimes are connected by the single parameter \(α\), yielding a continuous transition from Gaussian to Gumbel distribution. For the spectral radius, we obtain the exact rates of convergence both in the fixed-\(α\) regime and at the boundaries \(α= 0\) and \(α= +\infty\). For the rightmost eigenvalue \(\max_{1\leq j\leq n}\Re Z_{j}\), we establish the convergence rates in the boundary regimes, while for \(α\in (0, +\infty)\) we show that the limiting distribution, though not available in closed form, still interpolates continuously between the normal and Gumbel laws.
Abstract
Let be an random matrix with entries in , where each row is independently and uniformly sampled from the set of all vectors in containing exactly ones, with for some fixed constant . A recent result of Tran states that the smallest singular value is bounded below by with high probability. In this note, we establish a complementary upper bound for , proving that \[ \forall \varepsilon >0 \qquad \mathbb{P}\left(s_n(M)\le \frac{\sqrt{d}}{\varepsilon^2 n}\right)\ge 1-C_p\left(\varepsilon+\frac{1}{\sqrt{d}}\right), \]where is a positive constant depending only on . This result confirms that the least singular value of dense random combinatorial matrices is typically of the order .
Abstract
Let be an signed random combinatorial matrix whose rows are independent and uniformly distributed over the set of -vectors with exactly zero coordinates. Despite the dependence induced by the row constraints, we prove that there exist constants such that for any , \begin{align*} \textbf{P}\left(s_{n}(M_n)\le {\varepsilon}{n^{-1/2}}\right)\le C\varepsilon+e^{-cn}. \end{align*} In particular, the probability that is singular is exponentially small. Our approach builds on the Combinatorial Least Common Denominator (CLCD) introduced by Tran and develops the method in the present constrained setting.
Abstract
We analyse the eigenvectors of the adjacency matrix of a random inhomogeneous graph constructed from a specified degree sequence. We assume that the empirical degree sequence has bounded mean and variance. We show that near the edges of the spectrum, the eigenvectors are semilocalized in the sense that their mass concentrates around a small set of resonant vertices. For the extremal eigenvalues, we establish localization around a single vertex. In order to obtain effective estimates in the presence of highly inhomogeneous degrees, we introduce a new economical pruning procedure that carefully extracts a forest from the original graph, whose adjacency matrix is compared to that of the original graph using a suitably constructed local coupling to random trees with independent edges.
Abstract
We present a new method for constructing Hadamard matrices that combines transformer neural networks with local search in the PatternBoost framework. Our approach is designed for extremely sparse combinatorial search problems and is particularly effective for Hadamard matrices of Goethals--Seidel type, where Fourier methods permit fast scoring and optimisation. For orders between and , it produces large numbers of inequivalent Hadamard matrices, and in harder cases it succeeds where local search from random initialisation fails. The largest example found by our method has order . In addition to these new constructions, our experiments reveal that the transformer can discover and exploit useful hidden symmetry in the search space.
Abstract
Let for some fixed , and let be an random matrix with entries in , where each row is independently and uniformly sampled from the set of all vectors in containing exactly ones. We show that the empirical spectral distribution of the appropriately rescaled matrix converges in probability to the circular law provided that . As a crucial element of the proof, we obtain quantitative lower bounds on the smallest singular value of the shifted matrices whenever and for some absolute positive constant .
Abstract
We prove that two half-space models in the KPZ universality class, exponential last-passage percolation and a family of Poisson-avoiding metrics generalizing colored TASEP, converge to a common scaling limit. This scaling limit is the directed landscape in half-space, a random directed metric in the half-plane indexed by a parameter which determines the strength of the boundary interaction. As part of our analysis, we characterize the half-space directed landscape in terms of the half-space KPZ fixed point, and prove convergence of geodesics. We also give an explicit construction of joint stationary measures (or horizons) in half-space for the log-gamma polymer, the KPZ equation, exponential and geometric last passage percolation, and the directed landscape itself.
Abstract
We study the direct and inverse spectral theory for a class of finite Hermitian banded matrices. Using the theory of matrix orthogonal polynomials, we provide an explicit procedure for reconstructing a banded matrix from a matrix-valued measure that encodes its spectral data. We establish necessary and sufficient conditions for a measure to be the spectral measure of a matrix in the examined class. We further analyze the connections between this spectral analysis, block tridiagonalization algorithms, and the Toda lattice evolution on banded matrices.
Abstract
We classify irreducible finite-dimensional modules of a collection of real Lie superalgebras that includes the simple ones, their classical variants, complex Lie superalgebras after restriction of scalars, and all real Lie algebras. Our strategy is to reduce this classification to determining the orbits of the parity and conjugation functors on irreducible modules of the complexifications of the aforementioned algebras. Then we provide explicit results for the computation of these orbits. For Lie superalgebras of basic type or of type , our classification applies to any highest-weight parametrization of irreducible complex modules with respect to an arbitrary Borel subalgebra. As a consequence, in the special case of real simple Lie algebras we obtain a new perspective on the classification of real simple modules and establish a conceptual connection with Kostant's cascade of strongly orthogonal roots.
Abstract
A color Lie algebra is a generalization of a Lie (super)algebra by an Abelian group . The underlying vector space and defining relations of the algebra are graded by , and the color Lie algebra admits graded Casimir elements. Furthermore, its loop algebra admits graded central extensions. We present a general method for constructing 2nd order graded Casimir elements and graded central extensions for a given color Lie algebra and its loop algebra, respectively. We also show that there exists a large class of color Lie algebras admitting such graded Casimir elements or central extensions by providing three examples, namely, for , and and for .
Abstract
In this paper, we study the long-time behavior of a stochastic heat equation with multiplicative noise and localized control. We begin by analyzing the uncontrolled dynamics and derive explicit decay rates for both mean-square and almost sure exponential stability. These estimates show that the two notions of stability may hold under different conditions on the parameters, reflecting the interplay between the drift and the multiplicative noise. We then introduce a finite-dimensional feedback control acting on a measurable subset of positive measure, built from finitely many Fourier modes of the solution. In particular, we show that the number of controlled modes determines the decay rate and allows for arbitrarily fast stabilization in the mean-square sense. As a consequence, almost sure exponential stability is recovered via a probabilistic argument, so that both notions of stability are achieved within the same framework and with the same decay rate. As an application, we provide a new proof of controllability for the stochastic heat equation based on an iterative construction of adapted controls in feedback form, avoiding the use of the adjoint equation.
Abstract
Consider critical Bernoulli percolation on for large; let be distinct points in . We prove that the probability that all lie in the same open cluster, rescaled by an appropriate power of , converges as to an explicit constant. This confirms a conjecture of Aizenman and Newman.
Abstract
Given a pair of finite posets , the function counting integer-valued order preserving extensions of an order preserving map is given by a piecewise polynomial in . We provide a criterion for the nonnegativity of the coefficients of these multivariate polynomials and apply it to show that marked order polytopes of skew shapes are Ehrhart positive in a multivariate sense. This extends recent results of Ferroni-Morales-Panova on order polytopes of skew shapes and proves conjectures on the Ehrhart positivity of skew Gelfand-Tsetlin polytopes and -generalized Pitman-Stanley polytopes due to Alexanderson-Alhajjar and Dugan-Hegarty-Morales-Raymond, respectively.
Abstract
We derive the Euler (hyperbolic) hydrodynamic limit for the directed exclusion process (DEP), a one-dimensional conservative interacting particle system that preserves particle-hole symmetry while breaking left-right symmetry. The proof relies on an explicit multi-process coupling, which guarantees a strong form of attractiveness and macroscopic stability for the particle system. Further open questions about DEP are briefly discussed.
Abstract
Consider the following nonlinear one-dimensional stochastic fractional heat equation where is the fractional Laplacian on for , and is a Gaussian noise that is white in time and behaves in space as a fractional Brownian motion with Hurst index satisfying . When , Hu and Wang ({\it Ann. Inst. Henri Poincaré Probab. Stat.} {\bf 58} (2022) 379-423) studied the well-posedness of the solution and its Hölder continuity, removing the technical condition that was previously assumed in Hu et al. ({\it Ann. Probab.} {\bf 45} (2017) 4561-4616). Their approach relied on working in a weighted space with a suitable power decay function. For the case , inspired by Hu and Wang, we investigate the well-posedness of the stochastic fractional heat equation without imposing the technical condition of , which was required in the earlier work of Liu and Mao ({\it Bull. Sci. Math.} {\bf181} (2022) 103207). In our analysis, precise estimates of the heat kernel associated with the fractional Laplacian play a crucial role.
Abstract
We study the -deformation of the Demazure product model from arXiv:2407.21653. Consider the longest element in written as a reduced word in simple transpositions. Independently delete each transposition with probability and apply the -Demazure product to the remaining ones. We show that the law of the resulting permutation converges as to a deterministic permuton, which coincides with the case studied in arXiv:2407.21653 for adjusted probability . This resolves Conjecture 1.13 from arXiv:2407.21653 and identifies the limiting permuton explicitly.
Abstract
In a recent work [Phys. Rev. E 109, L042102 (2024)], interesting dimensional crossovers [from two- to one-dimensional (2D to 1D) scaling] were found in the growth of Kardar-Parisi-Zhang (KPZ) interfaces on rectangular substrates, with lateral sizes . Here, we extend this study to other universality classes for interface growth -- specifically, the Edwards-Wilkinson (EW), the Mullins-Herring (MH), and the Villain-Lai Das Sarma (VLDS) classes. From extensive simulations, we demonstrate that, in all systems with sufficiently large aspect ratio , the roughness scales with time in the growth regime as for and for , where in most cases. For the VLDS class, this crossover is also observed in the height distribution (HD), which approaches its characteristic probability density function for the 2D case at short times () and then crosses over to the asymptotic 1D HD. Dimensional crossovers are also found in the steady state regime, both in the roughness scaling as well as in the VLDS HD, which interpolate between the 2D and 1D ones as increases. The particular case , with , is also discussed in detail and reveals interesting features of the investigated systems. For instance, there exist a `special' exponent such that the temporal crossover cannot be observed for . Moreover, this leads the saturation roughness to display a nonuniversal scaling: , with .
Abstract
Quasisymmetric functions in superspace were introduced as a natural extension of classical quasisymmetric functions involving both commuting and anticommuting variables. In this paper, we first provide a characterization of the algebra of quasisymmetric functions in superspace as an algebra of invariants under a quasisymmetrizing action of the symmetric group. Furthermore, we complete the superspace analogue of the classical hierarchy of combinatorial Hopf algebras by introducing the algebra of quasisymmetric functions in noncommuting variables in superspace. We endow this algebra with a Hopf superalgebra structure and thoroughly investigate its -basis and monomial basis, which are indexed by set supercompositions. By restricting to the minimal elements of the underlying poset, we construct the Hopf superalgebra of superpermutations, serving as the superspace analogue of the Malvenuto--Reutenauer algebra. We provide explicit product and coproduct formulas for these bases in terms of super-shuffles and global descents. Finally, via an abelianization morphism, we apply these noncommutative structures to derive a product formula for fundamental quasisymmetric functions in superspace.
Abstract
This article derives full asymptotic expansions for integrals of the form \[ \int_{0}^{1}f(u)(1+q\cdot u^{n})^{w/n}du \] as , with parameters real and , or positive for . We relate the coefficients of the asymptotic expansions to Nielsen's generalized polylogarithms. For , we obtain an expansion in terms of multiple zeta values, which in this setting, reduce to ordinary zeta values. A key point is that for , the integrals typically produce alternating multiple zeta values; we formulate a precise symmetry constraint on the relevant coefficient sequence under which all coefficients reduce to polynomials in ordinary zeta values. We also translate this symmetry into a statement about a binomial transform, and we verify the condition for several classical Appell-type families, like Euler, Bernoulli, Genocchi, and Hermite. Finally, we obtain precise results about the convergence of norms of random variables.
Abstract
We study transport dynamics of free fermions subject to the external monitoring of a conserved charge over an extensive region. Focusing on bipartition protocols, we consider monitoring the total particle number over half of the system, and study the profiles of local charges and currents at hydrodynamic scales. While the Lindbladian describing the averaged dynamics is non-local, we show that the profiles can be understood in terms of localized impurities. We present a general framework based on the generalized hydrodynamics (GHD) picture, allowing for a hybrid numerical-analytic solution of the quench dynamics at hydrodynamic scales. We illustrate our approach for domain-wall initial states, showing that monitoring leads to discontinuities in the profiles that become more pronounced as the rate increases and that lead to the absence of transport in the Zeno limit of infinite monitoring rates. Our GHD framework could be naturally extended to interacting systems, paving the way for a systematic study of transport of integrable models subject to extensive-charge measurements.
Abstract
We study the problem of describing the set of real functionals on the quotient of the ring of symmetric functions that are nonnegative on the images of certain modified Hall-Littlewood symmetric functions. This question is equivalent to the problem, posed in [Adv Math 395, p.108087 (2022)], of describing the set of coadjoint-invariant measures for unitary groups over a finite field in the infinite-dimensional setting. Our main results constitute partial progress towards this problem. Firstly, we show that the desired set of functionals is very large, in the sense that it contains explicit families of examples depending on infinitely many parameters. Secondly, we provide an analogue of Kerov's mixing construction that produces new sought after functionals from known old ones. This construction depends on an explicit "-twisted action" of on itself and the resulting dual map that makes into a comodule. Finally, our third main result explains the relation between the -twisted comultiplication and the usual comultiplication on .
Abstract
In this paper, we study estimation of parameters in a two-parameter Potts model with colors and coupling matrix . We characterize concrete sufficient conditions for existence of the pseudo-likelihood estimator of the Potts model, in terms of the local magnetic fields, and give sufficient conditions for the validity of the above characterization. We then provide sufficient criteria for estimation of both parameters at the optimal rate . In particular, if is the scaled adjacency matrix of a graph , then we show that joint estimation is possible if either has bounded degree or is irregular. In contrast, we give an example of a graph sequence which is approximately regular and dense, where no consistent estimator exists. We also show that one-parameter estimation at the optimal rate holds under much milder conditions when the other parameter is known. Along the way, we develop a concentration result for mean-field Potts models using the framework of nonlinear large deviations. Compared to the Ising case, our results for the Potts case require a novel analysis across multiple colors.
Abstract
The enumeration of lozenge tilings of hexagons with holes has been studied intensively in recent years. Researchers tried to find shapes and positions of holes in hexagonal regions so that the number of lozenge tilings of the resulting regions is given by a simple product formula. In the present work, we consider new regions that are hybrids of regions studied by the first author (hexagons with intrusions) and Ciucu (F-cored hexagons). Then, we show that the tiling generating functions of these new regions under a certain weight are given by simple product formulas. To give a proof, we present shuffling theorems for lozenge tilings of hexagons with intrusions, which give simple relations between the tiling generating functions of two related hexagonal regions with intrusions.
Abstract
P. Di Francesco first introduced the "Aztec triangle" in his study of the relationship between the twenty-vertex model and domino tilings. He conjectured an exact formula for the number of tilings of the Aztec triangle, and it has since been proved by several authors. In an attempt to prove the conjecture, M. Ciucu showed that the tiling number of the Aztec triangle divides the tiling number of a new region called the "cruciform region," a superposition of two Aztec rectangles. Ciucu proved that the number of domino tilings of a cruciform region is given by a simple product formula. In this paper, we generalize Ciucu's tiling formula by providing a generating-function formula for the cruciform region.
Abstract
Understanding the realization of thermal equilibrium through the thermalization process in a many-body system is a fundamental and complex scientific question, bridging thermodynamics and classical dynamics and connecting to a host of physical phenomena, such as mechanical instabilities in a thermal environment. In this work, based on the harmonic lattice model, we investigate the thermalization process in both velocity and coordinate spaces, by examining microscopic dynamics on the atomic level. We show the distinct relaxation rates of the transverse and longitudinal components of the velocity, reveal the power law governing the nonlinear proliferation of dominant frequencies, and observe the concurrent rapid proliferations of frequencies and topological defects. We also show that the lattice system's persistent out-of-plane deformations exhibit two-stage fluctuation behaviors, characterized by distinct power laws of fractional exponents and associated with the broken up-down symmetry. This work demonstrates the rich dynamics underlying the thermalization process, and advances our understanding on the dynamical adaptations of many-body systems to external disturbances.
Abstract
In the present paper, we study the equilibrium fluctuations of a particle system in infinite volume with two conserved quantities and long-range dependence. More specifically, the model of interest is the so-called ABC model, in which three types of particles (A, B and C) exchange their locations between and at a rate that depends on the type of particles involved and is proportional to for . After rigorously identifying the normal modes associated to the conserved quantities (the density of particles of types and , say), we prove that their fluctuations converge to independent fractional stochastic partial differential equations (SPDEs), which are either Gaussian or the Stochastic Burgers equation, and whose nature is determined by the microscopic range of dependence and the strength of the asymmetry.
Abstract
We consider a class of large-scale interacting systems with one conservation law satisfying the ``degree-preserving property'', and study the classification of their invariant measures and their hydrodynamic limits. Under a few basic conditions, we show that if the generator of the process preserves the degree of polynomials of the state variables up to two, then the marginals of any product invariant measure of the process must belong to one of six specific distributions. This classification result is essentially a consequence of a known result in statistics on univariate natural exponential families due to C.N. Morris, which we apply here for the first time in the context of microscopic stochastic systems. In particular, we introduce a new model whose invariant measure is given by the generalized hyperbolic secant distribution. Additionally, under the same conditions, we show that, regardless of the specific model, the hydrodynamic equation is always the classical heat equation, with a diffusion coefficient that depends on the model. Our proof is based on deriving uniform bounds on second-moments of state variables, whose proof is achieved by relating a correlation function to a one-dimensional random walk whose jump rates are model-dependent and obtaining sharp bounds on its occupation times on specific domains.
Abstract
We consider certain noncolliding interacting particle systems driven by Brownian noise. A key example is drifted Brownian motions conditioned not to intersect and related models of eigenvalues of Hermitian random matrices. We establish limit theorems for the extremal particle. We find: (i) the scaling limit of the largest eigenvalue of Brownian motion over Hermitian, positive-definite matrices, (ii) Airy process limit for the largest eigenvalue of Dyson's Brownian motion for GUE started from generic initial conditions, and (iii) a Fredholm determinant formula for the running maximum of the top path among noncolliding Brownian bridges and, as a byproduct, a new formula for the law of largest eigenvalue in the Laguerre Orthogonal Ensemble as well as for a related point-to-line last passage percolation model.
Abstract
We show that the increments of the KPZ fixed point started from arbitrary initial data are \emph{mutually} absolutely continuous with respect to Brownian motion with diffusion parameter on compacts, extending the one-sided Brownian absolute continuity relation of the KPZ fixed point established in \cite{sarkar2021brownian}. We also show that additive Brownian motion is absolutely continuous with respect to the centred Airy sheet on compacts, but it is not mutually absolutely continuous globally. As applications, we show that with probability strictly between zero and one, there exist record times of the KPZ fixed point away from any reference point, obtain a characterisation for the hitting probabilities of the graph of the KPZ fixed point to be positive in terms of a certain thermal capacity in the sense of \cite{watson1978corrigendum, watson1978thermal} and compute essential suprema of Hausdorff dimensions of these random intersections. Finally, we compute essential suprema of Hausdorff dimensions of images of subsets in the plane under the Airy sheet and give a condition for the positivity of their Lebesgue measure in terms of Bessel-Riesz capacity.
Abstract
We use random matrix theory for the Circular Unitary Ensemble (CUE) to study moments of derivatives of the Riemann zeta function shifted a small distance from the critical line. The corresponding CUE moments are studied in the limit of large matrix size in two regimes: when the spectral parameter is (1) suitably far inside the unit disc, and (2) at a small distance from the unit circle. In case (1), we obtain an asymptotic formula as a combinatorial sum over contingency tables, while in case (2) we obtain a sum over certain determinants with multiplicative coefficients given by Kostka numbers. The latter result is also valid exactly on the unit circle. Then, we consider the analogous problem for mean values of derivatives of the zeta function with suitable shifts. Assuming the Lindelöf hypothesis, we show that this mean value gives rise to the same sum over contingency tables obtained in the CUE. For sufficiently low-order moments, we establish this result unconditionally.
Abstract
We consider Erdős-Rényi-type random hypergraphs that are non-uniform, in the sense that hyperedges of different sizes may coexist, and inhomogeneous, in that connection probabilities may depend on the hyperedge size. All parameters are allowed to scale with the hypergraph size. We study the random adjacency matrix whose -entry counts the number of hyperedges containing both vertices and , and characterize its expected limiting spectral distribution in terms of the connection probabilities and the hyperedge sizes. We provide a Pastur-type condition, in the sense of Chatterjee (2005), under which the matrix can be Gaussianized, as well as a more restrictive but simpler sufficient condition in terms of the generalized average degree of the model. As a second main result, based on such a Gaussianization, we characterize the limiting spectral distributions under non-sparse conditions as semicircle laws with an explicit parametric variance. The latter can be expressed as a convex combination of the variances arising in the uniform cases, with coefficients determined by the trade-off between the different sources of inhomogeneity.
Abstract
For first passage percolation (FPP) on Euclidean lattices with , it is expected that the variance of the first passage time between two points grows sublinearly in the distance with a universal exponent strictly smaller than . Following Kesten's upper bound (Ann. Appl. Probab., 1993) on the variance, Benjamini, Kalai and Schramm (Ann. Probab., 2003) used hypercontractivity to obtain an improvement of a factor of when passage times take two values with equal probability. This was later extended to more general classes of passage time distributions. However, unlike in exactly solvable planar models in last passage percolation where the variance is known to be , the best known upper bound for the variance of passage times has remained in all non-trivial variants of FPP. For a class of rotationally invariant Riemannian FPP on the plane, we show that the variance is for some . Our argument uses fluctuation estimates for passage times and geodesics derived in Basu, Sidoravicius and Sly (2023) together with a multi-scale argument to establish that the geodesic exhibits disorder chaos, i.e., upon resampling a small fraction of the underlying randomness, the updated geodesic has on average a small overlap with the original one; this, established at a large number of scales, leads to a polynomial improvement of the variance bound.
Abstract
We obtain bivariate asymptotics for one part monotone Hurwitz numbers in high genus (i.e. as both the size and the genus go to infinity). To do so, we start with a linear recurrence for these numbers obtained by Do and Chaudhuri. Then, we apply a recent method developped by Elvey-Price, Fang, Wallner and the second author to extract asymptotics from such recurrences.
Abstract
This work establishes a large deviation principle for the spectral measure of the Lax matrix associated to the periodic Toda chain of particles, subject to a generalised Gibbs measure. This large deviation principle is governed by a rate function which can be regarded as a generalisation of the free energy of the system. Such a large deviation principle is proven both for the model when the momentum is constrained to be zero and when it is allowed to fluctuate. Moreover, the large deviation principle is proven directly on the level of the representation of the generalised Gibbs partition function given in terms of the variables realising the classical separation of variables, \textit{i.e.} rectify the equations of motion. As such, this work paves the way towards the computation of the thermodynamic limit of dynamical correlation functions in the Toda chain subject to generalised Gibbs ensemble statistics.
March 2026
Abstract
In this paper, we introduce and study random Schur measures whose parameters are sampled from the Circular Unitary Ensemble. We show that Schur measures with CUE disorder exhibit behavior reminiscent of spin glasses.
Abstract
Fully inhomogeneous spin Hall-Littlewood symmetric rational functions arise as partition functions of certain path configurations in the higher spin six vertex models. They are multiparameter generalizations of the classical Hall-Littlewood symmetric polynomials. We establish two new generalizations of the classical Littlewood identity, where we express a weighted sum of 's over all partitions as a product of the Littlewood kernel and another simple product in one case, and a product of the Littlewood kernel and a Pfaffian in the other case. As a corollary we obtain a novel Littlewood identity for Hall-Littlewood symmetric polynomials. We also elaborate on the newly established connection between the fully inhomogeneous spin Hall-Littlewood symmetric rational functions and the modified Robbins polynomials, the latter being multivariate generating functions for alternating sign matrices. This connection allowed us to discover the two generalizations of the Littlewood identity and we provide a bijection between the underlying combinatorial models in the case where is strictly decreasing.
Abstract
The conjectured three generic local bulk statistics amongst all non-Hermitian random matrix symmetry classes have recently been extended to three generic local edge statistics. We study analytically and numerically complex spacing ratios and nearest-neighbour (NN) spacing distributions that characterise such local statistics. We choose the three simplest representatives of these universality classes, given by the Gaussian ensembles of complex Ginibre, complex symmetric and complex self-dual matrices, denoted by class A, AI and AII. In the first part, we analytically study the complex spacing ratio in class A, at finite matrix size . Introducing a conditional point process, we simplify existing expressions and show why an uncontrolled approximation introduced earlier converges well in the large- limit in the bulk. When specifying to the elliptic Ginibre ensemble, we present a parameter-dependent surmise for the complex spacing ratio, interpolating to that of the Gaussian unitary ensemble (GUE), where such a surmise is very accurate. In the second numerical part, we compare complex spacing ratios, its moments, and NN spacing distributions for all three ensembles with that of uncorrelated points, the two-dimensional (2D) Poisson process, both in the bulk and at the edge. The varying degree of repulsion within these different edge universality classes can be well understood in terms of an effective 2D Coulomb gas description, at different values of inverse temperature . We find indications that the complex spacing ratio does not fully unfold the local statistics at the edge. Finally we verify that for small argument, in all three symmetry classes the NN spacing distributions in the bulk and at the edge are consistent with the universal cubic repulsion.
Abstract
We study a family of periodically weighted Aztec diamond dimer models near their turning points. We establish that, asymptotically, as , their fluctuations there, scaled by , are described by a marked GUE-corners process. This limiting point process is constructed by assigning a Bernoulli mark independently to each particle in a realization of the GUE-corners process. The Bernoulli parameters associated with the random marks reflect the periodicity of the model in the limit. To prove this result we use a double-contour integral representation of the inverse Kasteleyn matrix on a higher-genus Riemann surface, which is well-suited for asymptotic analysis.
Abstract
Two-dimensional Coulomb gases on an annulus at a special inverse temperature are studied by using the orthogonal polynomial method borrowed from the theory of random matrices. The correlation functions among the Coulomb gas molecules are written in determinant forms and their asymptotic forms in the thermodynamic limit are evaluated. When the Coulomb gas system has a continuous rotational symmetry, the corresponding orthogonal polynomials are monomials, and one can see a universal behavior of the correlation functions in a thin annulus limit. In a system with a discrete rotational symmetry, the corresponding orthogonal polynomials are not in general monomials, and a breakdown of the universality is observed.
Abstract
In this article, we study a model of random permutations, which we call random standardized permutations, based on a sequence of i.i.d. random variables. This model generalizes others, such as the riffle-shuffle and the major-index-biased permutations. We first establish an exact result on the joint distribution of the number of cycles of given lengths, involving the notion of primitive words. From this result, we obtain various convergence results, most of which are proved using the method of moments. First we prove that the number of small cycles may have either a Poisson limit distribution, or a limit distribution given by a countable sum of independent geometric distributions. Then we establish a limit distribution for large cycles, which is the Poisson-Dirichlet process. Finally we prove a central limit theorem for the total number of cycles.
Abstract
Matrix valued (asymmetric) truncated Toeplitz operators are generally not complex symmetric. In this paper, we define a new conjugation with unique properties and study its relation to matrix valued asymmetric truncated Toeplitz operators. We also explore the connections between matrix valued asymmetric truncated Toeplitz operators and Hankel operators with matrix symbols.
Abstract
Recent work of Bornemann has uncovered hitherto hidden integrable structures relating to the asymptotic expansion of quantities at the soft edge of Gaussian and Laguerre random matrix ensembles. These quantities are spacing distributions and the eigenvalue density, and the findings cover the cases of the three symmetry classes orthogonal, unitary and symplectic. In this work we give a different viewpoint on these results in the case of the soft edge scaled density, and in the Laguerre case we initiate an analogous study at the hard edge. Our tool is the scalar differential equation satisfied by the latter, known from earlier work. Unlike integral representations, these differential equations in soft edge scaling variables isolate the function of which is the expansion variable. Moreover, they give information on the correction terms which supplements the findings from the work of Bornemann. In the case of the Gaussian ensemble, we can demonstrate analogous features for Dyson index , which suggests a broader class of models, namely the classical ensembles, with asymptotic expansions exhibiting integrable features. For the Laguerre ensembles at the hard edge, we give the explicit form of the correction at second order for unitary symmetry, and at first order in the orthogonal and symplectic cases. Various differential relations are demonstrated.
Abstract
On the lattice, vertices are assigned random weights . The point-to-point last passage percolation (LPP) time between and is the maximum total weight among all upward/right-directed paths connecting the two. Point-to-line LPP time is the maximum of these maximal total weights over . Asymptotic distributions and fluctuations of these LPP times have been studied for i.i.d. weights. The current study deals with identically distributed but not necessarily independent weights, and maximizes LPP times in the sense of convex dominance. In particular, maximal expected LPP times are identified, in the class of all weight couplings with a given marginal distribution. For the case of mean- exponentially distributed weights, there is a coupling for which is the shifted exponential variable , such that for all couplings and all convex non-decreasing functions for which these expectations are well defined. In contrast to , with variance and mean diverging to like , converges to for the commonly studied i.i.d. weights.
Abstract
It is shown that a generalization of the Painlevé-II equation (P-II) to a system of coupled equations with symmetry breaking terms is integrable. A Lax pair for this system is used to relate the asymptotic behavior of the solutions at different infinities via an asymptotically exact WKB approach. The analysis relies on an exact solution of the quantum mechanical Demkov-Osherov model (DOM). An application to the problem of unstable vacuum decay during a second order phase transition provides precise scaling of the number of excitations, including subdominant contributions.
Abstract
The Schur -, -multiple zeta functions were defined by Nakasuji and Takeda inspired by the tableau representation of Schur -, -functions. While a product of two Schur -functions expands as a linear combination of Schur -functions, we obtain a similar expansion formula for the Schur -multiple zeta functions by taking summation over the symmetric group permutating all the variables. We also introduce a expansion formula of skew Schur -multiple zeta functions by taking summation over the symmetric group. Furthermore, this skew type formula can be refined by restricting the symmetric group to its specific subgroup.
Abstract
This article has two main goals. First, it provides a self-contained exposition of the tangent plane method for the dimer model - a technique for analyzing arctic curves and limit shapes introduced by R. Kenyon and I. Prause (2020). Second, it extends this method to multiply connected domains through a nontrivial computation of the frozen boundary for the Aztec diamond with a hole. This computation yields the first explicit parametrization in terms of elliptic functions of a family of arctic curves of a multiply-connected region indexed by the height change (hole height). We also derive and visualize the corresponding limit height functions.
Abstract
We prove determinantal-Pfaffian formulae that simultaneously generalise the Pfaffian minor summation formula of Ishikawa and Wakayama and Byun's recent minor summation formula. These formulae are based on factorisation formulae for the determinant of the sum of a skew-symmetric matrix and a rank-1 matrix. Applications include a Cauchy-type identity for skew Schur functions.
Abstract
We present computational results on principal specializations of Schubert polynomials, which count reduced pipe dreams and reduced bumpless pipe dreams (RBPD). We find the first counterexample, at , to the Merzon-Smirnov conjecture (arXiv:1410.6857) that the maximum of over is attained at a layered permutation. The simulations suggest that equals the maximal layered permutations' constant from Morales-Pak-Panova (arXiv:1805.04341). We also explore the random permutation drawn from the distribution proportional to , revealing permuton-like asymptotics similar to those for Grothendieck polynomials by Morales-Panova-Petrov-Yeliussizov (arXiv:2407.21653). We implement and compare three recurrences for : the descent formula (Macdonald), transition formula (Lascoux--Schutzenberger), and cotransition formula (Knutson). For sampling uniformly random RBPDs (whose count is ), we show that reducedness breaks the sublattice property of the ASM lattice, preventing monotone CFTP and causing false coalescence. We develop an efficient MCMC sampler with macroscopic "droop" updates for connectivity and fast mixing. Our code computes up to and samples random RBPDs up to on a personal computer ( on a cluster).
Abstract
We investigate the stochastic dynamics of a quasiparticle within a gas of hard rods, focusing on the evolution of its mean, variance, and autocorrelation for two choices of initial states: (i) one with long-range (LR) correlations and (ii) the other without it. We derive analytical results for the phase space density correlations in the former case to complement the known results for the latter case. These results enable us to obtain expressions for the mean, variance, and autocorrelation of a quasiparticle, which are applicable to both initial states. The LR correlations introduce a diffusive-scale correction to the mean Euler generalized hydrodynamic (GHD) equations, modifying the standard local equilibrium form, and our findings reveal that the form of the correction term depends on the LR correlations present in the initial state.
Abstract
We consider the elliptic Ginibre ensembles in the real, complex and symplectic symmetry classes. As the matrix size tends to infinity, we derive the asymptotic behaviour of the upper tail large deviation probabilities for both the spectral radius and the rightmost eigenvalue. More generally, we obtain asymptotic formulas for the probability that an eigenvalue is found in a prescribed region outside the support of the elliptic law, thereby providing a unified framework in which the results for the spectral radius and the rightmost eigenvalue appear as special cases. The key ingredient of our analysis is the precise asymptotic behaviour of the associated one-point functions, which is of independent interest.
Abstract
We revisit the transfer-matrix approach to directed polymers in random media and show that a single ensemble of random transfer-matrix products provides a unified realization of the canonical one-point fluctuation laws in dimensions. For a fixed disorder realization, the polymer partition function is obtained as a contraction of the same product matrix , and different contractions reproduce the standard KPZ subclasses: Tracy-Widom GUE (point-to-point), GOE (point-to-line), GSE (half-space point-to-point), and Baik-Rains (stationary line-to-point). In each case, we observe free-energy fluctuation growth and convergence of standardized distributions with low-order cumulants close to the corresponding universal benchmarks. Viewing geometry-dependent subclasses as projections of a single matrix-product ensemble naturally suggests additional observables intrinsic to . As an example, we examine the leading eigenvalue whose logarithm exhibits scaling, while its standardized statistics remain distinct from the canonical Tracy-Widom laws within the accessible range. This transfer-matrix perspective thus organizes known KPZ one-point subclasses within a finite-dimensional matrix framework and highlights matrix-level fluctuation observables beyond geometry-selected universality classes.
Abstract
We determine the phase diagram of the Anderson tight-binding model on regular random graphs with Gaussian disorder distribution and sufficiently large degree. In particular, we prove that if the degree is fixed and the number of vertices goes to infinity, the spectrum asymptotically consists of a finite delocalized interval surrounded by two unbounded localized components. Our argument uses a recent description of the spectrum of the tight-binding model on the Bethe lattice (Aggarwal--Lopatto, 2025). By viewing the Bethe lattice as the local limit of a random regular graph, and proving several preliminary estimates, we transfer this characterization of the spectrum of the limiting model to its finite approximants.
Abstract
We study the fixed-time spatial covariance of the KPZ equation with flat initial profile. Using Malliavin calculus and a Clark-Ocone representation, we show that as , is governed by a boundary-layer regime near the initial time and satisfies as , where , is the flat stochastic heat equation solution, and is the one-dimensional heat kernel. In sharp contrast with the narrow-wedge regime, where Gu-Pu (2025, Theorem 1.1) proved that for each fixed , as , the flat initial profile exhibits Gaussian decay, yielding, to the best of our knowledge, the first exact spatial covariance asymptotic for the KPZ equation under flat initial data. We also establish an explicit closed-form formula for the second moment of the continuum directed random polymer partition function.
Abstract
We generalize well-known bijections between alternative tableaux and permutations to bijections between rhombic alternative tableaux (RAT) and assemblées of permutations. We show how these various bijections are connected. As a consequence, we find a refined enumeration formula for RAT. One of our bijections carries many statistics from RAT to assemblées; notably, it sends the number of free cells to the number of crossings, which answers a question of Mandelshtam and Viennot. We also find an -to- map from marked Laguerre histories to assemblées, answering a question of Corteel and Nunge.
Abstract
Duality transformations reveal unexpected equivalences between seemingly distinct models. We introduce an out-of-equilibrium generalisation of matrix product operators to implement duality transformations in one-dimensional boundary-driven Markov processes on lattices. In contrast to local dualities associated with generalised symmetries, here the duality operator intertwines two Markov processes via generalised exchange relations and realises the out-of-equilibrium duality globally. We construct these operators exactly for the symmetric simple exclusion process with distinct out-of-equilibrium boundaries. In this case, out-of-equilibrium boundaries are dual to equilibrium boundaries satisfying Liggett's condition, implying that the Gibbs-Boltzmann measure captures out-of-equilibrium physics when leveraging the duality operator. We illustrate this principle through physical applications.
Abstract
The asymmetric simple exclusion process (ASEP) is a fundamental stochastic model describing asymmetric many-particle diffusion with hard-core interactions on a one-dimensional lattice, and has been widely applied in the study of nonequilibrium transport phenomena. Motivated by the modeling of proton transport along oxygen networks in proton-conducting solid oxides, we extend the ASEP to a model defined on a one-dimensional backbone lattice with tree-like network branches. We derive the exact stationary distribution of this network ASEP and investigate its transport properties. By considering two representative network geometries for which physical quantities can be expressed in terms of certain hypergeometric series, we elucidate how the network geometry influences transport properties.
Abstract
It is shown that for a class of state integral models on shaped pseudo 3-manifolds, including the edge formulation of Teichmüller TQFT, the Boltzmann weight assigned to a tetrahedron solves the tetrahedron equation. The dihedral angles of the tetrahedron play the role of spectral parameters.
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
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.
Abstract
We investigate the infinite-dimensional limit of nonequilibrium surface growth by numerically integrating stochastic growth equations on a fully connected graph. In particular, we study the Edwards-Wilkinson (EW), Kardar-Parisi-Zhang (KPZ), and tensionless KPZ (TKPZ) equations. Using a network discretization adapted to the all-to-all interaction topology, we analyze the global roughness, height-fluctuation statistics, time power spectra, and two-time correlations. For the EW equation, we obtain an exact expression for the roughness that matches the numerical simulations and shows that the interface becomes flat as . We also compute analytically the time power spectrum, show that height fluctuations are Gaussian, and derive an explicit expression for the two-time height autocorrelation function, indicating that the aging behavior is trivial. For the KPZ equation, finite-size and strong-coupling effects can cause deviations from EW behavior at moderate system sizes , often accompanied by numerical instabilities; however, these differences disappear as increases. In the large- limit, KPZ dynamics converges to EW behavior, as the four observables analyzed exhibit identical scaling properties. Overall, our results indicate that on a fully connected graph the KPZ nonlinearity is irrelevant as , leading to EW-like dynamics with asymptotically flat interfaces. These findings are interpreted in the context of the upper critical dimension of the KPZ universality class.
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
For an arbitrary diffusion process with time-homogeneous drift and variance parameters and , let be times the total time spends in the strip .The limit as is the full halfline version of the local time of at zero, and can be thought of as the time spends along the straight line . We prove that is either infinite with probability 1 or distributed as a mixture of an exponential and a unit point mass at zero, and we give formulae for the parameters of this distribution in terms of , , , , and the starting point . The special case ofa Brownian motion is studied in more detail, leading in particular to a full process with continuous sample paths and exponentially distributed marginals. This construction leads to new families of bivariate and multivariate exponential distributions. Truncated versions of such `total relative time' variables are also studied. A relation is pointed out to a second order asymptotics problem in statistical estimation theory, recently investigated in Hjort and Fenstad (1992a, 1992b).
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
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 consider a card guessing game with complete feedback. An ordered deck of cards labeled up to is riffle-shuffled exactly one time. Given a value , the riffle shuffle is assumed to be unbalanced, such that the cut is expected to happen at position . The goal of the game is to maximize the number of correct guesses of the cards: one after another a single card is drawn from the top, and shown to the guesser until no cards remain. We provide a detailed analysis of the optimal guessing strategy and study the distribution of the number of correct guesses.
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 study the large-scale behaviour of a class of driven diffusive systems modelled by a Stochastic Partial Differential Equation, the Stochastic Burgers Equation (SBE) with general nonlinearity, at the critical dimension and in infinite volume. Our main result shows that, under a logarithmically superdiffusive space-time scaling, it is given by the same explicit Gaussian Fixed point obtained in [G. Cannizzaro, Q. Moulard, & F. Toninelli, arxiv.org/abs/2501.00344, 2025] for the quadratic SBE, but with suitably renormalised coefficients, thereby rigorously justifying and partly correcting the classical Physics derivation of the SBE in [H. van Beijeren, R. Kutner, & H. Spohn, Phys. Rev. Lett., 1986] based on Spohn's theory of nonlinear fluctuating hydrodynamics. Besides, ours is the first universality-type result for out-of-equilibrium systems and the first extension of [M. Hairer, J. Quastel, Forum of Mathematics, Pi, Vol. 6, 2018, e3], to the critical dimension and beyond weak coupling. The major challenge in our work is the mild growth condition on the nonlinearity which renders even the well-posedness of the microscopic equation non-trivial. Additional key novelties include the derivation of fine estimates on the non-quadratic part of the generator as well as a new approximation for the resolvent associated to the solution of the quadratic SBE.
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 introduce a nonparametric model for time-evolving, unobserved probability distributions from discrete-time data consisting of unlabelled partitions. The latent process is a two-parameter Poisson-Dirichlet diffusion, and observations arise via exchangeable sampling. Applications include social and genetic data where only aggregate clustering summaries are observed. To address the intractable likelihood, we develop a tractable inferential framework that avoids label enumeration and direct simulation of the latent state. We exploit a duality between the diffusion and a pure-death process on partitions, together with coagulation operators that encode the effect of new data. These yield closed-form, recursive updates for forward and backward inference. We compute exact posterior distributions of the latent state at arbitrary times and predictive distributions of future or interpolated partitions. This enables online and offline inference and forecasting with full uncertainty quantification, bypassing MCMC and sequential Monte Carlo. Compared to particle filtering, our method achieves higher accuracy, lower variance, and substantial computational gains. We illustrate the methodology with synthetic experiments and a social network application, recovering interpretable patterns in time-varying heterozygosity.
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
We study the asymptotic behavior of random dimer coverings of the one-periodic Aztec diamond in random environment. We investigate quenched limit theorems for the height function and we extend annealed limit theorems that were recently studied in [arXiv:2507.08560]. We consider more general choices of random edge weights (independence is not assumed) and we distinguish two cases where the random edge weights satisfy the Central Limit Theorem (CLT) under different scalings. For both cases, we prove convergence to the Gaussian Free Field for the quenched fluctuations. For the annealed version, it had been shown in [arXiv:2507.08560], that Gaussian Free Field fluctuations can be dominated by the much larger fluctuations of the random environment. To access quenched fluctuations we analyze the Schur process with random parameters in a way that allows to prove the annealed CLT for the height function for non i.i.d. weights. We consider specific examples where we determine the asymptotic fluctuations.
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.