Source author record

Mathew D. Penrose

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

Catalog footprint

Research overview

22works
8topics
4close collaborators

Published work

22 publications

preprint2023arXiv

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

Let X1,X2,…X_1,X_2, \ldots be independent identically distributed random points in a convex polytopal domain A⊂RdA \subset \mathbb{R}^d. Define the largest nearest neighbour link LnL_n to be the smallest rr such that every point of Xn:={X1,…,Xn}\mathcal X_n:=\{X_1,\ldots,X_n\} has another such point within distance rr. We obtain a strong law of large numbers for LnL_n in the large-nn limit. A related threshold, the connectivity threshold MnM_n, is the smallest rr such that the random geometric graph G(Xn,r)G(\mathcal X_n, r) is connected. We show that as n→∞n \to \infty, almost surely nLnd/log⁡nnL_n^d/\log n tends to a limit that depends on the geometry of AA, and nMnd/log⁡nnM_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,ϕ)G(λ,s,ϕ), obtained by placing a Poisson(λs2λs^2) number of vertices uniformly at random in a square of side ss, with edges placed between each pair x,yx,y of vertices with probability ϕ(∥x−y∥)ϕ(\|x-y\|), where ϕ:R+→[0,1]ϕ: {\bf R}_+ \to [0,1] is a finite-range connection function. This paper is concerned with the asymptotic behaviour of the graph G(λ,s,ϕ)G(λ,s,ϕ) in the large-ss 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(ϕ)λ_c(ϕ).

preprint2022arXiv

Random Euclidean coverage from within

Let X1,X2,…X_1,X_2, \ldots be independent random uniform points in a bounded domain A⊂RdA \subset \mathbb{R}^d with smooth boundary. Define the coverage threshold RnR_n to be the smallest rr such that AA is covered by the balls of radius rr centred on X1,…,XnX_1,\ldots,X_n. We obtain the limiting distribution of RnR_n and also a strong law of large numbers for RnR_n in the large-nn limit. For example, if AA has volume 1 and perimeter ∣∂A∣|\partial A|, if d=3d=3 then Pr⁡[nπRn3−log⁡n−2log⁡(log⁡n)≤x]\Pr[nπR_n^3 - \log n - 2 \log (\log n) \leq x] converges to exp⁡(−2−4π5/3∣∂A∣e−2x/3)\exp(-2^{-4}π^{5/3} |\partial A| e^{-2 x/3}) and (nπRn3)/(log⁡n)→1(n πR_n^3)/(\log n) \to 1 almost surely, and if d=2d=2 then Pr⁡[nπRn2−log⁡n−log⁡(log⁡n)≤x]\Pr[n πR_n^2 - \log n - \log (\log n) \leq x] converges to exp⁡(−e−x−∣∂A∣π−1/2e−x/2)\exp(- e^{-x}- |\partial A|π^{-1/2} e^{-x/2}). We give similar results for general dd, and also for the case where AA 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 AA be uniform.

preprint2020arXiv

Leaves on the line and in the plane

The Dead Leaves Model (DLM) provides a random tessellation of dd-space, representing the visible portions of fallen leaves on the ground when d=2d=2. For d=1d=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=2d=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,rn)G(n,r_n), nn vertices are placed randomly in Euclidean dd-space and edges are added between any pair of vertices distant at most rnr_n from each other. We establish strong laws of large numbers (LLNs) for a large class of graph parameters, evaluated for G(n,rn)G(n,r_n) in the thermodynamic limit with nrnd=nr_n^d = const., and also in the dense limit with nrnd→∞n r_n^d \to \infty, rn→0r_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−εd-\varepsilon, under thermodynamic scaling of the distance parameter.

preprint2016arXiv

Connectivity of soft random geometric graphs

