Source author record

Zakhar Kabluchko

Zakhar Kabluchko appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

68works
17topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

68 published item(s)

preprint2026arXiv

First passage times for decoupled random walks

Motivated by a connection to the infinite Ginibre point process, decoupled random walks were introduced in a recent article Alsmeyer, Iksanov and Kabluchko (2025). The decoupled random walk is a sequence of independent random variables, in which the $n$th variable has the same distribution as the position at time $n$ of a standard random walk with nonnegative increments. We prove distributional convergence in the Skorokhod space equipped with the $J_1$-topology of the running maxima and the first passage times of decoupled random walks. We show that there exist five different regimes, in which distinct limit theorems arise. Rather different functional limit theorems for the number of visits of decoupled standard random walk to the interval $[0,t]$ as $t\to\infty$ were earlier obtained in the aforementioned paper Alsmeyer, Iksanov and Kabluchko (2025). While the limit processes for the first passage times are inverse extremal-like processes, the limit processes for the number of visits are stationary Gaussian.

preprint2026arXiv

Zero distribution of multiplicative Hermite and Laguerre polynomials

It is well-known that, as $n\to\infty$, the zero distribution of the $n$-th Hermite polynomial converges to the semicircular law (the free normal distribution), while the zero distribution of the associated Laguerre polynomials converges to the Marchenko--Pastur law (the free Poisson distribution). In this paper, we establish multiplicative analogues of these results. We define the multiplicative Hermite and Laguerre polynomials by \begin{align*} H_n^*(x;s) &:= e^{-\frac 12 s ((x\partial_x)^2 - n x \partial_x) } (x-1)^n = \sum_{j=0}^n (-1)^{n-j} \binom nj e^{-\frac 12 s (j^2 - nj)} x^j, \\ L_n^*(x; b,c) &:= (x\partial_x + b)^c (x-1)^n = \sum_{j=0}^n (-1)^{n-j} \binom nj (j+b)^c x^j, \end{align*} where $n\in \mathbb N_0$, $\partial_x$ denotes the differentiation operator w.r.t. $x$, and $s\in \mathbb R$, $b\in \mathbb C$, $c\in \mathbb N_0$ are parameters. In the Hermite case, we show that, as $n\to\infty$, the zero distribution of $H_n^*(x;s/n)$ converges weakly to the free multiplicative normal distribution on the positive half-line (when $s>0$) or to the free unitary normal distribution on the unit circle $\{|z| = 1\}$ (when $s<0$). In the Laguerre case, we show that the zero distribution of $L_n^*(x; nβ, \lfloor n γ\rfloor)$ converges to the free multiplicative Poisson distribution on the positive half-line (when $γ>0$ and $β\in \mathbb R\backslash[0,1]$) or on the unit circle (when $γ>0$ and $β\in -\frac 12 + \sqrt{-1} \, \mathbb R$). All these results are obtained by essentially the same method, which treats the Hermite/Laguerre cases and the unitary/positive settings in a unified way.

preprint2023arXiv

Mod-$φ$ convergence of Stirling distributions and limit theorems for zeros of their generating functions

We study mod-$φ$ convergence of several probability distributions on the set of positive integers that involve Stirling numbers of both kinds and, as a consequence, derive various limit theorems for these distributions. We also derive closely related limit theorems for the distribution of zeros of the corresponding generating functions. For example, we identify the asymptotic distribution of zeros for the generating polynomial of the number of occupied boxes when $n$ balls are allocated equiprobably and independently among $θ$ boxes in the regime when $θ$ grows linearly with $n$.

preprint2023arXiv

Sectional Voronoi tessellations: Characterization and high-dimensional limits

The intersections of beta-Voronoi, beta-prime-Voronoi and Gaussian-Voronoi tessellations in $\mathbb{R}^d$ with $\ell$-dimensional affine subspaces, $1\leq \ell\leq d-1$, are shown to be random tessellations of the same type but with different model parameters. In particular, the intersection of a classical Poisson-Voronoi tessellation with an affine subspace is shown to have the same distribution as a certain beta-Voronoi tessellation. The geometric properties of the typical cell and, more generally, typical $k$-faces, of the sectional Poisson-Voronoi tessellation are studied in detail. It is proved that in high dimensions, that is as $d\to\infty$, the intersection of the $d$-dimensional Poison-Voronoi tessellation with an affine subspace of fixed dimension $\ell$ converges to the $\ell$-dimensional Gaussian-Voronoi tessellation.

preprint2022arXiv

Beta-star polytopes and hyperbolic stochastic geometry

Motivated by problems of hyperbolic stochastic geometry we introduce and study the class of beta-star polytopes. A beta-star polytope is defined as the convex hull of an inhomogeneous Poisson processes on the complement of the unit ball in $\mathbb{R}^d$ with density proportional to $(||x||^2-1)^β$, where $||x|| > 1$ and $β>d/2$. Explicit formulas for various geometric and combinatorial functionals associated with beta-star polytopes are provided, including the expected number of $k$-dimensional faces, the expected external angle sums and the expected intrinsic volumes. Beta-star polytopes are relevant in the context of hyperbolic stochastic geometry, since they are tightly connected to the typical cell of a Poisson-Voronoi tessellation as well as the zero cell of a Poisson hyperplane tessellation in hyperbolic space. The general results for beta-star polytopes are used to provide explicit formulas for the expected $f$- vector of the typical hyperbolic Poisson-Voronoi cell and the hyperbolic Poisson zero cell. Their asymptotics for large intensities and their monotonicity behaviour is discussed as well. Finally, stochastic geometry in the de Sitter half-space is studied as the hyperbolic analogue to recent investigations about random cones generated by random points on half-spheres in spherical or conical stochastic geometry.

preprint2022arXiv

Face numbers of high-dimensional Poisson zero cells

Let $\mathcal Z_d$ be the zero cell of a $d$-dimensional, isotropic and stationary Poisson hyperplane tessellation. We study the asymptotic behavior of the expected number of $k$-dimensional faces of $\mathcal Z_d$, as $d\to\infty$. For example, we show that the expected number of hyperfaces of $\mathcal Z_d$ is asymptotically equivalent to $\sqrt{2π/3}\, d^{3/2}$, as $d\to\infty$. We also prove that the expected solid angle of a random cone spanned by $d$ random vectors that are independent and uniformly distributed on the unit upper half-sphere in $\mathbb R^{d}$ is asymptotic to $\sqrt 3 π^{-d}$, as $d\to\infty$.

preprint2022arXiv

Lah distribution: Stirling numbers, records on compositions, and convex hulls of high-dimensional random walks

Let $ξ_1,ξ_2,\ldots$ be a sequence of independent copies of a random vector in $\mathbb R^d$ having an absolutely continuous distribution. Consider a random walk $S_i:=ξ_1+\cdots+ξ_i$, and let $C_{n,d}:=\text{conv}(0,S_1,S_2,\ldots,S_n)$ be the convex hull of the first $n+1$ points it has visited. The polytope $C_{n,d}$ is called $k$-neighborly if for every indices $0\leq i_0 <\cdots < i_k\leq n$ the convex hull of the $k+1$ points $S_{i_0},\ldots, S_{i_k}$ is a $k$-dimensional face of $C_{n,d}$. We study the probability that $C_{n,d}$ is $k$-neighborly in various high-dimensional asymptotic regimes, i.e. when $n$, $d$, and possibly also $k$ diverge to $\infty$. There is an explicit formula for the expected number of $k$-dimensional faces of $C_{n,d}$ which involves Stirling numbers of both kinds. Motivated by this formula, we introduce a distribution, called the Lah distribution, and study its properties. In particular, we provide a combinatorial interpretation of the Lah distribution in terms of random compositions and records, and explicitly compute its factorial moments. Limit theorems which we prove for the Lah distribution imply neighborliness properties of $C_{n,d}$. This yields a new class of random polytopes exhibiting phase transitions parallel to those discovered by Vershik and Sporyshev, Donoho and Tanner for random projections of regular simplices and crosspolytopes.

preprint2022arXiv

On weak Fraisse limits

