Leonid Petrov — arXiv Feed

A combinatorial formula for Wilson loop expectations on compact surfaces

2026-03-09T11:42:40-04:00

We give an almost purely combinatorial expression for Wilson loop expectations of the Yang-Mills holonomy process with values in the unitary group on a compact oriented surface, possibly with boundary and arbitrary boundary conditions. Our main result computes the non-normalized expectation of products of traces of holonomies along an arbitrary family of immersed curves with transverse self-intersections and no triple points. It is expressed as a sum over assignments of highest weights of the unitary group to the connected components of the complement of the curves. Each term is a product of a Gaussian exponential factor, dimensions of unitary representations, and local contributions at the intersection points given by the sine or cosine of an angle determined by the surrounding highest weights. As an application, we obtain a new and short proof of the Makeenko-Migdal equations on arbitrary compact surfaces.

Central Limits via Dilated Categories

2026-03-09T07:38:32-04:00

The Central Limit Theorem (CLT) establishes that sufficiently large sequences of independent and identically distributed random variables converge in probability to a normal distribution. This makes the CLT a fundamental building block of statistical reasoning and, by extension, in reasoning about computing systems that are based on statistical inference such as probabilistic programing languages, programs with optimisation, and machine learning components. However, there is no general theory of CLT-like results currently, which forces practitioners to redo proofs without having a good handle on the essential ingredients of CLT-type results. In this paper, we introduce dilated seminorm-enriched category theory as a unifying framework for central limits, and we establish an abstract central limit theorem within that framework. We illustrate how a strengthened version of the classical CLT and the law of large numbers can be obtained as instances of our framework. Moreover, we derive from our framework a novel central limit theorem for symplectic manifolds, the CLT for observables, which finds applications in statistical mechanics.

A generalization of Kadell's orthogonality ex-conjecture

2026-03-09T03:31:01-04:00

In 2000, Kadell gave an orthogonality conjecture for a symmetric function generalization of the Zeilberger--Bressoud $q$ -Dyson constant term identity. The non-zero part of Kadell's conjecture is a constant term identity indexed by a weak composition $v$ . This conjecture was first proved by Károlyi, Lascoux and Warnaar in 2015. They further formulated a closed-form expression for the above constant term when all parts of the composition $v$ are distinct. In 2021, Zhou obtained a recursion for this constant term for an arbitrary composition $v$ . In this paper, by categorizing the variables into two parts, we generalize Zhou's result.

A class of d-dimensional directed polymers in a Gaussian environment

2026-03-06T13:59:52-05:00

We introduce and analyze a broad class of continuous directed polymers in $\mathbb{R}^d$ driven by Gaussian environments that are white in time and spatially correlated, under Dalang's condition. Using an Itô-renormalized stochastic-heat-equation representation, we establish structural properties of the partition function, including positivity, stationarity, scaling, homogeneity, and a Chapman--Kolmogorov relation. On finite time intervals, we prove Brownian-type pathwise behavior, namely Hölder continuity and identification of the quadratic variation. We then obtain a sharp measure-theoretic dichotomy: the quenched polymer measure is singular with respect to Wiener measure if and only if $\widehat f(\mathbb{R}^d)=\infty$ (equivalently, the noise is non-trace-class), and it is equivalent otherwise. Finally, in dimension $d\ge 3$ , we prove diffusive behavior at large times in the high-temperature regime. This extends the Alberts--Khanin--Quastel framework from the $1+1$ white-noise setting to higher-dimensional Gaussian environments with general spatial covariance.

Massive holomorphicity of near-critical dimers and sine-Gordon model

2026-03-06T12:33:49-05:00

In this paper, we consider the near-critical dimer model in the setup of isoradial superpositions with Temperleyan boundary conditions. We show that the centered height function converges as the mesh size tends to zero to a limiting field which agrees with the (electromagnetically tilted) sine-Gordon model, whose derivative correlations are described by Grassmann variables (or equivalently determinants involving a massive Dirac operator). This answers a longstanding question in the field. A crucial part of the work is to develop a notion of discrete massive holomorphic functions and the tools to study such functions, in particular finding an exact discrete form of the massive Cauchy--Riemann equations, which is satisfied by the inverse Kasteleyn matrix. In comparison with previous studies, a key novelty of this part of our work is that the mass is not only allowed to be non-constant but can be complex-valued.

