Portrait of Matthias Irlbeck

Matthias Irlbeck

I am a postdoctoral researcher in mathematics at the University of Hamburg working together with Jack Hanson. Before that, I completed my PhD at the University of Groningen with Tobias Müller. My research is in discrete probability theory, with a focus on percolation, random graphs, and stochastic geometry.

Research

Details and animationClose details: On the shape of the typical Poisson-Voronoi cell in high dimensions

Model

A Poisson-Voronoi tessellation assigns each location in Rd\mathbb{R}^d to its nearest site in a homogeneous Poisson point process of intensity λ>0\lambda>0. To obtain the typical cell VV, add a site at the origin and take its Voronoi cell. This describes the tessellation as seen from a typical site; sampling a location in space would favour larger cells. We study how the shape of this random convex polytope changes as the dimension grows.

Main result

Let rr and RR be the radii of the largest inscribed and smallest enclosing balls, allowing both centres to vary. Write DD for the diameter and w‾\overline{w} for the mean width, the distance between parallel supporting hyperplanes averaged uniformly over directions in Sd−1S^{d-1}.

Among other geometric properties, we prove that

rR→d→∞P12,\frac{r}{R}\xrightarrow[d\to\infty]{\mathbb{P}}\frac12, DR→d→∞P2,\frac{D}{R}\xrightarrow[d\to\infty]{\mathbb{P}}2, w‾R→d→∞P2.\frac{\overline{w}}{R}\xrightarrow[d\to\infty]{\mathbb{P}}2.

The joint distribution of these ratios is independent of Poisson intensity, so the limits hold for any positive sequence λ=λ(d)\lambda=\lambda(d).

Illustration

Poisson-Voronoi cell shape

All cells that intersect a fixed ball are shown here.

Interpretation

The diameter and mean width are each asymptotic to 2R2R, while the largest inscribed ball has only about half the enclosing ball’s radius. Uniformly over all vertices, their distances from the origin divided by RR tend to one in probability.

Yet the whole cell remains far from spherical: with probability tending to one, its Hausdorff distance from every ball is at least a fixed positive fraction of DD.

The facets also vary greatly in size. For every fixed ε>0\varepsilon>0, the proportion, by number, of facets with diameter larger than εD\varepsilon D tends to zero in probability. Nevertheless, with probability tending to one, some facets have diameter at least cDcD for a fixed c>0c>0. The preprint develops this picture further for faces of other dimensions.

Details and animationClose details: Thresholds for colouring the random Borsuk graph

Model

Choose nn independent uniform points on the unit sphere SdS^d in Rd+1\mathbb{R}^{d+1}. For a small positive parameter α=α(n)\alpha=\alpha(n), join two points when their angular distance exceeds π−α\pi-\alpha, so edges connect nearly antipodal points. Increasing α\alpha adds edges and makes proper colouring harder: adjacent vertices must receive different colours. For fixed d≥2d\geq2, we study when kk colours cease to suffice as nn tends to infinity.

Main result

For every 2≤k≤d2\leq k\leq d, we locate the loss of kk-colourability between positive constant multiples of n−1/dn^{-1/d}, where the expected degree is of constant order. The later transition to requiring d+2d+2 colours occurs at order (log⁡n/n)1/d(\log n/n)^{1/d}, where the expected degree is logarithmic, as previously shown by Kahle and Martinez-Figueroa.

For two colours, we identify the sharp threshold

α2(n)=c2(d) n−1/d.\alpha_2(n)=c_2(d)\,n^{-1/d}.

Below this threshold by any fixed positive proportion, the graph is two-colourable with probability tending to one. Above it by any fixed positive proportion, it contains an odd cycle with probability tending to one.

For 3≤k≤d+13\leq k\leq d+1, we construct threshold sequences αk(n)\alpha_k(n) and a common density-one set of sample sizes. For every fixed positive relative margin, kk colours suffice below the threshold and fail above it with probability tending to one as nn grows within this set. The proportion of excluded sizes tends to zero.

Illustration

Random Borsuk graph

Example for a random Borsuk graph with n=32n=32 on S2S^2 and chromatic number 33.

Interpretation

The constant comes from continuum AB percolation: two independent Poisson point clouds in Rd\mathbb{R}^d, each of intensity λ\lambda, with edges between opposite types at Euclidean distance at most one. If λAB(d)\lambda_{\mathrm{AB}}(d) is its critical intensity and σd\sigma_d is the surface area of SdS^d, then

c2(d)=(σd λAB(d))1/d.c_2(d)=\bigl(\sigma_d\,\lambda_{\mathrm{AB}}(d)\bigr)^{1/d}.

Applying the antipodal map x↦−xx\mapsto-x to the points in one of two small antipodal caps and rescaling gives a local approximation by this model. The proof joins local connections above its percolation threshold to create a global odd cycle on the sphere, the obstruction to two-colouring.

