Source author record

Mathew D. Penrose

Mathew D. Penrose 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

22works
8topics
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

22 published item(s)

preprint2023arXiv

Largest nearest-neighbour link and connectivity threshold in a polytopal random sample

Let $X_1,X_2, \ldots $ be independent identically distributed random points in a convex polytopal domain $A \subset \mathbb{R}^d$. Define the largest nearest neighbour link $L_n$ to be the smallest $r$ such that every point of $\mathcal X_n:=\{X_1,\ldots,X_n\}$ has another such point within distance $r$. We obtain a strong law of large numbers for $L_n$ in the large-$n$ limit. A related threshold, the connectivity threshold $M_n$, is the smallest $r$ such that the random geometric graph $G(\mathcal X_n, r)$ is connected. We show that as $n \to \infty$, almost surely $nL_n^d/\log n$ tends to a limit that depends on the geometry of $A$, and $nM_n^d/\log n$ tends to the same limit.

preprint2022arXiv

Giant component of the soft random geometric graph

Consider a 2-dimensional soft random geometric graph $G(λ,s,ϕ)$, obtained by placing a Poisson($λs^2$) number of vertices uniformly at random in a square of side $s$, with edges placed between each pair $x,y$ of vertices with probability $ϕ(\|x-y\|)$, where $ϕ: {\bf R}_+ \to [0,1]$ is a finite-range connection function. This paper is concerned with the asymptotic behaviour of the graph $G(λ,s,ϕ)$ in the large-$s$ limit with $(λ,ϕ)$ fixed. We prove that the proportion of vertices in the largest component converges in probability to the percolation probability for the corresponding random connection model, which is a random graph defined similarly for a Poisson process on the whole plane. We do not cover the case where $λ$ equals the critical value $λ_c(ϕ)$.

preprint2022arXiv

Random Euclidean coverage from within

Let $X_1,X_2, \ldots $ be independent random uniform points in a bounded domain $A \subset \mathbb{R}^d$ with smooth boundary. Define the coverage threshold $R_n$ to be the smallest $r$ such that $A$ is covered by the balls of radius $r$ centred on $X_1,\ldots,X_n$. We obtain the limiting distribution of $R_n$ and also a strong law of large numbers for $R_n$ in the large-$n$ limit. For example, if $A$ has volume 1 and perimeter $|\partial A|$, if $d=3$ then $\Pr[nπR_n^3 - \log n - 2 \log (\log n) \leq x]$ converges to $\exp(-2^{-4}π^{5/3} |\partial A| e^{-2 x/3})$ and $(n πR_n^3)/(\log n) \to 1$ almost surely, and if $d=2$ then $\Pr[n πR_n^2 - \log n - \log (\log n) \leq x]$ converges to $\exp(- e^{-x}- |\partial A|π^{-1/2} e^{-x/2})$. We give similar results for general $d$, and also for the case where $A$ is a polytope. We also generalize to allow for multiple coverage. The analysis relies on classical results by Hall and by Janson, along with a careful treatment of boundary effects. For the strong laws of large numbers, we can relax the requirement that the underlying density on $A$ be uniform.

preprint2020arXiv

Leaves on the line and in the plane

The Dead Leaves Model (DLM) provides a random tessellation of $d$-space, representing the visible portions of fallen leaves on the ground when $d=2$. For $d=1$, we establish formulae for the intensity, two-point correlations, and asymptotic covariances for the point process of cell boundaries, along with a functional CLT. For $d=2$ we establish analogous results for the random surface measure of cell boundaries, and also determine the intensity of cells in a more general setting than in earlier work of Cowan and Tsang. We introduce a general notion of Dead Leaves Random Measures and give formulae for means, asymptotic variances and functional CLTs for these measures; this has applications to various other quantities associated with the DLM.

preprint2020arXiv

Limit theory of combinatorial optimization for random geometric graphs