Gaussian free field convergence of the six-vertex model with $-1\leq\Delta\leq-\frac12$

2026-03-06T08:31:41-05:00

We study the isotropic six-vertex model on $\mathbb{Z}^2$ with spectral parameter $\Delta\in[-1,-1/2]$ , that is, with weights $\mathbf{a}=\mathbf{b}=1$ and $\mathbf{c}\in[\sqrt{3},2]$ . We show that the associated height function converges, in the scaling limit, to a properly scaled full-plane Gaussian free field. The result extends to anisotropic weights $\mathbf{a}\neq\mathbf{b}$ upon using a suitable embedding of the lattice.

Fluctuations for the Sherrington--Kirkpatrick spin glass model near the critical temperature

2026-03-05T14:44:45-05:00

We consider the Sherrington--Kirkpatrick spin glass model with zero external field and at inverse temperature $\beta>0$ . Let $F_N(\beta)$ be the corresponding log-partition function. Under the assumption that $c_N:=N^{1/3}(1-\beta_N^2)$ is bounded away from $0$ , we prove that Var $(F_N(\beta_N)) = - \frac{1}{2} \log (1-\beta_N^2) -{\beta_N^2}/{2} + O( c_N^{-3/2}).$ As a consequence, we obtain Var $(F_N(1-c N^{-1/3})) = \frac16\log N + O(1)$ for any fixed constant $c\in(0,\infty)$ . We also prove a Gaussian central limit theorem for the centered and scaled $F_N(\beta_N)$ .

BBP Phase Transition for a Doubly Sparse Deformed Model

2026-03-05T00:30:43-05:00

We prove the equivalent of the Baik, Ben Arous, Péché (2004) phenomenon for a novel, doubly sparse model where both the Wigner noise matrix and signal vector(s) are sparse. Specifically, we consider a deformed sub-Gaussian sparse Wigner ensemble with a fixed number of sub-Gaussian spike vectors of the same-order sparsity added. We show that spike vectors with signals greater than one are correlated with the top eigenvectors of the deformed ensemble and that each spike vector of signal greater than one induces an outlier eigenvalue. Notably, our results hold in the supercritical sparsity regime for the Wigner matrix ( $q \gg \frac{\log n}{n}$ ) and for any sparse spike vector with an unbounded number of entries ( $np\to \infty$ ). No further relationship between the sparsities of the noise matrix ( $q$ ) and spike vectors ( $p$ ) is necessary. This generalizes the work of Benaych-Georges and Nadakuditi (2010) and Péché (2005).

Comparison theorems for the extreme eigenvalues of a random symmetric matrix

2026-03-04T13:30:04-05:00

This paper establishes a comparison theorem for the maximum eigenvalue of a sum of independent random symmetric matrices. The theorem states that the maximum eigenvalue of the matrix sum is dominated by the maximum eigenvalue of a Gaussian random matrix that inherits its statistics from the sum, and it strengthens previous results of this type. Corollaries address the minimum eigenvalue and the spectral norm. The comparison methodology is powerful because of the vast arsenal of tools for treating Gaussian random matrices. As applications, the paper improves on existing eigenvalue bounds for random matrices arising in spectral graph theory, quantum information theory, high-dimensional statistics, and numerical linear algebra. In particular, these techniques deliver the first complete proof that a sparse random dimension reduction map has the injectivity properties conjectured by Nelson & Nguyen in 2013.

Log Log Fluctuations of the Stochastic Heat Flow

2026-03-03T13:38:41-05:00

We study the stochastic heat flow with constant initial data and analyze its spatial average on the scale of $\varepsilon\ll1$ . We prove that the logarithm of the averaged process satisfies a pointwise central limit theorem: After being centered by $-\tfrac{1}{2}\log\log \varepsilon^{-1}$ and scaled down by $\sqrt{\log\log \varepsilon^{-1}}$ , it converges in distribution to a standard Gaussian.

Periodic KPZ fixed point with general initial conditions

2026-03-02T10:14:27-05:00

