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)
July 2026
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.
Abstract
We analyse the eigenvectors of the adjacency matrix of a random inhomogeneous graph constructed from a specified degree sequence. We assume that the empirical degree sequence has bounded mean and variance. We show that near the edges of the spectrum, the eigenvectors are semilocalized in the sense that their mass concentrates around a small set of resonant vertices. For the extremal eigenvalues, we establish localization around a single vertex. In order to obtain effective estimates in the presence of highly inhomogeneous degrees, we introduce a new economical pruning procedure that carefully extracts a forest from the original graph, whose adjacency matrix is compared to that of the original graph using a suitably constructed local coupling to random trees with independent edges.
Abstract
We present a new method for constructing Hadamard matrices that combines transformer neural networks with local search in the PatternBoost framework. Our approach is designed for extremely sparse combinatorial search problems and is particularly effective for Hadamard matrices of Goethals--Seidel type, where Fourier methods permit fast scoring and optimisation. For orders between and , it produces large numbers of inequivalent Hadamard matrices, and in harder cases it succeeds where local search from random initialisation fails. The largest example found by our method has order . In addition to these new constructions, our experiments reveal that the transformer can discover and exploit useful hidden symmetry in the search space.
Abstract
Let for some fixed , and let be an random matrix with entries in , where each row is independently and uniformly sampled from the set of all vectors in containing exactly ones. We show that the empirical spectral distribution of the appropriately rescaled matrix converges in probability to the circular law provided that . As a crucial element of the proof, we obtain quantitative lower bounds on the smallest singular value of the shifted matrices whenever and for some absolute positive constant .
Abstract
We prove that two half-space models in the KPZ universality class, exponential last-passage percolation and a family of Poisson-avoiding metrics generalizing colored TASEP, converge to a common scaling limit. This scaling limit is the directed landscape in half-space, a random directed metric in the half-plane indexed by a parameter which determines the strength of the boundary interaction. As part of our analysis, we characterize the half-space directed landscape in terms of the half-space KPZ fixed point, and prove convergence of geodesics. We also give an explicit construction of joint stationary measures (or horizons) in half-space for the log-gamma polymer, the KPZ equation, exponential and geometric last passage percolation, and the directed landscape itself.
Abstract
We study the direct and inverse spectral theory for a class of finite Hermitian banded matrices. Using the theory of matrix orthogonal polynomials, we provide an explicit procedure for reconstructing a banded matrix from a matrix-valued measure that encodes its spectral data. We establish necessary and sufficient conditions for a measure to be the spectral measure of a matrix in the examined class. We further analyze the connections between this spectral analysis, block tridiagonalization algorithms, and the Toda lattice evolution on banded matrices.
Abstract
We classify irreducible finite-dimensional modules of a collection of real Lie superalgebras that includes the simple ones, their classical variants, complex Lie superalgebras after restriction of scalars, and all real Lie algebras. Our strategy is to reduce this classification to determining the orbits of the parity and conjugation functors on irreducible modules of the complexifications of the aforementioned algebras. Then we provide explicit results for the computation of these orbits. For Lie superalgebras of basic type or of type , our classification applies to any highest-weight parametrization of irreducible complex modules with respect to an arbitrary Borel subalgebra. As a consequence, in the special case of real simple Lie algebras we obtain a new perspective on the classification of real simple modules and establish a conceptual connection with Kostant's cascade of strongly orthogonal roots.
Abstract
A color Lie algebra is a generalization of a Lie (super)algebra by an Abelian group . The underlying vector space and defining relations of the algebra are graded by , and the color Lie algebra admits graded Casimir elements. Furthermore, its loop algebra admits graded central extensions. We present a general method for constructing 2nd order graded Casimir elements and graded central extensions for a given color Lie algebra and its loop algebra, respectively. We also show that there exists a large class of color Lie algebras admitting such graded Casimir elements or central extensions by providing three examples, namely, for , and and for .
Abstract
In this paper, we study the long-time behavior of a stochastic heat equation with multiplicative noise and localized control. We begin by analyzing the uncontrolled dynamics and derive explicit decay rates for both mean-square and almost sure exponential stability. These estimates show that the two notions of stability may hold under different conditions on the parameters, reflecting the interplay between the drift and the multiplicative noise. We then introduce a finite-dimensional feedback control acting on a measurable subset of positive measure, built from finitely many Fourier modes of the solution. In particular, we show that the number of controlled modes determines the decay rate and allows for arbitrarily fast stabilization in the mean-square sense. As a consequence, almost sure exponential stability is recovered via a probabilistic argument, so that both notions of stability are achieved within the same framework and with the same decay rate. As an application, we provide a new proof of controllability for the stochastic heat equation based on an iterative construction of adapted controls in feedback form, avoiding the use of the adjoint equation.
Abstract
Consider critical Bernoulli percolation on for large; let be distinct points in . We prove that the probability that all lie in the same open cluster, rescaled by an appropriate power of , converges as to an explicit constant. This confirms a conjecture of Aizenman and Newman.
Abstract
Given a pair of finite posets , the function counting integer-valued order preserving extensions of an order preserving map is given by a piecewise polynomial in . We provide a criterion for the nonnegativity of the coefficients of these multivariate polynomials and apply it to show that marked order polytopes of skew shapes are Ehrhart positive in a multivariate sense. This extends recent results of Ferroni-Morales-Panova on order polytopes of skew shapes and proves conjectures on the Ehrhart positivity of skew Gelfand-Tsetlin polytopes and -generalized Pitman-Stanley polytopes due to Alexanderson-Alhajjar and Dugan-Hegarty-Morales-Raymond, respectively.
Abstract
We derive the Euler (hyperbolic) hydrodynamic limit for the directed exclusion process (DEP), a one-dimensional conservative interacting particle system that preserves particle-hole symmetry while breaking left-right symmetry. The proof relies on an explicit multi-process coupling, which guarantees a strong form of attractiveness and macroscopic stability for the particle system. Further open questions about DEP are briefly discussed.
Abstract
Consider the following nonlinear one-dimensional stochastic fractional heat equation where is the fractional Laplacian on for , and is a Gaussian noise that is white in time and behaves in space as a fractional Brownian motion with Hurst index satisfying . When , Hu and Wang ({\it Ann. Inst. Henri Poincaré Probab. Stat.} {\bf 58} (2022) 379-423) studied the well-posedness of the solution and its Hölder continuity, removing the technical condition that was previously assumed in Hu et al. ({\it Ann. Probab.} {\bf 45} (2017) 4561-4616). Their approach relied on working in a weighted space with a suitable power decay function. For the case , inspired by Hu and Wang, we investigate the well-posedness of the stochastic fractional heat equation without imposing the technical condition of , which was required in the earlier work of Liu and Mao ({\it Bull. Sci. Math.} {\bf181} (2022) 103207). In our analysis, precise estimates of the heat kernel associated with the fractional Laplacian play a crucial role.
Abstract
We study the -deformation of the Demazure product model from arXiv:2407.21653. Consider the longest element in written as a reduced word in simple transpositions. Independently delete each transposition with probability and apply the -Demazure product to the remaining ones. We show that the law of the resulting permutation converges as to a deterministic permuton, which coincides with the case studied in arXiv:2407.21653 for adjusted probability . This resolves Conjecture 1.13 from arXiv:2407.21653 and identifies the limiting permuton explicitly.
Abstract
In a recent work [Phys. Rev. E 109, L042102 (2024)], interesting dimensional crossovers [from two- to one-dimensional (2D to 1D) scaling] were found in the growth of Kardar-Parisi-Zhang (KPZ) interfaces on rectangular substrates, with lateral sizes . Here, we extend this study to other universality classes for interface growth -- specifically, the Edwards-Wilkinson (EW), the Mullins-Herring (MH), and the Villain-Lai Das Sarma (VLDS) classes. From extensive simulations, we demonstrate that, in all systems with sufficiently large aspect ratio , the roughness scales with time in the growth regime as for and for , where in most cases. For the VLDS class, this crossover is also observed in the height distribution (HD), which approaches its characteristic probability density function for the 2D case at short times () and then crosses over to the asymptotic 1D HD. Dimensional crossovers are also found in the steady state regime, both in the roughness scaling as well as in the VLDS HD, which interpolate between the 2D and 1D ones as increases. The particular case , with , is also discussed in detail and reveals interesting features of the investigated systems. For instance, there exist a `special' exponent such that the temporal crossover cannot be observed for . Moreover, this leads the saturation roughness to display a nonuniversal scaling: , with .
Abstract
Quasisymmetric functions in superspace were introduced as a natural extension of classical quasisymmetric functions involving both commuting and anticommuting variables. In this paper, we first provide a characterization of the algebra of quasisymmetric functions in superspace as an algebra of invariants under a quasisymmetrizing action of the symmetric group. Furthermore, we complete the superspace analogue of the classical hierarchy of combinatorial Hopf algebras by introducing the algebra of quasisymmetric functions in noncommuting variables in superspace. We endow this algebra with a Hopf superalgebra structure and thoroughly investigate its -basis and monomial basis, which are indexed by set supercompositions. By restricting to the minimal elements of the underlying poset, we construct the Hopf superalgebra of superpermutations, serving as the superspace analogue of the Malvenuto--Reutenauer algebra. We provide explicit product and coproduct formulas for these bases in terms of super-shuffles and global descents. Finally, via an abelianization morphism, we apply these noncommutative structures to derive a product formula for fundamental quasisymmetric functions in superspace.
Abstract
This article derives full asymptotic expansions for integrals of the form \[ \int_{0}^{1}f(u)(1+q\cdot u^{n})^{w/n}du \] as , with parameters real and , or positive for . We relate the coefficients of the asymptotic expansions to Nielsen's generalized polylogarithms. For , we obtain an expansion in terms of multiple zeta values, which in this setting, reduce to ordinary zeta values. A key point is that for , the integrals typically produce alternating multiple zeta values; we formulate a precise symmetry constraint on the relevant coefficient sequence under which all coefficients reduce to polynomials in ordinary zeta values. We also translate this symmetry into a statement about a binomial transform, and we verify the condition for several classical Appell-type families, like Euler, Bernoulli, Genocchi, and Hermite. Finally, we obtain precise results about the convergence of norms of random variables.
Abstract
We study transport dynamics of free fermions subject to the external monitoring of a conserved charge over an extensive region. Focusing on bipartition protocols, we consider monitoring the total particle number over half of the system, and study the profiles of local charges and currents at hydrodynamic scales. While the Lindbladian describing the averaged dynamics is non-local, we show that the profiles can be understood in terms of localized impurities. We present a general framework based on the generalized hydrodynamics (GHD) picture, allowing for a hybrid numerical-analytic solution of the quench dynamics at hydrodynamic scales. We illustrate our approach for domain-wall initial states, showing that monitoring leads to discontinuities in the profiles that become more pronounced as the rate increases and that lead to the absence of transport in the Zeno limit of infinite monitoring rates. Our GHD framework could be naturally extended to interacting systems, paving the way for a systematic study of transport of integrable models subject to extensive-charge measurements.
Abstract
We study the problem of describing the set of real functionals on the quotient of the ring of symmetric functions that are nonnegative on the images of certain modified Hall-Littlewood symmetric functions. This question is equivalent to the problem, posed in [Adv Math 395, p.108087 (2022)], of describing the set of coadjoint-invariant measures for unitary groups over a finite field in the infinite-dimensional setting. Our main results constitute partial progress towards this problem. Firstly, we show that the desired set of functionals is very large, in the sense that it contains explicit families of examples depending on infinitely many parameters. Secondly, we provide an analogue of Kerov's mixing construction that produces new sought after functionals from known old ones. This construction depends on an explicit "-twisted action" of on itself and the resulting dual map that makes into a comodule. Finally, our third main result explains the relation between the -twisted comultiplication and the usual comultiplication on .
Abstract
In this paper, we study estimation of parameters in a two-parameter Potts model with colors and coupling matrix . We characterize concrete sufficient conditions for existence of the pseudo-likelihood estimator of the Potts model, in terms of the local magnetic fields, and give sufficient conditions for the validity of the above characterization. We then provide sufficient criteria for estimation of both parameters at the optimal rate . In particular, if is the scaled adjacency matrix of a graph , then we show that joint estimation is possible if either has bounded degree or is irregular. In contrast, we give an example of a graph sequence which is approximately regular and dense, where no consistent estimator exists. We also show that one-parameter estimation at the optimal rate holds under much milder conditions when the other parameter is known. Along the way, we develop a concentration result for mean-field Potts models using the framework of nonlinear large deviations. Compared to the Ising case, our results for the Potts case require a novel analysis across multiple colors.
Abstract
The enumeration of lozenge tilings of hexagons with holes has been studied intensively in recent years. Researchers tried to find shapes and positions of holes in hexagonal regions so that the number of lozenge tilings of the resulting regions is given by a simple product formula. In the present work, we consider new regions that are hybrids of regions studied by the first author (hexagons with intrusions) and Ciucu (F-cored hexagons). Then, we show that the tiling generating functions of these new regions under a certain weight are given by simple product formulas. To give a proof, we present shuffling theorems for lozenge tilings of hexagons with intrusions, which give simple relations between the tiling generating functions of two related hexagonal regions with intrusions.
Abstract
P. Di Francesco first introduced the "Aztec triangle" in his study of the relationship between the twenty-vertex model and domino tilings. He conjectured an exact formula for the number of tilings of the Aztec triangle, and it has since been proved by several authors. In an attempt to prove the conjecture, M. Ciucu showed that the tiling number of the Aztec triangle divides the tiling number of a new region called the "cruciform region," a superposition of two Aztec rectangles. Ciucu proved that the number of domino tilings of a cruciform region is given by a simple product formula. In this paper, we generalize Ciucu's tiling formula by providing a generating-function formula for the cruciform region.
Abstract
Understanding the realization of thermal equilibrium through the thermalization process in a many-body system is a fundamental and complex scientific question, bridging thermodynamics and classical dynamics and connecting to a host of physical phenomena, such as mechanical instabilities in a thermal environment. In this work, based on the harmonic lattice model, we investigate the thermalization process in both velocity and coordinate spaces, by examining microscopic dynamics on the atomic level. We show the distinct relaxation rates of the transverse and longitudinal components of the velocity, reveal the power law governing the nonlinear proliferation of dominant frequencies, and observe the concurrent rapid proliferations of frequencies and topological defects. We also show that the lattice system's persistent out-of-plane deformations exhibit two-stage fluctuation behaviors, characterized by distinct power laws of fractional exponents and associated with the broken up-down symmetry. This work demonstrates the rich dynamics underlying the thermalization process, and advances our understanding on the dynamical adaptations of many-body systems to external disturbances.
Abstract
In the present paper, we study the equilibrium fluctuations of a particle system in infinite volume with two conserved quantities and long-range dependence. More specifically, the model of interest is the so-called ABC model, in which three types of particles (A, B and C) exchange their locations between and at a rate that depends on the type of particles involved and is proportional to for . After rigorously identifying the normal modes associated to the conserved quantities (the density of particles of types and , say), we prove that their fluctuations converge to independent fractional stochastic partial differential equations (SPDEs), which are either Gaussian or the Stochastic Burgers equation, and whose nature is determined by the microscopic range of dependence and the strength of the asymmetry.
Abstract
We consider a class of large-scale interacting systems with one conservation law satisfying the ``degree-preserving property'', and study the classification of their invariant measures and their hydrodynamic limits. Under a few basic conditions, we show that if the generator of the process preserves the degree of polynomials of the state variables up to two, then the marginals of any product invariant measure of the process must belong to one of six specific distributions. This classification result is essentially a consequence of a known result in statistics on univariate natural exponential families due to C.N. Morris, which we apply here for the first time in the context of microscopic stochastic systems. In particular, we introduce a new model whose invariant measure is given by the generalized hyperbolic secant distribution. Additionally, under the same conditions, we show that, regardless of the specific model, the hydrodynamic equation is always the classical heat equation, with a diffusion coefficient that depends on the model. Our proof is based on deriving uniform bounds on second-moments of state variables, whose proof is achieved by relating a correlation function to a one-dimensional random walk whose jump rates are model-dependent and obtaining sharp bounds on its occupation times on specific domains.
Abstract
We consider certain noncolliding interacting particle systems driven by Brownian noise. A key example is drifted Brownian motions conditioned not to intersect and related models of eigenvalues of Hermitian random matrices. We establish limit theorems for the extremal particle. We find: (i) the scaling limit of the largest eigenvalue of Brownian motion over Hermitian, positive-definite matrices, (ii) Airy process limit for the largest eigenvalue of Dyson's Brownian motion for GUE started from generic initial conditions, and (iii) a Fredholm determinant formula for the running maximum of the top path among noncolliding Brownian bridges and, as a byproduct, a new formula for the law of largest eigenvalue in the Laguerre Orthogonal Ensemble as well as for a related point-to-line last passage percolation model.
Abstract
We show that the increments of the KPZ fixed point started from arbitrary initial data are \emph{mutually} absolutely continuous with respect to Brownian motion with diffusion parameter on compacts, extending the one-sided Brownian absolute continuity relation of the KPZ fixed point established in \cite{sarkar2021brownian}. We also show that additive Brownian motion is absolutely continuous with respect to the centred Airy sheet on compacts, but it is not mutually absolutely continuous globally. As applications, we show that with probability strictly between zero and one, there exist record times of the KPZ fixed point away from any reference point, obtain a characterisation for the hitting probabilities of the graph of the KPZ fixed point to be positive in terms of a certain thermal capacity in the sense of \cite{watson1978corrigendum, watson1978thermal} and compute essential suprema of Hausdorff dimensions of these random intersections. Finally, we compute essential suprema of Hausdorff dimensions of images of subsets in the plane under the Airy sheet and give a condition for the positivity of their Lebesgue measure in terms of Bessel-Riesz capacity.
Abstract
We use random matrix theory for the Circular Unitary Ensemble (CUE) to study moments of derivatives of the Riemann zeta function shifted a small distance from the critical line. The corresponding CUE moments are studied in the limit of large matrix size in two regimes: when the spectral parameter is (1) suitably far inside the unit disc, and (2) at a small distance from the unit circle. In case (1), we obtain an asymptotic formula as a combinatorial sum over contingency tables, while in case (2) we obtain a sum over certain determinants with multiplicative coefficients given by Kostka numbers. The latter result is also valid exactly on the unit circle. Then, we consider the analogous problem for mean values of derivatives of the zeta function with suitable shifts. Assuming the Lindelöf hypothesis, we show that this mean value gives rise to the same sum over contingency tables obtained in the CUE. For sufficiently low-order moments, we establish this result unconditionally.
Abstract
We consider Erdős-Rényi-type random hypergraphs that are non-uniform, in the sense that hyperedges of different sizes may coexist, and inhomogeneous, in that connection probabilities may depend on the hyperedge size. All parameters are allowed to scale with the hypergraph size. We study the random adjacency matrix whose -entry counts the number of hyperedges containing both vertices and , and characterize its expected limiting spectral distribution in terms of the connection probabilities and the hyperedge sizes. We provide a Pastur-type condition, in the sense of Chatterjee (2005), under which the matrix can be Gaussianized, as well as a more restrictive but simpler sufficient condition in terms of the generalized average degree of the model. As a second main result, based on such a Gaussianization, we characterize the limiting spectral distributions under non-sparse conditions as semicircle laws with an explicit parametric variance. The latter can be expressed as a convex combination of the variances arising in the uniform cases, with coefficients determined by the trade-off between the different sources of inhomogeneity.
Abstract
For first passage percolation (FPP) on Euclidean lattices with , it is expected that the variance of the first passage time between two points grows sublinearly in the distance with a universal exponent strictly smaller than . Following Kesten's upper bound (Ann. Appl. Probab., 1993) on the variance, Benjamini, Kalai and Schramm (Ann. Probab., 2003) used hypercontractivity to obtain an improvement of a factor of when passage times take two values with equal probability. This was later extended to more general classes of passage time distributions. However, unlike in exactly solvable planar models in last passage percolation where the variance is known to be , the best known upper bound for the variance of passage times has remained in all non-trivial variants of FPP. For a class of rotationally invariant Riemannian FPP on the plane, we show that the variance is for some . Our argument uses fluctuation estimates for passage times and geodesics derived in Basu, Sidoravicius and Sly (2023) together with a multi-scale argument to establish that the geodesic exhibits disorder chaos, i.e., upon resampling a small fraction of the underlying randomness, the updated geodesic has on average a small overlap with the original one; this, established at a large number of scales, leads to a polynomial improvement of the variance bound.
Abstract
We obtain bivariate asymptotics for one part monotone Hurwitz numbers in high genus (i.e. as both the size and the genus go to infinity). To do so, we start with a linear recurrence for these numbers obtained by Do and Chaudhuri. Then, we apply a recent method developped by Elvey-Price, Fang, Wallner and the second author to extract asymptotics from such recurrences.
Abstract
This work establishes a large deviation principle for the spectral measure of the Lax matrix associated to the periodic Toda chain of particles, subject to a generalised Gibbs measure. This large deviation principle is governed by a rate function which can be regarded as a generalisation of the free energy of the system. Such a large deviation principle is proven both for the model when the momentum is constrained to be zero and when it is allowed to fluctuate. Moreover, the large deviation principle is proven directly on the level of the representation of the generalised Gibbs partition function given in terms of the variables realising the classical separation of variables, \textit{i.e.} rectify the equations of motion. As such, this work paves the way towards the computation of the thermodynamic limit of dynamical correlation functions in the Toda chain subject to generalised Gibbs ensemble statistics.
March 2026
Abstract
In this paper, we introduce and study random Schur measures whose parameters are sampled from the Circular Unitary Ensemble. We show that Schur measures with CUE disorder exhibit behavior reminiscent of spin glasses.
Abstract
Fully inhomogeneous spin Hall-Littlewood symmetric rational functions arise as partition functions of certain path configurations in the higher spin six vertex models. They are multiparameter generalizations of the classical Hall-Littlewood symmetric polynomials. We establish two new generalizations of the classical Littlewood identity, where we express a weighted sum of 's over all partitions as a product of the Littlewood kernel and another simple product in one case, and a product of the Littlewood kernel and a Pfaffian in the other case. As a corollary we obtain a novel Littlewood identity for Hall-Littlewood symmetric polynomials. We also elaborate on the newly established connection between the fully inhomogeneous spin Hall-Littlewood symmetric rational functions and the modified Robbins polynomials, the latter being multivariate generating functions for alternating sign matrices. This connection allowed us to discover the two generalizations of the Littlewood identity and we provide a bijection between the underlying combinatorial models in the case where is strictly decreasing.
Abstract
The conjectured three generic local bulk statistics amongst all non-Hermitian random matrix symmetry classes have recently been extended to three generic local edge statistics. We study analytically and numerically complex spacing ratios and nearest-neighbour (NN) spacing distributions that characterise such local statistics. We choose the three simplest representatives of these universality classes, given by the Gaussian ensembles of complex Ginibre, complex symmetric and complex self-dual matrices, denoted by class A, AI and AII. In the first part, we analytically study the complex spacing ratio in class A, at finite matrix size . Introducing a conditional point process, we simplify existing expressions and show why an uncontrolled approximation introduced earlier converges well in the large- limit in the bulk. When specifying to the elliptic Ginibre ensemble, we present a parameter-dependent surmise for the complex spacing ratio, interpolating to that of the Gaussian unitary ensemble (GUE), where such a surmise is very accurate. In the second numerical part, we compare complex spacing ratios, its moments, and NN spacing distributions for all three ensembles with that of uncorrelated points, the two-dimensional (2D) Poisson process, both in the bulk and at the edge. The varying degree of repulsion within these different edge universality classes can be well understood in terms of an effective 2D Coulomb gas description, at different values of inverse temperature . We find indications that the complex spacing ratio does not fully unfold the local statistics at the edge. Finally we verify that for small argument, in all three symmetry classes the NN spacing distributions in the bulk and at the edge are consistent with the universal cubic repulsion.
Abstract
We study a family of periodically weighted Aztec diamond dimer models near their turning points. We establish that, asymptotically, as , their fluctuations there, scaled by , are described by a marked GUE-corners process. This limiting point process is constructed by assigning a Bernoulli mark independently to each particle in a realization of the GUE-corners process. The Bernoulli parameters associated with the random marks reflect the periodicity of the model in the limit. To prove this result we use a double-contour integral representation of the inverse Kasteleyn matrix on a higher-genus Riemann surface, which is well-suited for asymptotic analysis.
Abstract
Two-dimensional Coulomb gases on an annulus at a special inverse temperature are studied by using the orthogonal polynomial method borrowed from the theory of random matrices. The correlation functions among the Coulomb gas molecules are written in determinant forms and their asymptotic forms in the thermodynamic limit are evaluated. When the Coulomb gas system has a continuous rotational symmetry, the corresponding orthogonal polynomials are monomials, and one can see a universal behavior of the correlation functions in a thin annulus limit. In a system with a discrete rotational symmetry, the corresponding orthogonal polynomials are not in general monomials, and a breakdown of the universality is observed.