Multivariate Normal Approximation for functionals of random polytopes

TL;DR

This paper establishes central limit theorems for various geometric functionals of random polytopes formed by Poisson point processes in convex bodies, including intrinsic volumes and the f-vector components.

Contribution

It provides the first multivariate normal approximation results for multiple geometric functionals of random polytopes, including the Will's functional and the oracle estimator.

Findings

Proves CLTs for intrinsic volumes and f-vector components.

Derives a multivariate CLT for several geometric functionals.

Establishes asymptotic normality for the volume estimator.

Abstract

Consider the random polytope, that is given by the convex hull of a Poisson point process on a smooth convex body in $\mathbb{R}^d$. We prove central limit theorems for continuous motion invariant valuations including the Will's functional and the intrinsic volumes of this random polytope. Additionally we derive a central limit theorem for the oracle estimator, that is an unbiased an minimal variance estimator for the volume of a convex set. Finally we obtain a multivariate limit theorem for the intrinsic volumes and the components of the $\mathbf{f}$-vector of the random polytope.

Equations

$$K_t:=\mathrm{conv}\left\{x:x\in {\eta }_t\cap K\right\}.$$

$$\tilde{\phi }(K_t):=\frac{\phi (K_t)-\text{mathds}E\left[\phi (K_t)\right]}{\sqrt{\text{mathds}V\left[\phi (K_t)\right]}}.$$

$$d_W(\tilde{\phi }(K_t),Z)=\mathcal{O}\left(t^{-\frac{1}{2}+\frac{1}{d+1}}\log (t{)}^{3+\frac{2}{d+1}}\right),$$

$$d_W(\tilde{W}(K_t),Z)=\mathcal{O}\left(t^{-\frac{1}{2}+\frac{1}{d+1}}\log (t{)}^{3+\frac{2}{d+1}}\right),$$

$${\hat{\vartheta }}_{oracle}(K_t)=V_d(K_t)+\frac{f_0(K_t)}{t},$$

$$\text{mathds}E\left[{\hat{\vartheta }}_{oracle}(K_t)\right]=V_d(K),$$

$$\text{mathds}V\left[{\hat{\vartheta }}_{oracle}(K_t)\right]=\frac{1}{t}\text{mathds}E\left[V_d(K\setminus K_t)\right]={\gamma }_d\Omega (K)(1+o(1))t^{-1-\frac{2}{d+1}},$$

$$d_W({\tilde{\vartheta }}_{oracle}(K_t),Z)=\mathcal{O}\left(t^{-\frac{1}{2}+\frac{1}{d+1}}\log (t{)}^{3d+\frac{2}{d+1}}\right),$$

$$F(K_t):=\left({\tilde{V}}_1(K_t),\dots ,{\tilde{V}}_d(K_t),{\tilde{f}}_0(K_t),\dots ,{\tilde{f}}_{d-1}(K_t)\right)\in \text{mathds}{R}^{2d}$$

$${\tilde{V}}_j(K_t):=\frac{V_j(K_t)-\text{mathds}E\left[V_j(K_t)\right]}{\sqrt{\text{mathds}V\left[V_j(K_t)\right]}}\text{and}{\tilde{f}}_k(K_t):=\frac{f_k(K_t)-\text{mathds}E\left[f_k(K_t)\right]}{\sqrt{\text{mathds}V\left[f_k(K_t)\right]}}.$$

$$d_3(F(K_t),N_{\Sigma (t)})=\mathcal{O}\left(t^{-\frac{1}{2}+\frac{1}{d+1}}\log (t{)}^{3d+\frac{2}{d+1}}\right)$$

$$x\in \eta \iff x\in \left\{y\in \text{mathds}{R}^d:\eta (\{y\})>0\right\}.$$

$$D_{x}F:=f(\eta +{\delta }_x)-f(\eta ),x\in \text{mathds}{R}^d,$$

$$\underset{\text{mathds}{R}^d}{\int }\text{mathds}E\left[(D_{x}F{)}^2\right]\mu (\mathrm{d}x)<\infty .$$

$$D_{x_1,x_2}^{2}F$$

$$=f(\eta +{\delta }_{x_1}+{\delta }_{x_2})-f(\eta +{\delta }_{x_1})-f(\eta +{\delta }_{x_2})+f(\eta ),x_1,x_2\in \text{mathds}{R}^d.$$

$$d_W(X,Y):=\underset{h\in \mathrm{Lip}(1)}{\sup }\left|\text{mathds}E\left[h(X)\right]-\text{mathds}E\left[h(Y)\right]\right|.$$

$${\tau }_1$$

$${\tau }_2$$

$${\tau }_3$$

$$d_W(F,Z)\le 2\sqrt{{\tau }_1}+\sqrt{{\tau }_2}+{\tau }_3.$$

$$\underset{1\le i_1\le i_2\le m}{\max }\underset{x\in \text{mathds}{R}^d}{\sup }\left|\frac{{\partial }^2}{\partial x_{i_1}\partial x_{i_2}}h(x)\right|\le 1\text{and}\underset{1\le i_1\le i_2\le i_3\le m}{\max }\underset{x\in \text{mathds}{R}^d}{\sup }\left|\frac{{\partial }^3}{\partial x_{i_1}\partial x_{i_2}\partial x_{i_3}}h(x)\right|\le 1.$$

$$d_3(X,Y):=\underset{h\in {\mathcal{H}}_m}{\sup }\left|\text{mathds}E\left[h(X)\right]-\text{mathds}E\left[h(Y)\right]\right|.$$

$${\gamma }_1$$

$${\gamma }_2$$

$${\gamma }_3$$

$$d_3(F,N_{\Sigma })\le \frac{m}{2}\sum _{i,j=1}^m\left|{\sigma }_{ij}-\mathrm{Cov}\left[F_i,F_j\right]\right|+m\sqrt{{\gamma }_1}+\frac{m}{2}\sqrt{{\gamma }_2}+\frac{m^2}{4}{\gamma }_3.$$

$$V_j(K)=\left(\genfrac{}{}{0ex}{}{d}{j}\right)\frac{{\kappa }_d}{{\kappa }_j{\kappa }_{d-j}}\underset{G(d,j)}{\int }{\Lambda }_j(K|\mathbb{L}){\nu }_j(\mathrm{d}\mathbb{L})$$

$$\phi (K)+\phi (L)=\phi (K\cup L)+\phi (K\cap L)$$

$$\phi (L)=\sum _{i=0}^d{c}_i{V}_i(L),$$

$${\mathcal{F}}_i^{\mathrm{Vis}}(P,x):=\left\{\mathfrak{F}\in {\mathcal{F}}_i(P):\forall z\in \mathfrak{F}:[x,z]\cap P=\left\{z\right\}\right\}.$$

$$\mathrm{link}(P,v):=\left\{\mathfrak{F}\in \mathcal{F}(P):v\notin \mathfrak{F}\text{ and }\exists \mathfrak{G}\in \mathcal{F}(P):\mathfrak{F}\subset \mathfrak{G}\ni v\right\},$$

$$f_i(P)=|{\mathcal{F}}_i(P)|.$$

Full text

Keywords. central limit theorem, multivariate limit theorem, intrinsic volumes, f-vector, random polytope, random convex hull, stochastic geometry, Poisson point process, oracle estimator, volume estimation

