Leonid Petrov — arXiv Feed

The Harish-Chandra isomorphism for supersymmetric spaces and ghost distributions

2026-06-09T15:57:05-04:00

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.

Large-time asymptotics of a new KdV soliton gas

2026-06-09T07:30:44-04:00

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$.

Dynamical models for the two-parameter Poisson-Dirichlet distribution and the Pitman-Yor process

2026-06-07T22:04:42-04:00

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 $\gamma$ are smaller than the corresponding one-parameter models in terms of quadratic forms or bilinear forms when $\gamma$ is positive. The well known Petrov's diffusion corresponds to $\gamma=0$ among the unlabelled diffusions. If $\gamma$ is the same as the stable parameter $\alpha$ 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.

Elliptic orbifold lines and integrable hierarchies

2026-06-07T21:21:53-04:00

We prove that the Gromov--Witten invariants of the elliptic orbifold lines $\mathbf{P}^1_{3,3,3}$ , $\mathbf{P}^1_{2,4,4}$ , and $\mathbf{P}^1_{2,3,6}$ 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 $\mathbf{P}^1$ 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.

The Gaussian structure of a perturbed KPZ

2026-06-05T11:02:13-04:00

We study the KPZ equation on a circle with an additive spatial perturbation $\partial_t h=\tfrac12\Delta h+\tfrac12|\nabla h|^2+\xi+ V$ , where $\xi$ is a spacetime white noise and $V$ is a smooth spatial function. When $V=0$ , 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 $p\in(1,\infty)$ , the corresponding Radon-Nikodym derivative belongs to $L^p$ , provided that $\int V^2$ 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.

The Rectangular Finite Free Heat Flow

2026-06-04T23:03:17-04:00

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.

Correlation inequalities for Schur positivity

2026-06-04T16:09:45-04:00

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.

Neural Networks and Schramm-Loewner Evolutions

2026-06-01T13:51:45-04:00

In this manuscript, we explore the application of neural networks to predict the natural parameter $\kappa \geq 0$ of Schramm-Loewner Evolution (SLE $_\kappa$ ) theory. SLE $_\kappa$ is a family of random fractal curves that has significant implications in Statistical Mechanics and Conformal Field Theory. This parameter $\kappa \geq 0$ 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 $\kappa \geq 0$ . In addition, there are three different statistical behaviors of the SLE curves as the parameter $\kappa$ changes in $[0, \infty).$ Leveraging the powerful pattern recognition capabilities of neural networks, this study aims to develop a predictive model that can estimate the $\kappa$ parameter with good accuracy.

Shape changing identities for permuted-basement nonsymmetric Macdonald polynomials

2026-06-01T11:43:41-04:00

Permuted-basement Macdonald polynomials $E_\alpha^\sigma(\mathbf{x};q,t)$ are nonsymmetric generalizations of symmetric Macdonald polynomials that form a basis for the polynomial ring $\mathbb{Q}(q,t)[\mathbf{x}]$ for each fixed $\sigma$ . 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 $E_\alpha^\sigma$ , $E_\alpha^{\sigma s_i}$ , $E_{s_i\alpha}^\sigma$ , and $E_{s_i\alpha}^{\sigma s_i}$ . 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 $E_\alpha^\sigma$ in the polynomials $\{E_{s_i\alpha}^\tau \}_\tau$ .

Universality of the cokernels of random p-adic matrices with inhomogeneously balanced columns

2026-06-01T08:37:15-04:00

In this paper, we prove universality of the distribution of the cokernels of a random $p$ -adic matrix with inhomogeneously balanced columns. More precisely, let $u \ge 0$ be an integer and $A(n)$ be a random $n \times (n+u)$ matrix over $\mathbb{Z}_p$ whose $i$ -th column is $\alpha_n(i)$ -balanced. We prove that if $\sum_{i=1}^{n+u} \exp(-\varepsilon \alpha_n(i)n) \to 0$ as $n \to \infty$ for every $\varepsilon>0$ , then the cokernels of $A(n)$ converge in distribution, as $n \to \infty$ , to the same limiting law as the cokernels of Haar-random $n \times (n+u)$ matrices over $\mathbb{Z}_p$ . This extends a universality theorem of Nguyen and Wood to random $p$ -adic matrices with inhomogeneously balanced columns.