We consider the periodic totally asymmetric simple exclusion process with a general initial condition that properly approximates a periodic upper-semicontinuous function. We find the large time limit of the rescaled space-time multipoint distribution of the height function in the relaxation time scale. The limiting functions form a consistent family of finite-dimensional distributions; thus, they define the periodic KPZ fixed point with a general upper-semicontinuous initial condition. The main technical novelty of the paper is a hitting expectation representation of the energy function and the characteristic function in the finite-time multipoint distribution formula obtained in arXiv:1912.10143. The representation of the characteristic function is partly inspired by the work of arXiv:1701.00018, arXiv:2509.03246, while the representation of the energy function is based on an extensive guess-and-check exploration.

Condensation in stochastic lattice gases with size-dependent stationary weights

2026-02-28T15:57:45-05:00

We consider stochastic lattice gases with stationary product weights and a polynomial perturbation vanishing with the system size that leads to condensation. If the density of particles exceeds a critical value the system phase separates into a bulk with homogeneous distribution of particles and a condensed phase. Depending on parameter values, the latter consists of a single macroscopic cluster or a diverging number of independent clusters on a smaller scale. We establish the condensation transition via the equivalence of ensembles and the main novelty is a derivation of the cluster size distribution using size-biased sampling, generalizing previous work on zero-range and inclusion processes. Simulations of zero-range processes illustrate our theoretical results on the condensate scale and size distribution.

The Airy line ensemble at the rough-smooth boundary

2026-02-27T09:52:09-05:00

We study the rough-smooth boundary in the two-periodic Aztec diamond, a random domino tiling model exhibiting three types of macroscopic regions. We show that the height function at this boundary converges to an independent sum of an Airy surface and an i.i.d. noise field with fluctuations governed by the full-plane smooth phase. Going further, we prove convergence of a family of Temperleyan backbone paths to the Airy line ensemble. This gives the first convergence result for a family of undirected paths converging to the Airy line ensemble, as well as Airy convergence at a noisy boundary.

Exact Anomalous Current Fluctuations in Quantum Many-Body Dynamics

2026-02-27T08:32:17-05:00

Fluctuations of integrated currents have attracted considerable interest over the past decades in the context of statistical mechanics. Recently, anomalous current fluctuations, characterized by the M-Wright function, were obtained exactly in a classical automaton [ $Ž$ . Krajnik et al., Phys. Rev. Lett. 128, 160601 (2022)], and previous studies have shown that the anomalous behavior can arise in a variety of classical systems. Despite the rapidly growing interest in such anomalous behaviors, which capture a universal aspect of one-dimensional many-body transport, the exact derivation of the M-Wright function in quantum many-body systems has remained elusive. In this Letter, we present the first exact microscopic derivation of the M-Wright function in quantum many-body dynamics by analyzing the integrated spin current in a one-dimensional Fermi-Hubbard model with infinitely strong repulsive interactions. Our results lay the groundwork for exploring anomalous integrated currents in a broad class of quantum many-body systems.

An exactly solvable evaporation-deposition PCA with long-distance interactions

2026-02-27T04:18:54-05:00

We consider a probabilistic cellular automaton (PCA) of evaporation-deposition on the one-dimensional lattice having $n$ sites with periodic boundary conditions, in which each site, during each epoch, can be in one of two states: $0$ and $1$ . Fix a positive integer $m\geqslant 2$ . There are two types of transitions at each discrete time, which are as follows: (i) the first site in every contiguous block of $m$ $0$ s becomes a $1$ with probability $p_1$ , and (ii) the first site in every contiguous block of $(m-1)$ $0$ s followed immediately by a $1$ also becomes a $1$ with probability $(1-p_2)$ . As in a PCA, all of these transitions occur simultaneously. We show that the resulting discrete-time Markov chain is ergodic, and we give an explicit formula for its limiting distribution, the partition function and the density. We also propose necessary and sufficient conditions for this Markov chain to be reversible. For $m=2$ , we provide a fully analytical expression for the free energy of this model.

Pfaffian point processes for coalescing particles via checkerboard duality

2026-02-26T06:23:29-05:00