MSC (2010). 52A22, 60D05, 60F05

1 Introduction

We denote by $\mathcal{K}^{d}$, $d\geq 2$, the collection of all compact convex sets in the Euclidean space $\mathds{R}^{d}$ and by $\mathcal{K}^{d}_{sm}$ the set of all smooth convex bodies in $\mathcal{K}^{d}$, which are all $K\in\mathcal{K}^{d}$ that have nonempty interior, boundary of differentiability class $\mathcal{C}^{2}$, and positive Gaussian curvature everywhere. Let $\eta_{t}$ be a Poisson point process on $\mathds{R}^{d}$ with intensity measure $\mu=t\Lambda_{d}$, where $t>0$ and $\Lambda_{d}$ denotes the $d$-dimensional Lebesgue measure on $\mathds{R}^{d}$, see e.g. [LP18] for more details on Poisson point processes. We fix $K\in\mathcal{K}^{d}_{sm}$ and investigate the random polytope $K_{t}\subset K$ defined as the convex hull of all points in $\eta_{t}\cap K$,

[TABLE]

The investigation of random convex hulls is one of the classical problems in stochastic geometry, see for instance the survey article [Hug13] and the introduction to stochastic geometry [HR16]. Functionals like the intrinsic volumes $V_{j}(K_{t})$ and the components $f_{j}(K_{t})$ of the $\mathbf{f}$-vector, see Section 3, of the random polytope $K_{t}$ have been studied prominently, see \citesCalkaYukich2014,Reitzner2010[Section 1]CalkaSchreiberYukich2013 and the references therein as well as the remarks and references on [LSY17, Theorem 5.5] for more details.

Central limit theorems for $V_{j}(K_{t})$ were proven in the special case that $K$ is the $d$-dimensional Euclidean unit ball, see [CSY13] and [Sch02]. Short proofs for the binomial case $K_{n}$, where $n$ i.i.d. uniformly distributed points in $K$ are considered instead of a Poisson point process, were derived in [TTW18]. Recently [LSY17] embedded the problem for both cases, the binomial and the Poisson case, in the theory of stabilizing functionals deriving central limit theorems for smooth convex bodies removing the logarithmic factors in the error of approximation improving the rate of convergence.

We extend the results of [TTW18] on the intrinsic volumes to the more general case of continuous and motion invariant valuations. This includes the total intrinsic volume functional (Wills functional) $W(K_{t})$ and a central limit theorem for the intrinsic volumes in the Poisson case similar to [TTW18, Theorem 1.1].

We also obtain a univariate central limit theorem for the oracle estimator $\hat{\vartheta}_{oracle}$, that was derived in the remarkable work of [BR16] on the estimation of the volume of a convex body.

Finally we investigate the components of the $\mathbf{f}$-vector $f_{j}(K_{t})$, $j\in\left\{{0,\ldots,d-1}\right\}$ defined as the number of $j$-dimensional faces of $K_{t}$ and derive a multivariate limit theorem on the random vector composed of the intrinsic volumes and the $\mathbf{f}$-vector.

For a continuous and motion invariant valuation $\varphi:\mathcal{K}^{d}\rightarrow\mathds{R}$ we define the random variable $\varphi(K_{t})$ and the corresponding standardization

[TABLE]

We prove a central limit theorem for $t\rightarrow\infty$ under some additional assumptions on the coefficients $c_{i}$ in the linear decomposition into the intrinsic volumes given by the remarkable theorem of Hadwiger stated here as Theorem 3.2, see [Che04, Kla95, Had57] for more details and different proofs. Given this restriction of our model to valuation functionals $\varphi(K_{t})$, that can be safely assumed to not loose variance compared to the intrinsic volume functionals $V_{j}(K_{t})$, we show the first main result of this paper:

We denote by $\operatorname{d}_{W}$ the Wasserstein distance, see (9) in Section 2.1 for the precise Definition.

Theorem 1.1** (Univariate Limit Theorem).**

Assume that $\varphi$ is a continuous and motion invariant valuation, with linear decomposition given by Hadwiger (13), such that $c_{i}c_{j}\geq 0$ for all $i,j\in\left\{{0,\ldots,d}\right\}$ and suppose that at least one index $k\in\left\{{1,\ldots,d}\right\}$ exists, such that $c_{k}\neq 0$. Then

[TABLE]

where $Z\overset{\mathclap{d}}{\sim}\mathcal{N}(0,1)$.

For $j\in\left\{{1,\ldots,d}\right\}$ we can directly obtain the central limit theorem for the standardized $j$-th intrinsic volume in the Poisson model by setting $\varphi(K_{t})=V_{j}(K_{t})$.

Further we can directly obtain a central limit theorem with rate of convergence for the total intrinsic volume, also known as the Wills functional, see [Wil73, Had75, McM75], setting the coefficients $c_{j}=1$, $j\in\left\{{0,\ldots,d}\right\}$ in Theorem 1.1:

Corollary 1.2** (Wills functional).**

Let $W(K_{t})$ denote the Wills functional, defined by $W(K_{t})=\sum_{j=0}^{d}V_{j}(K_{t})$ and denote by $\tilde{W}(K_{t})$ the corresponding standardization. Then

[TABLE]

where $Z\overset{\mathclap{d}}{\sim}\mathcal{N}(0,1)$.

In the remarkable work [BR16] on the estimation of the volume of a convex body $V_{d}(K)$ given the intensity $t>0$ is known, the oracle estimator

[TABLE]

is derived. This estimator is unbiased,

[TABLE]

and of minimal variance among all unbiased estimators (UMVU), see [BR16, eq. (3.2), Theorem 3.2]. Additionally the variance can be obtained by combining [BR16, Theorem 3.2] with [Rei05, Lemma 1] yielding

[TABLE]

for $t\rightarrow\infty$, where the constant $\gamma_{d}$ only depends on the dimension and is known explicitly and $\Omega(K)$ denotes the affine surface area of $K$. This enables us to prove the following univariate limit theorem for the oracle estimator:

Theorem 1.3** (Oracle estimator).**

Let $\hat{\vartheta}_{oracle}$ be the oracle estimator for a smooth convex body $K\in\mathcal{K}^{d}_{sm}$ and denote by $\tilde{\vartheta}_{oracle}$ the corresponding standardization. Then

[TABLE]

where $Z\overset{\mathclap{d}}{\sim}\mathcal{N}(0,1)$.

Finally we provide a multivariate limit theorem for the intrinsic volumes and the components of the $\mathbf{f}$-vector of the random polytope $K_{t}$, that were considered in [LSY17, Theorem 5.5] in the univariate case.

We denote by $\operatorname{d}_{3}$ the distance given by (10) in Section 2.2.

Theorem 1.4** (Multivariate Limit Theorem).**

Let

[TABLE]

be the vector of the standardized intrinsic volumes $\tilde{V}_{j}$ and standardized number of $k$-dimensional faces $\tilde{f}_{k}(K_{t})$, i.e.

[TABLE]

We denote by $F_{i}:=F_{i}(K_{t})$, $i=1,\ldots,2d$ the $i$-th component of the multivariate functional $F(K_{t})$. Define $\Sigma(t)=(\sigma_{ij}(t))_{i,j\in\left\{{1,\ldots,2d}\right\}}$ as the covariance matrix of $F(K_{t})$, i.e. $\sigma_{ij}(t):={\mathrm{Cov}}\left[F_{i},F_{j}\right]$ and $\sigma_{ii}(t)=1$ for all $i\in\left\{{1,\ldots,2d}\right\}$ and all $t>0$. Then