Stationary fluctuations for an exclusion process with mass and energy conservation

2026-06-01T05:02:39-04:00

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.

Brown measure convergence for the spectrum of polynomials in Ginibre matrices

2026-06-01T00:18:29-04:00

Fix a multivariate polynomial $\mathfrak{p}$ in $n$ non-commuting variables of arbitrary degree, and consider $n$ independent $N\times N$ complex Ginibre matrices $X_1^N,\cdots,X_n^N$ . We prove that the empirical spectral distribution of $P^N=\mathfrak{p}(X_1^N,\cdots,X_n^N)$ converges as $N$ tends to infinity to the so-called Brown measure of $\mathfrak{p}$ 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 $\mathfrak{p}$ of any degree. The main step in the proof is a least singular value lower bound for $P^N-z$ for almost all complex shifts $z$ , 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 $\mathbb{C}$ and finite moments of all orders.

Tandem Exclusion Process

2026-05-31T22:15:47-04:00

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, $110\to101$ at rate $\alpha$ (leader advancement) and $101\to011$ at rate $\beta$ (follower recovery). We prove that the stationary measure on the dynamically active sector is the Gibbs measure $\pi \propto q^{H(\eta)}$ , where $q=\beta/\alpha$ and $H(\eta)$ 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 $\rho \mapsto1-\rho$ , with its maximum occurring at densities strictly larger than $1/2$ , in contrast to the symmetric current $\rho(1-\rho)$ of the totally asymmetric simple exclusion process (TASEP). For $q>1$ , cooperative dynamics enhances the current above the TASEP value and generates strong spatial clustering; in the limit $q\to\infty$ , the current approaches $J\to\alpha \rho$ , 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.

q-Selberg extensions of Gasper's two stage q-beta integrals

2026-05-31T06:19:39-04:00

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 $BC$ .

Sharp behavior of the free energy for the two-dimensional directed polymer model

2026-05-28T20:58:44-04:00

We consider the directed polymer model on $\mathbb{Z}^d$ , in an i.i.d.\ random environment $\omega=(\omega_{n,x})_{n\geq 0,x\in\mathbb Z^d}$ , focusing on the critical dimension $d=2$ . 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 `` $\log$ -energy'': this property quantifies the regularity of the polymer measures at diffusive scales and we show that it propagates along open paths. Writing $\mathfrak{f}(\beta)$ for the quenched free energy, and setting $\lambda(\beta):=\log \mathbb E[e^{\beta \omega_{1,0}}]$ and $\sigma(\beta)^2:=e^{\lambda(2\beta)-2\lambda(\beta)}-1$ , our lower bound combined with Theorem 2.8 of Berger, Caravenna, and Turchi (2025) gives $-\mathfrak{f}(β) \asymp \exp{\Big(- \fracπ{σ^2(β)}\Big)},\quad \text{ as $\beta \downarrow 0$.}$

Gaussian Multiplicative Chaos for i.i.d. matrices

2026-05-28T10:02:26-04:00

We consider $N\times N$ matrices $X$ with independent, identically distributed entries, and prove that the sequence of measures $\frac{ | \det (X-z)|^\gamma}{\mathbb{E}[ | \det (X-z)|^\gamma]}$ converge to the Gaussian Multiplicative Chaos in the full subcritical regime $\gamma \in (0, 2 \sqrt{2})$ as $N \to \infty$ . 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 $K$ -point function of $| \det (X-z)|$ at any collection of mesoscopically separated points $z_i$ . Our methods are analytic and probabilistic in nature, relying in part on the dynamical approach based on Dyson Brownian motion.

On certain combinatorial expressions of TASEP transition probabilities

2026-05-28T03:38:48-04:00

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.

Cutoff profiles for conjugacy invariant random walks on symmetric groups

2026-05-27T13:27:47-04:00

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.

Non-stationary current fluctuations in 1D boundary-driven diffusive systems via Macroscopic Fluctuation Theory

2026-05-26T12:52:27-04:00

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.

Sum of the k Largest Eigenvalues of Symmetric Matrices: Theory and Applications

2026-05-26T04:47:25-04:00