In the random geometric graph $G(n,r_n)$, $n$ vertices are placed randomly in Euclidean $d$-space and edges are added between any pair of vertices distant at most $r_n$ from each other. We establish strong laws of large numbers (LLNs) for a large class of graph parameters, evaluated for $G(n,r_n)$ in the thermodynamic limit with $nr_n^d =$ const., and also in the dense limit with $n r_n^d \to \infty$, $r_n \to 0$. Examples include domination number, independence number, clique-covering number, eternal domination number and triangle packing number. The general theory is based on certain subadditivity and superadditivity properties, and also yields LLNs for other functionals such as the minimum weight for the travelling salesman, spanning tree, matching, bipartite matching and bipartite travelling salesman problems, for a general class of weight functions with at most polynomial growth of order $d-\varepsilon$, under thermodynamic scaling of the distance parameter.

preprint2016arXiv

Connectivity of soft random geometric graphs

Consider a graph on $n$ uniform random points in the unit square, each pair being connected by an edge with probability $p$ if the inter-point distance is at most $r$. We show that as $n\to\infty$ the probability of full connectivity is governed by that of having no isolated vertices, itself governed by a Poisson approximation for the number of isolated vertices, uniformly over all choices of $p,r$. We determine the asymptotic probability of connectivity for all $(p_n,r_n)$ subject to $r_n=O(n^{-\varepsilon})$, some $\varepsilon >0$. We generalize the first result to higher dimensions and to a larger class of connection probability functions.

preprint2016arXiv

Percolation of even sites for enhanced random sequential adsorption

Consider random sequential adsorption on a chequerboard lattice with arrivals at rate $1$ on light squares and at rate $λ$ on dark squares. Ultimately, each square is either occupied, or blocked by an occupied neighbour. Colour the occupied dark squares and blocked light sites {\em black}, and the remaining squares {\em white}. Independently at each meeting-point of four squares, allow diagonal connections between black squares with probability $p$; otherwise allow diagonal connections between white squares. We show that there is a critical surface of pairs $(λ, p)$, containing the pair $(1,0.5)$, such that for $(λ, p)$ lying above (respectively, below) the critical surface the black (resp. white) phase percolates, and on the critical surface neither phase percolates.

preprint2015arXiv

The strong giant in a random digraph

Consider a random directed graph on $n$ vertices with independent identically distributed outdegrees with distribution $F$ having mean $μ$, and destinations of arcs selected uniformly at random. We show that if $μ>1$ then for large $n$ there is very likely to be a unique giant strong component with proportionate size given as the product of two branching process survival probabilities, one with offspring distribution $F$ and the other with Poisson offspring distribution with mean $μ$. If $μ\leq 1$ there is very likely to be no giant strong component. We also extend this to allow for $F$ varying with $n$.

preprint2014arXiv

Continuum AB percolation and AB random geometric graphs

Consider a bipartite random geometric graph on the union of two independent homogeneous Poisson point processes in $d$-space, with distance parameter $r$ and intensities $λ,μ$. We show for $d \geq 2$ that if $λ$ is supercritical for the one-type random geometric graph with distance parameter $2r$, there exists $μ$ such that $(λ,μ)$ is supercritical (this was previously known for $d=2$). For $d=2$ we also consider the restriction of this graph to points in the unit square. Taking $μ= τλ$ for fixed $τ$, we give a strong law of large numbers as $λ\to \infty$, for the connectivity threshold of this graph.

preprint2014arXiv

Moments and central limit theorems for some multivariate Poisson functionals

This paper deals with Poisson processes on an arbitrary measurable space. Using a direct approach, we derive formulae for moments and cumulants of a vector of multiple Wiener-Itô integrals with respect to the compensated Poisson process. Second, a multivariate central limit theorem is shown for a vector whose components admit a finite chaos expansion of the type of a Poisson U-statistic. The approach is based on recent results of Peccati et al.\ combining Malliavin calculus and Stein's method, and also yields Berry-Esseen type bounds. As applications, moment formulae and central limit theorems for general geometric functionals of intersection processes associated with a stationary Poisson process of $k$-dimensional flats in $\R^d$ are discussed.

preprint2013arXiv

Limit theory for point processes in manifolds

