Leonid Petrov — arXiv Feed

Structure preserving properties of higher order moment closures for TASEP

2026-04-17T06:29:03-04:00

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.

Kardar-Parisi-Zhang physics in optically-confined continuous polariton condensates

2026-04-16T10:54:29-04:00

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.

Fluctuations for the Toda lattice

2026-04-15T14:59:53-04:00

In this paper we consider the Toda lattice $(\mathbf{p}(t);\mathbf{q}(t))$ at thermal equilibrium, meaning that its variables $(p_j)$ and $(e^{q_j - q_{j+1}})$ 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 $q_0$ 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.

Shocks, instability, and the twenty networks of infinite geodesics in the Directed Landscape

2026-04-14T12:57:03-04:00

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.

On the large genus of Hurwitz numbers

2026-04-14T11:31:06-04:00

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.

From Gaussian to Gumbel: extreme eigenvalues of complex Ginibre products with exact rates

2026-04-14T00:59:07-04:00

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.

An upper bound on the smallest singular value of dense random combinatorial matrices

2026-04-13T23:12:31-04:00

Let $M$ be an $n\times n$ random matrix with entries in $\{0, 1\}$ , where each row is independently and uniformly sampled from the set of all vectors in $\{0, 1\}^n$ containing exactly $d$ ones, with $d=pn$ for some fixed constant $p\in (0,1/2]$ . A recent result of Tran states that the smallest singular value $s_n(M)$ is bounded below by $c_p n^{-1/2}$ with high probability. In this note, we establish a complementary upper bound for $s_n(M)$ , 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 $C_p$ is a positive constant depending only on $p$ . This result confirms that the least singular value $s_n(M)$ of dense random combinatorial matrices is typically of the order $n^{-1/2}$ .

The smallest singular value of signed random combinatorial matrices

2026-04-13T13:34:04-04:00

Let $M_n$ be an $n\times n$ signed random combinatorial matrix whose rows are independent and uniformly distributed over the set of $\{-1,0,1\}$ -vectors with exactly $n/2$ zero coordinates. Despite the dependence induced by the row constraints, we prove that there exist constants $C,c > 0$ such that for any $\varepsilon\ge0$ , \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 $M_n$ 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.

Semilocalization for inhomogeneous random graphs

2026-04-13T12:27:34-04:00

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.

Generating Hadamard matrices with transformers

2026-04-13T03:18:49-04:00

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 $100$ and $250$ , 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 $244$ . In addition to these new constructions, our experiments reveal that the transformer can discover and exploit useful hidden symmetry in the search space.

The circular law for sparse random combinatorial matrices

2026-04-11T23:59:07-04:00

Let $\log^{2+\varepsilon} n \le d \le n/2$ for some fixed $\varepsilon \in (0,1)$ , and let $M_n$ be an $n\times n$ random matrix with entries in ${0,1}$ , where each row is independently and uniformly sampled from the set of all vectors in ${0,1}^n$ containing exactly $d$ ones. We show that the empirical spectral distribution of the appropriately rescaled matrix $M_n$ converges in probability to the circular law provided that $d=o(n)$ . As a crucial element of the proof, we obtain quantitative lower bounds on the smallest singular value of the shifted matrices $M_n-zI_n$ whenever $|z|\le \sqrt d \log\log d$ and $C\log n \le d \le n/2$ for some absolute positive constant $C$ .

The directed landscape in half-space

2026-04-11T00:21:33-04:00

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.

Banded Hermitian Matrices, Matrix Orthogonal Polynomials, and the Toda Lattice

2026-04-10T16:22:16-04:00

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.

Classification of irreducible real modules of real Lie superalgebras

2026-04-10T11:16:28-04:00

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 $\mathbf Q(n)$ , 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.

Graded Casimir elements and central extensions of color Lie algebras

2026-04-09T23:08:29-04:00

A color Lie algebra is a generalization of a Lie (super)algebra by an Abelian group $\Gamma$ . The underlying vector space and defining relations of the algebra are graded by $\Gamma$ , 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, $\mathfrak{sl}(2)$ for $\Gamma= \mathbb{Z}_3^2$ , and $\mathfrak{q}(n)$ and $\mathfrak{osp}(m|2n)$ for $\Gamma= \mathbb{Z}_2^2$ .

Stability for the stochastic heat equation with multiplicative noise via finite-dimensional feedback

2026-04-09T14:15:50-04:00

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.

Convergence of k-point functions in high dimensional percolation

2026-04-09T12:58:03-04:00

Consider critical Bernoulli percolation on $\mathbb{Z}^d$ for $d$ large; let $y_0, \dots, y_{k-1}$ be $k$ distinct points in $\mathbb{R}^d$ . We prove that the probability that $\{\lfloor n y_i\rfloor\}_{i=0}^{k-1}$ all lie in the same open cluster, rescaled by an appropriate power of $n$ , converges as $n \to \infty$ to an explicit constant. This confirms a conjecture of Aizenman and Newman.