Using the natural action of $S_\infty$ we show that a countable hereditary class $\mathcal C$ of finitely generated structures has the joint embedding property (JEP) and the weak amalgamation property (WAP) if and only if there is a structure $M$ whose isomorphism type is comeager in the space of all countable, infinitely generated structures with age in $\mathcal C$. In this case, $M$ is the weak Fraïssé limit of $\mathcal C$. This applies in particular to countable structures with generic automorphisms and recovers a result by Kechris and Rosendal [Proc. Lond. Math. Soc., 2007].

preprint2022arXiv

Random cones in high dimensions I: Donoho-Tanner and Cover-Efron cones

Two models of random cones in high dimensions are considered, together with their duals. The Donoho-Tanner random cone $D_{n,d}$ can be defined as the positive hull of $n$ independent $d$-dimensional Gaussian random vectors. The Cover-Efron random cone $C_{n,d}$ is essentially defined as the same positive hull, conditioned on the event that it is not the whole space. We consider expectations of various combinatorial and geometric functionals of these random cones and prove that they satisfy limit theorems, as $d$ and $n$ tend to infinity in a suitably coordinated way. This includes, for example, large deviation principles and central as well as non-central limit theorems for the expected number of $k$-faces and the $k$-th conic intrinsic volumes, as $n$, $d$ and possibly also $k$ tend to infinity simultaneously. Furthermore, we determine the precise high-dimensional asymptotic behaviour of the expected statistical dimension for both models of random cones, uncovering thereby another high-dimensional phase transition. As an application, limit theorems for the number of $k$-faces of high-dimensional polytopes generated by random Gale diagrams are discussed as well.

preprint2021arXiv

Angle sums of Schläfli orthoschemes

We consider the simplices $$ K_n^A=\{x\in\mathbb{R}^{n+1}:x_1\ge x_2\ge \ldots\ge x_{n+1},x_1-x_{n+1}\le 1,x_1+\ldots+x_{n+1}=0\} $$ and $$ K_n^B=\{x\in\mathbb{R}^n:1\ge x_1\ge x_2\ge \ldots\ge x_n\ge 0\}, $$ which are called the Schläfli orthoschemes of types $A$ and $B$, respectively. We describe the tangent cones at their $j$-faces and compute explicitly the sum of the conic intrinsic volumes of these tangent cones at all $j$-faces of $K_n^A$ and $K_n^B$. This setting contains sums of external and internal angles of $K_n^A$ and $K_n^B$ as special cases. The sums are evaluated in terms of Stirling numbers of both kinds. We generalize these results to finite products of Schläfli orthoschemes of type $A$ and $B$ and, as a probabilistic consequence, derive formulas for the expected number of $j$-faces of the Minkowski sums of the convex hulls of a finite number of Gaussian random walks and random bridges. Furthermore, we evaluate the analogous angle sums for the tangent cones of Weyl chambers of types $A$ and $B$ and finite products thereof.

preprint2021arXiv

Fluctuations of the Magnetization for Ising models on Erdős-Rényi Random Graphs -- the Regimes of Low Temperature and External Magnetic Field

We continue our analysis of Ising models on the (directed) Erdős-Rényi random graph $G(N,p)$. We prove a quenched Central Limit Theorem for the magnetization and describe the fluctuations of the log-partition function. In the current note we consider the low temperature regime $β>1$ and the case when an external magnetic field is present. In both cases, we assume that $p=p(N)$ satisfies $p^3N \to \infty$.

preprint2021arXiv

Phase transition for the volume of high-dimensional random polytopes

The beta polytope $P_{n,d}^β$ is the convex hull of $n$ i.i.d. random points distributed in the unit ball of $\mathbb{R}^d$ according to a density proportional to $(1-\lVert{x}\rVert^2)^β$ if $β>-1$ (in particular, $β=0$ corresponds to the uniform distribution in the ball), or uniformly on the unit sphere if $β=-1$. We show that the expected normalized volumes of high-dimensional beta polytopes exhibit a phase transition and we describe its shape. We derive analogous results for the intrinsic volumes of beta polytopes and, when $β=0$, their number of vertices.

preprint2021arXiv

The $β$-Delaunay tessellation II: The Gaussian limit tessellation

We study the weak convergence of $β$- and $β'$-Delaunay tessellations in $\mathbb{R}^{d-1}$ that were introduced in part I of this paper, as $β\to\infty$. The limiting stationary simplicial random tessellation, which is called the Gaussian-Delaunay tessellation, is characterized in terms of a space-time paraboloid hull process in $\mathbb{R}^{d-1}\times\mathbb{R}$. The latter object has previously appeared in the analysis of the number of shocks in the solution of the inviscid Burgers' equation and the description of the local asymptotic geometry of Gaussian random polytopes. In this paper it is used to define a new stationary random simplicial tessellation in $\mathbb{R}^{d-1}$. As for the $β$- and $β'$-Delaunay tessellation, the distribution of volume-power weighted typical cells in the Gaussian-Delaunay tessellation is explicitly identified, establishing thereby a new bridge to Gaussian random simplices. Also major geometric characteristics of these cells such as volume moments, expected angle sums and also the cell intensities of the Gaussian-Delaunay tessellation are investigated.

preprint2021arXiv

The $β$-Delaunay tessellation III: Kendall's problem and limit theorems in high dimensions

The $β$-Delaunay tessellation in $\mathbb{R}^{d-1}$ is a generalization of the classical Poisson-Delaunay tessellation. As a first result of this paper we show that the shape of a weighted typical cell of a $β$-Delaunay tessellation, conditioned on having large volume, is close to the shape of a regular simplex in $\mathbb{R}^{d-1}$. This generalizes earlier results of Hug and Schneider about the typical (non-weighted) Poisson-Delaunay simplex. Second, the asymptotic behaviour of the volume of weighted typical cells in high-dimensional $β$-Delaunay tessellation is analysed, as $d\to\infty$. In particular, various high dimensional limit theorems, such as quantitative central limit theorems as well as moderate and large deviation principles, are derived.

preprint2021arXiv

The typical cell of a Voronoi tessellation on the sphere

The typical cell of a Voronoi tessellation generated by $n+1$ uniformly distributed random points on the $d$-dimensional unit sphere $\mathbb S^d$ is studied. Its $f$-vector is identified in distribution with the $f$-vector of a beta' polytope generated by $n$ random points in $\mathbb R^d$. Explicit formulae for the expected $f$-vector are provided for any $d$ and the low-dimensional cases $d\in\{2,3,4\}$ are studied separately. This implies an explicit formula for the total number of $k$-dimensional faces in the spherical Voronoi tessellation as well.

preprint2020arXiv

An identity for the coefficients of characteristic polynomials of hyperplane arrangements

Consider a finite collection of affine hyperplanes in $\mathbb R^d$. The hyperplanes dissect $\mathbb R^d$ into finitely many polyhedral chambers. For a point $x\in \mathbb R^d$ and a chamber $P$ the metric projection of $x$ onto $P$ is the unique point $y\in P$ minimizing the Euclidean distance to $x$. The metric projection is contained in the relative interior of a uniquely defined face of $P$ whose dimension is denoted by $\text{dim}(x,P)$. We prove that for every given $k\in \{0,\ldots, d\}$, the number of chambers $P$ for which $\text{dim}(x,P) = k$ does not depend on the choice of $x$, with an exception of some Lebesgue null set. Moreover, this number is equal to the absolute value of the $k$-th coefficient of the characteristic polynomial of the hyperplane arrangement. In a special case of reflection arrangements, this proves a conjecture of Drton and Klivans [A geometric interpretation of the characteristic polynomial of reflection arrangements, Proc. Amer. Math. Soc., 138(8): 2873-2887, 2010].

preprint2020arXiv

Angle sums of random polytopes