Let $Y_i,i\geq1$, be i.i.d. random variables having values in an $m$-dimensional manifold $\mathcal {M}\subset \mathbb{R}^d$ and consider sums $\sum_{i=1}^nξ(n^{1/m}Y_i,\{n^{1/m}Y_j\}_{j=1}^n)$, where $ξ$ is a real valued function defined on pairs $(y,\mathcal {Y})$, with $y\in \mathbb{R}^d$ and $\mathcal {Y}\subset \mathbb{R}^d$ locally finite. Subject to $ξ$ satisfying a weak spatial dependence and continuity condition, we show that such sums satisfy weak laws of large numbers, variance asymptotics and central limit theorems. We show that the limit behavior is controlled by the value of $ξ$ on homogeneous Poisson point processes on $m$-dimensional hyperplanes tangent to $\mathcal {M}$. We apply the general results to establish the limit theory of dimension and volume content estimators, Rényi and Shannon entropy estimators and clique counts in the Vietoris-Rips complex on $\{Y_i\}_{i=1}^n$.

preprint2012arXiv

Random parking, Euclidean functionals, and rubber elasticity

We study subadditive functions of the random parking model previously analyzed by the second author. In particular, we consider local functions $S$ of subsets of $\mathbb{R}^d$ and of point sets that are (almost) subadditive in their first variable. Denoting by $ξ$ the random parking measure in $\mathbb{R}^d$, and by $ξ^R$ the random parking measure in the cube $Q_R=(-R,R)^d$, we show, under some natural assumptions on $S$, that there exists a constant $\bar{S}\in \mathbb{R}$ such that % $$ \lim_{R\to +\infty} \frac{S(Q_R,ξ)}{|Q_R|}\,=\,\lim_{R\to +\infty}\frac{S(Q_R,ξ^R)}{|Q_R|}\,=\,\bar{S} $$ % almost surely. If $ζ\mapsto S(Q_R,ζ)$ is the counting measure of $ζ$ in $Q_R$, then we retrieve the result by the second author on the existence of the jamming limit. The present work generalizes this result to a wide class of (almost) subadditive functions. In particular, classical Euclidean optimization problems as well as the discrete model for rubber previously studied by Alicandro, Cicalese, and the first author enter this class of functions. In the case of rubber elasticity, this yields an approximation result for the continuous energy density associated with the discrete model at the thermodynamic limit, as well as a generalization to stochastic networks generated on bounded sets.

preprint2012arXiv

Rank deficiency in sparse random GF[2] matrices

Let $M$ be a random $m \times n$ matrix with binary entries and i.i.d. rows. The weight (i.e., number of ones) of a row has a specified probability distribution, with the row chosen uniformly at random given its weight. Let $N(n,m)$ denote the number of left null vectors in ${0,1}^m$ for $M$ (including the zero vector), where addition is mod 2. We take $n, m \to \infty$, with $m/n \to α> 0$, while the weight distribution may vary with $n$ but converges weakly to a limiting distribution on ${3, 4, 5, ...}$; let $W$ denote a variable with this limiting distribution. Identifying $M$ with a hypergraph on $n$ vertices, we define the 2-core of $M$ as the terminal state of an iterative algorithm that deletes every row incident to a column of degree 1. We identify two thresholds $α^*$ and $\underlineα$, and describe them analytically in terms of the distribution of $W$. Threshold $α^*$ marks the infimum of values of $α$ at which $n^{-1} \log{\mathbb{E} [N(n,m)}]$ converges to a positive limit, while $\underlineα$ marks the infimum of values of $α$ at which there is a 2-core of non-negligible size compared to $n$ having more rows than non-empty columns. We have $1/2 \leq α^* \leq \underlineα \leq 1$, and typically these inequalities are strict; for example when $W = 3$ almost surely, numerics give $α^* = 0.88949 ...$ and $\underlineα = 0.91793 ...$ (previous work on this model has mainly been concerned with such cases where $W$ is non-random). The threshold of values of $α$ for which $N(n,m) \geq 2$ in probability lies in $[α^*,\underlineα]$ and is conjectured to equal $\underlineα$. The random row weight setting gives rise to interesting new phenomena not present in the non-random case that has been the focus of previous work.

preprint2011arXiv

Local central limit theorems in stochastic geometry

We give a general local central limit theorem for the sum of two independent random variables, one of which satisfies a central limit theorem while the other satisfies a local central limit theorem with the same order variance. We apply this result to various quantities arising in stochastic geometry, including: size of the largest component for percolation on a box; number of components, number of edges, or number of isolated points, for random geometric graphs; covered volume for germ-grain coverage models; number of accepted points for finite-input random sequential adsorption; sum of nearest-neighbour distances for a random sample from a continuous multidimensional distribution.