This paper establishes new upper bounds for the sum of the $k$ 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 $k$ for almost all graphs, thereby taking a significant step toward its complete resolution.

Mixing times for Glauber dynamics of lozenge tilings of the hexagon

2026-05-26T03:53:29-04:00

We prove that the continuous-time, single-flip Glauber dynamics for lozenge tilings of the size- $N$ hexagon mix in time $N^{2+o(1)}$ . This was predicted to hold on fairly general domains of diameter $N$ (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).

Effect of slow bonds on current fluctuations in the symmetric simple exclusion process

2026-05-25T13:33:26-04:00

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.

Classification of the eternal solutions and multiple coalescing shocks in the KPZ fixed point

2026-05-25T13:10:48-04:00

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.

On the maximal dimension of an irreducible representation of the symmetric group

2026-05-25T12:14:41-04:00

We prove that the maximal dimension $d_N$ of an irreducible representation of the symmetric group $S_N$ satisfies $d_N=\sqrt{N!} \, e^{-(\mathfrak{d}+o(1))\sqrt{N} }, \quad N\to \infty,$ for some constant $\mathfrak{d}>0$ . This answers a question raised by Vershik--Kerov in 1985.

A Borodin-Okounkov-Geronimo-Case identity for tilted Toeplitz minors

2026-05-24T06:14:14-04:00

We prove a Fredholm determinantal identity for the tilted Toeplitz minor $D_{N}^{ξ,θ}(\varphi):= \det\bigl[(θ_{i}ξ_{j}\varphi)_{i-j}\bigr]_{i,j=1}^{N},$ generalizing the Borodin-Okounkov-Geronimo-Case (BOGC) identity to oblique splittings of the Hardy space. The tilts $\xi_{j},\theta_{i}$ enter only through an oblique projection that multiplies the trace-class kernel $K$ inside the Fredholm determinant; the BOGC operator $A=I-K$ constructed from $\varphi$ 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 $\theta_i=1$ ) admits a bialternant form recovering Schur and Grothendieck polynomials as special cases. A Cauchy-Binet expansion realizes $D_{N}^{\xi,\theta}$ 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.

Hurwitz numbers of a fixed partition (m, 1^{n-m}) via enumeration of unrooted hypermaps

2026-05-24T05:05:59-04:00

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.

Constructing Maximal Bumpless Pipedreams for Double Grothendieck Polynomials

2026-05-23T06:52:24-04:00

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 $\mathsf{rajcode}(\cdot)$ 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 $\mathsf{rajcode}(w)$ and column weight $\mathsf{rajcode}(w^{-1})$ . 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.

A Matrix Model for Higher-Genus Fuss--Catalan Numbers

2026-05-22T17:29:31-04:00

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.

Geodesics and Wandering Exponents in Brochette First-Passage Percolation

2026-05-22T10:59:24-04:00

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.

Invariant measures for half-space geometric LPP: classification and the one force--one solution principle

2026-05-21T13:42:24-04:00

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 $\infty$ , the recentered solution at time $0$ 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.

A random walk approach to high-dimensional critical phenomena

2026-05-20T13:27:43-04:00

We present a "black box" proof of mean-field near-critical behaviour for a family of functions on $\mathbb Z^d$ ( ${d>2}$ ) 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 $|x|^{-d+2+\varepsilon}\exp[-c|x|/\xi]$ , for any $\varepsilon>0$ . 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, $|\varphi|^4$ ), 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.

The Martin boundary of the Directed Landscape

2026-05-20T12:29:06-04:00

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.

The non-Hermitian minor process

2026-05-20T10:53:03-04:00

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.

The extra slow Tamari lattice

2026-05-20T04:44:26-04:00

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.

Eigenvalue Distribution of p-adic Random Matrices Among Algebraic Extensions, with an Analogue for p-adic Random Polynomials

2026-05-19T15:56:04-04:00

