number of BSs is a Poisson distributed stochastic variable with mean value $\lambda A$. The BSs are uniformly and independently distributed over the area. Each of them is equipped with an array of $M$ antennas and serves $K$ single-antenna UEs, which are selected at random from a (potentially) very large set of UEs in the cell since we treat $K$ as an optimization variable. We assume that each UE connects to its closest BS such that the coverage area of a BS is its Poisson-Voronoi cell; see Fig. \[figureVoronoi\] for an illustration. The UEs are assumed to be uniformly distributed in the Poisson-Voronoi cell of their serving BS. Note that the geographic locations of UEs and BSs are correlated under this model, such that small cells serve more UEs per $\textrm{km}^2$ than larger cells. We interpret this as a judicious network deployment where the BSs are matched to a heterogenous distribution of UEs.[^3]

The translation invariance of PPPs allows us to perform statistical performance analysis for a *typical UE* (located for example in the origin), which is statistically representative for any other UE in the network [@Baccelli2008a]. Assume that this typical UE has the arbitrary index $k$ and is connected 