\section{Introduction} \label{sec:crackle:intro}

This paper treats the homology of simplicial complexes built via deterministic rules from a  random set of vertices. In particular, it shows that, depending on the randomness that generates the vertices, the homology of these complexes can either become trivial as the sample size grows, or can contain more and more complex structures.

The motivation for these results comes from applications of topological tools for pattern analysis, object identification, and especially for the analysis of data sets.  Typically, one starts with a collection of points and forms some simplicial complexes associated to these, and then takes their homology.  For example, the $0$-dimensional homology of such complexes can be interpreted as a version of clustering.  The basic philosophy behind this attempt is that topology has an essentially qualitative nature and should therefore be robust with respect to small perturbations.
Some recent references are \cite{aronshtam_vanishing_2010,babson_fundamental_2011,cohen_homotopical_2010,meshulam_homological_2009,pippenger_topological_2006} with two reviews, from different aspects, in \cite{adler_persistent_2010} and \cite{ghrist_barcodes:_2008}. Many of these papers find their {\it raison d'\^etre} in essentially
statistical problems, in which data generates the structures.

An important example occurs in the following manifold learning problem. Let $\mathcal{M}$ be an unknown manifold embedded in a Euclidean space, and suppose that we are a given a set of independent and identically distributed ($\mathrm{i.i.d.}$) random samples $\mathcal{X}_n = \set{X_1,\ldots,X_n}$ from the manifold. In order to recover the homology of $\mathcal{M}$, we consider  the homology of
\begin{eqnarray}
\label{U:equn}
U \ =\  \bigcup_{k=1}^n B_{\epsilon}(X_k),
\end{eqnarray}
where $B_\epsilon(X)$ is the Euclidean ball, in the ambient space,   of radius $\epsilon$ about the point $X$. The belief, or hope,  is that, for large enough $n$, the homology of $U$ will be equivalent to that of $\mathcal{M}$. A confounding issue arises
when the sample points do not necessarily lie on the manifold, but rather are perturbed from it by a random amount.
When this happens, it will follow from our results that the precise distribution behind the randomness plays a qualitatively important role. It is known that if the perturbations come from a bounded or strongly concentrated distribution, then they do not lead to much spurious homology, and the above line of attack, appropriately applied, works. For example,
it was shown in \cite{niyogi_topological_2011}  that for Gaussian noise it is  possible to clean the data and recover the underlying topology of $\mathcal{M}$ in a way that is essentially independent on the ambient dimension.
Both \cite{niyogi_finding_2008,niyogi_topological_2011} contain results of the form that, given  a nice enough $\mathcal{M}$, and any $\delta>0$,  there are explicit conditions on $n$ and $\epsilon$ such that the homology of $U$ is equal to the homology of $\mathcal{M}$ with a probability of at least $(1-\delta)$.  However, for  other distributions no such results exist, nor, in view of the results of this paper,  are they to be expected.

Figure \ref{fig:mfld} provides an illustrative example of what happens when sampling points from an
annulus and perturbing them with additional noise before reconstructing the annulus as in  \eqref{U:equn}.
In particular, it shows that if the additional noise is in some sense large then  sample points can appear basically anywhere, introducing extraneous homology elements.

\begin{figure}[hc!]
\label{fig:mfld}
\centering
\subfigure[]
{
  \centering
    {\includegraphics[scale=0.2]{ml1}}
    \label{fig:mfld_orig}
}
\subfigure[]
{
  \centering
    {\includegraphics[scale=0.2]{ml2}}
    \label{fig:mfld_no_noise}
}
\subfigure[]
{
  \centering
    {\includegraphics[scale=0.2]{ml3}}
    \label{fig:mfld_noise_low}
}
\subfigure[]
{
  \centering
    {\includegraphics[scale=0.2]{ml4}}
    \label{fig:mfld_noise_high}
}

\caption{(a) The original space  $\mathcal{M}$ (an annulus) that we wish to recover from random samples. (b) With the appropriate choice of radius, we can easily recover the homology of the original space from random samples from $\mathcal{M}$. (c) In the presence of bounded noise, homology recovery is undamaged. (d) In the presence of unbounded noise, many extraneous homology elements appear, and significantly interfere with homology recovery.}
\label{fig:complexes:simp_comp}
\end{figure}

In order to  be able, eventually, to extend the work in \cite{niyogi_topological_2011} beyond Gaussian noise, and make more concrete statements about the probabilistic features of the homology this extension generates, it is necessary to first focus on the behaviour of samples generated by pure noise, with no underlying manifold. In this case, thinking of the above setup,  the manifold $\mathcal{M}$ is simply the point at the origin, and the homology that we shall be trying to recapture is trivial. Nevertheless, we shall see that differing noise models can make this task extremely delicate, regardless of sample size.

\subsection{Some sample results}
To start being more concrete, let $$\mathcal{X}_n = \set{X_1,\ldots,X_n}$$ be a set of $n$ $\mathrm{i.i.d.}$ random samples in $\mathbb{R}^d$, from a common density function $f$. Recall that the
abstract simplicial complex  $\check{C}(\mathcal{X},\varepsilon)$  constructed according to the following rules is called the \v{C}ech  complex associated to $\mathcal{X}$ and $\varepsilon$:
\begin{enumerate}
\item The $0$-simplices  of  $\check{C}(\mathcal{X},\varepsilon)$   are the points in $\mathcal{X}$,
\item An $n$-simplex $\sigma=[x_{i_0},\ldots,x_{i_n}]$ is in $\check{C}(\mathcal{X},\varepsilon)$ if $\bigcap_{k=0}^{n}
B_{x_{i_k}}(\varepsilon) \ne  \emptyset$,
\end{enumerate}
An important result, known as the `nerve theorem', links \v{C}ech  complexes
 and the neighborhood set $U$ of \eqref{U:equn},  establishing that they are  homotopy equivalent  (cf.\ \cite{borsuk_imbedding_1948}).  In particular,  they have the same  Betti numbers, measures of homology that we shall concentrate on in what follows.

 If the sample distribution has a compact support $S$, then it is easy to show that,   for large enough $n$,
 \begin{eqnarray*}
\check{C}(\mathcal{X},\varepsilon) \ \simeq \ \bigcup_{k=1}^n B_\varepsilon(X_k) \ \approx\  \tube(S,\varepsilon)  \ \triangleq \ \big\{x\in \mathbb{R}^d: \min_{y\in S} \|x-y\|\,\leq \, \varepsilon \big\},
\end{eqnarray*}}\def\mathcal{C}{\mathcal{C}
where $\simeq$ denotes homotopy equivalence and $\norm{\cdot}$ is the standard $L^2$ norm in $\mathbb{R}^d$. Thus, there is not much to study in this case. However, when the support of the distribution is unbounded, interesting phenomena occur.

To study these phenomena, we shall consider three representative examples of probability densities. These are the  \emph{power-law}, \emph{exponential}, and the \emph{standard Gaussian} distributions,
whose density functions are given, respectively, by
\begin{align}
\label{eq:crackle:dist_pl}
f_{\mathrm{p}}(x) &\triangleq \frac{c_{\mathrm{p}}}{1+\norm{x}^\alpha},\\
\label{eq:crackle:dist_exp}
f_{\mathrm{e}}(x) &\triangleq c_{\mathrm{e}} e^{-\norm{x}},\\
\label{eq:crackle:dist_norm}
f_{\mathrm{g}}(x) &\triangleq c_{\mathrm{g}} e^{-\norm{x}^2/2},
\end{align}
where $\alpha > d$ and $c_{\mathrm{p}},c_{\mathrm{e}},c_{\mathrm{g}}$ are  appropriate normalization constants that will not be of concern to us.

For large samples from any of these distributions we shall show that  there exists a `core' - a region in which the density of points is very high and so
placing unit balls around them completely covers the region. Consequently, the \v{C}ech  complex inside the core is contractible. The size of the core obviously grows to infinity as the sample size $n$ goes to infinity, but its exact size will depend on the underlying distribution. For the three examples above, if we
denote the radius of the core by $R_n^{\mathrm{c}}$, we shall prove in Section \ref{sec:crackle:core} that
\[
R_n^{\mathrm{c}} \sim \begin{cases} (n/\log n)^{1/\alpha} & f(x) \propto \frac{1}{1+\norm{x}^\alpha} ,\\
\log n & f(x) \propto e^{-\norm{x}}, \\
\sqrt{2\log n}  & f(x) \propto e^{-\norm{x}^2/2}.
 \end{cases}
\]
Note that in all three cases we have tacitly assumed that the  cores are balls, a natural consequence of the spherical symmetry of the probability densities.

Beyond the core, the topology is more varied.
For fixed $n$,  there may be additional isolated components, but no longer enough placed densely  enough to connect with one another and to form  a contractible set. Indeed, we shall show that the individual components will typically have enough homology to be, individually, non-contractible.
Thus, in this region, the topology of the \v{C}ech  complex is highly nontrivial, and many homology elements of different orders appear.
We call this phenomenon `crackling', akin to the well known phenomenon caused by noise interference in audio signals and commonly referred to as crackling.

As for core size, the exact crackling behaviour depends on the choice of distribution.  It turns out that  Gaussian samples do not lead to crackling, but the other two cases do. To describe this, with some imprecision of notation we shall write $[a,b)$ not only for an interval on the real line, but also for the annulus
 \begin{eqnarray*}
 [a,b) \ \triangleq \  \big\{ x\in \mathbb{R}^d: a\leq \|x\| < b\big\}.
 \end{eqnarray*}}\def\mathcal{C}{\mathcal{C}
 In Sections \ref{sec:crackle:pl} and \ref{sec:crackle:exp} we  shall show that the exterior of the core can be divided into disjoint spherical annuli at radii
$$R_n^c  \ll  R_{d-1,n} \ll R_{d-2,n} \ll \cdots \ll R_{0,n}$$
 (defined differently for each of the two crackling distributions) with different types of crackling (i.e.\ of homology) dominating in different regions.

 In $[R_{0,n},\infty)$ there are mostly disconnected points, and no structures with nontrivial homology.
 In $[R_{1,n},R_{0,n})$ connectivity is a bit higher, and a finite number of $1$-cycles appear. In $[R_{2,n},R_{1,n})$ we have a finite number of $2$-cycles, while the number of $1$-cycles grows to infinity as $n\to\infty$. In general, in $[R_{k,n},R_{k-1,n})$, as $n\to\infty$ we have a finite number of $k$-cycles, infinitely many $l$-cycles for $l<k$, and no cycles of dimension $l>k$. In other words, the crackle starts with a pure dust at $R_{n,0}$ and as we get closer to the core, higher dimensional homology  gradually appears.  See Figure \ref{fig:crackle:layers} in the following section
for more details.

As we already mentioned, the Gaussian distribution is fundamentally different than the other two, and does not lead to crackling. In Section \ref{sec:crackle:norm} we show that, for the Gaussian distribution, there are hardly any points located outside the core. Thus, as $n\to \infty$, the union of balls around the sample points becomes a giant contractible ball of radius of order $\sqrt{2\log n}$.

It is now possible to understand a little better how the results of this paper relate to  the noisy manifold learning problem discussed above. For example, if the distribution of the noise is Gaussian, our results imply that if the manifold is well behaved, and the sample size is moderate, noise  outliers should not significantly interfere with homology recovery, since Gaussian noise does not introduce artificial homology elements with large samples. However, there is a delicate counterbalance here between `moderate' and `large'. Once the sample size is large,  the core is also large, and the reconstructed manifold will have the topology of
$\mathcal{M}\oplus B_{O(\sqrt{2\log n})}(0)$, where $\oplus$ is Minkowski addition. As $n$ grows, the core will eventually envelope any compact manifold, and thus the homology of $\mathcal{M}$ will be hidden by that of the core.

 On the other hand, if the distribution of the noise is power-law or exponential, then noise outliers will typically generate extraneous homology elements that, for almost any sample size, will complicate the estimation of the original manifold.
 Furthermore, increasing the sample size in no way solves this problem.
 Note that this issue is in addition to the fact that increasing the sample size will, as in the Gaussian case,
 create the problem of a large core concealing the  topology
 of $\mathcal{M}$.

 Thus, from a practical point of view, the message of this paper is that outliers cause problems in manifold estimation when  noise is present, a fact well known to all practitioners who have worked in the area. What is qualitatively new here is a quantification of how this happens, and how it relates to the distribution of the noise. We do not attempt {to solve this problem here, but unfortunately it follows from the results of this paper that algorithms for handling outliers will probably involve knowing at least the tail behaviour of the error distribution, despite the fact that in practical situations one does not generally want to take as known prior knowledge.

 \subsection{On persistence intervals}
 While the above discussion has concentrated on the persistence of noise induced crackle as sample sizes grow, and the regions in $\mathbb{R}^d$ in which different types of homology appear, the proofs below also yield information about the more classical persistence diagrams of topological data analysis (cf.\ \cite{carlsson_topology_2009,edelsbrunner_persistent_2008,edelsbrunner_book,ghrist_barcodes:_2008}).

  For example, in the two cases for which crackle persists -- the power-law and exponential cases  -- estimates of the type appearing in Section \ref{proofs:sec} indicate that, with high probability, there exist extremely long bars in the bar code representation of persistent homology. Up to lower order corrections, preliminary calculations show that  bar lengths for the $k$-th homology can be as large as  $O(n^{\alpha_k})$ for the power-law case, and  $\beta_k (\log \log n)$ for the exponential case, for appropriate $\alpha_k$ and $\beta_k$. More detailed studies of these phenomena will appear in a later publication.

\subsection{Poisson Processes}
\label{poisson:subsec}
Although we have described everything so far in terms of a random sample   $\mathcal{X}$ of $n$ points taken from a density $f$, there is another way to approach the results of this paper, and that is to replace the points of $\mathcal{X}$ with the points of a
$d$-dimensional Poisson process $\mathcal{P}_n$ whose intensity function is given by $\lambda_n = n f$. In this case the number of points is no longer fixed, but has mean $n$.

All the results of this paper stated for $\mathcal{X}$ hold, without any change, if we replace $\mathcal{X}$ by $\mathcal{P}$.

\subsection{Disclaimers}
Before starting the paper in earnest, and so as not to be accused of myopia, we note that the subject of manifold learning is obviously much broader that that described above, and  algorithms for `estimating' an underlying manifold from
a finite sample abound in the statistics and computer science literatures.
Very few of them, however, take an algebraic point of view that we or the literature quoted above take.
Furthermore, we note that other important results about the homology of Rips and \v{C}ech  complexes for various distributions can be found in the papers \cite{kahle_random_2011,kahle_limit_2010} and  \cite{bobrowski_distance_2011} .
However, the methods and emphases of these papers are rather different.

\section{Results}
\label{sec:results}
In this section we shall present all our main results, along with some  discussion, more technical than that of the Introduction. Recall from Section \ref{poisson:subsec} that although we present all results for the point set $\mathcal{X}$, they also hold if we replace the points of $\mathcal{X}$ by the points of an appropriate Poisson process.
All proofs are deferred to Section \ref{proofs:sec}.

\subsection{The Core of Distributions with Unbounded Support} \label{sec:crackle:core}

We start by examining the core of the power-law, exponential and Gaussian distributions. These distributions are spherically symmetric and the samples are concentrated near the origin.
By `core' we refer to a centered ball $B_{R_n} \triangleq B_{R_n}(0) \subset \mathbb{R}^d$ containing a very large number of points from the sample $\mathcal{X}_n$, such that
\[
    B_{R_n} \subset \bigcup_{X\in \mathcal{X}_n \cap B_{R_n}} B_1(X).
\]
i.e.\! the unit balls around the sample points completely cover $B_{R_n}$.
In this case the homology of $\bigcup_{X\in\mathcal{X}_n\cap B_{R_n}} B_1(X)$, or equivalently, of $\check{C}(\mathcal{X}_n \cap B_{R_n},1)$, is trivial.
Obviously, as $n\to\infty$, the radius $R_n$ grows as well.

Let $\set{R_n}_{n=1}^\infty$ be an increasing sequence of positive numbers.
Define by $C_n$ the event that $B_{R_n}$ is covered, i.e.
\[
    C_n \triangleq \set{B_{R_n} \subset \bigcup_{X\in \mathcal{X}_n \cap B_{R_n}} B_1(X)}.
\]
We wish to find the largest possible value of $R_n$ such that
$\prob{C_n} \to 1$.
The following theorem presents lower bounds for this value.

\begin{theorem}\label{thm:crackle:core} Let $\epsilon >0$, and define
\[
R_n^{\mathrm{c}} \triangleq \begin{cases} \param{\frac{\delta_{\mathrm{p}} n}{\param{\log n - e^{-\epsilon} \log \log n}}-1}^{1/\alpha} & f = f_{\mathrm{p}}, \\
\log n - \log\log\log n -\delta_{\mathrm{e}}-\epsilon & f = f_{\mathrm{e}}, \\
\sqrt{2\param{\log n -\log\log\log n -\delta_{\mathrm{g}}-\epsilon}}  & f = f_{\mathrm{g}},
 \end{cases}
\]
where the three distributions are given by \eqref{eq:crackle:dist_pl}--\eqref{eq:crackle:dist_norm}, and
\begin{align*}
\delta_{\mathrm{p}} &=  c_{\mathrm{p}}\alpha 2^{-d} d^{-(1+d/2)}, \\
\delta_{\mathrm{e}} &= (1+d/2)\log d +d\log 2-\log c_{\mathrm{e}} ,\\
\delta_{\mathrm{g}} &= (1+d/2)\log d + (d-1)\log 2 -\log c_{\mathrm{g}}.
\end{align*}
If $R_n \le R_n^{\mathrm{c}}$, then
\[
\prob{C_n} \to 1.
\]
\end{theorem}

Theorem \ref{thm:crackle:core} implies that the core size has a completely different order of magnitude for each of the three distributions. The heavy-tailed, power-law distribution has the largest core, while the core of the Gaussian distribution is the smallest.
In the following sections we shall study the behaviour of the \v{C}ech  complex outside the core.

\subsection{How Power-Law Noise Crackles} \label{sec:crackle:pl}

In this section we explore the crackling phenomenon in the power-law distribution $f=f_p$.
Let $B_{R_n}\subset \mathbb{R}^d$ be the centered ball with radius $R_n$, and let
\[
    \check{C}_n \triangleq \check{C}(\mathcal{X}_n \cap (B_{R_n})^c,1),
\]
be the \v{C}ech  complex constructed from sample points outside $B_{R_n}$.
We wish to study
\[
    \beta_{k,n} \triangleq \beta_k(\check{C}_n),
\]
the $k$-th Betti number of $\check{C}_n$.

Note that the minimum number of points required to form a $k$-dimensional cycle ($k\ge 1$) is $k+2$. For $k\ge 1$ and $\mathcal{Y}\subset \mathbb{R}^d$, denote
\[
    T_k(\mathcal{Y}) \triangleq \ind\set{\abs{\mathcal{Y}} = k+2,\ \beta_k(\check{C}(\mathcal{Y},1)) = 1},
\]
i.e.\! $T_k$ takes the value $1$ if $\check{C}(\mathcal{Y},1)$ is a minimal $k$-dimensional cycle, and $0$ otherwise. This indicator function will be used to define the limits of the Betti numbers.

\begin{theorem}\label{thm:crackle:mean_bk_pl}
If $\lim_{\ninf} n R_n^{-\alpha} = 0$, then
\begin{align*}
    \lim_{\ninf} \param{ n R_n^{d-\alpha} }^{-1}\mean{\beta_{0,n}} &= \mu_{\mathrm{p},0} , \\
     \lim_{\ninf}\param{n^{k+2} R_n^{d- \alpha(k+2)}}^{-1}\mean{\beta_{k,n}} &= \mu_{\mathrm{p},k} ,\quad 1 \le k \le d-1
\end{align*}
where
\begin{align}
\label{eq:crackle:mplo}
\mu_{\mathrm{p},0} &\triangleq \frac{s_{d-1} c_{\mathrm{p}}}{\alpha-d}, \\
\label{eq:crackle:mpl}
\mu_{\mathrm{p},k} &\triangleq \frac{s_{d-1}c_{\mathrm{p}}^{k+2}}{(\alpha(k+2)-d)(k+2)!}\int_{(\mathbb{R}^d)^{k+1}} T_k(0,\mathbf{y})d\mathbf{y}, \quad 1\le k\le d-1,
\end{align}
and where $s_{d-1}$ is the surface area of the $(d-1)$-dimensional unit sphere in $\mathbb{R}^d$.
\end{theorem}

Next, we define the following values, which will serve as critical radii for the crackle,
\begin{align*}
\Ro^\eps &\triangleq n^{\param{\frac{1}{\alpha-d}+\epsilon}}, \\
R_{0,n} &\triangleq R_{0,n}^0 ,\\
\Rk^\eps &\triangleq n^{\param{\frac{1}{\alpha-d/(k+2)}+\epsilon}}, \quad (k\ge 1) \\
R_{k,n} &\triangleq R_{k,n}^0.
\end{align*}
The following is a straightforward corollary of Theorem \ref{thm:crackle:mean_bk_pl}, and summarizes the behaviour of $\mean{\beta_{k,n}}$ in the power-law case.

\begin{corollary}\label{cor:crackle:mean_bk_pl}
For $k\ge 0$ and $\epsilon >0$,
\[
    \lim_{\ninf}\mean{\beta_{k,n}} = \begin{cases}
    0 & R_n = \Rk^\eps, \\
    \mu_{\mathrm{p},k} & R_n = R_{k,n} ,\\
    \infty & R_n = R_{k,n}^{-\epsilon},
    \end{cases}
\]
\end{corollary}

Theorem \ref{thm:crackle:mean_bk_pl} and Corollary \ref{cor:crackle:mean_bk_pl} reveal that the crackling behaviour is organized into separate `layers', see Figure \ref{fig:crackle:layers}. Dividing $\mathbb{R}^d$ into a sequence of annuli at radii
\[
    \Ro^\eps \gg R_{0,n} \gg R_{1,n}^\epsilon \gg R_{1,n} \gg \cdots \gg R_{d-1,n}^\epsilon \gg R_{d-1,n} \gg R_n^{\mathrm{c}},
\]
we observe a different behaviour of the Betti numbers in each annulus.
We shall briefly review the behaviour in each annulus, in a decreasing order of radii values. The following description is mainly qualitative, and refers to expected values only.
\begin{itemize}
\item $[\Ro^\eps,\infty)$ - there are hardly any points ($\beta_k\sim 0$, $0\le k \le d-1$).
\item $[R_{0,n},\Ro^\eps)$ - points start to appear, and $\beta_0\sim \mu_{\mathrm{p},0}$. The points are very few and scattered, so no cycles are generated ($\beta_k \sim 0$, $1\le k \le d-1$).
\item $[R_{1,n}^\epsilon,R_{0,n})$ - the number of components grows to infinity, but no cycles are formed yet ($\beta_0 \sim \infty$, and $\beta_k = 0$, $1 \le k \le d-1$).
\item $[R_{1,n},R_{1,n}^\epsilon)$ - a finite number of $1$-dimensional cycles show up, among the infinite number of components ($\beta_0 \sim \infty$, $\beta_1\sim\mu_{\mathrm{p},1}$, and $\beta_k = 0$, $1 \le k \le d-1$).
\item $[R_{2,n}^\epsilon,R_{1,n})$ - we have $\beta_0\sim\infty$, $\beta_1\sim\infty$, and $\beta_k\sim 0$ for $k\ge 1$.
\end{itemize}
This process goes on, until the $(d-1)$-dimensional cycles appear -
\begin{itemize}
\item $[R_{d-1},R_{d-1}^\epsilon)$ - we have $\beta_{d-1}\sim \mu_{\mathrm{p},d-1}$ and $\beta_k\sim\infty$ for $0\le k \le d-2$.
\item $[R_n^{\mathrm{c}},R_{d-1})$ - just before we reach the core, the complex exhibits the most intricate structure, with $\beta_k \sim \infty$ for $0\le k \le d-1$.
\end{itemize}

Note that there is a very fast phase transition as we move from the contractible core to the first crackle layer. At this point we do not know exactly where and how this phase transition takes place. A reasonable conjecture would be that the transition occurs at $R_n = n^{1/\alpha}$ (since at this radius the term $n R_n^{-\alpha}$ that appears in Theorem \ref{thm:crackle:mean_bk_pl} changes its limit, affecting the limiting Betti numbers). However, this remains for future work.

\begin{figure}[h]
\centering
  \includegraphics[scale=0.6]{crackle}
\caption[Crackle layers]{The layered behaviour of crackle. Inside the core ($B_{R_n^{\mathrm{c}}}$) the complex consists of a single component and no cycles. The exterior of the core is divided into separate annuli. Going from right to left, we see how the Betti numbers grow. In each annulus we present the Betti number that was most recently changed.
}
\label{fig:crackle:layers}
\end{figure}

\subsection{How Exponential Noise Crackles} \label{sec:crackle:exp}

In this section we focus on the exponential density function $f=f_{\mathrm{e}}$.
The results in this section are very similar to the those for the power law distribution, and we shall describe them briefly. Differences lie in the specific values of the $R_{k,n}$ and in the terms in the limit formulae.

\begin{theorem}\label{thm:crackle:mean_bk_exp}
If $\lim_{\ninf} ne^{-R_n} = 0$, then
\begin{align*}
    \lim_{\ninf} \param{n R_n^{d-1} e^{-R_n}}^{-1}\mean{\beta_{0,n}} &= \mu_{\mathrm{e},0} , \\
     \lim_{\ninf}\param{n^{k+2} R_n^{d-1} e^{-(k+2)R_n}}^{-1} \mean{\beta_{k,n}} &= \mu_{\mathrm{e},k} ,\ k\ge 1
\end{align*}
where
\begin{align}
\label{eq:crackle:mexpo}
\mu_{\mathrm{e},0} &\triangleq s_{d-1}c_{\mathrm{e}},\\
\label{eq:crackle:mexp}
\mu_{\mathrm{e},k} &\triangleq \frac{s_{d-1}c_{\mathrm{e}}^{k+2}}{(k+2)!} \int_0^\infty \int_{(\mathbb{R}^d)^{k+1}}T_k(0,\mathbf{y}) e^{-\param{(k+2)\rho + \sum_{i=1}^{k+1}y_i^1}} \prod_{i=1}^{k+1} \ind\set{y_i^1 > -\rho} d\mathbf{y} d\rho,
\end{align}
and where $y_i^1$ is the first coordinate of $y_i\in\mathbb{R}^d$.
\end{theorem}

Next, define
\begin{align*}
\Ro^\eps &\triangleq  \log n + \param{d-1+\epsilon}\log \log n, \\
R_{0,n} &\triangleq R_{0,n}^0 ,\\
\Rk^\eps &\triangleq \log n + \param{\frac{d-1}{k+2}+\epsilon}\log \log n, \quad (k\ge 1)\\
R_{k,n} &\triangleq R_{k,n}^0.
\end{align*}
From Theorem \ref{thm:crackle:mean_bk_exp} we can conclude the following.

\begin{corollary}\label{cor:crackle:mean_bk_exp}
For $k\ge 0$ and $\epsilon >0$,
\[
    \lim_{\ninf}\mean{\beta_{k,n}} = \begin{cases}
    0 & R_n = \Rk^\eps, \\
    \mu_{\mathrm{e},k} & R_n = R_{k,n} ,\\
    \infty & R_n = R_{k,n}^{-\epsilon},
    \end{cases}
\]
\end{corollary}

As in the power-law case, Theorem \ref{thm:crackle:mean_bk_exp} implies the same `layered' behaviour, the only difference being in the values of $R_{k,n}$. From examining the values of $R_n^{\mathrm{c}}$, and $R_{k,n}$ it is reasonable to guess that the phase transition in the exponential case occurs at $R_n = \log n$.

\subsection{Gaussian Noise Does Not Crackle} \label{sec:crackle:norm}

Simplicial complexes built over vertices sampled from
the standard Gaussian distribution exhibit a completely different behaviour to that we saw in the power-law and  exponential cases. Define
\[
\Ro^\eps \triangleq \sqrt{2\log n + (d-2+\epsilon)\log\log n},
\]
then
\begin{theorem}\label{thm:crackle:mean_bk_norm}
If $f=f_{\mathrm{g}}$, $\epsilon > 0$, and $R_n = \Ro^\eps$, then for $0\le k \le d-1$
\[
\lim_{\ninf} \mean{\beta_{k,n}} = 0.
\]
\end{theorem}
Note that in the Gaussian case $\lim_{\ninf}\param{ \Ro^\eps - R_n^{\mathrm{c}}} = 0$. This implies that as $n\to\infty$ we have the core which is contractible, and outside the core there is hardly anything. In other words, the ball placed around every new point we add to the sample immediately connects to the core, and thus, the Gaussian noise \textit{does not crackle}.

\section{Proofs}
\label{proofs:sec}
We now turn to proofs, starting with the proof of the main result  of Section \ref{sec:crackle:core}.

\subsection{The Core}

\begin{proof}[Proof of Theorem \ref{thm:crackle:core}]
The proof covers all three distributions, except for specific calculations near the end.
Take a grid on $\mathbb{R}^d$ of size $g = \frac{1}{2\sqrt{d}}$.
Let ${\cal{Q}}_n$ be the collection of cubes in this grid that are contained in  $B_{R_n}$. Let $\tilde{C}_n$ be the following event
\[
    \tilde{C}_n \triangleq \set{\forall Q\in {\cal{Q}}_n : Q\cap \mathcal{X}_n \ne \emptyset},
\]
i.e.\! $\tilde{C}_n$ is the event that every cube in ${\cal{Q}}_n$ contains at least one point from $\mathcal{X}_n$. Recall the definition of $C_n$,
\[
    C_n \triangleq \set{B_{R_n} \subset \bigcup_{X\in \mathcal{X}_n \cap B_{R_n}} B_1(X)}.
\]
Then it is easy to show that $\tilde{C}_n \subset C_n$. The complementary event $\tilde{C}_n^c$ is the event that at least one cube is empty. Thus,
\[
\probx{\tilde{C}_n^c} \le \sum_{Q\in {\cal{Q}}_n} \prob{Q \cap \mathcal{X}_n = \emptyset} = \sum_{Q\in {\cal{Q}}_n}(1-p(Q))^n \le \sum_{Q\in {\cal{Q}}_n}e^{-np(Q)}
\]
where
\[
p(Q) = \int_Q f(z)dz \ge  g^d f(R_n).
\]
In addition, the number of cubes that are contained in $B_{R_n}$ is less than $\param{2{{R_n}/{g}}}^d$.
Therefore,
\begin{equation}\label{eq:crackle:empty_box}
\probx{\tilde{C}_n^c} \le (2 g^{-1})^d R_n^d e^{-n g^d f(R_n) }.
\end{equation}
Now, choose any $\epsilon > 0$ and set
\[
R_n = R_n^{\mathrm{c}} \triangleq \begin{cases} \param{\frac{\delta_{\mathrm{p}} n}{\param{\log n - e^{-\epsilon} \log \log n}}-1}^{1/\alpha} & f = f_{\mathrm{p}}, \\
\log n - \log\log\log n -\delta_{\mathrm{e}}-\epsilon & f= f_{\mathrm{e}}, \\
\sqrt{2\param{\log n -\log\log\log n -\delta_{\mathrm{g}}-\epsilon}}  & f = f_{\mathrm{g}},
 \end{cases}
\]
where
\begin{align*}
\delta_{\mathrm{p}} &=  c_{\mathrm{p}}\alpha 2^{-d} d^{-(1+d/2)}, \\
\delta_{\mathrm{e}} &= \log d -\log c_{\mathrm{e}} - \log g^d, \\
\delta_{\mathrm{g}} &= \log(d/2) -\log c_{\mathrm{g}} - \log g^d.
\end{align*}
It is easy to verify that in all cases we have
\[
R_n^d e^{-n g^d f(R_n) }  \to 0.
\]
Thus, from \eqref{eq:crackle:empty_box} we conclude that $\probx{\tilde{C}_n} \to 1$. Since $\prob{C_n} \ge \probx{\tilde{C}_n}$ we now have that for $R_n = R_n^{\mathrm{c}}$, in each of the distributions,
\[
\prob{C_n} \to 1,
\]
which completes the proof.
\end{proof}

\subsection{Crackle - Notation and General Lemmas}

For $R_n > 0$, set
\[
    \cX_{n,R_n} \triangleq \mathcal{X}_n \cap (B_{R_n})^c,
\]
i.e.\! $\cX_{n,R_n}$ consists of the points of $\mathcal{X}_n$ located outside the ball $B_{R_n}$. Next, recall the definition of $T_k$,
\[
    T_k(\mathcal{Y}) \triangleq \ind\set{\abs{\mathcal{Y}} = k+2,\ \beta_k(\check{C}(\mathcal{Y},1)) = 1},
\]
for $\mathcal{Y} \subset \mathbb{R}^d$, and write
\begin{align*}
    S_{0,n} & \triangleq \abs{\cX_{n,R_n}}, \\
    \hat{S}_{0,n} &\triangleq \#\set{X \in \cX_{n,R_n} : X \textrm{ is a connected component of } \check{C}(\mathcal{X}_n,1)}\\
    S_{k,n} &\triangleq \sum_{\mathcal{Y} \subset \cX_{n,R_n}} T_k(\mathcal{Y}),\\
    \widehat{S}_{k,n} &\triangleq \sum_{\mathcal{Y} \subset \cX_{n,R_n}} T_k(\mathcal{Y})\ind\set{\check{C}(\mathcal{Y},1) \textrm{ is a connected component of } \check{C}(\mathcal{X}_n,1)},\\
    L_{k,n} &\triangleq \sum_{\mathcal{Y} \subset \cX_{n,R_n}} \ind\set{\abs{\mathcal{Y}} = k+3,\ \check{C}(\mathcal{Y},1) \textrm{ is connected}},
\end{align*}
where $k\ge 1$. Observe that
\begin{align}
\label{eq:crackle:betti_0_ineq}
\hat{S}_{0,n} \le &\beta_{0,n} \le S_{0,n} \\
\label{eq:crackle:betti_k_ineq}
    \widehat{S}_{k,n} \le &\beta_{k,n} \le \widehat{S}_{k,n} + L_{k,n},\quad k\ge 1
\end{align}
We will evaluate the limits of $\mean{S_{k,n}}$, $\meanx{\widehat{S}_{k,n}}$ and $\mean{L_{k,n}}$ and deduce from these the limit of $\mean{\beta_{k,n}}$.

In addition,  set
\begin{align*}
    \mathrm{e}_1 &\triangleq (1,0,\ldots,0) \in \mathbb{R}^d ,\\
    f(r) &\triangleq f(r \mathrm{e}_1), \ \ r \in \mathbb{R}, \\
    U(\mathbf{x}) &\triangleq \bigcup_{i=1}^{k} B_2(x_i), \ \ \mathbf{x}\in (\mathbb{R}^d)^k, \\
    p(\mathbf{x}) &\triangleq \int_{U(\mathbf{x})}f(z)dz,\ \ \mathbf{x}\in (\mathbb{R}^d)^k.
\end{align*}
The following two lemmas are purely technical, but will considerably simplify our computations later.

\begin{lemma}\label{lem:crackle:symmetry_zero}
Let  $f:\mathbb{R}^d\to\mathbb{R}$ be a spherically symmetric probability density.
Then,
\begin{align*}
\mean{S_{0,n}} &=  s_{d-1}n \int_{R_n}^\infty r^{d-1}f(r)dr ,\\
\meanx{\hat{S}_{0,n}} &=  s_{d-1}n \int_{R_n}^\infty r^{d-1}f(r)(1-np(r\mathrm{e}_1))^{n-1}dr,
\end{align*}
where $s_{d-1}$ is the volume of the $d-1$ dimensional unit sphere.
\end{lemma}

\begin{proof}
$S_{0,n}$ is simply a sum of Bernoulli variables, therefore
\[
\mean{S_{0,n}}  = n \prob{\norm{X} > R_n} = n\int_{\mathbb{R}^d} f(x)\ind\set{\norm{x}>R_n}dx.
\]
Writing the integral in polar coordinates yields
\[
\mean{S_{0,n}} = n \int_{R_n}^\infty \int_{S^{d-1}} f(r\theta)r^{d-1}J(\theta)d\theta dr,
\]
where $J(\theta) = \abs{\frac{\partial x}{\partial \theta}}$.
Since $f$ is spherically symmetric, $f(r\theta) = f(r)$, and therefore
\[
\mean{S_{0,n}} = s_{d-1} n \int_{R_n}^\infty r^{d-1} f(r) dr.
\]
The proof for $\hat{S}_{0,n}$ is similar, using the fact that the probability that a point $x\in\mathbb{R}^d$ is disconnected from the rest of the complex $\check{C}(\mathcal{X}_n,1)$ is $(1-p(x))^{n-1}$.

\end{proof}

\begin{lemma}\label{lem:crackle:symmetry}
Let $f:\mathbb{R}^d\to\mathbb{R}$ be a spherically symmetric probability density.
Then, for $k\ge 1$,
\begin{align*}
\mean{S_{k,n}} &=  s_{d-1} \binom{n}{k+2} \int_{R_n}^\infty r^{d-1}f(r)G_k(r)dr ,\\
\meanx{\widehat{S}_{k,n}} &= s_{d-1} \binom{n}{k+2}\int_{R_n}^\infty r^{d-1}f(r)\hat{G}_k(r)dr,
\end{align*}
where $s_{d-1}$ is the volume of the $d-1$ dimensional sphere, and where
\begin{align*}
    G_k(r) &\triangleq \int_{(\mathbb{R}^d)^{k+1}} f(\norm{r\mathrm{e}_1+\mathbf{y}}) T_k(0,\mathbf{y}) \prod_{i=1}^{k+1} \ind\set{\norm{r\mathrm{e}_1+y_i} > R_n} d\mathbf{y}, \\
\hat{G}_k(r) &\triangleq \int_{(\mathbb{R}^d)^{k+1}} f(\norm{r\mathrm{e}_1+\mathbf{y}}) T_k(0,\mathbf{y}) \prod_{i=1}^{k+1} \ind\set{\norm{r\mathrm{e}_1+y_i} > R_n} \\
&\times (1-p(r\mathrm{e}_1, r\mathrm{e}_1+\mathbf{y}))^{n-k-2}d\mathbf{y}.
\end{align*}

\end{lemma}

\begin{proof}
The proof is in the same spirit of the proof of Lemma \ref{lem:crackle:symmetry_zero}, but technically more complicated.
Thinking of $S_{k,n}$ as a sum of Bernoulli variables, we have that
\[
\mean{S_{k,n}} =  \binom{n}{k+2}\int_{(\mathbb{R}^d)^{k+2}}f(\mathbf{x}) T_k(\mathbf{x}) \prod_{i=1}^{k+2} \ind\set{\norm{x_i} > R_n} d\mathbf{x}.
\]
Let $I_k$ denote the integral above. Then, using the change of variables
\begin{align*}
    x_1 &\to x, \qquad x_i \to x+ y_{i-1} ,\ \  (i>1),
\end{align*}
yields
\begin{align*}
I_k &= \int_{\norm{x}\ge R_n}  \int_{(\mathbb{R}^d)^{k+1}} f(x)f(x+\mathbf{y}) T_k(x,x+\mathbf{y}) \prod_{i=1}^{k+1} \ind\set{\norm{x+y_i} > R_n}d\mathbf{y} dx \\
&=\int_{\norm{x}\ge R_n}  \int_{(\mathbb{R}^d)^{k+1}} f(x)f(x+\mathbf{y}) T_k(0,\mathbf{y}) \prod_{i=1}^{k+1} \ind\set{\norm{x+y_i} > R_n}d\mathbf{y} dx.
\end{align*}
Moving to polar coordinates yields
\begin{align*}
I_k &= \int_{R_n}^\infty \int_{S^{d-1}}  \int_{(\mathbb{R}^d)^{k+1}} f(r\theta)f(r\theta+\mathbf{y}) T_k(0,\mathbf{y}) \\
&\times\prod_{i=1}^{k+1} \ind\set{\norm{r\theta+y_i} > R_n} r^{d-1}J(\theta)d\mathbf{y} d\theta dr\\
&= \int_{R_n}^\infty r^{d-1}f(r) \int_{S^{d-1}}  J(\theta)\int_{(\mathbb{R}^d)^{k+1}} f(\norm{r\theta+\mathbf{y}}) T_k(0,\mathbf{y}) \\
&\times \prod_{i=1}^{k+1} \ind\set{\norm{r\theta+y_i} > R_n} d\mathbf{y} d\theta dr,
\end{align*}
where $J(\theta) = \abs{\frac{\partial x}{\partial \theta}}$, and $f(x) = f(\norm{x})$ by the spherical symmetry assumption.
Set
\[
    G_k(r,\theta) \triangleq \int_{(\mathbb{R}^d)^{k+1}} f(\norm{r\theta+\mathbf{y}}) T_k(0,\mathbf{y}) \prod_{i=1}^{k+1} \ind\set{\norm{r\theta+y_i} > R_n} d\mathbf{y}.
\]
Since $T_k$ is rotation invariant, it is easy to show that for every $\theta \in S^{d-1}$
\[
    G_k(r,\theta) = G_k(r,\mathrm{e}_1) \triangleq G_k(r).
\]
Thus,
\begin{equation}\label{eq:crackle:I_k}
    I_k = s_{d-1}\int_{R_n}^\infty r^{d-1}f(r)G_k(r)dr.
\end{equation}
This completes the proof for $S_{k,n}$. The proof for $\widehat{S}_{k,n}$ is similar.

\end{proof}

In what follows, we shall use the following elementary limits:
\begin{enumerate}
\item For every $k > 0$,
\begin{equation} \label{eq:binom_limit}
    \lim_{\ninf} n^{-k} \binom{n}{k}  = \frac{1}{k!}
\end{equation}
\item For every sequence $a_n\to 0$ and $k\ge0$,
\begin{equation}\label{eq:power_exp}
    \lim_{\ninf} \frac{(1-a_n)^{n-k}}{e^{-na_n}} = 1
\end{equation}
\end{enumerate}

\subsection{Crackle - The Power Law Distribution}

In this section we prove the results in Section \ref{sec:crackle:pl}. First, we need a few lemmas.

\begin{lemma}\label{lem:crackle:pl_dust}
If $f=f_{\mathrm{p}}$, and $R_n\to\infty$, then
\[
\lim_{\ninf} \param{n R_n^{d-\alpha}}^{-1} \mean{S_{0,n}} = \mu_{\mathrm{p},0},
\]
where $\mu_{\mathrm{p},0}$ is defined in \eqref{eq:crackle:mplo}.

If, in addition, $nR_n^{-\alpha}\to 0$, then
\[
\lim_{\ninf} \param{ n R_n^{d-\alpha} }^{-1}\meanx{\hat{S}_{0,n}} = \mu_{\mathrm{p},0}.
\]
\end{lemma}

\begin{proof}
From Lemma \ref{lem:crackle:symmetry_zero} we have that
\[
    \mean{S_{0,n}} = s_{d-1}n \int_{R_n}^\infty r^{d-1}f(r)dr.
\]
Making the change of variables $r\to R_n \rho$ yields
\begin{align*}
    \mean{S_{0,n}} &= s_{d-1}n  \int_1^\infty \frac{c_{\mathrm{p}}(R_n\rho)^{d-1}}{1+ (R_n\rho)^\alpha}R_n d\rho  \\
    &= s_{d-1}c_{\mathrm{p}} n R_n^{d-\alpha} \int_1^\infty \frac{\rho^{d-1}}{R_n^{-\alpha}+ \rho^\alpha} d\rho.
\end{align*}
Applying the dominated convergence theorem to the previous integral gives
\[
\lim_{\ninf} \param{nR_n^{d-\alpha}}^{-1}\mean{S_{0,n}} = s_{d-1}c_{\mathrm{p}}\int_1^\infty \rho^{d-1-\alpha}d\rho = \frac{s_{d-1}c_{\mathrm{p}}}{\alpha-d} = \mu_{\mathrm{p},0}.
\]
This proves the first part of the lemma.

Next, from Lemma \ref{lem:crackle:symmetry_zero} we have that
\[
    \meanx{\hat{S}_{0,n}} = s_{d-1}n \int_{R_n}^\infty r^{d-1}f(r)(1-p(r\mathrm{e}_1))^{n-1}dr.
\]
The power term is bounded by $1$ and therefore will not affect the conditions needed for dominated convergence. Thus, using \eqref{eq:power_exp}, we only need to evaluate its limit.
\[
p(r\mathrm{e}_1) = \int_{B_2(r\mathrm{e}_1)}f(z)dz = \int_{B_2(0)}\frac{c_{\mathrm{p}}}{1+\norm{r\mathrm{e}_1+z}}dz,
\]
and after the change of variables $r\to R_n\rho$ we have,
\[
p(R_n\rho\mathrm{e}_1)  =c_{\mathrm{p}} R_n^{-\alpha}\int_{B_2(0)} \frac{1}{R_n^{-\alpha}+\norm{\rho\mathrm{e}_1+R_n^{-1}z}^\alpha}dz.
\]
If $nR_n^{-\alpha}\to 0$, then, by dominated convergence, we have
\[
\lim_{\ninf} np(R_n\rho\mathrm{e}_1) =0.
\]
Thus,
\[
        \lim_{\ninf} (1-p(R_n\rho\mathrm{e}_1))^{n-1} = \lim_{\ninf} e^{-np(R_n\rho\mathrm{e}_1)} =1,
\]
and therefore we have
\[
    \lim_{\ninf}\param{nR_n^{d-\alpha}}^{-1}\meanx{\hat{S}_{0,n}} = \lim_{\ninf}\param{nR_n^{d-\alpha}}^{-1}\mean{S_{0,n}} = \mu_{\mathrm{p},0}.
\]
This completes the proof of the second part of the lemma.
\end{proof}

\begin{lemma}\label{lem:crackle:pl_holes}
If $f=f_{\mathrm{p}}$, and $R_n\to\infty$ then
\[
\lim_{\ninf} \param{n^{k+2} R_n^{d- \alpha(k+2)}}^{-1} \mean{S_{k,n}} = \mu_{\mathrm{p},k},
\]
where $\mu_{\mathrm{p},k}$ is defined in \eqref{eq:crackle:mpl}.
If, in addition,$n R_n^{-\alpha} \to 0$, then
\[
\lim_{\ninf} \param{n^{k+2} R_n^{d-\alpha(k+2)}}^{-1} \meanx{\widehat{S}_{k,n}} = \mu_{\mathrm{p},k}.
\]
\end{lemma}

\begin{proof}
The proof is in the spirit of the proof of Lemma \ref{lem:crackle:pl_dust}, but technically more complicated.
From Lemma \ref{lem:crackle:symmetry} we have that
\[
\mean{S_{k,n}} = \binom{n}{k+2} I_k,
\]
where
\[
    I_k = s_{d-1}\int_{R_n}^\infty r^{d-1}f(r)G_k(r)dr.
\]
Making the change of variables $r \to R_n\rho$ yields
\begin{align*}
I_k &= s_{d-1}R_n\int_1^\infty (R_n\rho)^{d-1} f(R_n\rho)  G_k(R_n\rho) d\rho \\
&= s_{d-1}c_{\mathrm{p}}^{k+2}(R_n)^{d-\alpha(k+2)}\int_1^\infty \int_{(\mathbb{R}^d)^{k+1}} \frac{\rho^{d-1}}{R_n^{-\alpha} + \rho^\alpha}\prod_{i=1}^{k+1}\frac{1}{R_n^{-\alpha} +\norm{\rho\mathrm{e}_1 + R_n^{-1}y_i}^\alpha} \\
&\ \ \times T_k(0,\mathbf{y}) \prod_{i=1}^{k+1} \ind\set{\norm{\rho\mathrm{e}_1+ R_n^{-1}y_i} > 1} d\mathbf{y}.
\end{align*}
Thus, using \eqref{eq:binom_limit},
\begin{align*}
&(n^{k+2} R_n^{d-\alpha (k+2)})^{-1}\mean{S_{k,n}} = \frac{s_{d-1}c_{\mathrm{p}}^{k+2}}{(k+2)!}\int_1^\infty \int_{(\mathbb{R}^d)^{k+1}} \frac{\rho^{d-1}}{R_n^{-\alpha} + \rho^\alpha}\\
&\quad \times T_k(0,\mathbf{y}) \prod_{i=1}^{k+1} \frac{1}{R_n^{-\alpha} +\norm{\rho\mathrm{e}_1 + R_n^{-1}y_i}^\alpha}\ind\set{\norm{\rho\mathrm{e}_1+ R_n^{-1}y_i} > 1} d\mathbf{y}.
\end{align*}
It is easy to show that the integrand is bounded by an integrable term, so the dominated convergence theorem applies, yielding
\begin{align*}
&\lim_{\ninf} (n^{k+2} R_n^{d-\alpha (k+2)})^{-1}\mean{S_{k,n}}\\
 & \qquad \qquad= \frac{s_{d-1}c_{\mathrm{p}}^{k+2}}{(k+2)!}\int_1^{\infty}\rho^{d-1-\alpha(k+2)}d\rho\int_{(\mathbb{R}^d)^{k+1}} T_k(0,\mathbf{y})d\mathbf{y} \\
 & \qquad \qquad= \frac{s_{d-1}c_{\mathrm{p}}^{k+2}}{(\alpha(k+2)-d)(k+2)!}\int_{(\mathbb{R}^d)^{k+1}} T_k(0,\mathbf{y})d\mathbf{y}\\
 & \qquad \qquad=\mu_{\mathrm{p},k}.
\end{align*}
This proves the first part of the lemma.

Next, the terms $G_k(r)$ and $\hat{G}_k(r)$ in Lemma \ref{lem:crackle:symmetry} differ only by the term $(1-p(r\mathrm{e}_1, r\mathrm{e}_1+\mathbf{y}))^{n-k-2}$, so dominated convergence still applies.
Now,
\[
p(r\mathrm{e}_1, r\mathrm{e}_1+\mathbf{y}) = \int_{U(r\mathrm{e}_1,r\mathrm{e}_1+\mathbf{y})} f(z)dz = \int_{U(0,\mathbf{y})}f(r\mathrm{e}_1+z)dz,
\]
and substituting $r\to R_n\rho$ yields,
\[
p(R_n\rho\mathrm{e}_1 , R_n\rho\mathrm{e}_1+\mathbf{y}) = c_{\mathrm{p}} R_n^{-\alpha}\int_{U(0,\mathbf{y})} \frac{1}{R_n^{-\alpha}+\norm{\rho\mathrm{e}_1+R_n^{-1}z}^\alpha}dz.
\]
If $nR_n^{-\alpha}\to 0$, then using the dominated convergence we have
\[
\lim_{\ninf} np(R_n\rho\mathrm{e}_1,R_n\rho\mathrm{e}_1 + \mathbf{y}) = 0.
\]
Thus,
\[
    \lim_{\ninf} e^{-np(R_n\rho\mathrm{e}_1, R_n\rho\mathrm{e}_1 + \mathbf{y})} =1,
\]
and therefore, using \eqref{eq:power_exp},
\begin{align*}
\lim_{\ninf} \param{n^{k+2} R_n^{d-\alpha (k+2)}}^{-1}\meanx{\widehat{S}_{k,n}}  &= \lim_{\ninf} \param{n^{k+2} R_n^{d-\alpha (k+2)}}^{-1}\meanx{S_{k,n}}  \\
&= \mu_{\mathrm{p},k}.
\end{align*}
This completes the proof of the second part of the lemma.
\end{proof}

\begin{lemma}\label{lem:crackle:pl_other}
If $f=f_{\mathrm{p}}$, and $R_n\to\infty$ then
\[
\lim_{\ninf} \param{n^{k+3} R_n^{d-\alpha(k+3)}}^{-1} \mean{L_{k,n}} = \hat{\mu}_{\mathrm{p},k},
\]
for some $\hat{\mu}_{\mathrm{p},k} > 0$.
\end{lemma}

\begin{proof}
The proof is very similar to the proof of Lemma \ref{lem:crackle:pl_holes}. We need only replace $T_k$ with an indicator function that tests whether a  sub-complex generated by $k+3$ points is connected. The exact value of $\hat{\mu}_{\mathrm{p},k}$ will not be needed anywhere.
\end{proof}

We can now prove Theorem \ref{thm:crackle:mean_bk_pl}.

\begin{proof}[Proof of Theorem \ref{thm:crackle:mean_bk_pl}]
To prove the limit for $\beta_{0,n}$ simply combine Lemma \ref{lem:crackle:pl_dust} with the inequality \eqref{eq:crackle:betti_0_ineq}.
To prove the limit for $\beta_{k,n}$, $k\ge 1 $,  combine Lemmas \ref{lem:crackle:pl_holes} and \ref{lem:crackle:pl_other} with the inequality  \eqref{eq:crackle:betti_k_ineq}.
\end{proof}

\subsection{Crackle - The Exponential Distribution}

In this section we wish to prove Theorem \ref{thm:crackle:mean_bk_exp}. We start with the following lemmas.

\begin{lemma}\label{lem:crackle:exp_dust}
If $f=f_{\mathrm{e}}$, and $R_n\to\infty$ then,
\[
\lim_{\ninf} \param{n R_n^{d-1} e^{-R_n}}^{-1} \mean{S_{0,n}} = \mu_{\mathrm{e},0},
\]
where $\mu_{\mathrm{e},0}$ is defined in \eqref{eq:crackle:mexpo}.

If, in addition, $ne^{-R_n}\to 0$ then,
\[
\lim_{\ninf} \param{n R_n^{d-1} e^{-R_n}}^{-1} \meanx{\hat{S}_{0,n}} = \mu_{\mathrm{e},0}.
\]
\end{lemma}

\begin{proof}
From Lemma \ref{lem:crackle:symmetry_zero} we have that
\[
    \mean{S_{0,n}} = s_{d-1}n \int_{R_n}^\infty r^{d-1}f(r)dr.
\]
Using the change of variables $r\to \rho + R_n$ yields
\begin{align*}
    \mean{S_{0,n}} &= s_{d-1}n  \int_0^\infty (\rho+R_n)^{d-1}c_{\mathrm{e}} e^{-(\rho+R_n)} d\rho  \\
    &= s_{d-1}c_{\mathrm{e}} n R_n^{d-1}e^{-R_n} \int_0^\infty \param{\frac{\rho}{R_n}+1}^{d-1}e^{-\rho} d\rho.
\end{align*}
Applying dominated convergence to the last integral yields,
\[
\lim_{\ninf} \param{n R_n^{d-1} e^{-R_n}}^{-1} \mean{S_{0,n}} = s_{d-1}c_{\mathrm{e}}\int_0^\infty e^{-\rho}d\rho = s_{d-1}c_{\mathrm{e}} = \mu_{\mathrm{e},0}.
\]
This proves the first part of the lemma.

Next, from Lemma \ref{lem:crackle:symmetry_zero} we have that
\[
    \meanx{\hat{S}_{0,n}} = s_{d-1}n \int_{R_n}^\infty r^{d-1}f(r)(1-p(r\mathrm{e}_1))^{n-1}dr.
\]
The power term will not affect the dominated convergence conditions. Thus, we only need to evaluate its limit.
\[
p(r\mathrm{e}_1) = \int_{B_2(r\mathrm{e}_1)}f(z)dz = \int_{B_2(0)}c_{\mathrm{e}} e^{-\norm{r\mathrm{e}_1+z}}dz,
\]
and after the change of variables $r\to \rho+R_n$ we have,
\[
p((\rho+R_n)\mathrm{e}_1) = \int_{B_2(0)} c_{\mathrm{e}} e^{-\norm{(\rho+R_n)\mathrm{e}_1 + z}}dz \le e^{-(R_n+\rho)} \int_{B_2(0)} c_{\mathrm{e}} e^{\norm{z}}dz.
\]
If $ne^{-R_n}\to 0$, then
\[
\lim_{\ninf} np((\rho+R_n)\mathrm{e}_1) =0.
\]
Thus,
\[
    \lim_{\ninf} e^{-np((\rho+R_n)\mathrm{e}_1)} =1,
\]
and therefore, using \eqref{eq:power_exp}, we have
\[
   \lim_{\ninf} \param{n R_n^{d-1} e^{-R_n}}^{-1} \meanx{\hat{S}_{0,n}} = \lim_{\ninf} \param{n R_n^{d-1} e^{-R_n}}^{-1} \meanx{S_{0,n}} = \mu_{\mathrm{e},0}.
\]
This completes the proof of the second part of the lemma.
\end{proof}

\begin{lemma}\label{lem:crackle:exp_holes}
If $f=f_{\mathrm{e}}$, and $R_n\to\infty$ then,
\[
\lim_{\ninf} \param{n^{k+2} R_n^{d-1} e^{-(k+2)R_n}}^{-1} \mean{S_{k,n}} = \mu_{\mathrm{e},k},
\]
where $\mu_{\mathrm{e},k}$ is defined in \eqref{eq:crackle:mexp}.

If, in addition, $ne^{-R_n}\to 0$ then,
\[
\lim_{\ninf} \param{n^{k+2} R_n^{d-1} e^{-(k+2)R_n}}^{-1} \meanx{\widehat{S}_{k,n}} = \mu_{\mathrm{e},k}.
\]
\end{lemma}

\begin{proof}
From Lemma \ref{lem:crackle:symmetry} we have that
\[
\mean{S_{k,n}} = \frac{n^{k+2}}{(k+2)!} I_k,
\]
where
\[
    I_k = s_{d-1}\int_{R_n}^\infty r^{d-1}f(r)G_k(r)dr.
\]
Making the change of variables $r \to \rho + R_n$ yields
\begin{align*}
    I_k &= s_{d-1}\int_0^\infty (\rho+R_n)^{d-1}f(\rho+R_n)G_k(\rho+R_n)d\rho \\
    &= s_{d-1}c_{\mathrm{e}}^{k+2} \int_0^\infty \int_{(\mathbb{R}^d)^{k+1}} (\rho+R_n)^{d-1} e^{-(\rho+R_n)} \prod_{i=1}^{k+1} e^{-\norm{(\rho+R_n)\mathrm{e}_1 + y_i}} \\
    &\quad\times T_k(0,\mathbf{y}) \prod_{i=1}^{k+1} \ind\set{\norm{(\rho+R_n)\mathrm{e}_1+y_i} > R_n} d\mathbf{y} d\rho \\
    &= s_{d-1}c_{\mathrm{e}}^{k+2}e^{-(k+2)R_n}R_n^{d-1} \int_0^\infty \int_{(\mathbb{R}^d)^{k+1}} \param{\frac{\rho}{R_n}+1}^{d-1} e^{-\rho} \\
    &\quad\times T_k(0,\mathbf{y}) \prod_{i=1}^{k+1} e^{-\norm{(\rho+R_n)\mathrm{e}_1 + y_i}} e^{R_n} \ind\set{\norm{(\rho+R_n)\mathrm{e}_1+y_i} > R_n} d\mathbf{y} d\rho.
\end{align*}
The last integral can be easily shown to satisfy the  conditions of the dominated convergence theorem.
In addition, it is easy to show that
\[
    \lim_{\ninf}  e^{-\norm{(\rho+R_n)\mathrm{e}_1 +y_i}}e^{R_n}  = e^{-\param{\rho +\iprod{\mathrm{e}_1,y_i}}} = e^{-(\rho + y_i^1)},
\]
where $y_i^1$ is the first coordinate of $y_i \in \mathbb{R}^d$, and also that
\[
    \lim_{\ninf}\ind\set{\norm{(\rho+R_n)\mathrm{e}_1+y_i} > R_n} = \ind\set{y_i^1 \ge -\rho}.
\]
Altogether, we have that
\begin{align*}
    &\lim_{\ninf} \param{n^{k+2} R_n^{{d-1}} e^{-(k+2)R_n}}^{-1}  \mean{S_{k,n}} \\
    &= \frac{s_{d-1}c_{\mathrm{e}}^{k+2}}{(k+2)!} \int_0^\infty \int_{(\mathbb{R}^d)^{k+1}} T_k(0,\mathbf{y}) e^{-\param{(k+2)\rho + \sum_{i=1}^{k+1}y_i^1}} \prod_{i=1}^{k+1} \ind\set{y_i^1 \ge -\rho} d\mathbf{y} d\rho,
\end{align*}
proving the first part of the lemma.

Next, as in the proof of Lemma \ref{lem:crackle:pl_holes}, we need to evaluate the term $p(r\mathrm{e}_1, r\mathrm{e}_1+\mathbf{y})$.
\[
p(r\mathrm{e}_1, r\mathrm{e}_1+\mathbf{y}) = \int_{U(0,\mathbf{y})}c_{\mathrm{e}} e^{-\norm{r\mathrm{e}_1+z}}dz \le  \int_{U(0,\mathbf{y})}c_{\mathrm{e}} e^{-(r-\norm{z})}dz.
\]
The change of variables $r\to\rho+R_n$ yields
\[
p((\rho+R_n)\mathrm{e}_1,(\rho+R_n)\mathrm{e}_1+\mathbf{y}) \le e^{-R_n}e^{-\rho} \int_{U(0,\mathbf{y})}c_{\mathrm{e}} e^{\norm{z}}dz.
\]
If $ne^{-R_n}\to 0$, then
\[
\lim_{\ninf} n p((\rho+R_n)\mathrm{e}_1,(\rho+R_n)\mathrm{e}_1+\mathbf{y}) = 0.
\]
Thus,
\[
    \lim_{\ninf} e^{-np(R_n\rho\mathrm{e}_1, R_n\rho\mathrm{e}_1 + \mathbf{y})} =1,
\]
and therefore,
\begin{align*}
&\lim_{\ninf} \param{n^{k+2} R_n^{{d-1}} e^{-(k+2)R_n}}^{-1}  \meanx{\widehat{S}_{k,n}} \\
&\qquad \qquad = \lim_{\ninf} \param{n^{k+2} R_n^{{d-1}} e^{-(k+2)R_n}}^{-1}  \mean{S_{k,n}}= \mu_{\mathrm{e},k}.
\end{align*}
This completes the proof.
\end{proof}

\begin{lemma}\label{lem:exp_other}
If $f=f_{\mathrm{e}}$, and $R_n\to\infty$ then
\[
\lim_{\ninf} \param{n^{k+3} R_n^{d-1} e^{-(k+3)R_n}}^{-1} \mean{L_{k,n}} = \hat{\mu}_{\mathrm{e},k}.
\]
where $\hat{\mu}_{\mathrm{e},k} > 0$.
\end{lemma}

\begin{proof}
As for the proof of Lemma \ref{lem:crackle:pl_other}, mimic now  the proof of Lemma \ref{lem:crackle:exp_holes}, replacing $T_k$ with an indicator function that tests whether a sub-complex generated by $k+3$ points is connected.
\end{proof}

\begin{proof}[Proof of Theorem \ref{thm:crackle:mean_bk_exp}]
The proof follows the same steps as the proof of Theorem \ref{thm:crackle:mean_bk_pl}.
\end{proof}

\subsection{Crackle - The Gaussian Distribution}

In this section we  prove Theorem \ref{thm:crackle:mean_bk_norm}.

\begin{proof}[Proof of Theorem \ref{thm:crackle:mean_bk_norm}]
From Lemma \ref{lem:crackle:symmetry_zero} we have that
\[
    \mean{S_{0,n}} = s_{d-1}n \int_{R_n}^\infty r^{d-1}f(r)dr.
\]
Making the change of variables $r \to (\rho^2 + R_n^2)^{1/2}$ which implies $dr = \frac{\rho}{(\rho^2+R_n^2)^{1/2}}d\rho$, we have
\begin{align*}
\mean{S_{0,n}} &=  {s_{d-1}c_{\mathrm{g}} n}e^{-R_n^2/2}\int_{0}^\infty (\rho^2+R_n^2)^{(d-2)/2}\rho e^{-\rho^2/2}d\rho \\
 &= s_{d-1}c_{\mathrm{g}} n e^{-R_n^2/2}R_n^{d-2}\int_{0}^\infty \param{\param{{\rho}/{R_n}}^2+1}^{(d-2)/2}\rho e^{-\rho^2/2}d\rho.
\end{align*}
The integrand is bounded, and applying dominated convergence we have

\[
\lim_{\ninf}\param{{n e^{-R_n^2/2} R_n^{d-2}}}^{-1}\mean{S_{0,n}}  = s_{d-1}c_{\mathrm{g}}.
\]
Taking $R_n = \Ro^\eps \triangleq \sqrt{2 \log n + \param{{d-2}+\epsilon} \log \log n}$, we have
\[
e^{-R_n^2/2} = n^{-1}(\log n)^{-(d-2+\epsilon)/2}
\]
and so
\[
\lim_{\ninf} {n e^{-R_n^2/2} R_n^{d-2}} = 0
\]
which implies that
\[
    \mean{S_{0,n}} \to 0.
\]

Finally, for every $0 \le k \le d-1$,
\[
\beta_{k,n} \le S_{0,n}.
\]
Therefore,
\[
\lim_{\ninf}\mean{\beta_{k,n}} = 0,
\]
completing the proof.
\end{proof}

\bibliographystyle{plain}