Ehrhart positivity for marked order polytopes

2026-04-09T11:55:15-04:00

Given a pair of finite posets $A \subseteq P$ , the function counting integer-valued order preserving extensions of an order preserving map $\lambda: A\rightarrow \mathbb{Z}$ is given by a piecewise polynomial in $\lambda$ . 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 $m$ -generalized Pitman-Stanley polytopes due to Alexanderson-Alhajjar and Dugan-Hegarty-Morales-Raymond, respectively.

Hydrodynamic limit of the directed exclusion process

2026-04-09T08:13:11-04:00

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.

Stochastic fractional heat equation with general rough noise

2026-04-08T21:37:31-04:00

Consider the following nonlinear one-dimensional stochastic fractional heat equation $\frac{\partial }{\partial t}u(t, x)= -(-Δ)^{α/2}u(t, x) +σ(t,x,u(t,x)) \dot{W}(t, x),$ where $-(-\Delta)^{\alpha/2}$ is the fractional Laplacian on $\mathbb R$ for $1<\alpha<2$ , and $\dot{W}$ is a Gaussian noise that is white in time and behaves in space as a fractional Brownian motion with Hurst index $H$ satisfying $\frac{3-\alpha}{4}<H<\frac12$ . When $\alpha=2$ , 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 $\sigma(0)=0$ 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 $\alpha \in (1,2)$ , inspired by Hu and Wang, we investigate the well-posedness of the stochastic fractional heat equation without imposing the technical condition of $\sigma(0)=0$ , 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 $-(-\Delta)^{\alpha/2}$ play a crucial role.

Random permutations from q-Demazure products

2026-04-07T20:05:14-04:00

We study the $q$ -deformation of the Demazure product model from arXiv:2407.21653. Consider the longest element $w_0$ in $S_n$ written as a reduced word in simple transpositions. Independently delete each transposition with probability $1-p$ and apply the $q$ -Demazure product to the remaining ones. We show that the law of the resulting permutation converges as $n \to \infty$ to a deterministic permuton, which coincides with the $q=0$ case studied in arXiv:2407.21653 for adjusted probability $p'=p(1-q)/(1-qp)$ . This resolves Conjecture 1.13 from arXiv:2407.21653 and identifies the limiting permuton explicitly.

Dimensional crossover in surface growth on rectangular substrates