[TABLE]

where $N_{\Sigma(t)}\overset{\mathclap{d}}{\sim}\mathcal{N}(0,\Sigma(t))$.

Note that $N_{\Sigma(t)}$ still depends on the intensity $t$. This gives rise to the following questions:

Open Problems:

The limit of the variances and co-variances are still unknown, thus we can not set $\sigma_{ij}$ to be the limit of the correlations (rescaled co-variances) $\sigma_{ij}(t)$. These limits would give rise to a limit theorem providing a fixed multivariate Gaussian distribution $\mathcal{N}(0,\Sigma)$ with fixed co-variance matrix $\Sigma$. In this case, the rate of the limit theorem will also contain the rate of the correlations on the right hand side of the bound, thus it would be beneficial to obtain the limit $\sigma_{ij}$ of $\sigma_{ij}(t)$ including an upper bound for $\lvert\sigma_{ij}(t)-\sigma_{ij}\rvert$, since the convergence of the correlations could be slower than the rate given in Theorem 1.4, slowing down the overall convergence. We should mention, that Calka, Schreiber and Yukich, were able to derive limits for the variance in the case that $K$ is the euclidean unit ball using the paraboloid growth process, see [CSY13], but up to our knowledge there are no results on the limit of the variance in a general (smooth) convex body and also no results on the co-variances of $F(K_{t})$.

Note that the Euler-Poincaré equations and more general the Dehn-Sommerville equations, see [Zie95, Chapter 8], imply linear dependencies on the components of the $\mathbf{f}$-vector. Thus the covariance matrix $\Sigma(t)$ is singular and therefore ${\mathrm{rank}}\left(\Sigma(t)\right)<2d$, which gives rise to the question what ${\mathrm{rank}}\left(\Sigma(t)\right)$ actually is and if this also applies to the limiting co-variance matrix $\Sigma$?

Remark 1.5**.**

Note that the univariate results can be derived for the $d_{3}$-distance using the multivariate result. Since the additional work that is needed to prove the univariate results alongside the multivariate limit theorem is small, we decided to directly proof the univariate results in the Wasserstein distance.

The paper is organized as follows. For the convenience of the reader, we repeat the relevant material on the Malliavin-Stein-Method for normal approximation of Poisson functionals in Section 2. In Section 3 we introduce some background material on convex geometry and corresponding estimations using floating bodies without proofs, to keep our presentation reasonably self-contained. The proofs of our main results are presented in Section 4, starting with the central limit theorem for valuations, Theorem 1.1, in Subsection 4.1. handling the intrinsic volumes. In Subsection 4.2 we prove the multivariate limit theorem, Theorem 1.4 by extending our proof to the components of the $\mathbf{f}$-vector. Finally we can combine the results derived in the proofs before to obtain the central limit theorem for the oracle estimator, Theorem 1.3.

2 Stochastic preliminaries

Let $\eta$ be a Poisson point process over the Euclidean space $(\mathds{R}^{d},\mathscr{B}^{d})$ with intensity measure $\mu$. One can think of $\eta$ as a random element in the space $\mathrm{N}_{\sigma}$ of all $\sigma$-finite counting measures $\chi$ on $\mathds{R}^{d}$, i.e. $\chi(B)\in\mathds{N}\cup\{\infty\}$ for all $B\in\mathscr{X}$, where the space $\mathrm{N}_{\sigma}$ is equipped with the $\sigma$-field $\mathscr{N}_{\sigma}$ generated by the mappings $\chi\rightarrow\chi(B)$, $B\in\mathscr{B}^{d}$. To simplify our notation we will often handle $\eta$ as a random set of points given by

[TABLE]

We call a random variable $F$ a Poisson functional if there exists a measurable map $f:\mathrm{N}_{\sigma}\rightarrow\mathds{R}$ such that $F=f(\eta)$ almost surely. The map $f$ is called the representative of $F$. We define the difference operator or so called “add-one-cost operator”:

Definition 2.1**.**

Let $F$ be a Poisson functional and $f$ its corresponding representative, then the first order difference operator is given by

[TABLE]

where $\delta_{x}$ denotes the Dirac measure with mass concentrated in $x$. We say that $F$ belongs to the domain of the difference operator, i.e. $F\in\operatorname{dom}(D)$, if ${\mathds{E}}\left[F^{2}\right]<\infty$ and

[TABLE]

The second order difference operator is obtained by iterating the definition:

[TABLE]

For a depper discussion of the underlying theory of Poisson point processes, Malliavin-Calculus, the Wiener-Itô-Chaos Expansion and the Malliavin-Stein Method see [PR16] and [LP18].

2.1 Malliavin-Stein-Method for the univariate case

We denote by $\operatorname{Lip}(1)$ the class of Lipschitz functions $h:\mathds{R}\rightarrow\mathds{R}$ with Lipschitz constant less or equal to one. Given two $\mathds{R}$-valued random variables $X,Y$, with ${\mathds{E}}\left\lvert X\right\rvert<\infty$ and ${\mathds{E}}\left\lvert Y\right\rvert<\infty$ the Wasserstein distance between the laws of $X$ and $Y$, written $\operatorname{d}_{W}(X,Y)$ is defined as

[TABLE]

Note that if a sequence $(X_{n})$ of random variables satisfies $\lim_{n\rightarrow\infty}\operatorname{d}_{W}(X_{n},Y)=0$, then it holds that $X_{n}$ converges to $Y$ in distribution, see [LP18, p. 219 and Proposition B.9]. Especially if $Y$ has standard Gaussian distribution, we obtain a central limit theorem by showing $\operatorname{d}_{W}(X_{n},Y)\rightarrow 0$, which we will achieve by rephrasing the bound given by [LPS16, Theorem 1.1] which is an extension of [Pec+10, Theorem 3.1], see also [LP18, Chapter 21.1, 21.2] for a slightly different form and proofs.

Theorem 2.2**.**

Let $F\in\operatorname{dom}(D)$ be a Poisson functional such that ${\mathds{E}}\left[F\right]=0$ and ${\mathds{V}}\left[F\right]=1$. Define

[TABLE]

and let $Z$ be a standard Gaussian random variable, then

[TABLE]

2.2 Malliavin-Stein-Method for the multivariate case

We denote by $\mathcal{H}_{m}$ the class of all $\mathcal{C}^{3}$-functions $h:\mathds{R}^{m}\rightarrow\mathds{R}$ such that all absolute values of the second and third partial derivatives are bounded by one, i.e.

[TABLE]

Given two $\mathds{R}^{m}$-valued random vectors $X,Y$ with ${\mathds{E}}\left\lVert X\right\rVert^{2}<\infty$ and ${\mathds{E}}\lVert Y\rVert^{2}<\infty$ the distance $\operatorname{d}_{3}$ between the laws of $X$ and $Y$, written $\operatorname{d}_{3}(X,Y)$ is defined as

[TABLE]