For two families of random polytopes we compute explicitly the expected sums of the conic intrinsic volumes and the Grassmann angles at all faces of any given dimension of the polytope under consideration. As special cases, we compute the expected sums of internal and external angles at all faces of any fixed dimension. The first family are the Gaussian polytopes defined as convex hulls of i.i.d. samples from a non-degenerate Gaussian distribution in $\mathbb R^d$. The second family are convex hulls of random walks with exchangeable increments satisfying certain mild general position assumption. The expected sums are expressed in terms of the angles of the regular simplices and the Stirling numbers, respectively. There are non-trivial analogies between these two settings. Further, we compute the angle sums for Gaussian projections of arbitrary polyhedral sets, of which the Gaussian polytopes are a special case. Also, we show that the expected Grassmann angle sums of a random polytope with a rotationally invariant law are invariant under affine transformations. Of independent interest may be also results on the faces of linear images of polyhedral sets. These results are well known but it seems that no detailed proofs can be found in the existing literature.

preprint2020arXiv

Angles of Random Simplices and Face Numbers of Random Polytopes

Pick $d+1$ points uniformly at random on the unit sphere in $\mathbb R^d$. What is the expected value of the angle sum of the simplex spanned by these points? Choose $n$ points uniformly at random in the $d$-dimensional ball. What is the expected number of faces of their convex hull? We answer these and some related questions of stochastic geometry. To this end, we compute expected internal angles of random simplices whose vertices are independent random points sampled from one of the following $d$-dimensional distributions: (i) the beta distribution with the density proportional to $(1-\|x\|^2)^β$, where $x$ is belongs to the unit ball in $\mathbb R^d$; (ii) the beta' distribution with the density proportional to $(1+\|x\|^2)^{-β}$, where $x\in\mathbb{R}^{d}$. These results imply explicit formulae for the expected face numbers of the following random polytopes: (a) the typical Poisson-Voronoi cell; (b) the zero cell of the Poisson hyperplane tessellation; (c) beta and beta' polytopes defined as convex hulls of i.i.d. samples from the corresponding distributions.

preprint2020arXiv

Beta polytopes and Poisson polyhedra: $f$-vectors and angles