2026-04-07T17:38:59-04:00

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 $L_y > L_x$ . 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 $\mathcal{R}=L_y/L_x$ , the roughness $W$ scales with time $t$ in the growth regime as $W \sim t^{\beta_{\text{2D}}}$ for $t \ll t_c$ and $W \sim t^{\beta_{\text{1D}}}$ for $t \gg t_c$ , where $t_c \sim L_x^{z_{2\text{D}}}$ 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 ( $t \ll t_c$ ) 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 $\mathcal{R}$ increases. The particular case $L_x = L_y^\delta$ , with $0 < \delta< 1$ , is also discussed in detail and reveals interesting features of the investigated systems. For instance, there exist a `special' exponent $\delta^* = z_{1\text{D}}/z_{2\text{D}}$ such that the temporal crossover cannot be observed for $\delta> \delta^*$ . Moreover, this leads the saturation roughness to display a nonuniversal scaling: $W_s \sim L_y^\Lambda$ , with $\Lambda= (1-\delta) \alpha_{1\text{D}} + \delta \alpha_{2\text{D}}$ .

On the quasisymmetric functions in superspace

2026-04-07T16:14:40-04:00

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

Asymptotic expansions of integrals and Nielsen's polylogarithms

2026-04-07T09:58:57-04:00

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 $n\rightarrow\infty$ , with parameters real $w\neq 0$ and $q\in(-1,1]$ , or positive $w$ for $q=-1$ . We relate the coefficients of the asymptotic expansions to Nielsen's generalized polylogarithms. For $q=-1$ , 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 $q=1$ , 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.

Generalized hydrodynamics of free fermions under extensive-charge monitoring

2026-04-07T09:14:45-04:00

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.

Hall-Littlewood-positive harmonic functionals on the algebra of symmetric functions

2026-04-06T12:26:59-04:00

We study the problem of describing the set of real functionals on the quotient $\textrm{Sym}/(p_2-1)$ 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 " $p_2$ -twisted action" of $\textrm{Sym}$ on itself and the resulting dual map that makes $\textrm{Sym}$ into a comodule. Finally, our third main result explains the relation between the $p_2$ -twisted comultiplication and the usual comultiplication on $\textrm{Sym}$ .

Joint Estimation in Potts Model

2026-04-06T08:42:05-04:00

In this paper, we study estimation of parameters in a two-parameter Potts model with $q$ colors and coupling matrix $A_N$ . 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 $\sqrt{N}$ . In particular, if $A_N$ is the scaled adjacency matrix of a graph $G_N$ , then we show that joint estimation is possible if either $G_N$ has bounded degree or is irregular. In contrast, we give an example of a graph sequence $G_N$ which is approximately regular and dense, where no consistent estimator exists. We also show that one-parameter estimation at the optimal rate $\sqrt{N}$ 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.

Lozenge Tilings of Hexagons with Intrusions II: Shuffling Phenomenon

2026-04-05T12:08:39-04:00

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.

Domino Tilings of Cruciform Regions

2026-04-05T11:23:15-04:00

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.

A molecular dynamics simulation of thermalization of crystalline lattice with harmonic interaction

2026-04-04T20:47:33-04:00

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.

Equilibrium fluctuations for a multi-species particle system with long jumps

2026-04-04T11:58:09-04:00

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 $x\in\mathbb{Z}$ and $x+z\in\mathbb{Z}$ at a rate that depends on the type of particles involved and is proportional to $|z|^{-\gamma-1}$ for $\gamma>0$ . After rigorously identifying the normal modes associated to the conserved quantities (the density of particles of types $A$ and $B$ , 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.

Classification of product invariant measures for degree-preserving conservative processes and their hydrodynamics

2026-04-03T22:24:37-04:00

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.

The extreme statistics of some noncolliding Brownian processes

2026-04-03T13:29:16-04:00

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.

The KPZ fixed point and Brownian motion share the same null sets

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

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

Higher order derivative moments of CUE characteristic polynomials and the Riemann zeta function

2026-04-03T10:04:33-04:00

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.

Semicircle laws with combined variance for non-uniform Erdős-Rényi hypergraphs

2026-04-02T06:33:14-04:00

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 $(u,v)$ -entry counts the number of hyperedges containing both vertices $u$ and $v$ , 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.

Rotationally invariant first passage percolation: Breaking the n/\log n variance barrier

2026-04-01T13:53:51-04:00

For first passage percolation (FPP) on Euclidean lattices $\mathbb{Z}^d$ with $d\ge 2$ , 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 $1$ . Following Kesten's $O(n)$ 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 $\log n$ 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 $\Theta(n^{2/3})$ , the best known upper bound for the variance of passage times has remained $O(n/\log n)$ in all non-trivial variants of FPP. For a class of rotationally invariant Riemannian FPP on the plane, we show that the variance is $O(n^{1-\varepsilon})$ for some $\varepsilon>0$ . 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.

High genus one part monotone Hurwitz numbers

2026-04-01T13:37:28-04:00

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.

Large deviations of the periodic Toda chain

2026-04-01T04:41:54-04:00

This work establishes a large deviation principle for the spectral measure of the Lax matrix associated to the periodic Toda chain of $N$ 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.

Disordered Schur Measures

2026-03-31T19:40:31-04:00

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.

Two Littlewood identities for fully inhomogeneous spin Hall-Littlewood symmetric rational functions

2026-03-31T10:57:55-04:00

Fully inhomogeneous spin Hall-Littlewood symmetric rational functions $F_\lambda$ arise as partition functions of certain path configurations in the $\mathfrak{sl}_2$ 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 $F_\lambda$ 's over all partitions $\lambda$ 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 $F_\lambda$ 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 $\lambda$ is strictly decreasing.

Three non-Hermitian random matrix universality classes of complex edge statistics: Spacing ratios and distributions

2026-03-30T09:59:38-04:00

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 $N$ . Introducing a conditional point process, we simplify existing expressions and show why an uncontrolled approximation introduced earlier converges well in the large- $N$ limit in the bulk. When specifying to the elliptic Ginibre ensemble, we present a parameter-dependent $N=3$ 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 $\beta$ . 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.

Marked GUE-corners process in doubly periodic dimer models

2026-03-29T19:28:42-04:00

We study a family of periodically weighted Aztec diamond dimer models near their turning points. We establish that, asymptotically, as $N\rightarrow\infty$ , their fluctuations there, scaled by $\sqrt{N}$ , 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.

Asymptotic correlation functions of Coulomb gases on an annulus

2026-03-27T06:47:14-04:00

Two-dimensional Coulomb gases on an annulus at a special inverse temperature $\beta= 2$ 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.

Cycle structure of random standardized permutations

2026-03-25T05:40:25-04:00

In this article, we study a model of random permutations, which we call random standardized permutations, based on a sequence of i.i.d. random variables. This model generalizes others, such as the riffle-shuffle and the major-index-biased permutations. We first establish an exact result on the joint distribution of the number of cycles of given lengths, involving the notion of primitive words. From this result, we obtain various convergence results, most of which are proved using the method of moments. First we prove that the number of small cycles may have either a Poisson limit distribution, or a limit distribution given by a countable sum of independent geometric distributions. Then we establish a limit distribution for large cycles, which is the Poisson-Dirichlet process. Finally we prove a central limit theorem for the total number of cycles.

On matrix valued (asymmetric) truncated Toeplitz operators

2026-03-24T09:45:18-04:00

Matrix valued (asymmetric) truncated Toeplitz operators are generally not complex symmetric. In this paper, we define a new conjugation with unique properties and study its relation to matrix valued asymmetric truncated Toeplitz operators. We also explore the connections between matrix valued asymmetric truncated Toeplitz operators and Hankel operators with matrix symbols.

Edge density expansions for the classical Gaussian and Laguerre ensembles

2026-03-24T05:11:12-04:00

Recent work of Bornemann has uncovered hitherto hidden integrable structures relating to the asymptotic expansion of quantities at the soft edge of Gaussian and Laguerre random matrix ensembles. These quantities are spacing distributions and the eigenvalue density, and the findings cover the cases of the three symmetry classes orthogonal, unitary and symplectic. In this work we give a different viewpoint on these results in the case of the soft edge scaled density, and in the Laguerre case we initiate an analogous study at the hard edge. Our tool is the scalar differential equation satisfied by the latter, known from earlier work. Unlike integral representations, these differential equations in soft edge scaling variables isolate the function of $N$ which is the expansion variable. Moreover, they give information on the correction terms which supplements the findings from the work of Bornemann. In the case of the Gaussian ensemble, we can demonstrate analogous features for Dyson index $\beta= 6$ , which suggests a broader class of models, namely the classical $\beta$ ensembles, with asymptotic expansions exhibiting integrable features. For the Laguerre ensembles at the hard edge, we give the explicit form of the correction at second order for unitary symmetry, and at first order in the orthogonal and symplectic cases. Various differential relations are demonstrated.

Convex upper bounds for last passage percolation (LPP) times in the Kardar-Parisi-Zhang (KPZ) universality class

2026-03-23T16:01:32-04:00

On the $Z^2$ lattice, vertices are assigned random weights $W(i,j)$ . The point-to-point last passage percolation (LPP) time $S_{M,N+1-M}$ between $(1,1)$ and $(M,N+1-M)$ is the maximum total weight among all upward/right-directed paths connecting the two. Point-to-line LPP time $R_N$ is the maximum of these maximal total weights over $M$ . Asymptotic distributions and fluctuations of these LPP times have been studied for i.i.d. weights. The current study deals with identically distributed but not necessarily independent weights, and maximizes LPP times in the sense of convex dominance. In particular, maximal expected LPP times are identified, in the class of all weight couplings with a given marginal distribution. For the case of mean- $1$ exponentially distributed weights, there is a coupling for which $R_N$ is the shifted exponential variable $R_N^* = N W(1,1) + \log(N!)$ , such that $E[\Psi(R_N)] \le E[\Psi(R_N^*)]$ for all couplings and all convex non-decreasing functions $\Psi$ for which these expectations are well defined. In contrast to ${{R_N^*} \over N}= W(1,1)+{{\log(N!)} \over N}$ , with variance $1$ and mean diverging to $\infty$ like $\log(N)$ , ${{R_N} \over N}$ converges to $2$ for the commonly studied i.i.d. weights.

Multivariable Painleve'-II equation: connection formulas for asymptotic solutions

2026-03-23T14:38:37-04:00

It is shown that a generalization of the Painlevé-II equation (P-II) to a system of coupled equations with symmetry breaking terms is integrable. A Lax pair for this system is used to relate the asymptotic behavior of the solutions at different infinities via an asymptotically exact WKB approach. The analysis relies on an exact solution of the quantum mechanical Demkov-Osherov model (DOM). An application to the problem of unstable vacuum decay during a second order phase transition provides precise scaling of the number of excitations, including subdominant contributions.

The Littlewood-Richardson rule for Schur P-, Q-multiple zeta functions

2026-03-22T22:00:43-04:00

The Schur $P$ -, $Q$ -multiple zeta functions were defined by Nakasuji and Takeda inspired by the tableau representation of Schur $P$ -, $Q$ -functions. While a product of two Schur $P$ -functions expands as a linear combination of Schur $P$ -functions, we obtain a similar expansion formula for the Schur $P$ -multiple zeta functions by taking summation over the symmetric group permutating all the variables. We also introduce a expansion formula of skew Schur $Q$ -multiple zeta functions by taking summation over the symmetric group. Furthermore, this skew type formula can be refined by restricting the symmetric group to its specific subgroup.