We study the distribution of eigenvalues of Haar-random matrices over $\mathbb{Z}_p$ among algebraic extensions of $\mathbb{Q}_p$ . Our results give $p$ -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, $p$ -adic eigenvalues are asymptotically evenly distributed among possible extension degrees. We also show that the maximal unramified extension $\mathbb{Q}_p^{\mathrm{un}}$ captures all but a bounded expected number of eigenvalues, and that the expected number of eigenvalues outside $\mathbb{Q}_p^{\mathrm{un}}$ 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 $\mathbb{Z}_p$ , using the correlation function formulas of Caruso (arXiv:2110.03942). These polynomial results are $p$ -adic analogues of the real-root counting results of Edelman-Kostlan, again with behavior different from the real setting.

Limit Profiles for Separation Distance

2026-05-18T16:17:46-04:00

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 $t_n+cw_n$ . 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.

Universal interface fluctuations in absorbing-state phase transitions

2026-05-17T22:59:54-04:00

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.

Coupled double Poisson brackets

2026-05-17T19:34:03-04:00

We introduce coupled double Poisson brackets on an associative algebra $A$ 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 $\operatorname{GL}_N$ -invariant Poisson structure on the representation scheme $\operatorname{Rep}_N(A)$ parametrizing $N$ -dimensional representations of $A$ , 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.

Multicritical Scaling Limit of Shifted Schur Measure

2026-05-15T04:54:12-04:00

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.

Hydrodynamic limits for TASEP with space-time discontinuities

2026-05-13T09:31:46-04:00

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.

Inhomogeneous q-Whittaker polynomials II: ring theorem and positive specializations

2026-05-13T08:26:20-04:00

We study inhomogeneous $q$ -Whittaker polynomials which extend both $q$ -Whittaker and stable Grothendieck polynomials. We prove that inhomogeneous $q$ -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 $q$ -Whittaker polynomials.

State-of-art minibatches via novel DPP kernels: discretization, wavelets, and rough objectives

2026-05-13T03:54:37-04:00

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.

Unimodality of q-Fibonomial coefficients for small cases

2026-05-12T19:44:53-04:00

Bergeron--Ceballos--Küstner introduced the $q$ -Fibonomial coefficients $ \qfibonom{m+n}{n}$, gave a combinatorial interpretation of the $q$ -Fibonomial coefficients via a weighted path-domino tiling model, and conjectured that these polynomials are unimodal. We prove the conjecture for $n\leq3$ . For the $n=2$ 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 $n=3$ case, we establish a more general unimodality result for certain products of $q$ -analogs and propose several related conjectures.

Central limit theorem for the homozygosity of the hierarchical Pitman-Yor process

2026-05-12T13:54:55-04:00

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.

The Geometry of Spectral Fluctuations: On Near-Optimal Conditions for Universal Gaussian CLTs, with Statistical Applications

2026-05-11T20:51:24-04:00

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.

Flocking with Multiple Types: Competition, Fluid Limits and Traveling Waves

2026-05-11T18:04:39-04:00

We study a class of interacting particle systems on $\mathbb{R}$ with two types. Particles evolve by independent jumps sampled from a fixed distribution, with type-dependent jump rates $v_+$ , $v_-$ 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 $v_+>v_-=0$ , 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.

Hydrodynamics and boundary-induced phase transitions in the n-species particle-exchange process

2026-05-08T07:43:27-04:00

The $n$ -species particle-exchange process (PEP( $n$ )) is an exclusion process in which particles of $n$ 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 $n$ coupled inviscid Burgers equations. This system of conservation laws is shown to admit Riemann invariants for arbitrary $n$ from which explicit solutions of the Riemann problem in terms of shock waves and rarefaction fans are obtained. We also introduce the open PEP( $n$ ), 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 $2n+1$ phases, with boundary-induced phase transitions analogous to those of the single-species asymmetric simple exclusion process.

Multilane Asymmetric Exclusion Process with stationary Bernoulli measure

2026-05-07T17:08:03-04:00

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.

\mathsf{GL}_N(\mathbb{C}) Brownian motion and stochastic PDE on entire functions

2026-05-07T11:34:40-04:00

We construct the full edge scaling limit of the singular values of Brownian motion on the general linear group $\mathsf{GL}_N(\mathbb{C})$ 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.

Asymmetric bilateral Bailey pairs and Rogers-Ramanujan type identities

2026-05-07T09:16:20-04:00

The theory of Bailey's transform provides a systematic method for deriving $q$ -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.