Note that if a sequence $(X_{n})$ of random vectors satisfies $\lim_{n\rightarrow\infty}\operatorname{d}_{3}(X_{n},Y)=0$, then it holds that $X_{n}$ converges to $Y$ in distribution, see [PZ10, Remark 2.16]. Especially if $Y$ is a $m$-dimensional centered Gaussian vector with covariance matrix $\Sigma\in\mathds{R}^{m\times m}$, we obtain a multivariate limit theorem by showing $\operatorname{d}_{3}(X_{n},Y)\rightarrow 0$. We will achieve this, similar to the univariate central limit theorem, by rephrasing the bound given by [SY18, Theorem 1.1] which extends [PZ10], to provide the multivariate analogon to the univariate result derived in [LPS16], that was stated here as Theorem 2.2 above.

Theorem 2.3**.**

Let $F=(F_{1},\ldots,F_{m})$, $m\geq 2$, be a vector of Poisson functionals $F_{1},\ldots,F_{m}\in\operatorname{dom}(D)$ with ${\mathds{E}}\left[F_{i}\right]=0$, $i\in\left\{{1,\ldots,m}\right\}$. Define

[TABLE]

and let $\Sigma=(\sigma_{ij})_{i,j\in\left\{{1,\ldots,m}\right\}}\in\mathds{R}^{m\times m}$ be positive semi-definite, then

[TABLE]

3 Geometric preliminaries

Fix $j\in\left\{{1,\ldots,d}\right\}$ and let $G(d,j)$ denote the Grassmannian of all $j$-dimensional linear subspaces of $\mathds{R}^{d}$ equipped with the uniquely determined Haar probability measure $\nu_{j}$, see [Sch14, Section 4.4]. For $k\in\mathds{N}$ the $k$-dimensional volume of the $k$-dimensional unit ball $\mathds{B}^{k}$ is denoted by $\kappa_{k}:=\pi^{\frac{k}{2}}\Gamma\left({1+\frac{k}{2}}\right)^{-1}$.

For a convex body $K\in\mathcal{K}^{d}$ the $j$-dimensional Lebesgue measure of the orthogonal projection of $K$ onto the linear subspace $\mathbb{L}\in G(d,j)$ is denoted by $\Lambda_{j}(K|\mathbb{L})$.

For $j\in\left\{{1,\ldots,d}\right\}$, the $j$-th intrinsic volume of $K$ is given by Kubota’s formula, see [SW08, p. 222]:

[TABLE]

and for $j=0$ the [math]-th intrinsic volume of $K$, $V_{0}(K)$ is the Euler characteristics of $K$, therefore we have $V_{0}(K)={\mathds{1}}\left\{{K\neq\emptyset}\right\}$. It is worth mentioning, that $V_{1}(K)$ is the mean with, $V_{d-1}(K)$ is the surface area up to multiplicative constants not depending on $K$ and $V_{d}(K)$ equals the $d$-dimensional Lebesgue-volume of $K$. The intrinsic volumes are crucial examples of continuous, motion invariant valuations:

Definition 3.1**.**

A real function on the space of convex bodies, $\varphi:\mathcal{K}^{d}\rightarrow\mathds{R}$ , is called a valuation, if and only if

[TABLE]

holds, whenever $K,L,K\cup L\in\mathcal{K}^{d}$. It is called continuous, if it is continuous according to the Hausdorff metric on $\mathcal{K}^{d}$, and it is called invariant under rigid motions if it is invariant under translations and rotations on $\mathds{R}^{d}$.

The theorem of Hadwiger [Had57, Kla95, Che04] states, that every continuous and motion invariant valuation $\varphi:\mathcal{K}^{d}\rightarrow\mathds{R}$ can be decomposed into a linear combination of intrinsic volumes:

Theorem 3.2** (Hadwiger).**

Let $\varphi$ be a continuous and motion invariant valuation. Then there existing coefficients $c_{i}\in\mathds{R}$, $i\in\left\{{0,\ldots,d}\right\}$, such that for all convex sets $L\in\mathcal{K}^{d}$ it holds, that

[TABLE]

where $V_{i}$ denotes the $i$-th dimensional intrinsic volume.

For further information on Hadwiger’s theorem, convex geometry and integral geometry we refer the reader to [Gru07, SW08, Sch14].

Let $P\in\mathcal{K}^{d}$ be a polytope and $i\in\left\{{0,\ldots,d}\right\}$. We denote by $\mathcal{F}_{i}(P)$ the set of all $i$-dimensional faces, $i$-faces for short, of $P$ and by $\mathcal{F}_{i}^{\operatorname{Vis}}(P,x)$ the restriction to those $i$-faces that can be seen from $x$, where we consider a face $\mathfrak{F}$ of $P$ to be seen by $x$ if all points $z\in\mathfrak{F}$ can be connected by a straight line $[z,x]$ to $x$ such that the intersection of this line with $P$ only contains the starting point $z$, i.e.

[TABLE]

The sets of all faces resp. all visible faces are denoted by $\mathcal{F}(P):=\bigcup_{i=0}^{d}\mathcal{F}_{i}(P)$ resp. $\mathcal{F}^{\operatorname{Vis}}(P,x):=\bigcup_{i=0}^{d}\mathcal{F}_{i}^{\operatorname{Vis}}(P,x)$.

For a vertex $v\in\mathcal{F}_{0}(P)$ the link of $v$ in $P$ is the set of all faces of $P$ that do not contain $v$ but are contained in a (higher dimensional) face that contains $v$, i.e.

[TABLE]

see [Zie95, Chapter 8.1] for a recent account of the theory.

The number of $i$-dimensional faces of $P$ will be denoted by $f_{i}(P)$, i.e.

[TABLE]

Note that the vector $(f_{-1}(P),f_{0}(P),\ldots,f_{d}(P))$ with $f_{-1}(P):=V_{0}(P)$ is the $\mathbf{f}$-vector of $P$, see [Zie95, Definition 8.16, p. 245] for more details.

3.1 Geometric estimations

We introduce the notion of the $\varepsilon$-floating body, following [Rei10, Section 2.2.3]. For a fixed $K\in\mathcal{K}^{d}$ and a closed halfspace $H$ we call the intersection $C=H\cap K$ a cap of $K$. If $C$ has volume $V_{d}(C)=\varepsilon$, we call $C$ an $\varepsilon$-cap of $K$. We define the function $v:K\rightarrow\mathds{R}$ by

[TABLE]

and the floating body with parameter $\varepsilon$, $\varepsilon$-floating body for short, as the level set

[TABLE]

which is convex, since it is the intersection of halfspaces. The wet part of $K$ is defined as $K(v\leq\varepsilon)$, where the name comes from the $3$-dimensional picture when $K$ is a box containing $\varepsilon$ units of water. Note that the $\varepsilon$-floating body is (up to its boundary) the remaining set of $K$, if all $\varepsilon$-caps are removed from $K$ and the wet part is the union of these caps. For the convenience of the reader, we will only use the notation for the floating body $K(v\geq\varepsilon)$ to prevent confusion with the wet part denoted by $K(v\leq\varepsilon)$. From now on, we will assume that the parameter $\varepsilon>0$ is sufficiently small. Thus we can use the following lemmas from [Vu05, Lemma 6.1-6.3]:

Lemma 3.3**.**

Let $C$ be an $\varepsilon$-cap of $K$, then there are two constants $c_{1},c_{2}\in(0,\infty)$ such that the diameter of $C$, ${\operatorname{diam}}(C)=\sup_{x,y\in C}\lVert x-y\rVert$, is bounded by