We study random polytopes of the form $[X_1,\ldots,X_n]$ defined as convex hulls of independent and identically distributed random points $X_1,\ldots,X_n$ in $\mathbb{R}^d$ with one of the following densities: $$ f_{d,β} (x) = c_{d,β} (1-\|x\|^2)^β, \qquad \|x\| < 1, \quad \text{(beta distribution, $β>-1$)} $$ or $$ \tilde f_{d,β} (x) = \tilde{c}_{d,β} (1+\|x\|^2)^{-β}, \qquad x\in\mathbb{R}^d, \quad \text{(beta' distribution, $β>d/2$)}. $$ This setting also includes the uniform distribution on the unit sphere and the standard normal distribution as limiting cases. We derive exact and asymptotic formulae for the expected number of $k$-faces of $[X_1,\ldots,X_n]$ for arbitrary $k\in\{0,1,\ldots,d-1\}$. We prove that for any such $k$ this expected number is strictly monotonically increasing with $n$. Also, we compute the expected internal and external angles of these polytopes at faces of every dimension and, more generally, the expected conic intrinsic volumes of their tangent cones. By passing to the large $n$ limit in the beta' case, we compute the expected $f$-vector of the convex hull of Poisson point processes with power-law intensity function. Using convex duality, we derive exact formulae for the expected number of $k$-faces of the zero cell for a class of isotropic Poisson hyperplane tessellations in $\mathbb R^d$. This family includes the zero cell of a classical stationary and isotropic Poisson hyperplane tessellation and the typical cell of a stationary Poisson--Voronoi tessellation as special cases. In addition, we prove precise limit theorems for this $f$-vector in the high-dimensional regime, as $d\to\infty$. Finally, we relate the $d$-dimensional beta and beta' distributions to the generalized Pareto distributions known in extreme-value theory.

preprint2020arXiv

Convex cones spanned by regular polytopes

We study three families of polyhedral cones whose sections are regular simplices, cubes, and crosspolytopes. We compute solid angles and conic intrinsic volumes of these cones. We show that several quantities appearing in stochastic geometry can be expressed through these conic intrinsic volumes. A list of such quantities includes internal and external solid angles of regular simplices and crosspolytopes, the probability that a (symmetric) Gaussian random polytope or the Gaussian zonotope contains a given point, the expected number of faces of the intersection of a regular polytope with a random linear subspace passing through its centre, and the expected number of faces of the projection of a regular polytope onto a random linear subspace.

preprint2020arXiv

Expected $f$-vector of the Poisson Zero Polytope and Random Convex Hulls in the Half-Sphere

We prove an explicit combinatorial formula for the expected number of faces of the zero polytope of the homogeneous and isotropic Poisson hyperplane tessellation in $\mathbb R^d$. The expected $f$-vector is expressed through the coefficients of the polynomial $$ (1+ (d-1)^2x^2) (1+(d-3)^2 x^2) (1+(d-5)^2 x^2) \ldots. $$ Also, we compute explicitly the expected $f$-vector and the expected volume of the spherical convex hull of $n$ random points sampled uniformly and independently from the $d$-dimensional half-sphere. In the case when $n=d+2$, we compute the probability that this spherical convex hull is a spherical simplex, thus solving an analogue of the Sylvester four-point problem on the half-sphere.

preprint2020arXiv

Faces in random great hypersphere tessellations

The concept of typical and weighted typical spherical faces for tessellations of the $d$-dimensional unit sphere, generated by $n$ independent random great hyperspheres distributed according to a non-degenerate directional distribution, is introduced and studied. Probabilistic interpretations for such spherical faces are given and their directional distributions are determined. Explicit formulas for the expected $f$-vector, the expected spherical Quermaßintegrals and the expected spherical intrinsic volumes are found in the isotropic case. Their limiting behaviour as $n\to\infty$ is discussed and compared to the corresponding notions and results in the Euclidean case. The expected statistical dimension and a problem related to intersection probabilities of spherical random polytopes is investigated.

preprint2020arXiv

Fluctuations for the partition function of Ising models on Erdös-Rényi random graphs

We analyze Ising/Curie-Weiss models on the Erdős-Rényi graph with $N$ vertices and edge probability $p=p(N)$ that were introduced by Bovier and Gayrard [J.\ Statist.\ Phys., 72(3-4):643--664, 1993] and investigated in two previous articles by the authors. We prove Central Limit Theorems for the partition function of the model and -- at other decay regimes of $p(N)$ -- for the logarithmic partition function. We find critical regimes for $p(N)$ at which the behavior of the fluctuations of the partition function changes.

preprint2020arXiv

Grassmann angles and absorption probabilities of Gaussian convex hulls

Let $M$ be an arbitrary subset in $\mathbb R^n$ with a conic (or positive) hull $C$. Consider its Gaussian image $AM$, where $A$ is a $k\times n$-matrix whose entries are independent standard Gaussian random variables. We show that the probability that the convex hull of $AM$ contains the origin in its interior coincides with the $k$-th Grassmann angle of $C$. Also, we prove that the expected Grassmann angles of $AC$ coincide with the corresponding Grassmann angles of $C$. Using the latter result, we show that the expected sum of $j$-th Grassmann angles at $\ell$-dimensional faces of a Gaussian simplex equals the analogous angle-sum for the regular simplex of the same dimension.

preprint2020arXiv

How long is the convex minorant of a one-dimensional random walk?

We prove distributional limit theorems for the length of the largest convex minorant of a one-dimensional random walk with independent identically distributed increments. Depending on the increment law, there are several regimes with different limit distributions for this length. Among other tools, a representation of the convex minorant of a random walk in terms of uniform random permutations is utilized.

preprint2020arXiv

Limit theorems for random points in a simplex

In this work the $\ell_q$-norms of points chosen uniformly at random in a centered regular simplex in high dimensions are studied. Berry-Esseen bounds in the regime $1\leq q < \infty$ are derived and complemented by a non-central limit theorem together with moderate and large deviations in the case where $q=\infty$. A comparison with corresponding results for $\ell_p^n$-balls is carried out as well.

preprint2020arXiv

Recursive Scheme for Angles of Random Simplices, and Applications to Random Polytopes

Consider a random simplex $[X_1,\ldots,X_n]$ defined as the convex hull of independent identically distributed random points $X_1,\ldots,X_n$ in $\mathbb{R}^{n-1}$ with the following beta density: $$ f_{n-1,β} (x) \propto (1-\|x\|^2)^β 1_{\{\|x\| < 1\}}, \qquad x\in\mathbb{R}^{n-1}, \quad β>-1. $$ Let $J_{n,k}(β)$ be the expected internal angle of the simplex $[X_1,\ldots,X_n]$ at its face $[X_1,\ldots,X_k]$. Define $\tilde J_{n,k}(β)$ analogously for i.i.d. random points distributed according to the beta' density $$ \tilde f_{n-1,β} (x) \propto (1+\|x\|^2)^{-β}, \qquad x\in\mathbb{R}^{n-1}, \quad β> \frac{n-1}{2}. $$ We derive formulae for $J_{n,k}(β)$ and $\tilde J_{n,k}(β)$ which make it possible to compute these quantities symbolically, in finitely many steps, for any integer or half-integer value of $β$. For $J_{n,1}(\pm 1/2)$ we even provide explicit formulae in terms of products of Gamma functions. We give applications of these results to two seemingly unrelated problems of stochastic geometry. (i) We compute the expected $f$-vectors of the typical Poisson-Voronoi cells in dimensions up to $10$. (ii) Consider the random polytope $K_{n,d} := [U_1,\ldots,U_n]$ where $U_1,\ldots,U_n$ are i.i.d. random points sampled uniformly inside some $d$-dimensional convex body $K$ with smooth boundary and unit volume. M. Reitzner proved the existence of the limit of the normalized expected $f$-vector of $K_{n,d}$: $$ \lim_{n\to\infty} n^{-{\frac{d-1}{d+1}}}\mathbb E \mathbf f(K_{n,d}) = \mathbf c_d \cdot Ω(K), $$ where $Ω(K)$ is the affine surface area of $K$, and $\mathbf c_d$ is an unknown vector not depending on $K$. We compute $\mathbf c_d$ explicitly in dimensions up to $d=10$ and also solve the analogous problem for random polytopes with vertices distributed uniformly on the sphere.

preprint2020arXiv

The maximum entropy principle and volumetric properties of Orlicz balls

We study the precise asymptotic volume of balls in Orlicz spaces and show that the volume of the intersection of two Orlicz balls undergoes a phase transition when the dimension of the ambient space tends to infinity. This generalizes a result of Schechtman and Schmuckenschläger [GAFA, Lecture notes in Math. 1469 (1991), 174--178] for $\ell_p^d$-balls. As another application, we determine the precise asymptotic volume ratio for $2$-concave Orlicz spaces $\ell_M^d$. Our method rests on ideas from statistical mechanics and large deviations theory, more precisely the maximum entropy or Gibbs principle for non-interacting particles, and presents a natural approach and fresh perspective to such geometric and volumetric questions. In particular, our approach explains how the $p$-generalized Gaussian distribution occurs in problems related to the geometry of $\ell_p^d$-balls, which are Orlicz balls when the Orlicz function is $M(t) = |t|^p$.

preprint2018arXiv

Distribution of complex algebraic numbers on the unit circle

For $-π\leqβ_1<β_2\leqπ$ denote by $Φ_{β_1,β_2}(Q)$ the number of algebraic numbers on the unit circle with arguments in $[β_1,β_2]$ of degree $2m$ and with elliptic height at most $Q$. We show that \[ Φ_{β_1,β_2}(Q)=Q^{m+1}\int\limits_{β_1}^{β_2}{p(t)}\,{\rm d}t+O\left(Q^m\,\log Q\right),\quad Q\to\infty, \] where $p(t)$ coincides up to a constant factor with the density of the roots of some random trigonometric polynomial. This density is calculated explicitly using the Edelman--Kostlan formula.

preprint2017arXiv

Adaptive Energy Saving Approximation for Random Stationary Processes

We consider a stationary process (with either discrete or continuous time) and find an adaptive approximating stationary process combining approximation quality and supplementary good properties that can be interpreted as additional smoothness or small expense of energy. The problem is solved in terms of the spectral characteristics of the approximated process by using classical analytic methods from prediction theory.

preprint2016arXiv

A central limit theorem and a law of the iterated logarithm for the Biggins martingale of the supercritical branching random walk

Let $(W_n(θ))_{n\in\mathbb N_0}$ be the Biggins martingale associated with a supercritical branching random walk and denote by $W_\infty(θ)$ its limit. Assuming essentially that the martingale $(W_n(2θ))_{n\in\mathbb N_0}$ is uniformly integrable and that $\text{Var} W_1(θ)$ is finite, we prove a functional central limit theorem for the tail process $(W_\infty(θ) - W_{n+r}(θ))_{r\in\mathbb N_0}$ and a law of the iterated logarithm for $W_\infty(θ)-W_n(θ)$, as $n\to\infty$.

preprint2016arXiv

A leader-election procedure using records

The study of the number of collisions in a Poisson-Dirichlet coalescent leads to the analysis of the following version of a stochastic leader-elec\-tion algorithm. Consider an infinite family of persons, labeled by $1,2,3,\ldots$, who generate iid random numbers from an arbitrary continuous distribution. Those persons who have generated a record value, that is, a value larger than the values of all previous persons, stay in the game, all others must leave. The remaining persons are relabeled by $1,2,3,\ldots$ maintaining their order in the first round, and the election procedure is repeated independently from the past and indefinitely. We prove limit theorems for a number of relevant functionals for this procedure, notably the number of rounds $T(M)$ until all persons among $1,\ldots,M$, except the first one, have left (as $M\to\infty$). For example, we show that the sequence $(T(M)-\log^{*}M)_{M\in\mathbb{N}}$, where $\log^{*}$ denotes the iterated logarithm, is tight, and study its weak subsequential limits. We further provide an appropriate and apparently new kind of normalization (based on tetrations) such that the original labels of persons who stay in the game until round $n$ converge (as $n\to\infty$) to some random non-Poissonian point process and study its properties. The results are applied to study subsequential distributional limits for the number of collisions in the Poisson-Dirichlet coalescent.

preprint2016arXiv

Ergodic decompositions of stationary max-stable processes in terms of their spectral functions

We revisit conservative/dissipative and positive/null decompositions of stationary max-stable processes. Originally, both decompositions were defined in an abstract way based on the underlying non-singular flow representation. We provide simple criteria which allow to tell whether a given spectral function belongs to the conservative/dissipative or positive/null part of the de Haan spectral representation. Specifically, we prove that a spectral function is null-recurrent iff it converges to $0$ in the Cesàro sense. For processes with locally bounded sample paths we show that a spectral function is dissipative iff it converges to $0$. Surprisingly, for such processes a spectral function is integrable a.s. iff it converges to $0$ a.s. Based on these results, we provide new criteria for ergodicity, mixing, and existence of a mixed moving maximum representation of a stationary max-stable process in terms of its spectral functions. In particular, we study a decomposition of max-stable processes which characterizes the mixing property.

preprint2016arXiv

Fractionally integrated inverse stable subordinators

A fractionally integrated inverse stable subordinator (FIISS) is the convolution of a power function and an inverse stable subordinator. We show that the FIISS is a scaling limit in the Skorokhod space of a renewal shot noise process with heavy-tailed, infinite mean `inter-shot' distribution and regularly varying response function. We prove local Hölder continuity of FIISS and a law of iterated logarithm for both small and large times.

preprint2016arXiv

Functional limit theorems for Galton-Watson processes with very active immigration

We prove weak convergence on the Skorokhod space of Galton-Watson processes with immigration, properly normalized, under the assumption that the tail of the immigration distribution has a logarithmic decay. The limits are extremal shot noise processes. By considering marginal distributions, we recover the results of Pakes [Adv. Appl. Probab., 11(1979), 31--62].

preprint2016arXiv

Inclusion-exclusion principles for convex hulls and the Euler relation

Consider $n$ points $X_1,\ldots,X_n$ in $\mathbb R^d$ and denote their convex hull by $Π$. We prove a number of inclusion-exclusion identities for the system of convex hulls $Π_I:=conv(X_i\colon i\in I)$, where $I$ ranges over all subsets of $\{1,\ldots,n\}$. For instance, denoting by $c_k(X)$ the number of $k$-element subcollections of $(X_1,\ldots,X_n)$ whose convex hull contains a point $X\in\mathbb R^d$, we prove that $$ c_1(X)-c_2(X)+c_3(X)-\ldots + (-1)^{n-1} c_n(X) = (-1)^{\dim Π} $$ for all $X$ in the relative interior of $Π$. This confirms a conjecture of R. Cowan [Adv. Appl. Probab., 39(3):630--644, 2007] who proved the above formula for almost all $X$. We establish similar results for the number of polytopes $Π_J$ containing a given polytope $Π_I$ as an $r$-dimensional face, thus proving another conjecture of R. Cowan [Discrete Comput. Geom., 43(2):209--220, 2010]. As a consequence, we derive inclusion-exclusion identities for the intrinsic volumes and the face numbers of the polytopes $Π_I$. The main tool in our proofs is a formula for the alternating sum of the face numbers of a convex polytope intersected by an affine subspace. This formula generalizes the classical Euler--Schläfli--Poincaré relation and is of independent interest.

preprint2016arXiv

Local universality for real roots of random trigonometric polynomials

Consider a random trigonometric polynomial $X_n: \mathbb R \to \mathbb R$ of the form $$ X_n(t) = \sum_{k=1}^n \left( ξ_k \sin (kt) + η_k \cos (kt)\right), $$ where $(ξ_1,η_1),(ξ_2,η_2),\ldots$ are independent identically distributed bivariate real random vectors with zero mean and unit covariance matrix. Let $(s_n)_{n\in\mathbb N}$ be any sequence of real numbers. We prove that as $n\to\infty$, the number of real zeros of $X_n$ in the interval $[s_n+a/n, s_n+ b/n]$ converges in distribution to the number of zeros in the interval $[a,b]$ of a stationary, zero-mean Gaussian process with correlation function $(\sin t)/t$. We also establish similar local universality results for the centered random vectors $(ξ_k,η_k)$ having an arbitrary covariance matrix or belonging to the domain of attraction of a two-dimensional $α$-stable law.

preprint2016arXiv

Mean width of regular polytopes and expected maxima of correlated Gaussian variables

An old conjecture states that among all simplices inscribed in the unit sphere the regular one has the maximal mean width. An equivalent formulation is that for any centered Gaussian vector $(ξ_1,\dots,ξ_n)$ satisfying $\mathbb Eξ_1^2= \dots =\mathbb Eξ_n^2=1$ one has $$ \mathbb E\,\max\{ξ_1,\dots,ξ_n\}\leq\sqrt{\frac{n}{n-1}}\, \mathbb E\,\max\{η_1,\dots,η_n\}, $$ where $η_1,η_2,\dots,$ are independent standard Gaussian variables. Using this probabilistic interpretation we derive an asymptotic version of the conjecture. We also show that the mean width of the regular simplex with $2n$ vertices is remarkably close to the mean width of the regular crosspolytope with the same number of vertices. Interpreted probabilistically, our result states that $$ 1\leq\frac{\mathbb E\,\max\{|η_1|,\dots,|η_n|\}}{\mathbb E\,\max\{η_1,\dots,η_{2n}\}} \leq\min\left\{\sqrt{\frac{2n}{2n-1}}, \, 1+\frac{C}{n\, \log n} \right\}, $$ where $C>0$ is an absolute constant. We also compute the higher moments of the projection length $W$ of the regular cube, simplex and crosspolytope onto a line with random direction, thus proving several formulas conjectured by S. Finch. Finally, we prove distributional limit theorems for the length of random projection as the dimension goes to $\infty$. In the case of the $n$-dimensional unit cube $Q_n$, we prove that $$ W_{Q_n} - \sqrt{\frac{2n}π} \overset{d}{\underset{n\to\infty}\longrightarrow} {\mathcal{N}} \left(0, \frac{π-3}π\right), $$ whereas for the simplex and the crosspolytope the limiting distributions are related to the Gumbel double exponential law.

preprint2016arXiv

Mode and Edgeworth expansion for the Ewens distribution and the Stirling numbers

We provide asymptotic expansions for the Stirling numbers of the first kind and, more generally, the Ewens (or Karamata-Stirling) distribution. Based on these expansions, we obtain some new results on the asymptotic properties of the mode and the maximum of the Stirling numbers and the Ewens distribution. For arbitrary $θ>0$ and for all sufficiently large $n\in\mathbb N$, the unique maximum of the Ewens probability mass function $$ \mathbb L_n(k) = \frac{θ^k}{θ(θ+1)\ldots(θ+n-1)} \genfrac{[}{]}{0pt}{}{n}{k}, \quad k=1,\ldots,n, $$ is attained at $k= \left\lfloor θ\log n + \frac{θΓ'(θ)}{Γ(θ)} - \frac 12\right\rfloor$ or $k=\left\lceil θ\log n + \frac{θΓ'(θ)}{Γ(θ)} + \frac 12\right\rceil$. We prove that the mode is $$ k=\left\lfloor θ\log n - \frac{θΓ'(θ)}{Γ(θ)}\right\rfloor $$ for a set of $n$'s of asymptotic density $1$, yet this formula is not true for infinitely many $n$'s.

preprint2016arXiv

Some extensions of linear approximation and prediction problems for stationary processes

Let $(B(t))_{t\in Θ}$ with $Θ={\mathbb Z}$ or $Θ={\mathbb R}$ be a wide sense stationary process with discrete or continuous time. The classical linear prediction problem consists of finding an element in $\overline{span}\{B(s),s\le t\}$ providing the best possible mean square approximation to the variable $B(τ)$ with $τ>t$. In this article we investigate this and some other similar problems where, in addition to prediction quality, optimization takes into account other features of the objects we search for. One of the most motivating examples of this kind is an approximation of a stationary process $B$ by a stationary differentiable process $X$ taking into account the kinetic energy that $X$ spends in its approximation efforts.

preprint2016arXiv

Stochastic integral representations and classification of sum- and max-infinitely divisible processes

Introduced is the notion of minimality for spectral representations of sum- and max-infinitely divisible processes and it is shown that the minimal spectral representation on a Borel space exists and is unique. This fact is used to show that a stationary, stochastically continuous, sum- or max-i.d. random process on $\mathbb{R}^d$ can be generated by a measure-preserving flow on a $σ$-finite Borel measure space and that this flow is unique. This development makes it possible to extend the classification program of Rosiński (Ann. Probab. 23 (1995) 1163-1187) with a unified treatment of both sum- and max-infinitely divisible processes. As a particular case, a characterization of stationary, stochastically continuous, union-infinitely divisible random measurable subsets of $\mathbb{R}^d$ is obtained. Introduced and classified are several new max-i.d. random field models including fields of Penrose type and fields associated to Poisson line processes.

preprint2016arXiv

Weak convergence of renewal shot noise processes in the case of slowly varying normalization

We investigate weak convergence of finite-dimensional distributions of a renewal shot noise process $(Y(t))_{t\geq 0}$ with deterministic response function $h$ and the shots occurring at the times $0 = S_0 < S_1 < S_2<\ldots$, where $(S_n)$ is a random walk with i.i.d.\ jumps. There has been an outbreak of recent activity around this topic. We are interested in one out of few cases which remained open: $h$ is regularly varying at $\infty$ of index $-1/2$ and the integral of $h^2$ is infinite. Assuming that $S_1$ has a moment of order $r>2$ we use a strong approximation argument to show that the random fluctuations of $Y(s)$ occur on the scale $s=t+g(t,u)$ for $u\in [0,1]$, as $t\to\infty$, and, on the level of finite-dimensional distributions, are well approximated by the sum of a Brownian motion and a Gaussian process with independent values (the two processes being independent). The scaling function $g$ above depends on the slowly varying factor of $h$. If, for instance, $\lim_{t\to\infty}t^{1/2}h(t)\in (0,\infty)$, then $g(t,u)=t^u$.

preprint2015arXiv

A characterization of the normal distribution using stationary max-stable processes

Consider the max-stable process $η(t) = \max_{i\in\mathbb N} U_i \rm{e}^{\langle X_i, t\rangle - κ(t)}$, $t\in\mathbb{R}^d$, where $\{U_i, i\in\mathbb{N}\}$ are points of the Poisson process with intensity $u^{-2}\rm{d} u$ on $(0,\infty)$, $X_i$, $i\in\mathbb{N}$, are independent copies of a random $d$-variate vector $X$ (that are independent of the Poisson process), and $κ: \mathbb{R}^d \to \mathbb{R}$ is a function. We show that the process $η$ is stationary if and only if $X$ has multivariate normal distribution and $κ(t)-κ(0)$ is the cumulant generating function of $X$. In this case, $η$ is a max-stable process introduced by R. L. Smith.

preprint2015arXiv

A functional central limit theorem for branching random walks, almost sure weak convergence, and applications to random trees

Let $W_{\infty}(β)$ be the limit of the Biggins martingale $W_n(β)$ associated to a supercritical branching random walk with mean number of offspring $m$. We prove a functional central limit theorem stating that as $n\to\infty$ the process $$ D_n(u):= m^{\frac 12 n} \left(W_{\infty}\left(\frac{u}{\sqrt n}\right) - W_{n}\left(\frac{u}{\sqrt n}\right) \right) $$ converges weakly, on a suitable space of analytic functions, to a Gaussian random analytic function with random variance. Using this result we prove central limit theorems for the total path length of random trees. In the setting of binary search trees, we recover a recent result of R. Neininger [Refined Quicksort Asymptotics, Rand. Struct. and Alg., to appear], but we also prove a similar theorem for uniform random recursive trees. Moreover, we replace weak convergence in Neininger's theorem by the almost sure weak (a.s.w.) convergence of probability transition kernels. In the case of binary search trees, our result states that $$ L\left\{\sqrt{\frac{n}{2\log n}} \left(EPL_{\infty} - \frac{EPL_n-2n\log n}{n}\right)\Bigg | G_{n}\right\} \to \{ω\mapsto N_{0,1}\}, \quad \text{a.s.w.},$$ where $EPL_n$ is the external path length of a binary search tree $X_n$ with $n$ vertices, $EPL_{\infty}$ is the limit of the Régnier martingale, and $L(\,\cdot\, |G_n)$ denotes the conditional distribution w.r.t. the $σ$-algebra $G_n$ generated by $X_1,\ldots,X_n$. A.s.w. convergence is stronger than weak and even stable convergence. We prove several basic properties of the a.s.w. convergence and study a number of further examples in which the a.s.w. convergence appears naturally. These include the classical central limit theorem for Galton-Watson processes and the Pólya urn.

preprint2015arXiv

Max-stable processes and stationary systems of Lévy particles

We study stationary max-stable processes $\{η(t)\colon t\in\mathbb R\}$ admitting a representation of the form $η(t)=\max_{i\in\mathbb N}(U_i+ Y_i(t))$, where $\sum_{i=1}^{\infty} δ_{U_i}$ is a Poisson point process on $\mathbb R$ with intensity ${\rm e}^{-u} {\rm d} u$, and $Y_1,Y_2,\ldots$ are i.i.d.\ copies of a process $\{Y(t)\colon t\in\mathbb R\}$ obtained by running a Lévy process for positive $t$ and a dual Lévy process for negative $t$. We give a general construction of such Lévy-Brown-Resnick processes, where the restrictions of $Y$ to the positive and negative half-axes are Lévy processes with random birth and killing times. We show that these max-stable processes appear as limits of suitably normalized pointwise maxima of the form $M_n(t)=\max_{i=1,\ldots,n} ξ_i(s_n+t)$, where $ξ_1,ξ_2,\ldots$ are i.i.d.\ Lévy processes and $s_n$ is a sequence such that $s_n\sim c \log n$ with $c>0$. Also, we consider maxima of the form $\max_{i=1,\ldots,n} Z_i(t/\log n)$, where $Z_1,Z_2,\ldots$ are i.i.d.\ Ornstein--Uhlenbeck processes driven by an $α$-stable noise with skewness parameter $β=-1$. After a linear normalization, we again obtain limiting max-stable processes of the above form. This gives a generalization of the results of Brown and Resnick [Extreme values of independent stochastic processes, J.\ Appl.\ Probab., 14 (1977), pp.\ 732--739] to the totally skewed $α$-stable case.

preprint2015arXiv

Maxima of independent, non-identically distributed Gaussian vectors

Let $X_{i,n},n\in \mathbb{N},1\leq i\leq n$, be a triangular array of independent $\mathbb{R}^d$-valued Gaussian random vectors with correlation matrices $Σ_{i,n}$. We give necessary conditions under which the row-wise maxima converge to some max-stable distribution which generalizes the class of Hüsler-Reiss distributions. In the bivariate case, the conditions will also be sufficient. Using these results, new models for bivariate extremes are derived explicitly. Moreover, we define a new class of stationary, max-stable processes as max-mixtures of Brown-Resnick processes. As an application, we show that these processes realize a large set of extremal correlation functions, a natural dependence measure for max-stable processes. This set includes all functions $ψ(\sqrt{γ(h)}),h\in \mathbb{R}^d$, where $ψ$ is a completely monotone function and $γ$ is an arbitrary variogram.

preprint2014arXiv

Asymptotic distribution of complex zeros of random analytic functions

Let $ξ_0,ξ_1,\ldots$ be independent identically distributed complex- valued random variables such that $\mathbb{E}\log(1+|ξ_0|)<\infty$. We consider random analytic functions of the form \[\mathbf{G}_n(z)=\sum_{k=0}^{\infty}ξ_kf_{k,n}z^k,\] where $f_{k,n}$ are deterministic complex coefficients. Let $μ_n$ be the random measure counting the complex zeros of $\mathbf{G}_n$ according to their multiplicities. Assuming essentially that $-\frac{1}{n}\log f_{[tn],n}\to u(t)$ as $n\to\infty$, where $u(t)$ is some function, we show that the measure $\frac{1}{n}μ_n$ converges in probability to some deterministic measure $μ$ which is characterized in terms of the Legendre-Fenchel transform of $u$. The limiting measure $μ$ does not depend on the distribution of the $ξ_k$'s. This result is applied to several ensembles of random analytic functions including the ensembles corresponding to the three two-dimensional geometries of constant curvature. As another application, we prove a random polynomial analogue of the circular law for random matrices.

preprint2014arXiv

Generalized Random Energy Model at Complex Temperatures

Motivated by the Lee--Yang approach to phase transitions, we study the partition function of the Generalized Random Energy Model (GREM) at complex inverse temperature $β$. We compute the limiting log-partition function and describe the fluctuations of the partition function. For the GREM with $d$ levels, in total, there are $\frac 12 (d+1)(d+2)$ phases, each of which can symbolically be encoded as $G^{d_1}F^{d_2}E^{d_3}$ with $d_1,d_2,d_3\in\mathbb{N}_0$ such that $d_1+d_2+d_3=d$. In phase $G^{d_1}F^{d_2}E^{d_3}$, the first $d_1$ levels (counting from the root of the GREM tree) are in the glassy phase (G), the next $d_2$ levels are dominated by fluctuations (F), and the last $d_3$ levels are dominated by the expectation (E). Only the phases of the form $G^{d_1}E^{d_3}$ intersect the real $β$ axis. We describe the limiting distribution of the zeros of the partition function in the complex $β$ plane (= Fisher zeros). It turns out that the complex zeros densely touch the positive real axis at $d$ points at which the GREM is known to undergo phase transitions. Our results confirm rigorously and considerably extend the replica-method predictions from the physics literature.

preprint2014arXiv

Intrinsic volumes of Sobolev balls with applications to Brownian convex hulls

A formula due to Sudakov relates the first intrinsic volume of a convex set in a Hilbert space to the maximum of the isonormal Gaussian process over this set. Using this formula we compute the first intrinsic volumes of infinite-dimensional convex compact sets including unit balls with respect to Sobolev-type seminorms and ellipsoids in the Hilbert space. We relate the distribution of the random one-dimensional projections of these sets to the distributions $S_1,S_2,C_1,C_2$ studied by Biane, Pitman, Yor [Bull.\ AMS 38 (2001)]. We show that the $k$-th intrinsic volume of the set of all functions on $[0,1]$ which have Lipschitz constant bounded by $1$ and which vanish at $0$ (respectively, which have vanishing integral) is given by $$ V_k = \frac{π^{k/2}}{Γ\left(\frac 32 k +1 \right)}, \text{ respectively } V_k = \frac{π^{(k+1)/2}}{2Γ\left(\frac 32 k +\frac 32\right)}. $$ This is related to the results of Gao and Vitale [Discrete Comput. Geom. 26 (2001), Elect. Comm. Probab. 8 (2003)] who considered a similar question for functions with a restriction on the total variation instead of the Lipschitz constant. Using the results of Gao and Vitale we give a new proof of the formula for the expected volume of the convex hull of the $d$-dimensional Brownian motion which is due to Eldan [Elect. J. Probab., to appear]. Additionally, we prove an analogue of Eldan's result for the Brownian bridge. Similarly, we show that the results on the intrinsic volumes of the Lipschitz balls can be translated into formulae for the expected volumes of zonoids (Aumann integrals) generated by the Brownian motion and the Brownian bridge. Also, these results have discrete versions for Gaussian random walks and bridges. Our proofs exploit Sudakov's and Tsirelson's theorems which establish a connection between the intrinsic volumes and the isonormal Gaussian process.

preprint2014arXiv

Limiting distribution for the maximal standardized increment of a random walk

Let $X_1,X_2,...$ be independent identically distributed random variables with $\mathbb E X_k=0$, $\mathrm{Var} X_k=1$. Suppose that $φ(t):=\log \mathbb E e^{t X_k}<\infty$ for all $t>-σ_0$ and some $σ_0>0$. Let $S_k=X_1+...+X_k$ and $S_0=0$. We are interested in the limiting distribution of the multiscale scan statistic $$ M_n=\max_{0\leq i <j\leq n} \frac{S_j-S_i}{\sqrt{j-i}}. $$ We prove that for an appropriate normalizing sequence $a_n$, the random variable $M_n^2-a_n$ converges to the Gumbel extreme-value law $\exp\{-e^{-c x}\}$. The behavior of $M_n$ depends strongly on the distribution of the $X_k$'s. We distinguish between four cases. In the superlogarithmic case we assume that $φ(t)<t^2/2$ for every $t>0$. In this case, we show that the main contribution to $M_n$ comes from the intervals $(i,j)$ having length $l:=j-i$ of order $a(\log n)^{p}$, $a>0$, where $p=q/(q-2)$ and $q\in{3,4,...}$ is the order of the first non-vanishing cumulant of $X_1$ (not counting the variance). In the logarithmic case we assume that the function $ψ(t):=2φ(t)/t^2$ attains its maximum $m_*>1$ at some unique point $t=t_*\in (0,\infty)$. In this case, we show that the main contribution to $M_n$ comes from the intervals $(i,j)$ of length $d_*\log n+a\sqrt{\log n}$, $a\in\mathbb R$, where $d_*=1/φ(t_*)>0$. In the sublogarithmic case we assume that the tail of $X_k$ is heavier than $\exp\{-x^{2-\varepsilon}\}$, for some $\varepsilon>0$. In this case, the main contribution to $M_n$ comes from the intervals of length $o(\log n)$ and in fact, under regularity assumptions, from the intervals of length $1$. In the remaining, fourth case, the $X_k$'s are Gaussian. This case has been studied earlier in the literature. The main contribution comes from intervals of length $a\log n$, $a>0$. We argue that our results cover most interesting distributions with light tails.

preprint2013arXiv

Roots of random polynomials whose coefficients have logarithmic tails

It has been shown by Ibragimov and Zaporozhets [In Prokhorov and Contemporary Probability Theory (2013) Springer] that the complex roots of a random polynomial $G_n(z)=\sum_{k=0}^nξ_kz^k$ with i.i.d. coefficients $ξ_0,\ldots,ξ_n$ concentrate a.s. near the unit circle as $n\to\infty$ if and only if ${\mathbb{E}\log_+}|ξ_0|<\infty$. We study the transition from concentration to deconcentration of roots by considering coefficients with tails behaving like $L({\log}|t|)({\log}|t|)^{-α}$ as $t\to\infty$, where $α\geq0$, and $L$ is a slowly varying function. Under this assumption, the structure of complex and real roots of $G_n$ is described in terms of the least concave majorant of the Poisson point process on $[0,1]\times (0,\infty)$ with intensity $αv^{-(α+1)}\,du\,dv$.

preprint2012arXiv

Complex Random Energy Model: Zeros and Fluctuations

The partition function of the random energy model at inverse temperature $β$ is a sum of random exponentials $Z_N(β)=\sum_{k=1}^N \exp(β\sqrt{n} X_k)$, where $X_1,X_2,...$ are independent real standard normal random variables (= random energies), and $n=\log N$. We study the large $N$ limit of the partition function viewed as an analytic function of the complex variable $β$. We identify the asymptotic structure of complex zeros of the partition function confirming and extending predictions made in the theoretical physics literature. We prove limit theorems for the random partition function at complex $β$, both on the logarithmic scale and on the level of limiting distributions. Our results cover also the case of the sums of independent identically distributed random exponentials with any given correlations between the real and imaginary parts of the random exponent.

preprint2012arXiv

Critical points of random polynomials with independent identically distributed roots

Let $X_1,X_2,...$ be independent identically distributed random variables with values in $\C$. Denote by $μ$ the probability distribution of $X_1$. Consider a random polynomial $P_n(z)=(z-X_1)...(z-X_n)$. We prove a conjecture of Pemantle and Rivin [arXiv:1109.5975] that the empirical measure $μ_n:=\frac 1{n-1}\sum_{P_n'(z)=0} δ_z$ counting the complex zeros of the derivative $P_n'$ converges in probability to $μ$, as $n\to\infty$.

preprint2012arXiv

Distribution of levels in high-dimensional random landscapes

We prove empirical central limit theorems for the distribution of levels of various random fields defined on high-dimensional discrete structures as the dimension of the structure goes to $\infty$. The random fields considered include costs of assignments, weights of Hamiltonian cycles and spanning trees, energies of directed polymers, locations of particles in the branching random walk, as well as energies in the Sherrington--Kirkpatrick and Edwards--Anderson models. The distribution of levels in all models listed above is shown to be essentially the same as in a stationary Gaussian process with regularly varying nonsummable covariance function. This type of behavior is different from the Brownian bridge-type limit known for independent or stationary weakly dependent sequences of random variables.

preprint2012arXiv

Estimation of Huesler-Reiss distributions and Brown-Resnick processes

Estimation of extreme-value parameters from observations in the max-domain of attraction (MDA) of a multivariate max-stable distribution commonly uses aggregated data such as block maxima. Since we expect that additional information is contained in the non-aggregated, single "large" observations, we introduce a new approach of inference based on a multivariate peaks-over-threshold method. We show that for any process in the MDA of the frequently used Hüsler-Reiss model or its spatial extension, the Brown-Resnick process, suitably defined conditional increments asymptotically follow a multivariate Gaussian distribution. This leads to computationally efficient estimates of the Hüsler-Reiss parameter matrix. Further, the results enable parametric inference for Brown-Resnick processes. A simulation study compares the performance of the new estimators to other commonly used methods. As an application, we fit a non-isotropic Brown-Resnick process to the extremes of 12 year data of daily wind speed measurements.

preprint2012arXiv

Random determinants, mixed volumes of ellipsoids, and zeros of Gaussian random fields

Consider a $d\times d$ matrix $M$ whose rows are independent centered non-degenerate Gaussian vectors $ξ_1,...,ξ_d$ with covariance matrices $Σ_1,...,Σ_d$. Denote by $\mathcal{E}_i$ the location-dispersion ellipsoid of $ξ_i:\mathcal{E}_i={\mathbf{x}\in\mathbb{R}^d : \mathbf{x}^\topΣ_i^{-1} \mathbf{x}\leqslant1}$. We show that $$ \mathbb{E}\,|\det M|=\frac{d!}{(2π)^{d/2}}V_d(\mathcal{E}_1,...,\mathcal{E}_d), $$ where $V_d(\cdot,...,\cdot)$ denotes the {\it mixed volume}. We also generalize this result to the case of rectangular matrices. As a direct corollary we get an analytic expression for the mixed volume of $d$ arbitrary ellipsoids in $\mathbb{R}^d$. As another application, we consider a smooth centered non-degenerate Gaussian random field $X=(X_1,...,X_k)^\top:\mathbb{R}^d\to\mathbb{R}^k$. Using Kac-Rice formula, we obtain the geometric interpretation of the intensity of zeros of $X$ in terms of the mixed volume of location-dispersion ellipsoids of the gradients of $X_i/\sqrt{\mathbf{Var} X_i}$. This relates zero sets of equations to mixed volumes in a way which is reminiscent of the well-known Bernstein theorem about the number of solutions of the typical system of algebraic equations.

preprint2012arXiv

Random Marked Sets

We aim to link random fields and marked point processes and therefore introduce a new class of stochastic processes which are defined on a random set in R^d. Unlike for random fields, the mark covariance function of a marked random set is in general not positive definite. This implies that in many situations the use of simple geostatistical methods appears to be questionable. Surprisingly, for a special class of processes based on Gaussian random fields, we do have positive definiteness for the corresponding mark covariance function and mark correlation function.

preprint2012arXiv

Universality for zeros of random analytic functions

Let $ξ_0,ξ_1,...$ be independent identically distributed (i.i.d.) random variables such that $\E \log (1+|ξ_0|)<\infty$. We consider random analytic functions of the form $$ G_n(z)=\sum_{k=0}^{\infty} ξ_k f_{k,n} z^k, $$ where $f_{k,n}$ are deterministic complex coefficients. Let $ν_n$ be the random measure assigning the same weight $1/n$ to each complex zero of $G_n$. Assuming essentially that $-\frac 1n \log f_{[tn], n}\to u(t)$ as $n\to\infty$, where $u(t)$ is some function, we show that the measure $ν_n$ converges weakly to some deterministic measure which is characterized in terms of the Legendre--Fenchel transform of $u$. The limiting measure is universal, that is it does not depend on the distribution of the $ξ_k$'s. This result is applied to several ensembles of random analytic functions including the ensembles corresponding to the three two-dimensional geometries of constant curvature. As another application, we prove a random polynomial analogue of the circular law for random matrices.

preprint2011arXiv

Functional limit theorems for sums of independent geometric Lévy processes

Let $ξ_i$, $i\in \mathbb {N}$, be independent copies of a Lévy process $\{ξ(t),t\geq0\}$. Motivated by the results obtained previously in the context of the random energy model, we prove functional limit theorems for the process \[Z_N(t)=\sum_{i=1}^N\mathrm{e}^{ξ_i(s_N+t)}\] as $N\to\infty$, where $s_N$ is a non-negative sequence converging to $+\infty$. The limiting process depends heavily on the growth rate of the sequence $s_N$. If $s_N$ grows slowly in the sense that $\liminf_{N\to\infty}\log N/s_N>λ_2$ for some critical value $λ_2>0$, then the limit is an Ornstein--Uhlenbeck process. However, if $λ:=\lim_{N\to\infty}\log N/s_N\in(0,λ_2)$, then the limit is a certain completely asymmetric $α$-stable process $\mathbb {Y}_{α;ξ}$.

preprint2011arXiv

Persistence of competing systems of branching random walks

We consider a system of independent branching random walks on $\R$ which start off a Poisson point process with intensity of the form $e_λ(du)=e^{-λu}du$, where $λ\in\R$ is chosen in such a way that the overall intensity of particles is preserved. Denote by $χ$ the cluster distribution and let $ϕ$ be the log-Laplace transform of the intensity of $χ$. If $λϕ'(λ)>0$, we show that the system is persistent (stable) meaning that the point process formed by the particles in the $n$-th generation converges as $n\to\infty$ to a non-trivial point process $Π_{e_λ}^χ$ with intensity $e_λ$. If $λϕ'(λ)<0$, then the branching population suffers local extinction meaning that the limiting point process is empty. We characterize (generally, non-stationary) point processes on $\R$ which are cluster-invariant with respect to the cluster distribution $χ$ as mixtures of the point processes $Π_{ce_λ}^χ$ over $c>0$ and $λ\in K_{\text{st}}$, where $K_{\text{st}}=\{λ\in\R: ϕ(λ)=0, λϕ'(λ)>0\}$.

preprint2010arXiv

Extremes of the standardized Gaussian noise

Let $\{ξ_n, n\in\Z^d\}$ be a $d$-dimensional array of i.i.d. Gaussian random variables and define $\SSS(A)=\sum_{n\in A} ξ_n$, where $A$ is a finite subset of $\Z^d$. We prove that the appropriately normalized maximum of $\SSS(A)/\sqrt{|A|}$, where $A$ ranges over all discrete cubes or rectangles contained in $\{1,\ldots,n\}^d$, converges in the weak sense to the Gumbel extreme-value distribution as $n\to\infty$. We also prove continuous-time counterparts of these results.

preprint2010arXiv

Limit laws for sums of independent random products: the lattice case

Let $\{V_{i,j}; (i,j)\in\N^2\}$ be a two-dimensional array of i.i.d.\ random variables. The limit laws of the sum of independent random products $$ Z_n=\sum_{i=1}^{N_n} \prod_{j=1}^{n} e^{V_{i,j}} $$ as $n,N_n\to\infty$ have been investigated by a number of authors. Depending on the growth rate of $N_n$, the random variable $Z_n$ obeys a central limit theorem, or has limiting $α$-stable distribution. The latter result is true for non-lattice $V_{i,j}$ only. Our aim is to study the lattice case. We prove that although the (suitably normalized) sequence $Z_n$ fails to converge in distribution, it is relatively compact in the weak topology, and describe its cluster set. This set is a topological circle consisting of semi-stable distributions.

preprint2010arXiv

Locally adaptive image denoising by a statistical multiresolution criterion

We demonstrate how one can choose the smoothing parameter in image denoising by a statistical multiresolution criterion, both globally and locally. Using inhomogeneous diffusion and total variation regularization as examples for localized regularization schemes, we present an efficient method for locally adaptive image denoising. As expected, the smoothing parameter serves as an edge detector in this framework. Numerical examples illustrate the usefulness of our approach. We also present an application in confocal microscopy.

preprint2010arXiv

Stationary systems of Gaussian processes

We describe all countable particle systems on $\mathbb{R}$ which have the following three properties: independence, Gaussianity and stationarity. More precisely, we consider particles on the real line starting at the points of a Poisson point process with intensity measure $\mathfrak{m}$ and moving independently of each other according to the law of some Gaussian process $ξ$. We classify all pairs $(\mathfrak{m},ξ)$ generating a stationary particle system, obtaining three families of examples. In the first, trivial family, the measure $\mathfrak{m}$ is arbitrary, whereas the process $ξ$ is stationary. In the second family, the measure $\mathfrak{m}$ is a multiple of the Lebesgue measure, and $ξ$ is essentially a Gaussian stationary increment process with linear drift. In the third, most interesting family, the measure $\mathfrak{m}$ has a density of the form $αe^{-λx}$, where $α>0$, $λ\in\mathbb{R}$, whereas the process $ξ$ is of the form $ξ(t)=W(t)-λσ^2(t)/2+c$, where $W$ is a zero-mean Gaussian process with stationary increments, $σ^2(t)=\operatorname {Var}W(t)$, and $c\in\mathbb{R}$.