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)
August 2026
Abstract
We study a general class of gradient interface models with Hamiltonian , , assuming essentially that the potential is even, on for some , and . We establish a Helffer-Sjöstrand representation for these models, and use it to prove that their scaling limits are Gaussian Free Fields (GFFs), and that their covariances decay at the same rate as the GFF. This extends results for strictly convex potentials to a large class of non-convex potentials and to arbitrary temperatures. Additionally, we prove Brascamp-Lieb and dimension-free Poincaré inequalities for the models. We obtain these results by representing the interface as a mixture of gradient interface models with strictly convex potentials, extending an idea by Biskup-Spohn who had considered mixtures of Gaussians at moderate inverse temperature . The construction of such a representation is one of the key new contributions of this work.
Abstract
This work provides explicit expressions, in all sectors of the Fock space: soliton, antisoliton and breather, for the fundamental building blocks -the form factors- of the generalised integral kernels which represent the exponential of the field operator in the 1+1 dimensional Sine-Gordon quantum field theory. The obtained form factors reproduce the various results obtained earlier in the literature using, in particular, the free field approach.
Abstract
We prove a hitting-time analogue of cutoff for the star transposition shuffle on the symmetric group S_n. Let tau be the first time at which every non-top card has been selected. We show that the shuffle is asymptotically mixed at time tau: more precisely, the total variation distance between the law of Y_tau and the uniform distribution on S_n is at most exp(-(log n)^(1/2+o(1))). Our proof compares the star transposition shuffle with the random transposition shuffle using simultaneous diagonalization of the two transition kernels, and then adapts the hitting-time strategy of Jain and Sawhney. This introduces a technique that can be applied to card shuffles that are not necessarily conjugacy invariant.
Abstract
We estimate the fractional moments of the normalized mass assigned by the Critical 2D Stochastic Heat Flow to small balls. Our results also cover the discrete case corresponding to the 2D directed polymer model and provide estimates that are uniform in all parameters. One key takeaway of our results is that the vanishing of the fractional moments is completely governed by the divergence of the second moment. We use a quite robust method, by refining the change of measure argument and introducing a novel coarse-graining procedure, reducing the proof to essentially second moment estimates (in fact, we also provide sharp second moment estimates for directed polymers, of independent interest).
Abstract
We study Toeplitz determinants for whose Fourier coefficients satisfy . This regime extends beyond and includes symbols with Fisher-Hartwig singularities. We develop an operator-theoretic approach based on the Baker-Campbell-Hausdorff formula that separates the quadratic term \[ \sum_{k=1}^{\infty}\min(k,n)f_kf_{-k} \] from the higher-order terms in the expansion of . We show that this quadratic term accounts for the possible growth with , while every fixed higher-order coefficient remains bounded. For symbols with bounded positive and negative Fourier parts, our estimates yield two-sided bounds for the determinant after removal of the quadratic contribution. For a broader admissible class, including Fisher-Hartwig-type symbols, we obtain uniform higher-order coefficient bounds and a central limit theorem for the associated CUE linear statistics. We also obtain bounds on mixed exponential moments for CUE-derived random fields beyond the characteristic polynomial.
Abstract
While the solution to the dimensional stochastic heat equation with multiplicative noise is closely related to the exponential of a Brownian motion, the two-dimensional picture exhibits an additional weak-to-strong disorder transition. In [CSZ '23], the critical two-dimensional stochastic heat flow (SHF) was constructed as the scaling limit of the partition function of dimensional directed polymers under the logarithmic intermediate-disorder scaling at criticality. The SHF is a random measure and, like many naturally occurring random measures, it is expected to exhibit rich intermittency. [CSZ '25] established that it is almost surely singular with respect to the Lebesgue measure. More recently, [GT '26] showed that the logarithm of the SHF averaged over small balls is asymptotically Gaussian, with both its mean and variance diverging as the ball radius tends to zero. In this paper we prove a sharp result quantifying the singularity of the support of the SHF as well as its intermittency. In particular, we show that, almost surely, for all small , up to a vanishing error, all the mass of the point-to-plane SHF in any domain is concentrated on balls of radius , each containing mass, thus precisely establishing its logarithmic fractal behavior. A key ingredient in the proof is a refined large-deviations theory, which allows access to conditional distributions, by taking advantage of the Gaussian-like behavior of the SHF at quasi-critical scales. A further useful observation that features prominently is that conditioning a Brownian motion on its endpoint being unusually large essentially induces a shift in the mean of its increments, and consequently, at small enough scales, their distributions do not alter significantly.
Abstract
We prove the conjecture of Huh, Matherne, Mészáros, and St.~Dizier that the normalization of every skew Schur polynomial in finitely many variables is Lorentzian. We first realize every nonzero skew Schur polynomial in finitely many variables as a specialization of a Schubert polynomial and prove that it is dually Lorentzian. The dual Jacobi--Trudi identity then identifies its normalization with the finite dual of a skew Schur polynomial obtained by rectangular complementation. As a consequence, skew Kostka numbers satisfy log-concavity inequalities along the root directions.
Abstract
We consider a class of two-dimensional continuous-time directed random polymers in a correlated Brownian environment. The partition functions of these polymers form a discrete approximation of the Critical Stochastic Heat Flow (SHF), constructed by Caravenna, Sun and Zygouras (2023). Using a recent moment based axiomatic characterization of the SHF by Tsai (2024), we prove that the point-to-point partition functions, viewed as random measures, converge to the critical SHF. The proof proceeds by first establishing convergence of all moments through the resolvent method introduced by Gu, Quastel and Tsai (2021) and based on Rajeev (1999) and Dimock and Rajeev (2004). This is followed by establishing tightness and verifying the remaining axioms in the characterization.
Abstract
We consider the problem of identifying the paths most likely to occur as geodesics in last passage percolation. Heuristics suggest that these modal paths should be the extreme corner paths, going straight between the corners of the cube of accessible vertices. We identify three mechanisms which favour such corner paths, and show that they are more likely to appear than all but a vanishingly small portion of paths. We make more specific comparisons in exponential last passage percolation, where moderate deviation estimates may be used to show that corner paths are nearly modal in a precise sense. Finally, we show a form of monotonicity of geodesic probabilities in a special case and conjecture that this holds in general.
Abstract
Following recent developments on matrix-valued orthogonal polynomials (MVOPs), we construct the matrix Szegő factorisation of the matrix weight function when it is supported on multiple intervals. We find that on the gaps the matrix function has unitary multiplicative jumps. In the next part, performing the Deift-Zhou steepest descent analysis we study the associated global parametrix which comes from the Riemann-Hilbert problem of the MVOPs. A solution is constructed visualising the rows of the parametrix as sections of vector bundles on a hyper-elliptic Riemann surface. An alternate view point to the global parametrix is also presented using a vanishing lemma. As an application we obtain strong asymptotics of MVOPs for multiple cuts which extends the works of Deaño, Kuijlaars, and Román.
Abstract
We prove, in the case , the extension to the -symmetric world of the original Macdonald positivity conjecture. This is achieved by giving a combinatorial interpretation of the Kostka coefficients in terms of standard fillings of the diagram associated to the -partition . This interpretation generalizes the one in the usual Macdonald case, which is given by a major index statistic on standard tableaux.
Abstract
We study rooted Gibbs-DLR measures in the directed polymer model on with an ergodic disorder distribution which satisfies an additional mild hypothesis. We prove that the set of extremal rooted Gibbs-DLR measures is closed and totally ordered, characterize extremality in terms of path coalescence, and show that each fully supported extremal rooted Gibbs measure canonically generates a globally consistent and coalescing family of extremal rooted Gibbs measures indexed by all lattice sites. These families are, moreover, totally ordered. Building on this structure, we prove strong existence and strong uniqueness of the associated Busemann process, together with an continuity theorem for the shift-covariant Busemann cocycles which is joint in the inverse temperature, the tilt parameter, and the random environment. This yields, as a corollary, in-probability continuity of the generated extremal Gibbs measures corresponding to directions of differentiability under bounded i.i.d. perturbations of the weights. In positive temperature, it shows in-probability convergence of the generated extremal Gibbs measures. At zero temperature, it also yields quenched subsequential large deviation principles for the corresponding positive-temperature rooted Gibbs-DLR measures on path space. The rate functions are determined by a zero-temperature Busemann cocycle and vanish precisely on the infinite geodesics generated by that cocycle.
Abstract
High-dimensional Principal Component Analysis (PCA) has become an essential tool in modern data analysis, offering dimensionality reduction and feature extraction. However, the presence of missing data introduces significant challenges, distorting the performance of PCA and complicating statistical inference. In this paper, we study the asymptotic behavior of PCA under a spiked population model with missing observations, leveraging recent advances in random matrix theory. We demonstrate that while the spiked sample eigenvalues exhibit asymptotic normality, the limiting parameters differ substantially from those in the complete data case, reflecting the non?trivial influence of the missing data mechanism. As an application of our results, we propose a test to evaluate the independent structure of a spiked population.
Abstract
We study the nested permutation classes . We relate the corresponding generating functions , to -difference equations. Solving these equations allows to show that the generation functions are not D-finite, and that , .
Abstract
In this paper, we investigate the monic orthogonal polynomials and the Hankel determinants associated with the generalized Freud weight \[w(x;T_{m};λ) = |x|^{2λ+1}\exp\biggl(-\sum_{k=1}^m t_k x^{2k}\biggr),\quad m \in \mathbb{Z}^+,\; t_{k} \in \mathbb{R} , x\in\mathbb{R}\setminus\{0\},\] where \(T_{m}=\{t_{1},t_{2},\cdots, t_{m}\}\), and \(λ>-1\).By employing ladder operators and compatibility conditions, we find that all members of the discrete Painlevé I hierarchy have a unified structure and the recurrence coefficient \(β_n\) of satisfies the -th member of the discrete Painlevé I hierarchy. Besides, we derive the second-order differential equation satisfied by , the partial derivatives of the recurrence coefficients \(β_n\) with respect to parameters \(t_1, t_2, \dots, t_{m-1}\) and the corresponding differential identities for . Based on the discrete Painlevé I hierarchy and the above differential identities, we obtain new partial differential equations satisfied by and .Using the discrete Painlevé I hierarchy and the asymptotic theory of linear difference equations, we derive the full asymptotic expansions of the recurrence coefficient , the nontrivial leading coefficient , and the Hankel determinant as , for general and . Notably, while the logarithmic term appears in the leading-order contributions, it is absent from the remainder terms in these expansions.We illustrate our results under the specific decic Freud weight .
Abstract
The coefficients that determine the Hilbert series of the superspace coinvariant ring are indexed by hook-shaped partitions. We give a manifestly positive combinatorial interpretation of these coefficients, together with several generating functions for them. Specializing this Hilbert series at , we show that its coefficients are differences of binomial coefficients. Consequently, this proves a conjecture of Sagan--Swanson (2024) that these coefficients are palindromic up to sign. More generally, for every we obtain closed-form expressions for the specialization.
Abstract
One hallmark of exactly solvable KPZ random growth models is product-form invariant measures. In the setting of exponential last-passage percolation (LPP), this corresponds to the independence of Busemann increments along any down-right path. However, this independence breaks down when multiple asymptotic directions are considered simultaneously, owing to the fact that jointly invariant measures are not jointly product-form. This paper shows that the failure of independence is one-sided: Busemann increments across arbitrary directions are negatively associated. As an application, we derive an exponential concentration inequality for sums of Busemann increments on the diffusive scale, even when the increments are not independent. While our argument relies on a Burke property that is special to exponential weights, all other proof ingredients$\unicode{x2014}$including hidden LPP monotonicities and braid relations for queueing maps$\unicode{x2014}$hold for arbitrary weights.
Abstract
We prove universality results for zeros of random -adic polynomials with independent, sufficiently non-concentrated coefficients. We show that when the degree goes to infinity, the limiting joint root statistics in arbitrary finite extensions of are universal and coincide with those of the Haar coefficient model, whose limiting root statistics were determined by Caruso (arXiv:2110.03942). In particular, our results extend the Haar-model formulas to a broad class of coefficient distributions and to general joint root statistics over finite extensions of . Our approach is based on a new method, which we call the resultant distribution method. Instead of studying the roots of a random polynomial directly, we analyze the distributions of the resultant valuations against suitable fixed test polynomials. We prove that these resultant distributions determine the limiting law of the polynomial and, together with suitable degree estimates, yield convergence of general root statistics. As a further application, we apply the same method to characteristic polynomials of random -adic matrices with independent entries. We prove that the distinguished factors of the characteristic polynomials converge to universal limiting distributions independent of the entry distribution, thereby showing that such a random matrix model is remarkably robust.
Abstract
We provide limit theory for the trace of the squared sample correlation matrix , constructed from observations of a -dimensional random vector with iid components. If the entries have finite fourth moment and and grow proportionally, it is known that satisfies a central limit theorem (CLT) and the centering and scaling sequences are universal in the sense that they do not depend on the entry distribution. Under symmetry and regular variation assumption with index and any growth rate of the dimension, we prove that the universal CLT remains valid for . For , we identify a critical dimension growth at which the fluctuations of become non-Gaussian. Moreover, if the dimension grows faster and we establish a non-universal CLT with norming sequences depending on the value of . Our findings are illustrated in a simulation study.
Abstract
We study strong asymptotics of orthogonal polynomials on an interval with varying weights, extending the framework of Totik's theorem on weighted orthogonal polynomials. Our results generalize this theorem in several directions. In particular, we allow the supports of the associated equilibrium measures to converge to a proper subinterval of the original interval of orthogonality or even collapse to a point. These extensions are motivated by applications to multiple orthogonality where such varying measures arise naturally.
Abstract
A finite density of mutually avoiding directed polymers in a quenched random medium is a minimal model of glassy line matter. The dilute theory, solved by replica Bethe ansatz, predicts an interaction free energy proportional to and disorder cumulants with distinct power-law dependences on the density , but direct numerical tests have been hindered by the combinatorially large many-polymer transfer matrix. We recast the problem as filling logarithmic eigenvalues of a single-polymer transfer-matrix product, obtaining the quenched free energy, its cumulants, and a disorder-induced linear spectral edge consistent with the replica prediction.
Abstract
We study partitions which attain the largest Littlewood-Richardson coefficient. More precisely, we prove that the largest is attained at partitions such that and is a disjoint union of squares. We conjecture that for , all largest must satisfy this property. We confirm this conjecture numerically, for .
Abstract
This paper provides formulas for Koornwinder polynomials in analogy with the creation formula, the alcove walk formula and the non-attacking fillings formula for the type Macdonald polynomials. We state the creation formula in terms of the divided-difference operators used in Schubert calculus, and we use a box-greedy reduced word to reformulate the alcove walk formula in terms of uncompressed set-valued tableaux. Then two types of compression, ``around-the-end compression'' and ``across-the--gap compression'', are used to derive a formula for Koornwinder polynomials in terms of compressed set-valued tableaux. Throughout we work in the full generality of relative Koornwinder polynomials, which are the analogues of the permuted basement Macdonald polynomials used in the type case.
Abstract
We consider the largest eigenvalues of minors of the classical Gaussian random matrices with a finite rank spike in the critical BBP regime. We show that they converge to a multi-level interlacing particle system whose lowest level is the Airy point process. Our construction and scaling limits extend in a natural way to non-classical .
Abstract
For ordinary permutations on letters, the distribution of the number of fixed points of random permutations is well known to approach the Poisson distribution in total variation distance as super-exponentially quickly. We use Stein's method to get related results for the number of fixed points of random permutations of multisets. Given a sequence of multisets on letters whose expected number of fixed points converges to a constant , we must have and the distribution of number of fixed points converges to the distribution as . If , then the rate of convergence in total variation distance can be as slow as .
July 2026
Abstract
We introduce and solve the Interacting Quantum Symmetric Exclusion Process (IQSEP), a family of models describing the stochastic quantum hopping of charged particles along the edges of a lattice, with hopping amplitudes that depend on the occupations of neighbouring sites. In the absence of interactions, they reduce to the standard quantum simple symmetric exclusion process, exhibiting coherent diffusive transport. For interactions of order one, they capture incoherent diffusive transport and its fluctuations, characterized by density-dependent diffusivity and mobility, making contact with the macroscopic fluctuation theory. By rescaling the interaction strength appropriately with the lattice mesh, we define a mesoscopic scaling regime that retains a finite coherence length in the continuous thermodynamic limit. This regime interpolates between coherent behavior at small length scales and incoherent behavior at large scales. The resulting scaling theory accounts for fluctuations of quantum coherences in interacting diffusive systems, going beyond the scope of standard fluctuating hydrodynamics.
Abstract
Kardar--Parisi--Zhang (KPZ) universality provides an example of macroscopic scaling generated by microscopic violation of detailed balance. While one- and two-dimensional realizations have been explored in driven condensates and growing interfaces, demonstrating KPZ scaling in three spatial dimensions remains a major challenge. Here we propose a three-dimensional exciton-polariton crystal as a platform for observation of 3D KPZ universality. Starting from a stochastic driven-dissipative Gross-Pitaevskii equation for a condensate formed in a three-dimensional photonic-crystal lower-polariton band, we eliminate the massive density and reservoir modes and obtain an effective -dimensional KPZ equation for the condensate phase. Numerical simulations of both the KPZ equation and the full driven-dissipative polariton model show an intermediate-asymptotic regime in which the first-order coherence obeys $-\ln |\gone(0,\Delta t)|\propto |\Delta t|^{2\beta}$ and $-\ln |\gone(\Delta r,0)|\propto |\Delta r|^{2\chi}$, with exponents consistent with the KPZ benchmarks , . Our results identify three-dimensional polariton crystals as a controllable quantum fluid route to higher-dimensional nonequilibrium universality.
Abstract
We introduce and study the perturbed -corners process, a deformation of the classical -corners process. The latter is a probability measure on interlacing arrays of real numbers that, for the classical values , describes the joint distribution of eigenvalues of principal submatrices of a Gaussian random matrix with the corresponding GOE/GUE/GSE symmetry. For classical , the perturbed process arises from adding a deterministic diagonal matrix to a GOE/GUE/GSE matrix. The eigenvalues of the resulting random matrix depend symmetrically on , but this symmetry does not extend to the whole corners process. We extend the construction to all via multivariate Bessel functions, and analyze the crystallization () of the resulting interlacing array, proving a law of large numbers and a central limit theorem in two regimes. For fixed perturbation, the eigenvalue repulsion dominates, and the array crystallizes on the same lattice of polynomial-derivative roots, with the same discrete Gaussian free field fluctuations, as in the unperturbed case treated by Gorin and Marcus (arXiv:1706.07393). New phenomena appear in the second regime, where the 's grow linearly in . Then the external source competes with the repulsion at leading order. Here the array freezes on a deformed lattice characterized by a coupled system of optimality equations. The fluctuations are governed by the same discrete Gaussian free field, now attached to the deformed lattice. The deformed lattice equations decouple in the special case of a single spike in the last coordinate. Then the deformed lattice is obtained explicitly by applying one shifted derivative followed by iterated ordinary derivatives.
Abstract
The Strong Szegő Limit Theorem is a theorem about the asymptotics of the determinants of large Toeplitz matrices. It can be reformulated as a probabilistic statement about eigenvalue statistics of random unitary matrices. We prove a multivariate generalization of the theorem in this latter form, replacing a single unitary with a system of independent random unitaries. It turns out that the standard proofs of the classical theorem do not generalize to this setting, and instead we must import tools from random matrix theory, free probability, and noncommutative function theory. In addition, we obtain results about the averaged determinants of random unitary pencils, as the size of the unitaries tends to infinity; and an auxiliary result giving a formula for the spectral radius of a matricial sum of free Haar unitaries.
Abstract
We give explicit formulas for the smallest non-trivial irreducible representation of quantum affine superalgebras in types D, (in the all-fermionic parity), F, (in the distinguished parity), and G, (in the distinguished parity), both in the Drinfeld-Jimbo and in the new Drinfeld realizations. We use this information to obtain an explicit expression for the corresponding -matrices.
Abstract
We construct a finite-energy Wigner counterexample and a concrete one-parameter family of Gaussian-regularised Bernoulli entry laws yielding counterexamples to the conjectural dimensional monotonicity of the quadratically penalised logarithmic energy for mean empirical spectral distributions by Chafaï, Dadoun, and Youssef.
Abstract
The Plancherel-random infinite Young tableau arises from applying the infinite Robinson-Schensted-Knuth (RSK) correspondence to a sequence of i.i.d. random variables distributed uniformly on the unit interval [0,1]. Building on the isomorphism of dynamical systems established in prior work, we provide an explicit geometric characterization of the inverse map. Our approach makes use of a previously unexplored structure, the jeu de taquin forest on the infinite tableau, where edges connect boxes according to local comparison rules. We prove that each tree in this forest almost surely extends toward infinity with a well-defined asymptotic direction, establishing a canonical bijection between trees and values in the i.i.d. input sequence. The ordering is recovered through tree lifetimes under iterated jeu de taquin transformations.
Abstract
The totally asymmetric simple exclusion process (TASEP) with open boundaries is a finite Markov chain describing particles hopping between adjacent sites on a one-dimensional lattice with left and right boundary transitions governed by parameters and . The multispecies TASEP is a higher-rank generalization in which particles have different species. Multiline queues were introduced by Ferrari and Martin (2007) to compute the stationary distribution of the multispecies TASEP on a circle. It has remained an open problem to find a combinatorial formula for the stationary distribution of the multispecies open-boundary TASEP. Using Kirillov--Reshetikhin crystals of type , we construct type multiline queues and a corresponding Ferrari--Martin pairing algorithm that projects them to TASEP configurations. This yields a combinatorial formula for the stationary distribution of the multispecies open-boundary TASEP for the specialization.
Abstract
We study the non-equilibrium density fluctuations of a one-dimensional two-species symmetric simple exclusion process in contact with critical slow boundary reservoirs. The boundary rates are of order , which leads to Robin boundary conditions at the macroscopic level. We prove that the coupled fluctuation field associated with the empirical densities of two species converges to a generalized Ornstein-Uhlenbeck process. The limiting process is characterized by a linear martingale problem whose coefficients reflect the combined effects of bulk diffusion, species conversion, and the boundary reservoirs.
Abstract
In this paper, by the method of comparing coefficients, we establish a new generating function of multivariable Al-Salam-Carlitz polynomials which contains both Rogers' and Bowman's symmetirc expansion formulas and the classical -Mehler formula as special cases. Some new -series identities related to the multivariable Al-Salam-Carlitz polynomials are also presented.
Abstract
The symmetric exclusion process (SEP) is a classical model of interacting particles on a graph, in which multiple particles execute random walks subject to the constraint that no two particles may occupy the same vertex. The version of the process in which the number of particles is fixed, is by now very well understood, through the study of reversible Markov chains. In this paper, we study the non-equilibrium version of the SEP, in which some vertices of the graph interact with external heat baths, held at different temperatures, that mediate transport through the graph. The resulting steady-state distribution is still the stationary distribution of a natural Markov chain, but the chain is no longer reversible. This means that even the steady-state distribution is very hard to describe, and indeed is the subject of a rich literature, which has focused almost exclusively on graphs which are finite subsets of the lattice . Very little is known about other important quantities, notably the mixing time. We derive various quantitative results about the non-equilibrium SEP for arbitrary graphs with arbitrary heat baths. Our main result is a bound on the mixing time in terms of the worst-case hitting time of a single particle to any of the heat baths; this bound is tight up to two logarithmic factors in the size of the graph. We then use this result to derive a simple algorithm, running in time roughly , that computes the steady-state joint occupation distribution for any set of vertices; since the steady-state distribution is very elusive, this provides a potentially useful tool to study its -wise correlations for small values of . Finally, we prove that the spectrum of the Markov chain associated with the non-reversible dynamics for the SEP is independent of the heat bath temperatures; thus in particular the spectrum is always real.
Abstract
Describing emergent macroscopic phenomena in low spatial dimensions is known to be notoriously challenging, primarily due to strong interactions that render perturbative approaches inapplicable. On the other hand, low-dimensional systems host a wealth of unorthodox phenomena. A prominent example is the emergence of superdiffusive transport in one-dimensional interacting systems, traditionally studied in the framework of nonlinear fluctuating hydrodynamics. After identifying and discussing internal inconsistencies in the previous formulations, in this work we develop a general and systematic approach for constructing effective field theories of one-dimensional hydrodynamic systems in the form of coupled stochastic Langevin-type equations compatible with the physical requirements of local equilibrium states such as the thermodynamic Maxwell relation and fluctuation-dissipation symmetry. We implemented a general numerical integration scheme and exemplified our construction on a simple model of two interacting hydrodynamic modes with a non-Gaussian stationary equilibrium measure.
Abstract
We develop criteria for collision and non-collision of reversible infinitely many interacting diffusion processes in the real line. The approach is potential-theoretic and is based on capacity estimates for symmetric Dirichlet forms on the configuration space. Our main results are model-independent in the sense that no determinantal or Pfaffian structure, prescribed interaction potential, or explicit labelled stochastic differential equation is required. The non-collision criterion involves only the second and third correlation functions of the reversible measure, whereas the collision criterion is based on a local lower bound for a finite-volume conditional density of the reversible measure. As an application, we identify the sharp collision threshold for the diffusion associated with the -symmetric Dirichlet form: the collision set is polar if and only if . This provides an infinite-particle counterpart of the classical collision threshold for the finite-particle Dyson Brownian motion with inverse temperature .
Abstract
KPZ-type extremal fluctuations have recently been proved for several models of random walks in space-time random environments (RWRE) in dimensions. A general moment criterion predicts the spatial scale at which this behavior should occur, but does not by itself guarantee non-trivial fluctuations at that scale. In this paper, we show that the averaging process provides an instance in which this criterion is not sharp: because of a degeneracy in the update mechanism, the actual KPZ limit appears only beyond the predicted scale. The relevant critical contribution arises from the interplay of two distinct fluctuation mechanisms, a phenomenon that appears to be rather special within the RWRE setting. Our proof builds on Dobrushin-type local-time limit theorems for zero-sum additive functionals, refined estimates for tilted -point motions, and the recent moment-based axiomatic characterization of Cole--Hopf solutions to the one-dimensional KPZ equation. We also identify the behavior on the two sides of the critical scale, thereby sharply separating the subcritical, critical, and supercritical regimes of the model.
Abstract
The atypical behaviour of random walks in time-dependent random environment was recently related to Kardar-Parisi-Zhang (KPZ) growth. While this is now well-understood in spatial dimension , further efforts are necessary to better understand these connections in dimensions . In this paper, we study this problem numerically for and , focusing on a discrete model with Dirichlet distributed transition probabilities. This model is a generalization of an integrable model in , and it has the advantage of admitting an explicit, product-form, stationary measure. We verify that the growth of the variance of the logarithm of point-to-point probabilities, namely from the origin to position in time , is compatible with KPZ growth in dimension and . In spatial dimension , we confirm the existence of a phase transition as the angle increases and we obtain a lower bound based on an exact second moment calculation. We find that in the weak disorder phase the point-to-point probability acquires a heavy tailed distribution, and that in the strong disorder phase the cumulants of its logarithm grow with time. Further, we show that for this model, we can compute exactly the sample to sample variance of the thermal average and that it is related to the extreme diffusion coefficient introduced recently.
Abstract
The KP equation is a prototypical -dimensional integrable PDE. Its soliton solutions are famously parametrized by the Sato Grassmannian. In seminal work, Kodama and Williams made the surprising discovery that the combinatorics of soliton solutions are intimately related to the combinatorics of the totally positive Grassmannian as pioneered by Postnikov. They introduced novel algorithmic methods inspired by polyhedral structures arising from tropical geometry. Soliton solutions to the 2D Toda lattice and the Davey--Stewartson equation, two closely related integrable systems with soliton solutions, are also classified by the Sato Grassmannian. Kodama suggested that the methods of his work with Williams could generalize to these two integrable equations. In this work, we show that this is indeed the case. We derive algorithms to produce contour plots from elements in the totally nonnegative Grassmannian in both cases. In the asymptotic setting, we recover and refine previous work of Biondini and Wang; as well as Biondini, Kireyev and Maruno.
Abstract
We study two classes of oscillatory Jacobi-type weight functions , . Other restrictions are imposed on and to guarantee the existence of the associated orthogonal polynomials. By using the ladder operators established in the recent literature for monic orthogonal polynomials associated with Jacobi-type weight functions and three compatibility conditions, we derive two coupled difference equations satisfied by the three-term recurrence coefficients. Compared with the existing results, these equations are structurally simpler and of lower order. Once the initial values are determined, the recurrence coefficients can be computed at any stage through the difference equations. The obtained expressions enable us to conjecture the symbolic forms for the recurrence coefficients.
Abstract
Determinantal point processes (DPPs) are widely used as probabilistic models for diverse random subsets, but their approximation error under model misspecification has not been fully characterized. We study the population-level approximation of a strictly positive target distribution p* by DPPs under the forward Kullback-Leibler divergence. Using information-geometric analysis and the standard quality-diversity decomposition of an L-ensemble kernel, in which the diagonal quality component Q encodes item-specific weights and the diversity component D controls repulsive interactions among items, we show that the quality component can be chosen uniquely to match all first-order inclusion probabilities of p*. The DPP approximation problem therefore reduces to the optimization of the diversity component. This reduction yields a global optimality result for attractively dependent targets: under conditions including weak positive association, the independent product distribution with the same first-order marginals, corresponding to D=I, is an optimal DPP approximation. For more general target distributions, positively correlated pairs yield lower bounds on the approximation error. On the repulsive side, a matching of disjoint negatively correlated pairs yields an upper bound on the approximation error, or equivalently a guaranteed improvement over the independent approximation. We further study local optimality around D=I and, more generally, around block-diagonal diversity matrices by analyzing perturbations between their blocks.
Abstract
Disordered systems such as spin glasses and polymers characteristically exhibit random energy landscapes with many macroscopically separated energetic valleys corresponding to near-ground states. This high complexity renders these systems extremely sensitive to perturbations of external parameters. For instance, the support of associated Gibbs measures may change macroscopically under such perturbations, a phenomenon known as chaos in the literature. In experiments, chaotic phenomena are typically studied via temperature perturbations. In this article, we initiate the rigorous study of temperature-chaotic properties of the continuum directed random polymer (CDRP), a canonical model in the KPZ universality class. The CDRP is driven by white noise and is parametrized by inverse temperature , and is known [Wu '26, Das-Zhu '24] to converge in the zero-temperature limit to the directed landscape constructed in [Dauvergne-Ortmann-Virág '22], the putative universal scaling limit of models in the KPZ universality class. The main result of this article considers the CDRP free energies coupled through the same white noise at a pair of inverse temperatures , and shows that they decouple in the limit , converging to a pair of independent directed landscapes. This is the first such "energetic de-correlation across temperatures" result. Our key estimate measures the "pivotality" or "influence" of spatially thin strips in models of last passage percolation. As a byproduct, the proof strategy also allows to show that the directed landscape is a two-dimensional black noise (in the sense of [Tsirelson-Vershik '98]), previously conjectured by Virág. This provides the third known example of a two-dimensional black noise after critical planar percolation [Schramm-Smirnov '11] and the Brownian web [Ellis-Feldheim '16].
Abstract
We construct the periodic directed landscape, which is the conjectured scaling limit for periodic models in the Kardar-Parisi-Zhang universality class. We establish the convergence of periodic exponential last passage percolation to the periodic directed landscape. Moreover, we confirm conjectures by Baik, Liao and Liu on the local structure of the periodic KPZ fixed point, and establish the convergence of periodic ASEPs to periodic KPZ fixed points, coupled according to the same periodic directed landscape. Our main tool to construct the periodic directed landscape (and prove convergence to it) is a technique for gluing full-space directed landscapes (and their prelimits), which is of independent interest.
Abstract
We study stationary measures of the log-gamma polymer on a finite diagonal strip and in a half-space. We establish the phase diagram for the stationary measure on the strip. To this end we develop a representation of {the Laplace transform of} this stationary measure by independent random variables. We also present an analytic approach that extends Barraquand's contour integral representation of the Laplace transform. We prove that, as the strip width tends to infinity, the stationary measure of the log-gamma polymer on the strip converges to the stationary measure of the half-space log-gamma polymer. Finally, we derive a contour integral formula for the Laplace transform of the stationary measure of the half-space log-gamma polymer.
Abstract
We give a Graham-positive expansion for the product of two double Schubert polynomials indexed by two inverse Grassmannian permutations. Surprisingly, the nonzero structure constants are double Schubert polynomials in two disjoint sets of equivariant variables. We also give a positive expansion for the product of two single Schubert polynomials indexed by a -avoiding permutation (e.g., a Grassmannian permutation) and an inverse Grassmannian permutation. Unexpectedly, the nonzero structure constants are Edelman--Greene coefficients.
Abstract
A random sequence sampled from a known continuous distribution is observed with the objective to choose an item with the overall rank one. A rejected item cannot be recalled and is immediately erased from the memory. Under this memory constraint, the choice problem is not amenable to recursive methods of optimal stopping and becomes a global optimisation task. We focus on a heavy-traffic form of the problem with infinitely many choice opportunities, which we state in terms of a planar Poisson process (PPP). Symmetries of the PPP are used to derive basic structural properties of the optimal stopping rule, including the balance at the boundary equation, and two key integral identities. Throughout, we make throrough comparison to the classic full-information counterpart of the problem, revisiting both discrete- and continuous-time models. The optimal value, stopping rule and other characteristics of the problem are determined analytically and approximated numerically with high precision.
Abstract
We study the class of functions counting accepting paths of non-deterministic log-space Turing machines and construct methods to prove containment in . We prove that a large number of classical combinatorial and number theoretic functions belong to this class: classical functions from enumerative combinatorics (multinomial coefficients, Catalan numbers, linear extensions of trees, Stirling numbers, etc), algebraic combinatorics (number of standard Young tableaux, etc), discrete geometry, number theoretic functions, representation theoretic multiplicities in a large class of cases. We show that -plethysm coefficients of bounded length outer partition can be counted by log-space polytime verifiers. We pose numerous questions and conjectures on containment and its generalizations, that suggest venues for conditionally disproving -completeness. While studying which combinatorial functions are in provides a formal way of (dis)proving the existence of combinatorial interpretations, the lower class serves as an analogue for functions computable in polynomial time.
Abstract
We study the vector Dyson equation with parameter in the complex upper half-plane , where and is a nonnegative matrix, not necessarily symmetric. This equation has a unique vector solution , for which we establish a complete measure decomposition and prove regularity. We then develop a graph-theoretic approach to the singularity and stability problem for non-symmetric matrices . The graph structure of identifies the possible degeneracies of the stability operator as approaches the real axis. In particular, for non-backtracking matrices, we prove square-root growth at regular edges, cubic-root growth at regular cusps, and complete stability estimates. We also obtain the corresponding estimates for symmetric matrices in the periodic setting.
Abstract
This article associates to each positive real number a dimer model on a snake graph with activity parameter whose structure is determined by the continued fraction expansion of . When is rational, the model is finite and gives rise to a probability measure on perfect matchings. For irrational , the model is infinite, and is defined as a limit over rational approximations to ; the main technical result of the paper shows that this limit is well defined and independent of the choice of rational approximants. denotes the odds that a -random perfect matching includes a distinguished edge. When is rational, , where is the algebraic -deformation introduced by Morier-Genoud and Ovsienko. This agreement, together with evidence from the irrational case, suggests a close connection between the probabilistic and algebraic constructions.
Abstract
We study energy transport in a finite one-dimensional unpinned harmonic chain with stochastic nearest-neighbor momentum exchanges and Langevin heat baths at its endpoints. Such systems are known to exhibit superdiffusive transport driven by long-wavelength acoustic modes, leading to fractional macroscopic behavior. While fractional heat equations have been rigorously derived for infinite chains, the corresponding boundary conditions for finite systems in contact with heat baths remain unclear due to the nonlocality of the fractional Laplacian. Under the superdiffusive time scaling , where is the system size, we prove that the averaged microscopic energy profile converges, as , to a temperature field solving a fractional heat equation on , with the generator given by a Neumann fractional Laplacian and additional nonlocal boundary terms induced by the heat baths. Our results provide a rigorous derivation of macroscopic boundary conditions for superdiffusive heat transport in open chains and introduce new boundary conditions for fractional Laplacians, that are motivated by a physical model.
Abstract
We give an algebraic branching program of polynomial size which computes Cayley determinant of right quantum matrices. This is a rare example of an efficient computation of a noncommutative determinant, and the first such example for quantum groups. We extend the results to the -Cayley determinant of -right quantum matrices, as well as to their multiparameter generalization. The proofs are entirely combinatorial, as we relate Cayley, Moore and Valiant determinants using bijections/involutions on words. We then employ the celebrated determinant construction of Mahajan and Vinay (SODA'97), to obtain the results.
Abstract
Non-stationary time evolution of interacting particle systems is in general a rather difficult topic however, exceptional examples are known where hidden processes within the model make the description manageable. These hidden processes often take the form of a finite number of interacting random walks and in some cases, rather than being hidden, were very explicitly revealed as second class particles associated with the model. The examples of asymmetric exclusion and exponential bricklayers model are known in this context, where such distributional structure was demonstrated in infinite volume. Here we find boundary mechanisms that save this remarkable structure in finite volume of the model with semi-open boundaries that let ordinary particles, but not the shocks, through. This finding also allows to characterise nontrivial two-species, non-reversible stationary distributions subject to our special boundary rates.
Abstract
The spiked mixture model (SMM) has been introduced as a probabilistic model that generalizes the single-spike (Wishart) model to a mixture model form. With applications ranging from imaging mass spectrometry in the life sciences to hyperspectral imaging in computer vision, it is crucial to understand under which circumstances its signals can be recovered from noisy measurements. The highly multiplexed nature of these measurement types furthermore necessitates such analysis to hold in highdimensional settings. In this paper, we prove that the extreme eigenvalues of the covariance matrix from the SMM exhibit a phase transition in high-dimensional regimes. We show that this phase transition, and thus signal recovery by extreme eigenvalues, depends on several interacting factors: the correlation between spikes (i.e., how similar in content underlying signals are), the energy parameters (i.e., the absolute strength of each underlying signal), and the mixture probabilities (i.e., how likely it is to encounter each underlying signal). This work provides sharp information-theoretic bounds on the parameters needed to detect one or more spikes from extreme eigenvalues of the SMM covariance matrix, and these guarantees could potentially impact any application of the SMM. Understanding this interplay could serve as a tool for driving experimental design in analytical chemistry and life sciences.
Abstract
We compute the dimension of the center of the Temperley--Lieb algebra over a field of characteristic zero for every nonzero value of the parameter . The proof uses the cellular filtration by cup number, together with known facts about the representation theory of the Temperley--Lieb algebra, especially the structure of its standard modules and their radicals. Dilation and compression maps compare the induced graded pieces of the center at levels and , giving an upper bound of one for each such piece. A deformation argument gives the matching lower bound, and hence . We also prove that every central element is fixed by the canonical anti-automorphism and by the natural diagram-reflection automorphism. Finally, we give a congruence criterion for the trivial-radical case and record a Gram-matrix computation for leading terms.
Abstract
The type D ASEP is an asymmetric two--species interacting particle system on , in which two separately conserved species hop, bind into a composite ``bound pair'', and split. The model, along with its reversible measures and orthogonal polynomial duality, was constructed using the representation theory of . The reversible measures and orthogonal polynomial duality are each a product of two copies of the single-species ASEP reversible measures and orthogonal polynomial duality. In this paper, we study the long-time asymptotics of the type D ASEP. In the fixed-- regime, using an exact current--decoupling identity, we prove that the asymptotic hydrodynamic limit and Tracy--Widom fluctuations decouple, as predicted from the duality. In the weak--asymmetry (Edwards--Wilkinson) regime, when , we prove that the two density fluctuation fields \underline{decouple}: each converges to a linear stochastic heat equation, with no cross--coupling in either the drift or the noise, the limiting noises having vanishing cross--correlation. More surprisingly, we then prove that the two limiting normal random variables are \underline{correlated} with a seemingly new correlation function. The correlation is exactly equal to , with the positive parts of the normal random variables having correlations expressed by the Bessel--Struve function. This paper, with the exception of the abstract and introduction, was written entirely by Claude Opus 4.8 and Fable 5. The proofs were then formalized in Lean, using Aristotle by Harmonic AI. The human author of this paper verified the proofs manually.
Abstract
We study the global fluctuations of the height function associated with the Airy line ensemble. Using its determinantal structure and a steepest-descent analysis of the extended Airy kernel, we prove that, after a suitable rescaling, the centered height function converges to an explicit pullback of the Gaussian free field. The convergence holds in the sense of joint moments of linear statistics against compactly supported continuous test functions.
Abstract
We prove Krattenthaler's conjecture from 2001 about the -phenomenon for lozenge tilings of semiregular hexagons. In a first step we reduce the problem to the case of regular hexagons. In a second step we further reduce the question to a special case already covered in the literature. This is achieved by vast application of the celebrated Zeilberger Algorithm and the Holonomic Ansatz.
Abstract
We provide a definitive analysis of the number of correct guesses in the complete-feedback card guessing game after an asymmetric single-shelf shuffle with parameter . We explicitly describe the optimal strategy that maximizes the expected number of correct guesses. We study the number of correct guesses, under an optimal strategy, using methods from analytic combinatorics. In addition, we find the explicit distribution for the number of correct guesses, and thus find the mean and the variance for the number of correct guesses. We show that the distribution of the number of correct guesses is log-concave. We also study the limiting behaviour of the number of correct guesses as the number of cards goes to infinity. In particular, we prove a (local) central limit theorem and a large deviation principle with an explicit rate function. We also prove phase transitions for the number of correct guesses near and . Prior to this work, only the optimal strategy and the expected number of correct guesses under the optimal strategy were known for the special case of .
Abstract
We study the weakly coupled elliptic Anderson model with spatial white noise on the four-dimensional torus, which provides a basic example of a critical SPDE requiring renormalization at arbitrarily high orders. With coupling where is sufficiently small, we prove that the Green's function of the corresponding random Schrödinger operator, suitably centered and rescaled, converges to a centered Gaussian random field with explicit covariance. The main difficulty is that, for such critical models, one must expand up to order , while the perturbative expansion contains factorially many pairings and a growing number of renormalization terms. To overcome this, we construct a truncated renormalized parametrix and prove sharp high-order bounds for its remainder. A central ingredient is a multiscale analysis based on a new version of Hepp trees, combined with new estimates for summations over permutations. These estimates reveal a precise balance between logarithmic losses from scale summation and factorial gains from the structure of primitive pairings. The methods developed here are intended as a first step toward a general theory for critical SPDEs with weak couplings.
Abstract
We study a -dimensional directed polymer in a spatially correlated random environment generated by power-law tail variables: , where the variables are i.i.d. and have a regularly varying right tail with exponent . The spatial covariance of the environment has long-range decay with Hurst parameter . We identify the limiting fluctuations of the log-partition function in the intermediate disorder regime and show that the critical tail exponent is . When , the model has the same scaling limits as the corresponding Gaussian spatially correlated polymer: if , the centered log-partition function converges to the logarithm of the solution of the stochastic heat equation driven by fractional spatial noise; if , its normalized fluctuation converges to a centered Gaussian law. In the regime , at the scale , the log-partition function still satisfies Gaussian fluctuation. The main ingredient is a truncation comparison argument adapted to long-range moving-average environments, together with an invariance principle for polynomial chaos. Due to the non-locality of the environments, we perform a far-near field analysis, as well as multiscale analysis, to prove that the truncated version does not change the log-partition function at the corresponding scales.
Abstract
Recently, there has been significant interest in applying the method of moments developed by Wood and others to study distributions of finite abelian groups that arise in number theory and combinatorics. When the moments do not grow too fast, they determine a unique distribution. We construct large families of distributions that have the same moments. These families include several distributions that arise naturally in the study of sandpile groups of families of random graphs. Wood determined the distribution of Sylow -subgroups of sandpile groups of Erdős--Rényi random graphs. This was extended by Mészáros to sandpile groups of random -regular graphs, who observed an interesting special case when is even and . We study Sylow -subgroups of sandpile groups of random bipartite graphs and similarly find a special case for . Although this distribution differs from that of Mészáros, we show that they have the same moments and fit into our broader construction. To compute the moments of the distributions we study, we apply combinatorial tools from the theory of Hall--Littlewood functions.
Abstract
We show that the lower edge of the appropriately scaled size Laguerre beta-ensemble with parameter converges to the process as when and . This completes the picture of the possible edge scaling limits of the Laguerre beta-ensemble with a fixed . When our proof establishes operator level convergence of the inverse of the scaled Dumitriu-Edelman tridiagonal matrix to the inverse of the stochastic Airy operator. Our methods allow us to prove similar operator level limits for the known soft edge scaling limits of the Laguerre and Gaussian beta-ensembles. For we give a different argument that relies on coupling and a result of Dumaz-Li-Valko for the transition between the hard and soft edge limits of the Laguerre beta-ensemble.
Abstract
We prove that the universal local point processes of random matrix theory are characterized by their loop equation hierarchies. More precisely, for every rational , the point process is the unique solution of the bulk loop equation hierarchy, and the point process is the unique solution of the edge loop equation hierarchy. These uniqueness results provide a direct route to universality: it suffices to verify the corresponding approximate loop equations for the ensemble. In many models, these equations follow from local laws and integration by parts.
Abstract
We study non-exchangeable interacting diffusions with pairwise interaction strengths encoded by a sequence of matrices. Under suitable structural and denseness conditions on these matrices, we prove a universal Central Limit Theorem for the global fluctuation field. As the number of particles becomes large, it converges in distribution to the unique solution of a stochastic partial differential equation (SPDE), the same Gaussian limit as in the exchangeable mean field case. The result applies, for instance, to scaled adjacency matrices of -regular graphs when . A spatial interaction model shows that the denseness threshold is sharp. The proof proceeds with an analysis in negative Sobolev spaces, building on sharp quantitative propagation of chaos results together with functional inequalities.
Abstract
Given any , let be a random matrix over , with all entries independent and -balanced (lying in each residue class mod with probability at most ). Wood proved that as the distribution of approaches Cohen and Lenstra's conjectured distribution of class groups. Given such that , we prove that the distribution of still approaches the Cohen--Lenstra distribution even if we weaken the hypothesis by allowing up to entries per column and up to entries per row of to not be -balanced. We also weaken the independence condition by allowing certain types of dependence between the entries of each column. In addition, we prove that, for any , the cokernels of random band matrices of width with -balanced entries in the band and arbitrary entries outside of it will also approach the Cohen--Lenstra distribution, which answers a question of Kang--Lee--Yu.
Abstract
We study invariant measures for soliton systems described by Mealy automata. Motivated by recently introduced soliton models associated with 2-letter, 3-state Mealy automata, we formulate the time evolution induced by Mealy automata on bi-infinite configuration spaces. We provide sufficient conditions for the invariance of Bernoulli product measures and derive a criterion for the invariance of two-sided space-homogeneous Markov distributions. We then apply these general results to three soliton models, which can be interpreted as variants of the box-ball system (BBS). For two of these models, BBS-S(2) and BBS-V(2), we prove that Bernoulli product measures are invariant. For the remaining model, BBS-C(2), we establish a more general result: the invariance of two-sided space-homogeneous Markov distributions, which include Bernoulli product measures as a special case. Furthermore, for all three models, we compute the phase shift associated with the interaction of two solitons, as well as the velocity of an isolated soliton. Although the latter has already been studied previously, both quantities constitute fundamental characteristics for understanding the generalized hydrodynamics of these systems. These results provide a foundation for the study of invariant measures, generalized Gibbs ensembles, and generalized hydrodynamic behavior in Mealy-automaton soliton systems.
Abstract
Let be an real-valued random matrix with independent, mean-zero, variance-one entries whose fourth moments are uniformly at most . Suppose that there exists such that the entries of satisfy We prove that there are constants , depending only on and , such that for every fixed invertible matrix and every , In the Gaussian case, we also show that the above estimate is sharp in the sense that
Abstract
We give a new proof of the celebrated theorem of Knutson and Tao that the spectra of triples of Hermitian matrices exactly correspond to positions of boundary rays of honeycombs. Most importantly, our proof gives new insights into why honeycombs are related to Hermitian matrices in the first place. Our proof is axiomatic: We distill four essential properties shared by honeycombs and spectra of Hermitian triples, and show that any two objects sharing these four properties must be equivalent. In this way, we argue that honeycombs are `model organisms' for Hermitian triples: they are families of objects satisfying the same defining properties, but in more obvious ways.
Abstract
Determinantal point processes have recently emerged as a kernel-based alternative to standard independent sampling for constructing efficient minibatches, coresets, and other compact representations of large-scale datasets. In particular, sampling mechanisms based on DPPs are believed to demonstrate better approximation properties compared to classical i.i.d. samplers, even at the scale of the exponent. One of the key strengths of DPP based samplers is that they can be deployed over very general spaces, in contrast to more classical sampling methods beyond i.i.d. which tend to work in very well-structured settings, principally Euclidean spaces. In this work, we establish explicit rate guarantees for determinantal sampling in spaces that extend far beyond known Euclidean setups, focusing on spectral kernels obtained from eigenspaces of naturally associated Laplacian and other Markov diffusion operators. This includes, in particular, Riemannian manifolds and weighted networks. In determinantal sampling from compact Riemannian manifolds, we establish sampling rates that automatically pick up the intrinsic dimensionality of the underlying manifold. In the setting of networks, we investigate DPP-based samplers on the celebrated k-nearest neighbour graphs, as well as weighted random geometric graphs, and demonstrate a similar improved dependence on the intrinsic dimensionality of the data. Overall, our approach achieves guarantees of that match known rates on Euclidean spaces of comparable dimension. In terms of techniques, we connect to the celebrated Weyl's Law for manifold spectra, and leverage tools from the theory of Markov diffusions and Dirichlet forms as well as certain ingredients from the theory of pseudodifferential operators, which could be of independent interest in this area.
Abstract
We prove equality conditions for the Ahlswede--Daykin (AD) inequality and the Fortuin--Kasteleyn--Ginibre (FKG) inequality. We then present a number of applications and special cases of these equality conditions. These include Björner's and Fishburn's inequalities for linear extensions of finite posets, the Lam--Postnikov--Pylyavskyy (LPP) and the Okounkov inequalities for Schur positivity of products of Schur functions. We conclude with equality conditions for the Ahlswede--Daykin--Schur (ADS) inequality recently introduced in Chan--Chen--Pak--Soskin (2026), which is an AD type extension of the LPP inequality.
Abstract
We prove a mesoscopic central limit theorem for linear eigenvalue statistics of correlated Hermitian random matrices. The class considered here includes Wigner and Wigner-type matrices, as well as models whose entry correlations decay polynomially in the distance between index pairs. The proof combines a multivariate cumulant expansion with multi-resolvent local laws and a detailed analysis of the resulting variance kernel on the operator-level.
Abstract
We study a bi-infinite system of interacting Brownian particles on the real line with singular asymmetric interactions mediated by the collision local times. Particles perform Brownian motions, and when neighboring particles collide, the associated local time is split in proportions and . We first develop well-posedness theory for the particle system, proving pathwise uniqueness and strong existence under natural growth assumptions on the initial configuration and local times. We also identify a family of stationary distributions for the infinite-dimensional process of gaps between successive particles: for every , the product measure with i.i.d. Exp gaps is invariant. Our main result concerns the equilibrium fluctuations of the associated particle-count (height) function in a weakly asymmetric regime. Taking so that , initializing the interparticle gaps with i.i.d. Exp random variables, and applying a microscopic Hopf-Cole transform to the diffusively rescaled count function, we prove convergence, as , to the multiplicative stochastic heat equation (SHE) with Brownian exponential initial data. Equivalently, the logarithm of the limit is the Hopf-Cole solution of the KPZ equation with two-sided Brownian initial data. The proof combines localization to finite particle subsystems via chains of collisions, Brownian last-passage percolation estimates, a key local time cancellation property, and a martingale problem for a scale-adapted mollification of the Hopf-Cole field, whose space-time regularity is tuned to match that of the limiting SHE. The resulting fluctuation theorem places these Brownian particle systems with asymmetric singular collision dynamics within the KPZ universality class.
Abstract
The goal of the minimal model program for the Hilbert scheme of points on a surface aims is to describe the (stable) base loci of all divisors, their associated birational models, and the maps between them. We answer all of these questions for the Hilbert scheme of points on the blowup of the affine plane at the origin. The birational models are Brill-Noether loci in a larger Hilbert scheme, and nested variants thereof, and the wall-crossing maps are described as explicit projections. We also establish several new facts about the homogeneous coordinate ring of this Hilbert scheme, including finding the minimal set of line bundles whose sections generate the ring.
Abstract
We prove that the number of lozenge tilings of a certain triangular region is given by the formula \[T_n=\prod_{\substack{1\leq a<b\leq 3n+2\\(a,b)\not=(n+1,2n+2)}}\left|1+ζ^a+ζ^b\right|^{1/3},\] where . This answers a question of Ciucu and Krattenthaler, both by finding the exact formula and by explaining why has many prime factors. The proof reduces the lozenge tiling enumeration problem to evaluating the determinant of the bipartite adjacency matrix of the dual graph of , and then evaluates this determinant by diagonalising .
Abstract
We present an embedded (in ) degeneration of Lusztig varieties (which generalize type Hessenberg varieties) to certain unions of Richardson varieties, giving a simultaneous reproof (and extension) of results of Anderson--Tymoczko, Harada--Horiguchi--Masuda--Park, and Kim. Although torus-equivariant, the degeneration is not Gröbner. In the case that the Lusztig variety is the permutahedral toric variety, this degeneration provides a subdivision of the permutahedron into Bruhat interval polytopes, and we prove a more general result showing equivariant degenerations of projective toric varieties produce subdivisions of the moment polytope (as was shown in the Gröbner case by Sturmfels). A Gröbner degeneration would result in a {\em regular} subdivision, and despite our degeneration not being Gröbner we show in types that our subdivisions of the permutahedron are indeed regular.
Abstract
Bergeron, Garsia, Haiman and Tesler conjectured in 1999 that, for all partitions , the polynomial has nonnegative integer coefficients, where is the Bergeron--Garsia nabla operator, which acts diagonally on the modified Macdonald basis, and is the monomial symmetric function. In this article, we prove this conjecture, and more generally that for all . We establish a recursion showing that has an expansion with coefficients in in the symmetric functions , where denotes the operator introduced by Haglund, Morse and Zabrocki. Combining this expansion with the compositional shuffle theorems of Carlsson--Mellit and Mellit, and with the Schur positivity of LLT polynomials, completes the proof. The same method, using the -positivity of column LLT polynomials after the substitution , also gives an -positive analogue.
June 2026
Abstract
In previous work \cite{GW}, we developed a theory of modulated \(2j-k\) bi-orthogonal polynomial systems \(\{P_n(z;r),Q_n(z;r)\}\) and \(j-2k\) bi-orthogonal polynomial systems \(\{R_n(z;r),S_n(z;r)\}\), which generalize the classical \(j-k\) Toeplitz systems. In the present paper, we further develop this theory in several directions. We derive simplified and unified recurrence relations for both families of polynomials, prove a more transparent Christoffel--Darboux formula, and give Riemann--Hilbert characterizations of the \(2j-k\) and \(j-2k\) systems.
Abstract
We establish the explicit leading order asymptotics, with a quantitative error bound, of tail probabilities of collision times for a class of integrable space-time inhomogeneous Markov chains, in discrete and continuous time. The corresponding process conditioned not to intersect arises in interacting particle systems with local push-block interactions thereby confirming a recent prediction. The generic discrete nature of the spatial inhomogeneities rules out powerful coupling-with-Brownian-motion techniques, so our proof strategy proceeds instead via a novel steepest-descent analysis combined with a Karlin--McGregor semigroup expansion in terms of dominant-index contributions.
Abstract
We study the macroscopic fluctuation theory (MFT) of a multispecies generalization of the symmetric simple exclusion process (mSSEP) on the infinite line, in which particles of species -- including, possibly, a vacancy species -- exchange positions at unit rate. Working with the full, redundant set of coarse-grained densities $\bfrho=\{\rho_0,\dots,\rho_N\}$ keeps the relabelling symmetry of the model manifest throughout. We first extend the argument of Derrida and Gerschenfeld to the multispecies setting, showing that the cumulant generating function of the multispecies current between two regions of an arbitrary graph depends on the boundary densities $ \bfrho_L,\bfrho_R$ and on the fugacities $\bflambda$ only through a single scalar variable . We then formulate the MFT saddle-point equations for the mSSEP on the infinite line and show that they define an integrable system: they are of Landau--Lifshitz type, and a Zakharov--Takhtajan gauge transformation recasts them in AKNS form. Solving the resulting linear scattering problem by the inverse scattering method, we recover the cumulant generating function for the multispecies current, as well as the initial and final density profiles conditioned on a prescribed current fluctuation. In particular, we show that coincides with the function obtained by Derrida and Gerschenfeld for the single-species SSEP, now derived for an arbitrary number of species directly from the integrable structure of the multispecies MFT equations.
Abstract
We prove that the differential of the log of the partition function for the -matrix model with quartic interactions converges in a certain double-scaling regime to the differential of the -function for the string equation. This confirms the convergence of the critical Ising model on random surfaces to the topological minimal model, which was stated in the works of Douglas and Shenker, Brézin and Kazakov, and Gross and Migdal. Our analysis is based on a steepest-descent analysis of a Riemann-Hilbert problem associated to a family of biorthogonal polynomials. New features in the matching problem in the construction of local parametrices appear.
Abstract
Let be a whole plane Gaussian free field, and let be a bounded domain in two dimensions. We study the asymptotics as of the Liouville quantum gravity (LQG) heat trace, defined as the integral over of the on-diagonal LQG heat kernel. Our main result is to show that the second term in the spectral expansion as of the expected heat trace is governed by a nontrivial exponent, given by the KPZ (Knizhnik--Polyakov--Zamolodchikov) relation. A similar but stronger (almost sure) result applies to the related notion of heat content. Along the way we obtain various results on the short-term behaviour of the heat kernel, notably solving a conjecture of \cite{BW} concerning its annealed asymptotics, and showing the finiteness of all moments of the properly rescaled heat kernel.
Abstract
For every , we prove global existence for the two-dimensional generalized parabolic Anderson model on the whole space with nonlinearity , driven by an enhanced noise . Here the noise has polynomially weighted spatial Besov--Hölder regularity , and is the corresponding renormalized second-order object. If is globally Lipschitz, the solution is unique. The proof combines a weight-compatible annular high-low decomposition with a paracontrolled transport representation. The final remainder is estimated simultaneously in weighted norm and in a higher order weighted parabolic Hölder norm, with strictly different polynomial weights. This weight gap absorbs the polynomial losses arising from the enhanced noise and the transport coefficient. Several key refinements of the earlier works allow the maximum-principle and Schauder estimates to yield a global a priori bound for a larger range of . Uniqueness is proved in a time-dependent exponentially weighted topology.
Abstract
The generalized Yang-Baxter equations are multi-site versions of the standard Yang-Baxter equation. When spectral parameters are included, such equations are expected to lead to integrable Hamiltonians with local interactions involving multiple degrees of freedom. In this work we characterize both the continuous and discrete symmetries of these equations required to establish an equivalence class of solutions. We find that the set of such symmetries depend on the number of sites on which the equation is supported. In several cases there are more symmetries than the standard Yang-Baxter equation, thus placing heavy constraints on the number of inequivalent solutions and the associated integrable models.
Abstract
We prove an off-diagonal expansion of the kernel of the Toeplitz operator whose symbol is the indicator function of a compact domain with smooth boundary in a complete symplectic manifold of bounded geometry. Using our approach, we extend two results to the non-compact setting: the first concerns the asymptotics of the trace of polynomials in this operator, and the second establishes a Weyl law for this Toeplitz operator.
Abstract
This paper is motivated by a recent result of Pitman and Yakubovich stating that a partially exchangeable, stationary (PES) random partition of is exchangeable. This echoes an earlier theorem of Kallenberg on the equivalence of spreadability (contractability) and exchangeability for infinite partitions. We revise the hierarchy of symmetries with a focus on partitions of finite sets , and ask about the extent to which these relaxed symmetry types differ from exchangeability, framing the question in terms of the geometry of the polytope of distributions defined by the symmetry constraints. We show that single-orbit exchangeable partitions remain extreme among spreadable and PES distributions, and that is the first case where non-exchangeable extreme partitions occur, causing the associated polytopes to deviate from a simplex structure.
Abstract
We derive a determinantal expression for the inhomogeneous asymmetric five vertex model in a rectangular geometry with arbitrary boundary conditions at the bottom and top. Standard non-intersecting lattice path, or free fermion, approaches are not applicable and the determinantal form thus is not immediate.
Abstract
We study a piecewise constant function with a finite number of discontinuities in any interval. We assume that the associated walk satisfying , pinned by , has finite length excursions over past minima. This is the continuous generalization of an initial ball configuration in the discrete Box Ball System introduced by Takahashi and Satsuma, where solitons of integer sizes are identified. We extend the slot decomposition developed by Ferrari, Nguyen, Rolla and Wang in the discrete setting to the continuous case. Each soliton of is represented by a point in two dimensional space, one coordinate for position and the other for the soliton height, mapping to a point configuration. We consider a distribution on walks given by a product measure on the decomposition of the path into excursions over past minima. Excursions are distributed as products of their solitons weights, which are determined by the soliton heights. We show that when the weight function is in the slot decomposition of is a Poisson process. This extends to the continuous case an approach of Ferrari and Gabrielli. As an example, we compute the intensity measure of the Poisson process associated to the asymmetric telegraph process introduced by Kac. In a forthcoming paper we discuss the dynamic properties.
Abstract
A fundamental identity for the number of vacillating tableaux was originally obtained through the representation theory of partition algebras. We extend this identity to arbitrary differential posets and show that it, together with several analogous identities, follows directly from the structural properties of differential posets. Specializing to Young's lattice and its Cartesian powers, we further obtain new bijective proofs via a simple deletion--insertion process.
Abstract
Motivated by the surge of interest in statistics of non-Hermitian random matrices as a framework for description of universal characteristics of dissipative chaotic quantum many-body systems, we address the problem of characterizing the parametric correlations of spectral densities. Considering parameter-dependent ensemble of complex Ginibre matrices we derive an explicit, closed-form expression for the parametric number covariance in the systems of symmetry class for eigenvalues in a circular domain containing on average a finite number of eigenvalues in the spectral bulk. This behavior is expected to be universal, as further supported by numerical evidence for the real Ginibre ensemble, non-Hermitian Bernoulli Wigner matrices and bi-unitarily invariant ensembles. We also discuss a relation between parametric correlations of spectral densities and the distribution of the so-called eigenvector non-orthogonality factor, which attracted considerable interest in recent years.
Abstract
In this work we consider an Aztec diamond model split into two unequal regions which are asymptotically fixed in size. Each region is weighted with a distinct two-periodic weighting. We refer to this model as the t-split two-periodic Aztec diamond, to signify its difference from the previous work title Split Two-Periodic Aztec Diamond, where the model was split into two equal regions. We derive an integral expression for the correlation kernel of the model and give a partial description of the scaling limit behavior, along with a conjecture for the remainder. We refer to the larger and smaller sides of the model as the dominant and non-dominant sides, and to the location of the weight change as the interface. The dominant side exhibits a limit shape that depends only on its own weighting and is identical to that of the two-periodic Aztec diamond, while the non-dominant side appears to have a novel limit shape that depends on both weightings and the location of the interface. Lastly, we consider the complete limit shape in the case where the dominant side two-periodic parameter goes to 0.
Abstract
We present a new classification of the class of all the hypergeometric orthogonal polynomial sequences. The classification uses properties of the coefficients and in the three-term recurrence relation satisfied by the orthogonal polynomial sequences. Such coefficients are rational functions of and are determined by a set of six parameters. The partial fraction decomposition of is a linear combination of five linearly independent functions of with coefficients , for , that are polynomials in our six parameters. For each value of a parameter , associated with the eigenvalues, we classify the according to the set of coefficients that are nonzero. There are certain particular values of that must be considered separately. Our classification yields a collection of 53 disjoint classes and it is different from the Askey scheme in several aspects. It is similar to the classifications proposed recently by Koornwinder.
Abstract
In an asymmetric shelf shuffle, a deck of cards is dealt sequentially from the bottom and assigned one of the shelves uniformly at random. The card is placed at the top of the assigned shelf with probability , and at the bottom of the assigned shelf with probability . Analysis of the shelf shuffle has gained much attention recently, and the case was first treated by Diaconis--Fulman--Holmes [Ann. Appl. Prob. 23 (2013), no. 4, 1692--1720]. In this paper, we extend the analysis of the shelf shuffle to general . In particular, we study the distribution of cycles, cycle lengths, number of descents, number of valleys, number of inversions, and the RSK shape of a permutation obtained from an asymmetric shelf shuffle. Our results extend the analysis of Diaconis--Fulman--Holmes to arbitrary . Furthermore, our analysis of the distribution of descents and inversions is new even for .
Abstract
A mechanical shuffler consists of shelves. A deck of cards, arranged in increasing order, is dealt from the bottom sequentially. Each card is assigned a shelf uniformly at random and placed on the top (bottom) of the existing pile with probability () independently. We refer to this as asymmetric shelf-shuffle. We find the law of the permutation induced by the asymmetric shelf-shuffle and show that the pair consisting of the number of descents and the number of valleys is a sufficient statistic. This generalizes a result of Diaconis, Fulman, and Holmes (Ann. Appl. Prob., 2013) corresponding to the case . For , Chen and Ottolini (ECP, 2025) established the cutoff in the total variation distance near . We establish the cutoff for the asymmetric shelf shuffle. Let be the uniform measure on the set of all permutations of . For a fixed and , we show that \[\TV\left(ν_{n, \lfloor cn^{3/2}\rfloor }^{(p)}, ν_n\right)=1-2Φ\left(-\frac{|2p-1|}{4\sqrt{3}c}\right)+O_{c, p}(n^{-1/2})\;.\] We also establish the cutoff in the separation distance near and in the relative entropy near . In both cases, we also obtain the cutoff profile explicitly.
Abstract
The Fomin local rules for Schensted column insertion can be seen as a two-lane box-ball system, in which a carrier moves particles forward or laterally. Running such two-lane dynamics in parallel on a periodic lattice gives rise to a two-dimensional generalization of the box-ball system, which we call the \emph{skew column RSK dynamics}. Equivalently, this is a deterministic dynamics on pairs of skew semistandard Young tableaux . We prove that this dynamics exhibits solitonic behavior and construct an explicit bijection that linearizes the time evolution. The resulting coordinates consist of two horizontally weak tableaux recording the asymptotic soliton data, integer riggings , and a weakly decreasing sequence of integers . A key feature of the construction is an explicit projection from the skew column RSK dynamics to the classical box-ball system; under this projection, the rigging is precisely the Kerov--Kirillov--Reshetikhin rigging of the associated box-ball configuration. Our proof uses two commuting affine crystal structures on pairs of skew tableaux and a novel connectivity theorem for distinguished subgraphs of tensor products of Kirillov--Reshetikhin crystals. We also derive Greene-type formulas for the soliton lengths in terms of last-passage percolation on the associated cylindrical environment. Finally, by taking generating functions in the linearizing coordinates, we obtain bijective proofs of Cauchy and Kawanaka--Littlewood-type identities for transformed Hall--Littlewood polynomials.
Abstract
In this paper, we establish a general comparison theorem for the order statistics of the edge eigenvectors for generalized Wigner matrices. Consequently, we derive the Gumbel law for the maximal edge eigenvector component and prove the universality of the Gaussian fluctuations of the order statistics in an intermediate regime close to the maximum. In addition, our comparison result also implies a quantitative first order estimate for moderately small order statistics.
Abstract
This work investigates the stability of outliers from the circular law, via the convergence of their associated diagonal overlaps between eigenvectors - also known as the squared eigenvalue condition numbers. We consider and compare two paradigmatic cases, namely: 1) the Complex Ginibre Ensemble conditioned on the existence of an outlier, and 2) the outlier induced by a rank-one Hermitian perturbation of a Complex Ginibre matrix. In both cases, we prove almost sure convergence towards a specific constant that only depends on the radius of the outlier and its status - either conditioned or induced. These results can be generalized to other complex integrable ensembles with the same techniques, and complement our understanding of eigenvalue stability in non-Hermitian ensembles.
Abstract
In this paper, we consider the Cauchy problem for the defocusing modified Korteweg-de Vries (mKdV) equation with a step-like initial data. Based on the Riemann-Hilbert problem associated with the mKdV equation, we derive the long-time asymptotic expansion of the solution to the defocusing mKdV equation in two transition regions using the nonlinear steepest descent method. It comes out that the leading term in the expansion is shown to match the corresponding background constants. The subleading term, however, decays at the order , and its coefficient is derived from the associated Painlevé \uppercase\expandafter{\romannumeral34\relax} model.
Abstract
We prove that the Opdam-Cherednik kernel, also known as the nonsymmetric Opdam hypergeometric function, can be written as the Laplace transform of a positive measure supported on the convex hull of the Weyl group orbit of its argument. As a consequence, the trigonometric Dunkl intertwining operator is positivity preserving. The main ingredient in the proof is a new formula for the Opdam-Cherednik kernel as a degeneration of nonsymmetric Macdonald polynomials. As a further application, we prove majorization inequalities for Macdonald polynomials and Heckman-Opdam hypergeometric functions associated with arbitrary root systems.
Abstract
Many-body integrable systems can be understood as a gas of quasiparticles. They propagate ballistically and drive large-scale transport. However, with the exception of the hard rods system, no tools have been available to numerically track such quasiparticles. Focusing on the Toda fluid, whose integrability relies on the availability of a Lax pair, we present a numerical scheme to track quasiparticle trajectories as determined by the time-dependent eigenvectors of the Lax matrix. Simulating the Toda fluid in thermal equilibrium, this tracking scheme is used to numerical confirm Brownian motion of a quasiparticle. Simulated is also the motion of a tagged particle. Our numerical results for the diffusion constant matches with a novel TBA prediction. We believe our numerical scheme can be extended to other classical many-particle models possessing a Lax matrix.
Abstract
The asymmetric simple exclusion process (ASEP) is an important particle model with deep connections to orthogonal polynomials. Motivated by this connection, Corteel and Williams introduced the Koornwinder moments at , which generalize the moments of Askey--Wilson polynomials. They showed that the partition function of the two-species ASEP is equal to for a one-row partition . In this paper, we investigate a conjecture of Rains on the positivity of the minimal numerator of the Koornwinder moment . We derive the first explicit formula for this moment, thereby obtaining a precise formulation of the conjecture by determining the minimal denominator of . We also propose a generalization of the conjecture for the more general Koornwinder moments indexed by two partitions at special parameter values. We prove the generalized Rains' conjecture in two special cases: and . For , we construct a lattice path model and obtain a combinatorial formula for in terms of non-intersecting lattice paths. For , we establish an explicit product formula for and give a combinatorial interpretation using lecture hall tableaux.
Abstract
We prove a local bijection between alternating sign matrices (ASMs) and descending plane partitions (DPPs) in the refined fibres of Behrend--Di Francesco--Zinn-Justin. The map is constructed by moving a boundary line through the ASM and DPP pictures. At each step only a bounded neighbourhood of that boundary is changed. In growth-diagram language, the boundary states are ordinary ASM row frontiers and shifted DPP path frontiers, and the elementary moves are reversible Fomin-type square rules. The proof checks that overlapping local moves are confluent, that every square is invertible, and that the completed boundary again satisfies the ordinary ASM and DPP coordinate inequalities. The only terminal choice is resolved by a radius-two fan rule on the visible height contour. The same local framework gives a statistic on refined ASMs whose generating function is the DPP degree enumerator. The statistic has a Striker--Fulmek quadratic core \[ Q_{\rm SF}(A)= \sum_{1\le r<s\le n}\ \sum_{1\le b<a\le n} (n-r+1)A_{r,a}A_{s,b}, \] plus a bounded fan-boundary correction and, outside the normal boundary gauge, an explicit compact-potential endpoint term. In the normal gauge this is \[ W_{\rm ASM}^{\rm SF}(A)=Q_{\rm SF}(A)+B_{\rm fan}(A). \] Thus the paper gives both a local ASM--DPP bijection and a statistic-preserving refined \(q\)-enumeration. The supplementary files reproduce selected local checks and the displayed arithmetic examples.
Abstract
A generalised inverse first passage problem asks whether, given a probability measure on , one can find a boundary such that the stopping time:\[τ:=\inf\left\{t:Λ\int_0^t ω(W_s-b(s))ds \geq U\right\}\] has distribution , where , and is a monotonic decreasing function. We construct a branching-selection particle system whose hydrodynamic limit is governed by a free boundary problem and connect this to the generalised inverse first passage problem. In the -particle system, particles move as independent Brownian motions, branch at a prescribed rate, and are removed at a rate proportional to their location relative to a position which is a function of the empirical distribution. We identify the limit of as the solution of the inverse first passage problem.
Abstract
A symmetric function is called Schur-positive if it admits an expansion in the Schur basis with nonnegative coefficients. In this paper, we study the Schur-positivity of symmetric functions naturally associated with set partitions, with respect to a descent set function that considers i as descent, if i and i+1 share a block in the partition. The Schur expansion involves hook-shaped Young diagrams, and the corresponding coefficients are given by Touchard-Riordan polynomials, which enumerate matchings by their number of crossings.
Abstract
Balanced words on a finite alphabet are those words in which every letter of the alphabet occurs the same number of times. The notion of descents and major index extends in a natural way to words. It is known that the bivariate generating polynomials for descents and major index over balanced words on the alphabet with occurrences each is palindromic, but a bijective proof has been missing even for balanced binary words. We give an explicit bijection proving this result. For permutations (which are also balanced), our bijection is different from the complementation map. We also show that for balanced binary words, this bijection simultaneously flips the ascent and comajor index as well.
Abstract
In this paper, we introduce and prove a set of rank recurrence relations for -Whittaker and Macdonald operators. We also derive an explicit expression for the -th power of the -Whittaker operators in terms of the -deformed binomial probability distribution, and we express the rank recurrence relations for the Macdonald operators in terms of the Cauchy determinant.
Abstract
We introduce a censored version of the stochastic six-vertex model. We show that for parameters , this model started from the initial condition is stochastically dominated at any time by the blocking measure. This is a partial analog of the censoring inequality for monotone spin systems. In particular, this result allows us to control the behavior of second-class particles. The proof uses parabolic Kazhdan--Lusztig -polynomials, whose appearance is explained using a connection between the stochastic six-vertex model and the Iwahori--Hecke algebras of symmetric groups. Furthermore, we find an intertwining relation for this process using normalized parabolic Kazhdan--Lusztig -polynomials as an intertwining kernel.
Abstract
We prove the Harish-Chandra isomorphism theorem for supersymmetric spaces, describing the polynomial algebra of eigenvalues of invariant differential operators. The polynomials obtained satisfy novel invariance conditions, which remain somewhat mysterious. We also prove the Harish-Chandra isomorphism for ghost distributions, which satisfy a `square root' of the invariance conditions coming from invariant differential operators. All proofs are algebraic, and rely on a rank-one reduction argument and the Chevalley restriction theorem.
Abstract
We study the large-time asymptotic behavior of a new KdV soliton gas. We first introduce a pure-soliton Riemann--Hilbert(RH) problem with \(2N\) poles and two different types of residue conditions. We show that, as \(N\to\infty\), this discrete problem converges to primitive-potential RH problem introduced by Dyachenko, Zakharov, and Zakharov, and the jump matrix of this soliton gas RH problem has two nonzero reflection coefficients. To analyze the large-time behavior, we apply the Deift--Zhou nonlinear steepest descent method together with an appropriate \(g\)-function mechanism. Through a sequence of transformations, the original RH problem is reduced to explicitly solvable model problems on an associated hyperelliptic Riemann surface. This allows us to derive an explicit leading-order asymptotic formula for the solution in terms of Jacobi elliptic function. The result provides a rigorous asymptotic description of a new KdV soliton gas and extends the available analysis beyond the previously studied case \(r_2\equiv 0\).
Abstract
In this paper, we introduce and study a family of diffusion processes associated with the Pitman-Yor process and the two-parameter Poisson-Dirichlet distribution. The diffusion coefficients indexed by a non-negative parameter are smaller than the corresponding one-parameter models in terms of quadratic forms or bilinear forms when is positive. The well known Petrov's diffusion corresponds to among the unlabelled diffusions. If is the same as the stable parameter in the Pitman-Yor process, we obtain both labelled and unlabelled reversible diffusion processes with the Pitman-Yor process and the two-parameter Poisson-Dirichlet distribution as the corresponding reversible measures. We construct these processes analytically through Dirichlet forms. In comparison with existing models in the literature, our models possess two fundamental new features. Firstly, our labelled model is characterized by an explicit generator, which is the first among all models studied so far. This makes it possible to establish the crucial integration-by-parts formula. Secondly, a novel foundational structure of the two-parameter distributions is the existence of a diversity index (a multiple of local time) for the positive stable parameter. By slowing down the diffusion in our model, the essential role of the diversity index is revealed in the evolution of the population. Additionally, we also obtain properties including ergodicity, path behaviour, and finite dimensional approximations.
Abstract
We prove that the Gromov--Witten invariants of the elliptic orbifold lines , , and satisfy a certain system of Hirota Quadratic (or Bilinear) Equations. Our result is the analogue of the so-called Toda conjecture in the Gromov-Witten theory of or more precisely its non-extended version. A new feature in our constructions is a certain bilinear operator whose principal symbol can be expressed in terms of elliptic theta functions.
Abstract
We study the KPZ equation on a circle with an additive spatial perturbation , where is a spacetime white noise and is a smooth spatial function. When , it is well-known that the unique invariant measure is the Brownian bridge. In the presence of the perturbation, we show that the equation admits a unique invariant measure that is absolutely continuous with respect to the Brownian bridge. We further prove the measure has a finite relative entropy with respect to the law of the bridge and that, for any , the corresponding Radon-Nikodym derivative belongs to , provided that is sufficiently small. The proof uses the discretization and mollification scheme of \cite{FQ}, together with an application of the log-Sobolev and spectral gap inequalities for the underlying Gaussian measure.
Abstract
We define and study the rectangular finite free heat flow, a dynamical system on polynomials that plays the role of the heat equation in the setting of rectangular finite free probability. We show several equivalent characterizations of the evolution (including PDE and gradient flow formulations), establish basic properties of the dynamics, and determine the asymptotic distributions of the polynomial roots in the long-time and high-degree limits. We also discuss connections with Calogero-Moser systems and Dunkl processes, and we show that the rectangular finite free heat flow describes the mean curvature expansion of a family of compact Lie group orbits.
Abstract
We generalize the Ahlswede--Daykin inequality (1978) to a Schur positive \emph{ADS inequality}, which also contains the Lam--Postnikov--Pylyavskyy inequality (2007) as a special case. We then present a number of further generalizations and applications. Notably, we resolve Mihalcea's conjecture on log-supermodularity of stable Grothendieck polynomials.
Abstract
In this manuscript, we explore the application of neural networks to predict the natural parameter of Schramm-Loewner Evolution (SLE) theory. SLE is a family of random fractal curves that has significant implications in Statistical Mechanics and Conformal Field Theory. This parameter plays an important role in the theory as there are models of Planar Statistical Physics that are proven to have SLE as scaling limits as well as others that are conjectured to have this limit for various choices of the parameter . In addition, there are three different statistical behaviors of the SLE curves as the parameter changes in Leveraging the powerful pattern recognition capabilities of neural networks, this study aims to develop a predictive model that can estimate the parameter with good accuracy.
Abstract
Permuted-basement Macdonald polynomials are nonsymmetric generalizations of symmetric Macdonald polynomials that form a basis for the polynomial ring for each fixed . There are combinatorial formulas for them as generating functions over composition-shaped non-attacking fillings. In this extended abstract, we bijectively prove identities for the relationship between , , , and . These identities correspond to two combinatorial operations on non-attacking fillings: (1) swapping adjacent entries in the basement, generalizing a result of Alexandersson (2019), and (2) swapping adjacent parts in the shape, which yields a straightening rule for expanding in the polynomials .
Abstract
In this paper, we prove universality of the distribution of the cokernels of a random -adic matrix with inhomogeneously balanced columns. More precisely, let be an integer and be a random matrix over whose -th column is -balanced. We prove that if as for every , then the cokernels of converge in distribution, as , to the same limiting law as the cokernels of Haar-random matrices over . This extends a universality theorem of Nguyen and Wood to random -adic matrices with inhomogeneously balanced columns.
Abstract
We introduce a novel exclusion process with two conservation laws, mass and energy, designed to mimic the essential features of continuous systems like interacting oscillators within the framework of interacting particle systems. This distinguishes our model from conventional multi-species processes where only particle numbers are conserved. As a basis for our fluctuation analysis, we first show that applying nonlinear fluctuating hydrodynamics (NFH) to this model reveals a wide variety of universality classes depending on the parameter choices. The main objective of this work is to study the stationary fluctuations of these conserved quantities. For a suitable choice of parameters, we rigorously show that the fluctuation fields converge to uncoupled stochastic Burgers equations (SBE) in the scaling limit. The proof relies on the second-order Boltzmann-Gibbs principle that we establish for this model, along with the spectral gap estimate and the equivalence of ensembles. Of independent interest is our general proof of the diagonalizability of the Jacobian matrix for the macroscopic current with distinct real eigenvalues. While this property is often taken as given in the physics literature, we establish it rigorously for multi-component systems even when the eigenvectors cannot be explicitly computed, offering a firm mathematical foundation for a broad class of models.
Abstract
Fix a multivariate polynomial in non-commuting variables of arbitrary degree, and consider independent complex Ginibre matrices . We prove that the empirical spectral distribution of converges as tends to infinity to the so-called Brown measure of evaluated at free circular variables. For polynomials of degree at most 2, the convergence was proven by Cook, Guionnet, and Husson \cite{cook2022spectrum}, and we prove that the convergence in fact holds for polynomials of any degree. The main step in the proof is a least singular value lower bound for for almost all complex shifts , and we prove this via a least singular value lower bound for a wide class of tensorized Ginibre matrices of finite type with a deterministic shift, which is of independent interest. We further show that the Brown measure convergence holds beyond Gaussians: the same convergence holds when the entry law has mean 0, variance 1, bounded density on and finite moments of all orders.
May 2026
Abstract
We introduce the \emph{tandem exclusion process} (TEP), a one-dimensional stochastic lattice model motivated by tandem running in ants. Particles evolve through two cooperative local transitions, at rate (leader advancement) and at rate (follower recovery). We prove that the stationary measure on the dynamically active sector is the Gibbs measure , where and counts neighboring occupied pairs, and derive exact closed-form expressions for the stationary current and spatial correlations using transfer-matrix methods. The current is asymmetric under particle--hole exchange , with its maximum occurring at densities strictly larger than , in contrast to the symmetric current of the totally asymmetric simple exclusion process (TASEP). For , cooperative dynamics enhances the current above the TASEP value and generates strong spatial clustering; in the limit , the current approaches , corresponding to nearly unconstrained collective transport. These results suggest that tandem coordination alone can substantially enhance collective transport efficiency at moderate and high densities, even without pheromone-mediated long-range communication.
Abstract
Gasper gave a two-step extension of the Askey-Roy formula for a beta-integral type defined on the unit circle. This involved increasing the number of parameters and adding a balancing condition. We derive a multidimensional generalization of the final formula in Gasper's extension. By considering a two-step degeneration involving parameters our generalization becomes Tarasov-Varchenko's formula -- itself a multidimensional generalization of the original Askey-Roy formula. The proof is done by Aomoto's method. This generalized integration by parts strategy makes use of a class of functions known as the fundamental invariants of type .
Abstract
We consider the directed polymer model on , in an i.i.d.\ random environment , focusing on the critical dimension . Our main contribution is to give a sharp lower bound on the free energy in the high-temperature regime. Our proof uses a percolation argument inspired by Lacoin (2010), for which we introduce a key property of bounded ``-energy'': this property quantifies the regularity of the polymer measures at diffusive scales and we show that it propagates along open paths. Writing for the quenched free energy, and setting and , our lower bound combined with Theorem 2.8 of Berger, Caravenna, and Turchi (2025) gives
Abstract
We consider matrices with independent, identically distributed entries, and prove that the sequence of measures converge to the Gaussian Multiplicative Chaos in the full subcritical regime as . Our result holds for both symmetry classes and in particular is new even for real Ginibre matrices, and is the first such convergence for any non-invariant ensemble of random matrices. We also establish the asymptotics for the -point function of at any collection of mesoscopically separated points . Our methods are analytic and probabilistic in nature, relying in part on the dynamical approach based on Dyson Brownian motion.
Abstract
We study combinatorial structures arising from finite-time transition probabilities of the Totally Asymmetric Simple Exclusion Process with open boundary conditions. While much of the existing combinatorial theory regarding the TASEP concerns the steady-state distribution, we focus instead on the transient dynamics. We first show that the enumeration of transition sequences between two configurations of the open TASEP is equivalent to the enumeration of standard Young tableaux of a family of non-classical shapes which have been of recent interest in the combinatorial literature. This extends to the open-boundary setting the correspondence between the TASEP with periodic boundaries and cylindric tableaux. We then introduce a family of tableau-like objects associated with Young diagrams in which repetitions of cells are allowed, subject to the partial order induced by the diagram. For each diagram, we collect the numbers of these objects into an exponential generating function. We prove that the entries of the homogeneous open TASEP transition matrix can be expressed as signed sums of such generating functions over suitable families of diagrams. This gives a combinatorial and order-theoretic interpretation of finite-time transition probabilities for the open TASEP, analogous to the combinatorial mappings known for steady-state probabilities.
Abstract
We prove asymptotic equivalents for finite-level representations of symmetric groups, that is, for Young diagrams having all but finitely many boxes on their first row. We deduce that random walks on symmetric groups generated by conjugacy classes with a macroscopic number of fixed points have a Poissonian cutoff profile. We also prove that the random involution walk exhibits cutoff and find its cutoff profile. Finally, we obtain numerics for the random transposition walk on a deck of 52 cards, giving concrete estimates on the question that originally motivated Diaconis and Shahshahani.
Abstract
While Macroscopic Fluctuation Theory (MFT) has been highly successful in analyzing non-equilibrium steady states, its application to non-steady-state processes remains limited. In this study, we apply MFT to the relaxation process of one-dimensional boundary-driven diffusive systems coupled to particle reservoirs at both ends. We exactly derive the current variance for systems with a constant diffusion coefficient and arbitrary mobility, as well as the cumulant generating function for the current in Reflective Brownian Motion (RBM). Our results demonstrate that non-steady current fluctuations during the approach to a steady state can be quantitatively described within the MFT framework.
Abstract
This paper establishes new upper bounds for the sum of the largest eigenvalues of symmetric matrices. When applied to the adjacency matrix of a graph, our results improve upon a related bound due to Mohar {\bf [On the sum of k largest eigenvalues of graphs and symmetric matrices, J. Combin. Theory Ser. B 99 (2009) 306--313]}. Furthermore, in the case of the Laplacian matrix, we prove that the well-known Brouwer's conjecture {\bf [Spectra of Graphs, Springer, New York, 2012]} holds for small values of for almost all graphs, thereby taking a significant step toward its complete resolution.
Abstract
We prove that the continuous-time, single-flip Glauber dynamics for lozenge tilings of the size- hexagon mix in time . This was predicted to hold on fairly general domains of diameter (on the basis of the ``Lifshitz law'' heuristic) but had previously only been established in domains such that the associated limit shape has no frozen facets. To access the hexagon, we introduce a multi-scale comparison argument between the height function of the random tiling evolving under the Glauber dynamics and the limit shape of a volume-tilted tiling (whose tilting parameter varies suitably in time).
Abstract
The symmetric simple exclusion process (SSEP) is a paradigmatic model of classical non-equilibrium dynamics. Exact results for large deviations of particle current in the SSEP have been obtained in various settings using integrability-based methods. In this Article, we discuss how these results are modified in the presence of localized slow bonds. We consider three conventional geometries: (a) a finite one-dimensional lattice weakly coupled to unequal reservoirs at its boundaries, (b) a semi-infinite one-dimensional lattice weakly coupled to a boundary reservoir, and (c) an infinite one-dimensional lattice with localized slow bonds near the origin. For each case, we present exact expressions for the large deviation function of current and validate them through rare-event simulations based on the cloning algorithm. In connection with our results, we present an elementary derivation of the exact large deviation function for the current in the semi-infinite SSEP, complementing recent results obtained through more elaborate techniques.
Abstract
We give a complete classification of the eternal solutions for the KPZ fixed point. Each of these is a (possibly infinite) patching together of the known eternal solutions, called Busemann functions. The resulting evolution of the KPZ fixed point exhibits a shock at each of the boundaries between the different Busemann functions. Moving forward in time, the shocks coalesce, while moving backwards in time, additional shocks can form. We describe several geometric properties of this tree of shocks.
Abstract
We prove that the maximal dimension of an irreducible representation of the symmetric group satisfies for some constant . This answers a question raised by Vershik--Kerov in 1985.
Abstract
We prove a Fredholm determinantal identity for the tilted Toeplitz minor generalizing the Borodin-Okounkov-Geronimo-Case (BOGC) identity to oblique splittings of the Hardy space. The tilts enter only through an oblique projection that multiplies the trace-class kernel inside the Fredholm determinant; the BOGC operator constructed from is unchanged. Baik-Liao-Liu (arXiv:2603.01964) and Liu-Tripathi (arXiv:2604.24747) have recently shown that the same tilted Toeplitz minor admits a contour Fredholm-determinantal representation, in connection with the periodic Totally Asymmetric Simple Exclusion Process (TASEP). In the periodic TASEP application of Baik-Liao-Liu, the formula plays an important role in identifying the periodic KPZ fixed point with general initial data. Our formula is a companion to their Fredholm determinant and readily reduces to the original BOGC identity. The one-sided tilted Toeplitz minor (that is, when all ) admits a bialternant form recovering Schur and Grothendieck polynomials as special cases. A Cauchy-Binet expansion realizes as a restricted sum over partitions of products of Jacobi-Trudi type determinants, generalizing Gessel's theorem. In the pure-shift setting this specializes to a skew Schur expansion. Finally, for finite Laurent exponential symbols, we record explicit resolvent-block flow identities and formulate the associated finite-dimensional closure problem. We also illustrate a possible asymptotic application leading to finite-rank perturbations of the Airy kernel.
Abstract
This manuscript studies a special case of the Hurwitz enumeration problem: for branched covers from genus g compact Riemann surface to the Riemann sphere, with three branch points, and require the branching data at one of the branch points to be of the partition (m 1^{n-m}), obtain a formula of Hurwitz number. The Hurwitz enumeration problem can be transformed into enumeration of a class of unrooted hypermaps. We first provide a enumeration formula for rooted hypermaps, thereby obtaining the weighted Hurwitz numbers. Next give the quantitative relationship between the enumeration of unrooted hypermaps and that of rooted hypermaps in the general setting. Finally, combining these two results, we obtain a formula for the unrooted hypermaps or the Hurwitz numbers in the special setting.
Abstract
Pipedreams and bumpless pipedreams are two combinatorial models that compute double Grothendieck polynomials. While studying matrix Schubert varieties, Pechenik, Speyer, and Weigandt defined a permutation statistic that captures the leading monomial of the top-degree components of a Grothendieck polynomial. Combinatorially, their result implies that there exists a unique pipedream (or bumpless pipedream) with row weight and column weight . A construction of such a pipedream was subsequently given by Chou and Yu. In this paper, we resolve the bumpless pipedream version of this problem by providing an explicit algorithm.
Abstract
The genus--g Fuss--Catalan (FC) number counts the number of ways to obtain a genus-g surface by identifying the edges of a pn--gon via p-valent hyperedges. For p=2 these are the genus--g Catalan numbers which are generated as the trace correlations in the Gaussian matrix model (GUE). Here we construct a simple two-matrix model which generates the higher-genus Fuss--Catalan numbers for any p as the coefficients of its 1/N-expansion. We obtain exact sum rules and an explicit formula for the higher-genus Fuss--Catalan numbers which generalises the Harer--Zagier formula to p>2. We discuss the relation of the higher-genus FC numbers to the intersection numbers and the Euler characteristic of the moduli space of spin-p curves.
Abstract
We study geodesics in the Brochette first-passage percolation model, where edges on the same axis-parallel line share a common random passage time, inducing long-range dependence. We focus on the maximal transversal deviation H n of geodesics from the origin to ne 1 . We prove existence of geodesics under mild assumptions and establish the order of magnitude of H n depending on the behavior of the passage-time distribution near its infimum. These results yield explicit wandering exponents in this dependent setting.
Abstract
We prove a complete characterization of the extremal invariant measures for half-space geometric last-passage percolation with an arbitrary boundary parameter. This is the first result of its kind for a model in the KPZ universality class that has boundary effects and an unbounded domain. A description of a class of invariant measures was previously given in a work of Barraquand and Corwin, where it was conjectured that these should comprise all extremal invariant measures. To complete the classification, we prove a one force--one solution principle: when started in the distant past from an arbitrary initial condition with a given asymptotic slope at , the recentered solution at time converges to a process which is distributed as the associated invariant measure with the specified slope. This limiting process is called the Busemann process, the first of its kind constructed for a half-space model. The Busemann process across all slopes is distributed as the joint invariant measure for geometric half-space LPP, recently constructed by Dauvergne and Zhang. There, it was conjectured that the constructed family of jointly invariant measures comprises all extremal jointly invariant measures; our analysis also confirms this conjecture. When the model has a strong (attractive) boundary, the collection of slopes for the invariant measures has a discontinuity, which does not arise in the full-space case. To handle this difficulty, we combine the control of the directions of semi-infinite geodesics with techniques from the theory of half-space Gibbsian line ensembles. Along the way, we classify the set of directions of semi-infinite geodesics for half-space geometric LPP, confirming a recent conjecture of Dauvergne and Zhang.
Abstract
We present a "black box" proof of mean-field near-critical behaviour for a family of functions on () satisfying a short list of assumptions. The functions represent two-point functions of a lattice statistical mechanical model in the subcritical or critical regimes, and are proved to have decay of the form , for any . The black box applies to several models for which commonplace methods can be used to verify the assumptions. Applications include models of self-avoiding walk, percolation, spins (Ising, XY, ), and lattice trees, all above their upper critical dimensions. The proof is based on random walk techniques, and provides a new, unified, probabilistic, and relatively simple proof of mean-field near-critical behaviour.
Abstract
In the directed landscape, the Martin boundary coincides with the horofunction boundary. We show that functions in this boundary are precisely the eternal solutions possessing a spatial growth rate, and that the minimal Martin boundary is given by the Busemann functions. Moreover, every eternal solution can be expressed as a max-plus convex combination of countably many Busemann functions. Horofunctions are exactly those eternal solutions that admit a representation in terms of at most two Busemann functions with a common growth rate. As a consequence of instability, not all horofunctions are Busemann functions, and the Martin boundary is strictly larger than its minimal part.
Abstract
We show that the log-determinant of leading principal minors of large non-Hermitian random matrices converges in distribution to a 2+1 dimensional Gaussian field, which is logarithmically correlated for the parabolic distance, reminiscent to the Edwards-Wilkinson universality class.
Abstract
We introduce the extra slow Tamari lattices, a new family of lattices defined on faithfully balanced tableaux. These tableaux arise naturally from the representation theory of type \( A \) quivers, and our construction extends the classical Tamari lattice and the slow Tamari lattice. We explicitly describe meets and joins in the extra slow Tamari lattices, and then prove that they are lattices. We then show that they are semidistributive, trim, polygonal, and congruence uniform. Their join-irreducible elements are described in terms of a three-color analogue of the positive roots of type \( A \), which leads to descriptions of their spines and congruence lattices. We also obtain several enumerative results for the extra slow Tamari lattices and their associated structures. Finally, we derive new structural and enumerative results for the slow Tamari lattices.
Abstract
We study the distribution of eigenvalues of Haar-random matrices over among algebraic extensions of . Our results give -adic analogues of the real-eigenvalue counting results of Edelman-Kostlan-Shub for the real Ginibre ensemble, but with a different degree behavior: while real eigenvalues form only a vanishing proportion in the real Ginibre ensemble, -adic eigenvalues are asymptotically evenly distributed among possible extension degrees. We also show that the maximal unramified extension captures all but a bounded expected number of eigenvalues, and that the expected number of eigenvalues outside has a finite positive limit with an explicit upper bound. The proof uses correlation function formulas from the author's previous joint work with Van Peski (arXiv:2601.06283), together with uniform estimates over varying finite extensions. We also prove analogous results for roots of random Haar polynomials over , using the correlation function formulas of Caruso (arXiv:2110.03942). These polynomial results are -adic analogues of the real-root counting results of Edelman-Kostlan, again with behavior different from the real setting.
Abstract
This paper studies limit profiles for the separation distance. A limit profile records the limiting shape of the distance to stationarity inside the cutoff window, at times of the form . We start with two famous card shuffles, a general setup for inverse riffle shuffles and random transpositions, and we determine their separation distance limit profiles. We then develop a spectral comparison technique and study continuity properties in the style of [Nes24; Nes25], adapted to separation distance. The comparison method is illustrated through random transpositions, as well as random walks on product groups and the hypercube.
Abstract
Despite similarities between models exhibiting absorbing phase transitions (APTs) and those showing Kardar-Parisi-Zhang (KPZ) growth, the relationship between these universal fluctuations has remained elusive. We numerically study (1+1)-dimensional interfaces of (2+1)-dimensional models showing APTs of directed percolation (DP) and compact directed percolation (CDP) classes with an active boundary, finding a universal crossover from short-time APT-governed fluctuations to long-time KPZ fluctuations. Upon rescaling time and length by the APT correlation time and length, the cumulants of the interface height distributions collapse onto a single scaling function. The fluctuation properties of the discrete Domany-Kinzel model and the continuum stochastic Fisher-Kolmogorov-Petrovsky-Piskunov (sFKPP) equation coincide, indicating that the KPZ growth parameters are determined solely by fundamental properties of the APT. For the CDP sFKPP equation, a dimensionless parameter tunes both the critical interface distribution and the KPZ parameters, with the interface properties of the biased voter model recovered in a limiting case. These results uncover a universal crossover in which KPZ fluctuations emerge from APT fluctuations at long times, linking paradigmatic universality classes of nonequilibrium scale-invariant phenomena.
Abstract
We introduce coupled double Poisson brackets on an associative algebra as pairs consisting of a generalized Van den Bergh's double Poisson bracket and a generalized Fairon--McCulloch's right double Poisson bracket subject to a cross-Jacobi identity. Each of Van den Bergh's double brackets, Fairon--McCulloch's right double brackets, and also Ginzburg--Schedler's wheeled Poisson brackets induces a -invariant Poisson structure on the representation scheme parametrizing -dimensional representations of , thereby satisfying the Kontsevich--Rosenberg principle. Wheeled Poisson brackets seem to be the most general such structures, and while their relation to Van den Bergh's double Poisson brackets is known, their relation to Fairon--McCulloch's right double Poisson brackets has remained open. We fill this gap and establish a bijection between pairs of coupled double Poisson brackets and wheeled Poisson brackets of Ginzburg and Schedler. On free polynomial algebras, we furthermore establish a one-to-one correspondence between linear coupled double Poisson brackets and a new algebraic structure that we call Poisson-left-pre-Lie algebras, and describe quadratic ones via solutions of the associative and classical Yang--Baxter equations satisfying a compatibility condition.
Abstract
We investigate the multicritical scaling limit of the shifted Schur measures. Under an appropriate scaling limit and specific conditions on the continuous parameters, we explicitly determine the limit shape of strict partitions distributed according to the shifted Schur measure. We then show that, under a multicritical condition, the edge scaling limit of the correlation function converges to a determinant of the higher-order Airy kernel. This rigorously demonstrates a transition from a Pfaffian point process to a determinantal distribution in the scaling limit.
Abstract
We develop a hydrodynamic theory for a height-dependent version of the totally asymmetric simple exclusion process in which the jump rate at a growth site is sampled from a macroscopic two-dimensional speed function evaluated at the spatial coordinate and the current height level. The speed function is allowed to have discontinuities along locally finitely many curves. Through the TASEP height-function representation, the process is coupled to an inhomogeneous directed last-passage percolation model whose exponential rates vary discontinuously in the two macroscopic LPP coordinates. Combining the law of large numbers for this last-passage model with an extension of the variational coupling method, we prove a hydrodynamic limit for the height function and for the associated particle density. The limiting current is characterised by a Lax-Oleinik type variational formula built from the discontinuous last-passage shape function. We then identify the first-order PDE structure selected by the microscopic dynamics. At points of differentiability of the limiting current and continuity of the sampled coefficient, the current solves a Hamilton-Jacobi equation whose Hamiltonian depends discontinuously on the spatial variable and on the value of the solution itself. At discontinuities, the variational formula leads to a natural envelope-based discontinuous viscosity formulation, and we prove that the limiting current satisfies this formulation. Finally, when the coefficient has only spatial discontinuities, we prove uniqueness of the Hamilton-Jacobi solution in the natural class of nondecreasing Lipschitz currents, and identify its spatial derivative as the maximal-current weak solution of the associated scalar conservation law with discontinuous flux.
Abstract
We study inhomogeneous -Whittaker polynomials which extend both -Whittaker and stable Grothendieck polynomials. We prove that inhomogeneous -Whittaker polynomials (in countably many variables) form a basis of certain commutative ring extending the ring of symmetric functions to a subring of its completion. We then describe positive specializations of that ring and relate them with a subset of Macdonald-positive specializations of the ring of symmetric functions. We also show some related probability distributions obtained from positive specializations of inhomogeneous -Whittaker polynomials.
Abstract
Determinantal point processes (DPPs) have emerged as a kernelized alternative to vanilla independent sampling for generating efficient minibatches, coresets and other parsimonious representations of large-scale datasets. While theoretical foundations and promising empirical performance have been demonstrated, there are two challenges for current proposals for DPP-based coresets or minibatches. The first is the need for families of DPPs with certain key variance reduction properties, usually constructed in a continuous setting, of which there are few known examples. The second is the need for an ad-hoc construction of a discrete DPP defined on a given dataset, that inherits such variance reduction. In this work, we contribute to the programme of establishing DPPs as a subsampling toolbox for ML by advancing on these two fronts. First, we propose new DPPs on the Euclidean space based on wavelets, with provably better accuracy guarantees than the best known rates. Second, we introduce a general method to convert such continuous DPPs, which are more amenable to proving analytical statements, into discrete kernels, which are pertinent for subsampling tasks such as minibatch and coreset constructions. This conversion mechanism simultaneously preserves the desired variance decay and reveals a low-rank decomposition of the discrete kernel, which makes sampling the corresponding DPP computationally inexpensive. En route, we enlarge the class of ML tasks amenable to improvements via DPP-based minibatches and coresets to include objective functions with arbitrarily low regularity, and rate guarantees that explicitly adapt to this regularity.
Abstract
Bergeron--Ceballos--Küstner introduced the -Fibonomial coefficients \( \qfibonom{m+n}{n}\), gave a combinatorial interpretation of the -Fibonomial coefficients via a weighted path-domino tiling model, and conjectured that these polynomials are unimodal. We prove the conjecture for . For the case, we give a combinatorial proof of both unimodality and symmetry by defining a nearly symmetric saturated chain decomposition on the set of tilings. For all three cases, we give an algebraic proof. Finally, for the case, we establish a more general unimodality result for certain products of -analogs and propose several related conjectures.
Abstract
The hierarchical Pitman-Yor process is a discrete random measure used as a prior in Bayesian nonparametrics. It is motivated by the study of groups of clustered data exhibiting power law behavior. Our focus in this paper is on the Gaussian behavior of a family of statistics, namely the power sum symmetric polynomials for the vector of weights of the process, as the concentration parameters tend to infinity. We establish a central limit theorem and obtain explicit representations for the asymptotic variance, with the latter clearly showing the impact of each component in the hierarchical structure. These results are crucial for understanding the asymptotic behavior of the sampling formulas associated with the process. In comparison with the known results for the hierarchical Dirichlet process, the results for the hierarchical Pitman-Yor process are mathematically more challenging and structurally more revealing of power law behavior.
Abstract
We study linear spectral statistics of high dimensional sample covariance matrices in a regime where the empirical spectral distribution remains governed by the classical sample covariance law but the fluctuation theory is nonclassical. Our starting point is a decomposition of the covariance of centered quadratic forms into a universal Gaussian part and a model dependent fourth order correction. This leads to an abstract framework, termed GHOST, for universal Gaussian central limit theorems under structured fourth order effects. Under this framework, we prove a Gaussian central limit theorem for linear spectral statistics, with explicit mean and covariance corrections determined by a bilinear fourth order kernel. Boundary examples show that the conditions are close to necessary for a broad universal Gaussian closure. We then develop a blockwise mixed radial model that verifies the abstract assumptions and makes the correction explicit. The correction splits into an entrywise fourth moment component and a lockwise energy fluctuation component. The latter may change the fluctuation scale, leading to a phase transition at the level of fluctuations. As an application, we study sphericity testing. Under the spherical null, the general correction collapses to a single scalar parameter, yielding a feasible data driven correction of John's test.
Abstract
We study a class of interacting particle systems on with two types. Particles evolve by independent jumps sampled from a fixed distribution, with type-dependent jump rates , and stochastic type switching driven by non-local order-based interactions. The switching rates depend on the empirical distribution through the proportion of opposite-type particles located ahead, leading to a nonlinear and discontinuous dependence on the empirical measure outside the standard Lipschitz McKean-Vlasov framework. Our first main result is a law of large numbers for the empirical measure process: we prove convergence, along subsequences, to a deterministic measure-valued process characterized by a McKean-Vlasov equation. The proof combines tightness in Wasserstein space with a martingale characterization of limit points. A uniqueness argument based on a Kolmogorov-Smirnov-type distance adapted to the ordering structure yields convergence of the full empirical measure sequence and, in turn, propagation of chaos on finite time intervals. We then study the long-time behavior of the limiting dynamics. Because the system has persistent drift, invariant distributions do not arise; instead, we analyze traveling waves, corresponding to stationary profiles in a moving frame. For exponential jump distributions, the associated non-local integro-differential system admits a local description. In the regime , this further reduces to a coupled system of non-linear ODEs, allowing a phase-plane analysis that yields a traveling wave as a heteroclinic orbit connecting two equilibria. We also identify the wave speed and mass partition, and derive tail asymptotics by spectral analysis of the linearized system.
Abstract
The -species particle-exchange process (PEP()) is an exclusion process in which particles of different species exchange positions on neighbouring sites with rates chosen such that the invariant measure on the discrete torus is a product measure. We address the large-scale hydrodynamic behaviour of this process which yields a system of coupled inviscid Burgers equations. This system of conservation laws is shown to admit Riemann invariants for arbitrary from which explicit solutions of the Riemann problem in terms of shock waves and rarefaction fans are obtained. We also introduce the open PEP(), in which particles are exchanged with boundary reservoirs. For a distinguished manifold of boundary rates, we prove that the invariant measure is the same product measure as in the periodic system. The hydrodynamic description in terms of Riemann invariants is used to derive the stationary phase diagram explicitly in terms of microscopic boundary rates. In the generic case, the steady state exhibits phases, with boundary-induced phase transitions analogous to those of the single-species asymmetric simple exclusion process.
Abstract
We consider an Asymmetric Exclusion Process evolving on parallel mutually interacting lanes with neighbouring nearest hoppings of hardcore particles. Number of particles on each lane is conserved. We find a choice of the hopping rates, for which the process has Bernouilli stationary product measure, and calculate the stationary particle currents as a function of average particle densities.
Abstract
We construct the full edge scaling limit of the singular values of Brownian motion on the general linear group starting from general conditions. We show that the limiting paths solve an infinite system of SDE with log-interaction and have a Gibbs resampling property with exponential Brownian bridges. Moreover, we show that the evolution of the limiting rescaled reverse characteristic polynomial solves a stochastic partial differential equation with a non-linear multiplicative noise and linear drift. From a special initial condition the resulting line ensemble coincides, in logarithmic coordinates, with a line ensemble constructed by Ahn which arises as a universal scaling limit of singular values of products of random matrices. We prove some analogous results on the evolution of limiting characteristic polynomials for two models whose stationary measures are given by the Hua-Pickrell and Bessel stochastic zeta functions respectively.
Abstract
The theory of Bailey's transform provides a systematic method for deriving -identities, the key factor of which is the Bailey pair. The concept of Bailey pair was first extended to bilateral version by Paule. In this paper, following Rogers' work on Fourier series, we derive two asymmetric bilateral Bailey pairs. By inserting them into the bilateral Bailey chains, we obtain several identities of Rogers-Ramanujan type, Andrews-Gordon type and also identities on false theta functions. Furthermore, based on the Bailey lattice due to Dousse, Jouhet and Konan, we get an asymmetric bilateral Bailey lemma which leads to identities on Appell-Lerch series. Moreover, by using the asymmetric bilateral Bailey lemmas due to Andrews and Warnaar, we get some identities on false theta functions and the generalized Hecke-type series.
Abstract
We study the behavior of the elastic polymer, a model of a directed polymer in a continuous Gaussian random environment that is independent in time and correlated in space, as the dimension of the environment is taken to infinity. We give an explicit asymptotic formula for the free energy, which is given in terms of the distribution of the inner product of two sampled configurations, which we also obtain an implicit formula for. From this, we provide an explicit characterization of both the low- and high-temperature phases of this model in terms of the spatial correlation function of the environment. We find asymptotics for the wandering exponent when the spatial correlation function is either an exponential or a power-law decay. Our results show that when the correlations are either suitably weak or short ranged, the model is asymptotically diffusive. On the other hand, for suitably strong long ranged correlations, the model is asymptotically superdiffusive. Moreover, we show that this transition coincides exactly with another transition where the model goes from being one-step replica symmetry breaking to full-step replica symmetry breaking. This rigorously confirms many of the findings of Mezard and Parisi [53] in the physics literature.
Abstract
We establish new transformation formulas involving theta functions for certain bilateral basic hypergeometric series. From these, we construct companion -series for a class of -series such that the asymptotic expansion of their quotient admits a simple closed form. This allows us to prove some conjectures of McIntosh on asymptotic transformations of -series.
Abstract
We provide multiple combinatorial expansion formulas - in terms of snake graphs, labelled posets, matrices, and -walks - for elements in generalized cluster algebras associated to arcs on punctured orbifolds and illustrate their equivalence. This work generalizes and unifies existing work on combinatorial expansion formulas from surfaces and unpunctured orbifolds.
Abstract
Andrews investigated parity conditions in the Rogers-Ramanujan-Gordon theorem. Under the conditions that even parts or odd parts appear an even number of times, Andrews discovered two Rogers-Ramanujan-Gordon type partition theorems and derived corresponding generating functions. In the Rogers-Ramanujan-Gordon theorem, there are two parameters and , where is the maximum number of consecutive parts and , and is the maximum number of parts equal to . Andrews' first theorem deals with the case , while the second theorem concerns the case where is even and is odd. These two partition identities have different infinite product forms on the right-hand side. In this paper, we consider the case and use Bailey's lemma to obtain an Andrews-Gordon type identity whose right-hand side coincides with that of Andrews' identity for the case . We were unable to find a suitable combinatorial interpretation of the infinite sum form of this expression in terms of partitions, but with the help of lattice paths, we provide an appropriate combinatorial interpretation.
Abstract
Let be a basic classical Lie superalgebra, a boundary admissible level of , where is a positive integer and is the dual Coxeter number of . In this paper, we classify the irreducible ordinary modules for the affine vertex operator superalgebra associated to any basic classical Lie superalgebra . More specifically, if is a basic classical Lie superalgebra of type I, we prove that has exactly inequivalent irreducible ordinary modules. If is a finite dimensional simple Lie algebra or a basic classical Lie superalgebra of type II, we prove that itself is the only irreducible ordinary -module.
Abstract
On a finite weighted graph, the dimer model is a probability measure on its dimer covers, that assigns to any cover a probability proportional to the product of the weights of its edges. For planar bipartite graphs, dimer correlations are encoded by the inverse of the so-called Kasteleyn matrix; for a large graph, typically taken as a finite domain in a periodic graph, this inverse matrix is known explicitly only for a handful of examples. In all previously known examples, the Newton polygon -- a convex lattice polygon that classifies periodic graphs up to local moves -- is either a triangle or a quadrilateral. Our main results are the following. For any (minimal) periodic planar bipartite graph, we construct an -dimensional family of finite subgraphs for which we obtain an explicit inverse Kasteleyn matrix; here is the number of sides of the Newton polygon. Their boundaries are formed by zig-zag paths and their overall shape is reminiscent of an astroid; we call them astroidal zig-zag graphs (AZ graphs). If the Newton polygon is the unit square then the corresponding AZ graph is the celebrated Aztec diamond with its size as the parameter. Our inverse Kasteleyn matrices are given by a double contour integral on the corresponding spectral curve for any Fock weighting of the graph. This includes, in particular, all periodic weightings. For periodic weightings, we asymptotically analyze the resulting inverse Kasteleyn matrices. We establish a phase separation in large AZ graphs into asymptotically frozen, rough (liquid), and smooth (gaseous) regions, and obtain an explicit parametrization of the `arctic curve'. We also compute the deterministic limit of the height function, known as the limit shape, and prove the convergence of the local dimer correlations to the translation-invariant Gibbs measure of the slope predicted by the limit shape.
Abstract
Wilson loops provide the central gauge-invariant probe of confinement in lattice gauge theory. This survey reviews the statistical-mechanical formulation of lattice gauge ensembles, the strong-coupling and duality mechanisms behind area laws, finite-temperature and continuum scaling diagnostics, and the mathematical status of Wilson-loop confinement.
Abstract
Based on the ballistic macroscopic fluctuation theory, the integration of the spin correlation function (spin conductivity) is analyzed for the spin-1/2 XXZ chain in the critical regime. In the time when the magnetization of an infinite spin chain fluctuates from an initial state with a wavelength as long as the infinite length , the equal-time two-point spin correlation function is scaled up to . In the state where the ballistic spin transport decays at high temperature , the diffusive transport remains on a large scale. We show that the spin conductivity is proportional to in the limit and its high temperature proportionality constant diverges in the case where one-quasiparticle magnetization is infinitely large. This analysis informs that the superdiffusive spin transport is driven by the -scaled long-range spin correlation and sheds a light on the dynamic scaling in spin transport at the isotropic point.
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.