Coalescing particles on a line merge when they meet. When every site is initially occupied, only finitely many particles survive at any positive time, and their positions form a Pfaffian point process: all correlation functions are determined by pairwise quantities arranged in antisymmetric matrices. Previous proofs of this structure relied on analytic methods specific to time-homogeneous dynamics. We identify the checkerboard duality as the structural reason for the Pfaffian: on a discrete planar graph, binary random choices create two complementary non-crossing forests, one tracing ancestral lineages backward and the other carrying the coalescing particles forward as domain boundaries. This duality converts the absence of particles in an interval into coalescence of ancestral lineages at its endpoints. A cancellative labeling then converts coalescence into annihilation, for which a companion paper provides a Pfaffian formula. The resulting empty-interval formula holds for any discrete graph with the checkerboard structure and arbitrary inhomogeneous edge probabilities, requiring no symmetry or specific distributions. This covers settings beyond existing methods, including totally asymmetric dynamics and position-dependent transition rules, and yields an explicit Pfaffian point process in each setting. For Brownian motion, the formula recovers the known Pfaffian point process and the empty-interval probabilities previously derived by PDE methods.

Numerical ranges of non-normal random matrices: elliptic Ginibre and non-Hermitian Wishart ensembles

2026-02-26T00:00:00-05:00

The numerical range of a non-normal matrix plays a central role as a descriptor of non-normal effects beyond spectral information. We study a class of fundamental non-Hermitian random matrix ensembles that interpolate between the Hermitian and non-Hermitian regimes. Our analysis focuses on the elliptic Ginibre ensemble and its chiral counterpart, as well as on non-Hermitian Wishart matrices. For each of these models, we explicitly characterise the geometry of the numerical range in the large-system limit. In particular, we show that for the elliptic Ginibre ensemble and its chiral version, the limiting numerical range is an ellipse, whereas for the non-Hermitian Wishart ensemble it is described by a non-elliptic envelope.

Multiple Poisson-Dirichlet diffusions on generalized Kingman simplices

2026-02-26T00:00:00-05:00

We construct a new class of infinite-dimensional diffusions taking values in a generalized Kingman simplex. Our model describes the temporal evolution of the relative frequencies of infinitely-many types which are "labeled" by an arbitrary finite number of marks or colors, but "unlabeled" within each mark. We start with a finite-dimensional construction which extends to Wright-Fisher diffusions a self-similarity property known for Dirichlet distributions, and corresponds to a multiple skew-product representation of the Wright-Fisher diffusion relative to the marks in the population. After ranking decreasingly the frequencies within each mark, we identify the limit in distribution of the resulting diffusion when the number of types for each mark goes to infinity, and describe its infinitesimal operator. The limiting process reduces to a diffusion in the Thoma simplex in the special case of only two marks, whereas the infinitely-many-neutral-alleles model is recovered when all frequencies have the same mark. The stationary measure of the limit diffusion is shown to be the recently introduced multiple Poisson-Dirichlet distribution, which extends Kingman's Poisson-Dirichlet distribution and is the de Finetti representing measure for a family of random partitions whose elements are marked.

Convergence of Half-Space Last Passage Percolation Away from the Boundary to the Directed Landscape

2026-02-26T00:00:00-05:00

In this note, we prove convergence of the half-space exponential last passage percolation (LPP) model, away from the boundary, to the directed landscape. Our approach couples the half-space and full-space LPP models and constructs two barrier events based on the monotonicity of last passage paths. Combining this coupling with moderate deviation estimates for both models and the known convergence of full-space LPP to the directed landscape, we establish the desired convergence.

The Quantum Symmetric Simple Exclusion Process in the Continuum and Free Processes

2026-02-26T00:00:00-05:00

The quantum symmetric simple exclusion process (QSSEP) is a recent extension of the symmetric simple exclusion process, designed to model quantum coherent fluctuating effects in noisy diffusive systems. It models stochastic nearest-neighbor fermionic hopping on a lattice, possibly driven out-of-equilibrium by boundary processes. We present a direct formulation in the continuum, and establish how this formulation captures the scaling limit of the discrete version. In the continuum, QSSEP emerges as a non-commutative process, driven by free increments, conditioned on the algebra of functions on the ambiant space to encode spatial correlations. We actually develop a more general framework dealing with conditioned orbits with free increments which may find applications beyond the present context. We view this construction as a preliminary step toward formulating a quantum extension of the macroscopic fluctuation theory.