[TABLE]

Lemma 3.4**.**

Let $x$ be a point on the boundary $\partial K$ of $K$ and $D(x,\varepsilon)$ the set of all points on the boundary which are of distance at most $\varepsilon$ to $x$. Then the convex hull of $D(x,\varepsilon)$ has volume at most $c_{3}\varepsilon^{d+1}$, where $c_{3}\in(0,\infty)$ is some constant not depending on $\varepsilon$.

Lemma 3.5**.**

Let $C$ be an $\varepsilon$-cap of $K$. The union of all $\varepsilon$-caps intersecting $C$ has volume at most $c_{4}\varepsilon$, where $c_{4}\in(0,\infty)$ is some constant not depending on $\varepsilon$.

4 Proofs of the main results

To shorten our notation we write $K^{x}_{t}$ resp. $K^{y}_{t}$ for the convex hull of $(\eta_{t}+\delta_{x})\cap K$ resp. $(\eta_{t}+\delta_{y})\cap K$ and $K^{xy}_{t}$ for the convex hull of $(\eta_{t}+\delta_{x}+\delta_{y})\cap K$. Further we will use $\mathbf{C}\in(0,\infty)$ to denote a constant, that can depend on the dimension and the convex set $K$ but is independent of the intensity of our Poisson point process $t$. For sake of brevity we will not mention this properties of $\mathbf{C}$ in the following, additionally the value of $\mathbf{C}$ may also differ from line to line. We will use $g(t)\ll f(t)$ to indicate that $g(t)$ is of order at most $f(t)$, i.e.

[TABLE]

where $c$ and $t_{0}$ are constants not depending on $t$. We will use $g(t)=\Theta(f(t))$ to indicate that $g(t)$ is of the same order of $f(t)$, i.e.

[TABLE]

For sufficiently large $t>0$, we define $\varepsilon_{t}:=c\tfrac{\log(t)}{t}$ with $c>0$ and denote by $K(v\geq\varepsilon_{t})$ the $\varepsilon_{t}$-floating body of $K$. Let $A(\varepsilon_{t},t):=\left\{{K(v\geq\varepsilon_{t})\subseteq K_{t}}\right\}$ be the event, that the $\varepsilon_{t}$-floating body is contained in the random polytope $K_{t}$. Recall the well known Lemma from \citesBaranyDalla1997Vu2005[Lemma 2.2]Reitzner2010 in a slightly modified version for the Poisson case:

Lemma 4.1**.**

For any $\beta\in(0,\infty)$ there exists a constant $c(\beta,d)\in(0,\infty)$ only depending on $\beta$ and the space dimension $d$, such that the event $A(\varepsilon_{t},t)$, that the $\varepsilon_{t}$-floating body is contained in the random polytope $K_{t}$ occurs with high probability. More precisely, the probability of the complementary event $A^{c}(\varepsilon_{t},t)$ has polynomial decay with exponent $-\beta$ for $t\rightarrow\infty$, i.e.

[TABLE]

whenever $t$ is sufficiently large.

Note that the parameter $\beta$ can be freely chosen in $(0,\infty)$, thus for $\beta=16d+1$, which is sufficiently large for all our purposes, we get $c(\beta,d)$ and therefore we can define $\varepsilon_{t}$ such that $K(v\geq\varepsilon_{t})\subseteq K_{t}$ with high probability according to Lemma 4.1.

We will use the following estimation of subsets of $G(d,j)$ from [BFV10, Lemma 1] to handle the projections arising from Kubota’s formula in our proof of Theorem 1.1:

Lemma 4.2**.**

For $z\in\mathbb{S}^{d-1}$ and $\mathbb{L}\in G(d,j)$ we define the angle $\sphericalangle(z,\mathbb{L})$ as the minimum of all angles $\sphericalangle(z,x)$, $x\in\mathbb{L}$. For sufficiently small $\alpha>0$ one has that

[TABLE]

4.1 Proof of Theorem 1.1: valuation functional

We first recall that the valuation functional $\varphi(K_{t})$ can be decomposed with Hadwiger (13) into the linear combination of intrinsic volumes, thus the variance ${\mathds{V}}\left[\varphi(K_{t})\right]$ can be rewritten into

[TABLE]

For $V_{i}$, $i\in\left\{{1,\ldots,d}\right\}$ we will use the variance bound from [LSY17, eq. 5.20, eq. 5.22, eq. 5.23], see also [CSY13, Corollary 7.1] and [Rei05]:

[TABLE]

Since $V_{0}(K_{t})$ is the Euler characteristics of $K_{t}$ we have $V_{0}(K_{t})={\mathds{1}}\left\{{K_{t}\neq\emptyset}\right\}$ and therefore $V_{0}(K_{t})$ is a Bernoulli distributed random variable with success probability ${\mathds{P}}\left(V_{0}(K_{t})=1\right)=1-e^{-t\Lambda_{d}(K)}$. The expectation is given by ${\mathds{E}}\left[V_{0}(K_{T})\right]=1-e^{-t\Lambda_{d}(K)}$ and the variance by ${\mathds{V}}\left[V_{0}(K_{t})\right]=(1-e^{-t\Lambda_{d}(K)})e^{-t\Lambda_{d}(K)}$, which can be bounded by

[TABLE]

Lemma 4.3**.**

For all $i,j\in\left\{{0,\ldots,d}\right\}$ the intrinsic volumes of $K_{t}$ are non negatively correlated and their covariances are bounded from above by the same order of magnitude as the variance, i.e.

[TABLE]

• Proof:

Since $D_{x}V_{j}(K_{t})\geq 0$ for all $x\in\mathds{R}^{d}$ it follows from the Harris-FKG inequality for Poisson processes, see [Las16, Theorem 11], that ${\mathds{E}}\left[V_{i}(K_{t})V_{j}(K_{t})\right]\geq{\mathds{E}}\left[V_{i}(K_{t})\right]{\mathds{E}}\left[V_{j}(K_{t})\right]$, which directly implies the lower bound on the covariances. The Cauchy-Schwarz inequality implies

[TABLE]

thus, using (14) and (15), the upper bound on the covariances is obtained. ∎

We are now in a position to bound the variance of our valuation functional with the following Lemma:

Lemma 4.4**.**

Under the assumptions of Theorem 1.1 the variance of the valuation functional is bounded by

[TABLE]

• Proof:

We assumed $c_{i}c_{j}\geq 0$ for all $i,j$ and that there exists at least one index $k\in\left\{{1,\ldots,d}\right\}$ such that $c_{k}\neq 0$. Thus Lemma 4.3 implies

[TABLE]

and

[TABLE]

which completes the proof. ∎

The crucial part in the proof of Theorem 1.1 is the application of the general bound given by Theorem 2.2, thus we need to investigate the moments occuring in $\tau_{1}$, $\tau_{2}$ and $\tau_{3}$. In the first step, we adapt and slightly extend the proof from [TTW18] for the binomial case, to work in the Poisson case, yielding upper bounds on the moments of the first and second order difference operators applied to the intrinsic volumes $V_{j}(K_{t})$ which will be used in the second step to derive the bounds for the valuation functional.

First order difference operator:

Fix $x\in K$ and $j\in\{{1,\ldots,d}\}$, then conditioned on the event $A(\varepsilon_{t},t)$, it follows that