preprint2010arXiv

Asymptotic normality of maximum likelihood estimator for cooperative sequential adsorption

We have shown in previous work that statistical inference for cooperative sequential adsorption model can be based on maximum likelihood estimation. In this paper we continue this research and establish asymptotic normality of the maximum likelihood estimator in thermodynamic limit. We also perform and discuss some numerical simulations of the model.

preprint2010arXiv

Martingale representation for Poisson processes with applications to minimal variance hedging

We consider a Poisson process $η$ on a measurable space $(\BY,\mathcal{Y})$ equipped with a partial ordering, assumed to be strict almost everwhwere with respect to the intensity measure $λ$ of $η$. We give a Clark-Ocone type formula providing an explicit representation of square integrable martingales (defined with respect to the natural filtration associated with $η$), which was previously known only in the special case, when $λ$ is the product of Lebesgue measure on $\R_+$ and a $σ$-finite measure on another space $\BX$. Our proof is new and based on only a few basic properties of Poisson processes and stochastic integrals. We also consider the more general case of an independent random measure in the sense of Itô of pure jump type and show that the Clark-Ocone type representation leads to an explicit version of the Kunita-Watanabe decomposition of square integrable martingales. We also find the explicit minimal variance hedge in a quite general financial market driven by an independent random measure.

preprint2010arXiv

Percolation and limit theory for the Poisson lilypond model

The lilypond model on a point process in $d$-space is a growth-maximal system of non-overlapping balls centred at the points. We establish central limit theorems for the total volume and the number of components of the lilypond model on a sequence of Poisson or binomial point processes on expanding windows. For the lilypond model over a homogeneous Poisson process, we give subexponentially decaying tail bounds for the size of the cluster at the origin. Finally, we consider the enhanced Poisson lilypond model where all the balls are enlarged by a fixed amount (the enhancement parameter), and show that for $d > 1 $ the critical value of this parameter, above which the enhanced model percolates, is strictly positive.

preprint2010arXiv

Strict inequalities of critical probabilities on Gilbert's continuum percolation graph

Any infinite graph has site and bond percolation critical probabilities satisfying $p_c^{site}\geq p_c^{bond}$. The strict version of this inequality holds for many, but not all, infinite graphs. In this paper, the class of graphs for which the strict inequality holds is extended to a continuum percolation model. In Gilbert's graph with supercritical density on the Euclidean plane, there is almost surely a unique infinite connected component. We show that on this component $p_c^{site} > p_c^{bond}$. This also holds in higher dimensions.

preprint2009arXiv

Limit theorems for random spatial drainage networks

Suppose that under the action of gravity, liquid drains through the unit $d$-cube via a minimal-length network of channels constrained to pass through random sites and to flow with nonnegative component in one of the canonical orthogonal basis directions of $\R^d$, $d \geq 2$. The resulting network is a version of the so-called minimal directed spanning tree. We give laws of large numbers and convergence in distribution results on the large-sample asymptotic behaviour of the total power-weighted edge-length of the network on uniform random points in $(0,1)^d$. The distributional results exhibit a weight-dependent phase transition between Gaussian and boundary-effect-derived distributions. These boundary contributions are characterized in terms of limits of the so-called on-line nearest-neighbour graph, a natural model of spatial network evolution, for which we also present some new results. Also, we give a convergence in distribution result for the length of the longest edge in the drainage network; when $d=2$, the limit is expressed in terms of Dickman-type variables.

preprint2007arXiv

Multivariate normal approximation in geometric probability

Consider a measure $μ_λ= \sum_x ξ_x δ_x$ where the sum is over points $x$ of a Poisson point process of intensity $λ$ on a bounded region in $d$-space, and $ξ_x$ is a functional determined by the Poisson points near to $x$, i.e. satisfying an exponential stabilization condition, along with a moments condition (examples include statistics for proximity graphs, germ-grain models and random sequential deposition models). A known general result says the $μ_λ$-measures (suitably scaled and centred) of disjoint sets in $R^d$ are asymptotically independent normals as $λ\to \infty$; here we give an $O(λ^{-1/(2d + ε)})$ bound on the rate of convergence. We illustrate our result with an explicit multivariate central limit theorem for the nearest-neighbour graph on Poisson points on a finite collection of disjoint intervals.