Scaling limits for some Mittag-Leffler queues

2026-02-26T00:00:00-05:00

In this paper, we consider five models of heavy-tailed queues involving Mittag-Leffler distributions that generalize the classical $M/M/1$ queues. These models are suitable modifications of previously defined models in such a way that the classical $M/M/1$ queue can be recovered by a suitable selection of parameters. We provide the distribution of inter-arrival and service times of both the original and modified queueing models. We then study the scaling limits of all the proposed models and we argue that the behaviour of the limiting processes can be used to characterise the traffic regime of the queues.

Martin Boundary and the Nonlinear Multiplicative Stochastic Heat Equation in Weak Disorder

2026-02-26T00:00:00-05:00

We study the nonlinear multiplicative stochastic heat equation on Dirichlet spaces with white in time noise under weak disorder. We show that positive invariant random fields with uniformly bounded second moments are in one-to-one correspondence with bounded positive harmonic functions. The proof combines a remote past pullback construction with a uniqueness argument that applies a contraction inspired by chaos expansion. As a consequence, the space of invariant measures inherits geometric structure from the Martin Boundary. We further establish a small-noise Gaussian fluctuation result within each harmonic sector and show that the long-time behavior of solutions is completely determined by the Martin boundary data of the initial condition. These results reveal a direct connection between stochastic PDE dynamics and boundary theory in potential analysis.

Asymptotics and zero distribution of geometric polynomials

2026-02-26T00:00:00-05:00

We obtain some results on the asymptotic behavior and zero distribution of the so-called geometric polynomials. The asymptotics is given both on compact subsets of $\C\setminus [-1,0]$ and on compact subsets of the interval $(-1,0)$ . The zeros of these polynomials are simple and lie in $(-1,0]$ ; moreover, the zeros of consecutive polynomials interlace. Its zero distribution is a measure whose density is similar to Cauchy weight. Some orthogonality properties of these polynomials are also proved.

A Phase Transition For Repeated K-Averages

2026-02-26T00:00:00-05:00

Let $x_1,\dots,x_{n}$ be a fixed sequence of real numbers. At each stage, pick $k$ integers $\{I_{i}\}_{1\leq i \leq k}$ uniformly at random without replacement and then for each $i \in \{1,2,\dots,k\}$ replace $x_{I_i}$ by $(x_{I_1}+x_{I_2}+\dots+x_{I_k})/k$ . It is easy to observe that all the co-ordinates converge to $(x_1+\dots+x_n)/n$ . In this article, we extend the result of \cite{chatterjee2019note} by establishing order of decay of the expected $L^{2}$ distance. Furthermore, we establish the mixing time to be in between $\frac{n}{k \log k}\log n$ and $\frac{n}{k-1}\log n$ .

A Characterization of the Macdonald Hypergeometric Series {}_r\Phi_s(x;q,t) and {}_r\Phi_s(x,y;q,t) via q-Difference Equations

2026-02-26T00:00:00-05:00

In two widely circulated manuscripts from the 1980s, I. G. Macdonald introduced certain multivariate hypergeometric series ${}_pF_q(x;\alpha)$ and ${}_pF_q(x,y;\alpha)$ and their $q$ -analogs ${}_r\Phi_s(x;q,t)$ and ${}_r\Phi_s(x,y;q,t)$ . These series are given by explicit expansions in Jack and Macdonald polynomials, and they generalize the hypergeometric functions of one and two matrix arguments from statistics. In a recent joint paper with Siddhartha Sahi, we constructed differential operators that characterize the Jack series ${}_pF_q$ thereby answering a question of Macdonald. In this paper we construct analogous $q$ -difference operators that characterize the Macdonald series ${}_r\Phi_s$ . More precisely, we construct three $q$ -difference operators $\mathcal A^{(x,y)}$ , $\mathcal B^{(x)}$ , $C^{(x)}$ . The equation $\mathcal A^{(x,y)}(f(x,y))=0$ characterizes ${}_r\Phi_s(x,y;q,t)$ , while the equations $\mathcal B^{(x)}(f(x))=0$ and $\mathcal C^{(x)}(f(x))=0$ each characterize ${}_r\Phi_s(x;q,t)$ . These characterizations are subject to certain symmetry, boundary and stability condition. In the special case of ${}_2\Phi_1(x;q,t)$ , our operator $\mathcal B^{(x)}$ was previously constructed by Kaneko in 1996.