[TABLE]

thus we can restrict the following to the case $x\in K\setminus K(v\geq\varepsilon_{t})$.

For $x\in K\setminus K(v\geq\varepsilon_{t})$ we define $z$ to be the closest point to $x$ on the boundary $\partial K$. Since $K$ is smooth $z$ is uniquely determined, if $\varepsilon_{t}$ is sufficiently small.

The visible region of $z$ is defined as the set of all points that can be connected to $z$ by a straight line in $K$ avoiding $K(v\geq\varepsilon_{t})$, i.e.

[TABLE]

Note that, given the sandwiching $K(v\geq\varepsilon_{t})\subset K_{t}\subset K$, a random point $x$ can influence the random polytope only in the visibility region.

We construct a full-dimensional spherical cap $C$ such that $K^{x}_{t}\setminus K_{t}\subseteq C$. The definition of the visible region, that was first used in [Vu05] and [BR10], is crucial in the following steps:

Let $y_{1},y_{2}\in\operatorname{Vis}_{z}(t)$, then there existing two $\varepsilon_{t}$-caps $C_{1}$ and $C_{2}$ such that the straight line $[y_{1},z]$ resp. $[y_{2},z]$ is contained in $C_{1}$ resp. $C_{2}$, thus

[TABLE]

Since the diameter of any $\varepsilon_{t}$-cap $C$ of $K$ can be bounded by $\mathbf{C}\varepsilon_{t}^{\frac{1}{d+1}}$, see Lemma 3.3, it follows directly that the diameter of the visibility region can be bounded by

[TABLE]

Let $D(z,\rho)$ be the set of all points on the boundary $\partial K$ which are of distance at most $\rho$ to $z$, i.e.

[TABLE]

and denote the cap, that is given by the convex hull of $D(z,\rho)$ by $C$, i.e.

[TABLE]

By construction, we have $K^{x}_{t}\setminus K_{t}\subseteq\operatorname{Vis}_{z}(t)\subseteq C$. It follows from Lemma 3.4, that the volume of $C$ is of order at most $\frac{\log(t)}{t}$.

Fix a linear subspace $\mathbb{L}\in G(d,j)$, then one has that the set-difference of the projection of $K^{x}_{t}$ and $K_{t}$ onto the subspace $\mathbb{L}$ is contained in the projection of $C$ onto $\mathbb{L}$,

[TABLE]

The $j$-dimensional volume of the projected cap $C|\mathbb{L}$ can be bounded in its order of magnitude by

[TABLE]

Depending on the angle between $z$ and $\mathbb{L}$, $\sphericalangle(z,\mathbb{L})$, the part $K^{x}_{t}\setminus K_{t}$ is not visible for the orthogonal projection on $\mathbb{L}$ since it is hidden behind $K(v\geq\varepsilon_{t})$, i.e.

[TABLE]

for sufficiently large $t$. To obtain a bound on the maximal angle $\sphericalangle(z,\mathbb{L})$ where the projection does not vanish we approximate $K$ by a ball $\mathds{B}^{d}(z_{c},r)$ with center $z_{c}$ and radius $r$ such that $\mathds{B}^{d}(z_{c},r)\subseteq K$ and $\mathds{B}^{d}(z_{c},r)\cap\partial K=\left\{{z}\right\}$. Indeed we approximate the boundary $\partial K$ of $K$ from the inside of $K$ by a ball, which is possible, since $K$ is sufficiently smooth.

We repeat the construction of the cap $C$ for $\mathds{B}^{d}(z_{c},r)$ with the corresponding $\varepsilon_{t}$-floating body $\mathds{B}^{d}(v\geq\varepsilon_{t})$ of the ball to obtain the cap $C_{\mathds{B}}$ and define $\alpha$ to be the central angle of $C_{\mathds{B}}$ in $\mathds{B}^{d}(z_{c},r)$. It follows from Lemma 3.4, that the volume of $C_{\mathds{B}}$ is of order at most $\frac{\log(t)}{t}$, since $\rho_{\mathds{B}}\ll({\frac{\log(t)}{t}})^{\frac{1}{d+1}}$, which yields

[TABLE]

Thus it follows from $\mathds{B}^{d}(v\geq\varepsilon_{t})\subseteq K(v\geq\varepsilon_{t})\subseteq K_{t}\subseteq K_{t}^{x}$ that $K_{t}^{x}|\mathbb{L}=K_{t}|\mathbb{L}$ if $\sphericalangle(z,\mathbb{L})$ is of larger order than $\alpha$, therefore we have

[TABLE]

Using (22) and (21) it follows that

[TABLE]

Finally we use Kubota’s formula (11) together with (18) and Lemma 4.2 to obtain

[TABLE]

where $c(d,j)=\binom{d}{j}\frac{\kappa_{d}}{\kappa_{j}\kappa_{d-j}}$ can be omitted since we are bounding $D_{x}V_{j}(K_{t})$ in its order of magnitude with respect to $t$.

Second order difference operator:

Fix $x,y\in K$ and $j\in\left\{{1,\ldots,d}\right\}$, similar to (18) and conditioned on the event $A(\varepsilon_{t},t)$, we have that

[TABLE]

To further restrict $x,y\in K\setminus K(v\geq\varepsilon_{t})$ we show the following Lemma:

Lemma 4.5**.**

Fix two convex bodies $P,K\in\mathcal{K}^{d}$ with $P\subset K$ and two points $x,y\in K\setminus P$. Denote by $P^{xy}$, $P^{x}$ and $P^{y}$ the convex hulls of $P\cup\left\{{x,y}\right\}$, $P\cup\left\{{x}\right\}$ resp. $P\cup\left\{{y}\right\}$. We define the visibility region of $x$ with respect to $P$ by

[TABLE]

If $\mathcal{V}_{x}(P)\cap\mathcal{V}_{y}(P)=\emptyset$, then

[TABLE]

and further it follows for all valuations $\psi:\mathcal{K}^{d}\rightarrow\mathds{R}$ that the second order difference operator of $\psi(P)$ vanishes, i.e.

[TABLE]

• Proof:

Using $P^{x}\subseteq\mathcal{V}_{x}(P)\cup P$ and $P^{y}\subseteq\mathcal{V}_{y}(P)\cup P$ it follows directly from $\mathcal{V}_{x}(P)\cap\mathcal{V}_{y}(P)=\emptyset$ that $P^{x}\cap P^{y}\subseteq(\mathcal{V}_{x}(P)\cap\mathcal{V}_{y}(P))\cup P=P$. Additionally the inclusion $P\subseteq P^{x}\cap P^{y}$, that follows directly from the definition of the convex hull, gives (23).

Again, it is immediate that $P^{x}\subseteq P^{xy}$ and $P^{y}\subseteq P^{xy}$, thus it remains to prove that $P^{xy}\subseteq P^{x}\cup P^{y}$. Assume $z\in P^{xy}\setminus(P^{x}\cup P^{y})$, then there existing $\lambda_{1},\lambda_{2}\in[0,1]$ and $u\in P^{x}$, $v\in P^{y}$ such that

[TABLE]