Consider a graph on nn uniform random points in the unit square, each pair being connected by an edge with probability pp if the inter-point distance is at most rr. We show that as n→∞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,rp,r. We determine the asymptotic probability of connectivity for all (pn,rn)(p_n,r_n) subject to rn=O(n−ε)r_n=O(n^{-\varepsilon}), some ε>0\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 11 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 pp; otherwise allow diagonal connections between white squares. We show that there is a critical surface of pairs (λ,p)(λ, p), containing the pair (1,0.5)(1,0.5), such that for (λ,p)(λ, 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 nn vertices with independent identically distributed outdegrees with distribution FF having mean μμ, and destinations of arcs selected uniformly at random. We show that if μ>1μ>1 then for large nn 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 FF and the other with Poisson offspring distribution with mean μμ. If μ≤1μ\leq 1 there is very likely to be no giant strong component. We also extend this to allow for FF varying with nn.

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 dd-space, with distance parameter rr and intensities λ,μλ,μ. We show for d≥2d \geq 2 that if λλ is supercritical for the one-type random geometric graph with distance parameter 2r2r, there exists μμ such that (λ,μ)(λ,μ) is supercritical (this was previously known for d=2d=2). For d=2d=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 kk-dimensional flats in Rd\R^d are discussed.

preprint2013arXiv

Limit theory for point processes in manifolds

Let Yi,i≥1Y_i,i\geq1, be i.i.d. random variables having values in an mm-dimensional manifold M⊂Rd\mathcal {M}\subset \mathbb{R}^d and consider sums ∑i=1nξ(n1/mYi,{n1/mYj}j=1n)\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,Y)(y,\mathcal {Y}), with y∈Rdy\in \mathbb{R}^d and Y⊂Rd\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 mm-dimensional hyperplanes tangent to M\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 {Yi}i=1n\{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 SS of subsets of Rd\mathbb{R}^d and of point sets that are (almost) subadditive in their first variable. Denoting by ξξ the random parking measure in Rd\mathbb{R}^d, and by ξRξ^R the random parking measure in the cube QR=(−R,R)dQ_R=(-R,R)^d, we show, under some natural assumptions on SS, that there exists a constant Sˉ∈R\bar{S}\in \mathbb{R} such that % lim⁡R→+∞S(QR,ξ)∣QR∣ = lim⁡R→+∞S(QR,ξR)∣QR∣ = Sˉ \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 ζ↦S(QR,ζ)ζ\mapsto S(Q_R,ζ) is the counting measure of ζζ in QRQ_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 MM be a random m×nm \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)N(n,m) denote the number of left null vectors in 0,1m{0,1}^m for MM (including the zero vector), where addition is mod 2. We take n,m→∞n, m \to \infty, with m/n→α>0m/n \to α> 0, while the weight distribution may vary with nn but converges weakly to a limiting distribution on 3,4,5,...{3, 4, 5, ...}; let WW denote a variable with this limiting distribution. Identifying MM with a hypergraph on nn vertices, we define the 2-core of MM 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 WW. Threshold α∗α^* marks the infimum of values of αα at which n−1log⁡E[N(n,m)]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 nn having more rows than non-empty columns. We have 1/2≤α∗≤α‾≤11/2 \leq α^* \leq \underlineα \leq 1, and typically these inequalities are strict; for example when W=3W = 3 almost surely, numerics give α∗=0.88949...α^* = 0.88949 ... and α‾=0.91793...\underlineα = 0.91793 ... (previous work on this model has mainly been concerned with such cases where WW is non-random). The threshold of values of αα for which N(n,m)≥2N(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

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+\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 dd-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>1d > 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 pcsite≥pcbondp_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 pcsite>pcbondp_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 dd-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 Rd\R^d, d≥2d \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(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=2d=2, the limit is expressed in terms of Dickman-type variables.

preprint2007arXiv

Multivariate normal approximation in geometric probability

Consider a measure μλ=∑xξxδxμ_λ= \sum_x ξ_x δ_x where the sum is over points xx of a Poisson point process of intensity λλ on a bounded region in dd-space, and ξxξ_x is a functional determined by the Poisson points near to xx, 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 RdR^d are asymptotically independent normals as λ→∞λ\to \infty; here we give an O(λ−1/(2d+ε))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.

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.