Determinant and Pfaffian formulas for particle annihilation

2026-02-26T00:00:00-05:00

When particles on a line collide, they may annihilate-both are destroyed. Computing exact annihilation probabilities has been difficult because collisions reduce the particle count, while determinantal methods require a fixed count throughout. The ghost particle method, introduced in a companion paper for coalescence, keeps destroyed particles walking as invisible ghosts that restore the missing dimension. We apply this method to annihilation: when two particles annihilate, both trajectories continue as invisible walkers, yielding an exact determinantal formula that specifies the number of annihilations, where survivors end up, and where ghosts end up. For complete annihilation (no survivors), the formula simplifies to a Pfaffian-an algebraic relative of the determinant built from pairwise quantities-connecting to Pfaffian point process theory. The annihilation formula also yields results about coalescence: pairwise coalescence can be reinterpreted as complete annihilation, producing a Pfaffian coalescence formula. These formulas are exact for any finite initial configuration and apply to discrete lattice paths, birth-death chains, and continuous diffusions including Brownian motion.

On Card guessing after a single shelf shuffle

2026-02-26T00:00:00-05:00

We consider a card guessing game with complete feedback. An ordered deck of $n$ cards labeled $1$ up to $n$ is shelf-shuffled exactly one time. One after the other a single card is drawn from the shuffled deck. The guesser makes has guess and the card is shown until no cards remain. We provide a distributional analysis of the number of correct guesses under the optimal strategy. We re-obtain the previously derived expectation and add a complete description of the distribution. We also obtain a central limit theorem for the number $n$ of cards tending to infinity. Furthermore, we discuss an unbalanced, biased shelf shuffle and show how to derive the extend our analysis, also adding the complete position matrix. Finally, a refined analysis of the number of correct guesses is carried out, distinguishing between pure luck guesses and certified correct guesses.

Combinatorial formulas for symmetric Macdonald polynomials by superizations

2026-02-26T00:00:00-05:00

In this paper, we derive new combinatorial formulas for symmetric Macdonald polynomials $P_{\lambda}(X;q,t)$ and integral Macdonald polynomials $J_{\lambda}(X;q,t)$ , in terms of several new statistics and the major index for a partition $\lambda$ . Compared to previous formulas, these new ones contain the fewest terms. Moreover, three existing formulas for symmetric Macdonald polynomials established by Corteel--Mandelshtam--Williams (2022), Corteel--Haglund--Mandelshtam--Mason--Williams (2022) and Mandelshtam (2025) are recovered. Our proof relies on a new statistic on super fillings, employing the superization formulas of Haglund--Haimain--Loehr (2005) and Ayyer--Mandelshtam--Martin (2023), together with our recent approach to modified Macdonald polynomials.

Integral formula for the propagator of the one-dimensional Hubbard model

2026-02-24T22:56:10-05:00

We present an exact integral formula for the multi-particle propagator of the one-dimensional Fermi--Hubbard model on an infinite lattice. The proof is based on the nested Bethe ansatz without relying on the string hypothesis. Our formula enables an explicit integral representation of the time evolution of arbitrary finite-particle wave functions and thereby provides a foundation for the exact analysis of nonequilibrium dynamics in the Hubbard model. It can further be applied to related open quantum models.

Numerical ranges of non-normal random matrices: elliptic Ginibre and non-Hermitian Wishart ensembles

2026-02-24T00:00:00-05:00

Multiple Poisson-Dirichlet diffusions on generalized Kingman simplices

2026-02-23T00:00:00-05:00

Conditional thinning and multiplicative statistics of Laguerre-type orthogonal polynomial ensembles

2026-02-20T05:06:31-05:00