where we can safely assume that $u$ and $v$ are chosen such that $\lambda_{1}$ and $\lambda_{2}$ are maximized. Note that $u\in P^{y}$ implies $z\in P^{y}$ resp. $v\in P^{x}$ implies $z\in P^{x}$, a contradiction, which leads to the remaining case $u\in P^{x}\setminus P^{y}$ and $v\in P^{y}\setminus P^{x}$. By construction it now follows that $[x,z]\cap P=\emptyset$ and $[y,z]\cap P=\emptyset$, which yields $z\in\mathcal{V}_{x}(P)$ and $z\in\mathcal{V}_{y}(P)$ contrary to $\mathcal{V}_{x}(P)\cap\mathcal{V}_{y}(P)=\emptyset$ which gives (24).

The second order difference operator of $\psi(P)$ is given by

[TABLE]

Using (23) and (24) we can rewrite the first term according to the valuation property (12) to

[TABLE]

which gives (25) when substituted in the representation of $D^{2}_{x,y}\psi(P)$. ∎

Since $\operatorname{Vis}_{x}(t)\cap\operatorname{Vis}_{y}(t)=\emptyset$ implies the conditions of Lemma 4.5 for $P=K_{t}\supseteq K(v\geq\varepsilon_{t})$ it follows that

[TABLE]

Taking $D_{x,y}^{2}V_{j}(K_{t})=D_{x}(D_{y}V_{j}(K_{t}))$ we obtain $\left\lvert D_{x,y}^{2}V_{j}(K_{t})\right\rvert\leq D_{x}V_{j}(K^{y}_{t})+D_{x}V_{j}(K_{t})$ where we immediately see, that the second term $D_{x}V_{j}(K_{t})$ is the first order difference operator that we have bounded before. Using $K(v\geq\varepsilon_{t})\subseteq K_{t}\subseteq K_{t}^{y}$ we can substitute $K_{t}$ with $K_{t}^{y}$ in the proof for the first order difference operator to obtain

[TABLE]

The results of the prior discussion can be summarized in the following lemma bounding the order of magnitude of the $p$-th absolute moment of the first and second order difference operator of the intrinsic volumes $V_{j}(K_{t})$.

Lemma 4.6**.**

Let $p\in\left\{{1,\ldots,8}\right\}$, $j\in\left\{{0,\ldots,d}\right\}$ and $x,y\in K$, then

[TABLE]

• Proof:

Let $j\neq 0$. On the event $A(\varepsilon_{t},t)$ we use the bounds derived before and on the complementary event $A^{c}(\varepsilon_{t},t)$ it is sufficient to use the estimations $D_{x}V_{j}(K_{t})\leq V_{j}(K)$ resp. $\left\lvert D_{x,y}^{2}V_{j}(K_{t})\right\rvert\leq 2V_{j}(K)$, thus

[TABLE]

and

[TABLE]

Since ${\mathds{P}}\left(A(\varepsilon_{t},t)\right)\leq 1$ and ${\mathds{P}}\left(A^{c}(\varepsilon_{t},t)\right)\leq t^{-\beta}$, with $\beta=16d+1>p$, see Lemma 4.1, our claim follows for $j\in\left\{{1,\ldots,d}\right\}$.

Let $j=0$. We use $V_{0}(K_{t})={\mathds{1}}\left\{{K_{t}\neq\emptyset}\right\}$ to derive

[TABLE]

and

[TABLE]

thus using ${\mathds{P}}\left(K_{t}=\emptyset\right)=e^{-t\Lambda_{d}(K)}$ the claim follows by bounding the exponential decay with $(\frac{\log(t)}{t})^{p}$. ∎

Our next objective is to prove corresponding bounds on the moments of the first and second order difference operator of the valuation functional we wish to study.

Lemma 4.7**.**

Let $p\in\left\{{1,\ldots,8}\right\}$ and $x,y\in K$, then

[TABLE]

• Proof:

Since $D_{x}$ is linear we obtain

[TABLE]

and the generalized Hölder inequality, yields

[TABLE]

Therefore we can use (26) to bound all moments on the right hand side yielding (28). The proof of (29) is similar using (27) instead of (26). ∎

Before we apply the previous lemma to the error bounds $\tau_{1}$, $\tau_{2}$ and $\tau_{3}$ we introduce two estimations for the domain of integration. From [Bár08, Theorem 6.3] it follows directly that

[TABLE]

Denote by $C(x)$ resp. $C(z)$ the caps constructed according to (20) for points $x,z\in K\setminus K(v\geq\varepsilon_{t})$, then for every fixed $x\in K\setminus K(v\geq\varepsilon_{t})$ we have

[TABLE]

Recall that the volumes of $C(x)$ and $C(z)$ are of order at most $\frac{\log(t)}{t}$, thus Lemma 3.5 yields

[TABLE]

for all $x\in K\setminus K(v\geq\varepsilon_{t})$.

Applying the previous results to the error bound $\tau_{1}$ yields

[TABLE]

In the same manner we can see that

[TABLE]

Combining these three bounds with Theorem 2.2 leads to

[TABLE]

for $Z\overset{\mathclap{d}}{\sim}\mathcal{N}(0,1)$, completing the proof of Theorem 1.1.

4.2 Proof of Theorem 1.4: multivariate functional

We start to investigate the moments of the first and second order difference operators applied to the components of the $\mathbf{f}$-vector by combining combinatorial results from [Rei05a] with the floating body and economic cap covering approach from [Rei10] and [TTW18].

First order difference operator:

Fix $x\in K$ and $j\in\left\{{0,\ldots,d-1}\right\}$, then conditioned on the event $A(\varepsilon_{t},t)$, it follows similar to (18) that

[TABLE]

thus we can restrict the following to the case $x\in K\setminus K(v\geq\varepsilon_{t})$.

Let $K_{t}$ be fixed and assume $x\not\in K_{t}$. Since the polytope $K_{t}$ is simplicial and all verticies are in general position almost surely, analysis similar to that in [Rei05a, Section 4] allows us to decompose $D_{x}f_{j}(K_{t})$ into the number of $j$-faces gained, denoted by $f_{j}^{+}$, and the number of $j$-faces lost, denoted by $f_{j}^{-}$:

[TABLE]

Every $j$-face gained in $K_{t}^{x}$ is the convex hull of $x$ and a $(j-1)$-face in $\operatorname{link}(K_{t}^{x},x)$. Additionally it can easily be seen that every $(j-1)$-face in $\operatorname{link}(K_{t}^{x},x)$ is also contained in $\mathcal{F}_{j-1}^{\operatorname{Vis}}(K_{t},x)$, thus

[TABLE]

On the other hand, the $j$-faces in $\mathcal{F}_{j}(K_{t})$ that are lost have to be visible from $x$, thus

[TABLE]

Note that not all visible $j$-faces are removed, to gain the exact amount of lost $j$-face, one has to calculate the number of $j$-faces that are visible and not contained in the link of $x$ in the new polytope $K_{t}^{x}$, i.e. $\left\lvert\mathcal{F}_{j}^{\operatorname{Vis}}(K_{t},x)\setminus\operatorname{link}(K_{t}^{x},x)\right\rvert$, as we will see later for the second order difference operator.

Let $z$ be the closest point to $x$ on the boundary of $\partial K$, then it follows immediately that every visible $i$-face $\mathfrak{F}\in\mathcal{F}_{i}^{\operatorname{Vis}}(K_{t},x)$ has to be as subset of $\operatorname{Vis}_{z}(t)$. Since $(i+1)$ pairwise distinct points are needed to form an $i$-face, the number of visible $i$-faces can be bound by the number of points in $\operatorname{Vis}_{z}(t)$, i.e.