Details and animationClose details: Non-vanishing uniqueness threshold for hyperbolic Poisson-Voronoi percolation in dimension at least three

Model

Start with a homogeneous Poisson point process of intensity λ>0\lambda>0 with respect to hyperbolic volume in Hd\mathbb{H}^d. Assign each location to its nearest site, and colour each resulting Voronoi cell black independently with probability pp. Points belong to the same black cluster when a continuous path through black cells joins them.

Unlike in Euclidean space, infinitely many unbounded black clusters can coexist. We study the uniqueness threshold pu(λ)p_u(\lambda): the infimum of pp for which exactly one unbounded black cluster exists with positive probability, equivalently almost surely in this model.

Main result

For every fixed dimension d≥3d\geq3, we prove

pu(λ)≥cd>0for every λ>0.\begin{gathered} p_u(\lambda)\geq c_d>0\\[4pt] \text{for every }\lambda>0. \end{gathered}

The constant depends only on the dimension. Thus a sufficiently small black-cell probability cannot produce a unique unbounded cluster, however sparse the tessellation becomes. This answers a question of Grebík and Recke.

Illustration

Hyperbolic Poisson-Voronoi tessellation

One finite cluster of black cells in the Poincaré half-space model

Interpretation

The ambient geometry makes a decisive difference. D’Achille, Grebík, Khezeli, Recke and Wilkens show that the uniqueness threshold tends to zero as the site intensity tends to zero on the Riemannian product H2×H2\mathbb{H}^2\times\mathbb{H}^2. Our result rules out this behaviour in Hd\mathbb{H}^d for d≥3d\geq3.

The proof bounds the probability that two distant points belong to the same black cluster. At low intensities, it extracts reduced paths from black-cell connections and bounds their expected counts for paths first reaching a distant region. The comparison uses paths with no guaranteed edge between nonconsecutive vertices in an auxiliary model whose edge choices are independent once point locations are fixed. A separate estimate handles larger intensities. Together these estimates rule out uniqueness at small pp.

Details and animationClose details: Poisson-Voronoi percolation in high dimensions

Model

Start with a homogeneous Poisson point process of intensity λ>0\lambda>0 in Rd\mathbb{R}^d and assign each location to its nearest site. Colour the resulting Voronoi cells black independently with probability pp. Cells sharing a facet are neighbours; connected black cells form clusters.

We study the critical probability pc(d)p_c(d), the infimum of pp for which an unbounded black cluster exists with positive probability. Rescaling space preserves cell adjacencies, so the threshold depends on the dimension but not on the Poisson intensity.

Main result

As d→∞d\to\infty, we determine its precise asymptotic behaviour:

pc(d)=ed 2d(1+o(1)).p_c(d)=\frac{e}{d\,2^d}\bigl(1+o(1)\bigr).

Writing mdm_d for the expected number of neighbours of a typical cell, sampled by its site, we also prove

md=d 2de(1+o(1)),pc(d) md⟶1.\begin{gathered} m_d=\frac{d\,2^d}{e}\bigl(1+o(1)\bigr),\\[4pt] p_c(d)\,m_d\longrightarrow1. \end{gathered}

Thus the critical black fraction is asymptotically the reciprocal of the mean number of neighbouring cells.

Illustration

Voronoi percolation cluster

One finite cluster of black cells for d=3d=3.

Interpretation

Exploring a black cluster means repeatedly discovering new black neighbours. A branching-process approximation treats those discoveries as independent and predicts survival when their mean exceeds one, suggesting the threshold p md≈1p\,m_d\approx1.

For every fixed ε∈(0,1)\varepsilon\in(0,1) and all sufficiently large dd, there is almost surely no unbounded black cluster at p=(1−ε)/mdp=(1-\varepsilon)/m_d, whereas an unbounded cluster exists with positive probability at p=(1+ε)/mdp=(1+\varepsilon)/m_d.

Making this prediction rigorous requires handling unbounded random degrees and spatial dependencies. For the upper bound, we explore only selected connections, couple these finite explorations to branching random walks, and join them through comparison with site percolation on Z2\mathbb{Z}^2. For the lower bound, we count paths in a more generous adjacency model that includes every actual Voronoi connection. The complete proof appears in Chapter 3 of my PhD thesis.

Short CV

  1. October 2026–present

    Postdoctoral researcher

    University of Hamburg

    Advisor: Jack Hanson
  2. 2022–2026

    PhD in Mathematics

    University of Groningen

  3. 2019–2021

    Master’s degree in (Financial) Mathematics

    LMU Munich

  4. 2016–2019

    Bachelor’s degree in (Financial) Mathematics

    LMU Munich

Contact

Feel free to reach out to me!

Email
Office
Room 908, Geomatikum
Bundesstraße 55, 20146 Hamburg, Germany