We study the local statistics of orthogonal polynomial ensembles near a hard edge, subject to a multiplicative deformation of the measure. Probabilistically, this deformation corresponds to a position-dependent conditional thinning of the particles. We prove that, under critical hard edge scaling and for a large class of potentials and deformation symbols, the correlation kernel of the conditional ensemble converges to a universal limit, which we identify as the conditional thinned Bessel point process. We derive an explicit expression for this limiting kernel in terms of the solution to a nonlocal integrable system depending on a parameter. For a special choice of the parameter, this system was recently identified in the study of multiplicative statistics of the Bessel point process. Our results establish that this system governs the full correlation structure of the conditional Bessel point process, extending the classical connection between the standard Bessel kernel and the Painlev\'e V equation.

Convergence of Half-Space Last Passage Percolation Away from the Boundary to the Directed Landscape

2026-02-20T00:00:00-05:00

Multispecies inhomogeneous t-PushTASEP with general capacity

2026-02-19T03:52:26-05:00

We study an $n$ -species $t$ -PushTASEP, an integrable long-range stochastic process, on a one-dimensional periodic lattice with inhomogeneities $x_1,\ldots,x_L$ and arbitrary capacity $l$ at each lattice site. The Markov matrix is identified with an alternating sum of commuting transfer matrices over all fundamental representations of $U_t(\widehat{sl}_{n+1})$ . Stationary probabilities are expressed in a matrix product form involving a fusion of quantized corner transfer matrices for the strange five-vertex model introduced by Okado, Scrimshaw, and the second author. The resulting partition function, which serves as the normalization factor of the stationary probabilities, is obtained from the $l=1$ case by a finite plethystic substitution of length $l$ .

Entrance laws for coalescing and annihilating Brownian motions

2026-02-18T09:55:43-05:00

Systems of instantaneously annihilating or coalescing Brownian motions on the line are considered. The extreme points of the set of entrance laws for this process are shown to be Pfaffian point processes at all times and their kernels are identified.

The Quantum Symmetric Simple Exclusion Process in the Continuum and Free Processes

2026-02-18T00:00:00-05:00

Scaling limits for some Mittag-Leffler queues

2026-02-18T00:00:00-05:00

Martin Boundary and the Nonlinear Multiplicative Stochastic Heat Equation in Weak Disorder

2026-02-18T00:00:00-05:00

Asymptotics and zero distribution of geometric polynomials

2026-02-17T00:00:00-05:00

A Phase Transition For Repeated K-Averages

2026-02-17T00:00:00-05:00

Estimating π with a Coin

2026-02-16T01:01:32-05:00

We describe a simple Monte Carlo method for estimating $\pi$ by tossing a coin. Although the underlying Catalan-number series identities appear implicitly in the probability theory literature, the probabilistic interpretation of $\frac\pi{4}$ presented here seems to be new.

A probabilistic interpretation for interpolation Macdonald polynomials

2026-02-13T16:55:56-05:00

Previous work of Ayyer, Martin, and Williams gave a probabilistic interpretation of the Macdonald polynomials $P_{\lambda}(x_1,\dots,x_n;1,t)$ at $q=1$ in terms of a Markov chain called the multispecies $t$ -Push TASEP, a Markov chain involving particles of types $\lambda_1,\dots,\lambda_n$ hopping around a ring. In particular, they showed that for each composition $\eta$ obtained by permuting the parts of $\lambda$ , the stationary probability of being in state $\eta$ is proportional to the ASEP polynomial $F_{\eta}(x_1,\dots,x_n; 1,t)$ , and the normalizing constant (or partition function) is $P_{\lambda}(x_1,\dots,x_n; 1,t)$ . There is an inhomogeneous generalization of Macdonald polynomials due to Knop and Sahi called interpolation Macdonald polynomials $P^*_{\lambda}(x_1,\dots,x_n;q,t)$ , as well as an inhomogeneous generalization of ASEP polynomials called interpolation ASEP polynomials $F^*_{\eta}(x_1,\dots,x_n;q,t)$ that we introduced in previous work. In this article we introduce a new Markov chain called the interpolation $t$ -Push TASEP, and show that its steady state probabilities and partition function are given by the interpolation ASEP polynomials and the interpolation Macdonald polynomial, evaluated at $q=1$ . This generalizes the previous result of Ayyer, Martin, and Williams.