[TABLE]

for all $i\in\left\{{0,\ldots,d-1}\right\}$. Combining these steps we obtain

[TABLE]

and further for $p\in\left\{{1,\ldots,8}\right\}$,

[TABLE]

The binomial theorem and the fact that $\eta\left({\operatorname{Vis}_{z}(t)}\right)$ is Poisson distributed with parameter $\mu(\operatorname{Vis}_{z}(t))$ yields

[TABLE]

where $\genfrac{\{}{\}}{0.0pt}{}{i}{k}$ denotes the Stirling numbers of the second kind. Recall that $\mu(\operatorname{Vis}_{z}(t))=t\Lambda_{d}(\operatorname{Vis}_{z}(t))\ll t\frac{\log(t)}{t}=\log(t)$ and $j\in\left\{{0,\ldots,d-1}\right\}$, thus

[TABLE]

Conditioned on the complementary event $A^{c}(\varepsilon_{t},t)$ we need to slightly modify the proof, replacing $\eta(\operatorname{Vis}_{z}(t))$ by $\eta(K)$, the number of all points in $K$, and using the Cauchy-Schwarz inequality to separate the expectation of the indicator, from the moments of $\eta(K)$:

[TABLE]

since $\beta=16d+1\geq 2pd$.

Second order difference operator:

Fix $x,y\in K$ and $j\in\left\{{0,\ldots,d-1}\right\}$, similar to the intrinsic volumes handled before we have

[TABLE]

We show the following Lemma, to obtain a restriction on $x,y$ for the components of the $\mathbf{f}$-vector similar to that derived from Lemma 4.5 for the intrinsic volumes:

Lemma 4.8**.**

Fix a $d$-dimensional polytope $P\subset K$ contained in a convex body $K\in\mathcal{K}^{d}$ and two points $x,y\in K\setminus P$. Let $\mathcal{V}_{x}(P)$ and $\mathcal{V}_{y}(P)$ denote the visibility regions of $x$ and $y$ with respect to $P$ be defined as in Lemma 4.5. If $\mathcal{V}_{x}(P)\cap\mathcal{V}_{y}(P)=\emptyset$, then

[TABLE]

for all $j\in\left\{{0,\ldots,d-1}\right\}$.

• Proof:

Denote by $P^{xy}$, $P^{x}$ and $P^{y}$ the convex hulls of $P\cup\left\{{x,y}\right\}$, $P\cup\left\{{x}\right\}$ resp. $P\cup\left\{{y}\right\}$, then we can decompose the number of $j$-faces of $P^{xy}$ into the number of $j$-faces of $P^{x}$ and the gained and lost $j$-faces obtained from adding $y$, i.e.

[TABLE]

Since the visible regions are disjoint we have $\mathcal{F}^{\operatorname{Vis}}(P^{x},y)=\mathcal{F}^{\operatorname{Vis}}(P,y)$ and $\operatorname{link}(P^{xy},y)=\operatorname{link}(P^{y},y)$, thus

[TABLE]

Similar to $f_{j}(P^{xy})$ we can decompose $f_{j}(P)$ by counting the $j$-faces in $P^{y}$ and subtracting the difference that arises from the addition of $y$ to $P$:

[TABLE]

yielding

[TABLE]

which is the desired conclusion. ∎

Since $\operatorname{Vis}_{x}(t)\cap\operatorname{Vis}_{y}(t)=\emptyset$ implies the conditions of Lemma 4.5 for $P=K_{t}\supseteq K(v\geq\varepsilon_{t})$ it follows that

[TABLE]

and similar to $D_{x,y}^{2}V_{j}(K_{t})$ we derive

[TABLE]

Having disposed of this preliminary steps, we can now summarize the results in the following Lemma bounding the order of magnitude of the $p$-th moment of the difference operator of the $\mathbf{f}$-vector components $f_{j}(K_{t})$.

Lemma 4.9**.**

Let $p\in\left\{{1,\ldots,8}\right\}$, $j\in\left\{{0,\ldots,d-1}\right\}$ and $x,y\in K$, then

[TABLE]

• Proof:

The proof is similar to the one of Lemma 4.6. ∎

We are left with the task of applying our estimations on the bound $\gamma_{1}$, $\gamma_{2}$ and $\gamma_{3}$ given by Theorem 2.3. Since we consider the multivariate functional given by (6) we have to distinguish the following three cases depending on the combination of functionals $F_{i}$ and $F_{j}$ using the corresponding bounds for the variance given by (14) for the intrinsic volumes and by

[TABLE]

for the components of the $\mathbf{f}$-vector, see [Rei05]. We denote by $\gamma_{1}(i,j)$, $\gamma_{2}(i,j)$ resp. $\gamma_{3}(i)$ the integral in $\gamma_{1}$, $\gamma_{2}$ resp. $\gamma_{3}$, then it follows from Lemma 4.6, 4.9 and the estimations on the domain of integration (30) and (31), that

[TABLE]

and

[TABLE]

It can easily be seen that the speed of convergence is dominated by the case $i,j\in\left\{{d+1,\ldots,2d}\right\}$ thus we can rewrite the bound in Theorem 2.3 using $\sigma_{ij}(t)={\mathrm{Cov}}\left[F_{i},F_{j}\right]$ to

[TABLE]

which completes the proof.

4.3 Proof of Theorem 1.3: oracle functional

Recall that the oracle estimator is given by $\hat{\vartheta}_{oracle}(K_{t})=V_{d}(K_{t})+t^{-1}f_{0}(K_{t})$, and its variance asymptotics is given by (4),

[TABLE]

for $t\rightarrow\infty$, where the constant $\gamma_{d}$ only depends on the dimension and is known explicitly and $\Omega(K)$ denotes the affine surface area of $K$. Rescaling of $\vartheta(K_{t})$ yields

[TABLE]

where the scaling $t^{\frac{1}{2}+\frac{1}{d+1}}$ resp. $t^{-\frac{1}{2}+\frac{1}{d+1}}$ corresponds with the asymptotic variance of $V_{d}(K_{t})$ resp. $f_{0}(K_{t})$, see (14) and (32). Therefore we can use the previous results to deduce bounds on the moments of the first and second order difference operator of the standardized oracle estimator $\tilde{\vartheta}(K_{t})$.

Lemma 4.10**.**

Let $p\in\left\{{1,\ldots,4}\right\}$ and $x,y\in K$, then

[TABLE]

• Proof:

Using the binominal theorem and the Cauchy–Schwarz inequality it follows directly with Lemma 4.6 and Lemma 4.9 that

[TABLE]

Since $j+d(p-j)\leq dp$ the desired conclusion follows. The proof for the second order difference operator is similar. ∎

Applying these estimations to the bound $\tau_{1}$, $\tau_{2}$ and $\tau_{3}$ in Theorem 2.2 yields

[TABLE]

thus

[TABLE]

for $Z\overset{\mathclap{d}}{\sim}\mathcal{N}(0,1)$, completing the proof of Theorem 1.3.

Acknowledgments. The author wishes to express his thanks to Matthias Schulte for helpful suggestions on the multivariate limit theorem and to Martina Juhnke-Kubitzke for several helpful comments concerning Lemma 4.5 and 4.8.