A Characterization of the Macdonald Hypergeometric Series {}_r\\Phi_s(x;q,t) and {}_r\\Phi_s(x,y;q,t) via q-Difference Equations

2026-02-13T00:00:00-05:00

Determinant and Pfaffian formulas for particle annihilation

2026-02-13T00:00:00-05:00

On Card guessing after a single shelf shuffle

2026-02-13T00:00:00-05:00

Combinatorial formulas for symmetric Macdonald polynomials by superizations

2026-02-13T00:00:00-05:00

Factorised stationary states for a long range misanthrope process

2026-02-11T08:44:15-05:00

The misanthrope process is an interacting particle system where particles move between neighbouring sites with hop rates depending only on the number of particles at the departure and arrival sites. Motivated by a discretised version of the Hammersley--Aldous--Diaconis process, we introduce a partially asymmetric long range misanthrope process (PALRMP) on a finite one-dimensional lattice with periodic boundary conditions where particles can move between sites that are not necessarily neighbours, as long as there are no particles in between the departure and arrival sites. In this model, each site $\ell$ has an inhomogeneous rate parameter $x_\ell$ associated to it, and the hop rate of a particle moving from site $k$ to site $\ell$ depends upon the parameter associated to the target site $x_\ell$ , the direction the particle moves, and the number of particles at sites $k$ and $\ell$ . We also consider the homogeneous PALRMP, where all the $x_\ell$ 's are 1. We find necessary and sufficient conditions on the hop rates under which the stationary distribution is of factorised form for both the PALRMP and the homogeneous PALRMP, as well as the extreme variants, namely the ones where the particle motion is totally asymmetric (TALRMP) and symmetric (SLRMP). As an illustrative example, we study in detail the discrete Hammersley--Aldous--Diaconis process.

Geodesic networks and the disjointness gap in the directed landscape

2026-02-10T10:03:23-05:00

The directed landscape is a random directed metric on the plane that arises as the scaling limit of metric models in the KPZ universality class. For a pair of points p, q, the disjointness gap G(p; q) measures the shortfall when we optimize length over pairs of disjoint paths from p to q versus optimizing over all pairs of paths. Any spatial marginal of G is simply the gap between the top two lines in an Airy line ensemble. In this paper, we show that when the start and end time are fixed, the disjointness gap fully encodes the set of exceptional geodesic networks. The correspondence uses simple features of the disjointness gap, e.g. zeroes, local minima. We give a similar correspondence relating semi-infinite geodesic networks to a Busemann gap function. The proofs are deterministic given a list of soft properties related to the coalescent geometry of the directed landscape.

Moments of CβE field partition function, \\mathsf{Sine}_β correlations and stochastic zeta

2026-02-09T09:41:40-05:00

We prove a conjecture of Fyodorov and Keating on the supercritical moments of the partition function of the C $\beta$ E field or equivalently the supercritical moments of moments of the characteristic polynomial of the C $\beta$ E ensemble for general $\beta>0$ and general real moment exponents. Moreover, we give the first expression for all correlation functions of the $\mathsf{Sine}_\beta$ point process for all $\beta>0$ . The main object behind both results is the Hua-Pickrell stochastic zeta function introduced by Li and Valk\'{o}.

Bulk heights of the KPZ line ensemble

2026-02-06T12:00:19-05:00

For $t > 0$ , let $\{\mathcal{H}^{(t)}_n, n \in \mathbb{N}\}$ be the KPZ line ensemble with parameter $t$ , satisfying the homogeneous $\mathbf{H}$ -Brownian Gibbs property with $\mathbf{H}(x) =e^x$ . We prove quantitative concentration estimates for the $n$ th line $\mathcal{H}^{(t)}_n$ which yield the asymptotics $\mathcal{H}^{(t)}_n = n \log n + o(n^{3/4 + \epsilon})$ as $n \to \infty$ . A key step in the proof is a general integration by parts formula for $\mathbf{H}$ -Brownian Gibbs line ensembles which yields the identity $\mathbb{E} \exp(\mathcal{H}^{(t)}_{n + 1}(x) - \mathcal{H}^{(t)}_n (x)) = n t^{-1}$ for any $n, t, x$ .