# Coverage and Capacity Scaling Laws in Downlink Ultra-Dense Cellular Networks

## Abstract

Driven by new types of wireless devices and the proliferation of bandwidth-intensive applications, data traffic and the corresponding network load are increasing dramatically. Network densification has been recognized as a promising and efficient way to provide higher network capacity and enhanced coverage. Most prior work on performance analysis of ultra-dense networks (UDNs) has focused on random spatial deployment with idealized singular path loss models and Rayleigh fading. In this paper, we consider a more precise and general model, which incorporates multi-slope path loss and general fading distributions. We derive the tail behavior and scaling laws for the coverage probability and the capacity considering strongest base station association in a Poisson field network. Our analytical results identify the regimes in which the signal-to-interference-plus-noise ratio (SINR) either asymptotically grows, saturates, or decreases with increasing network density. We establish general results on when UDNs lead to worse or even zero SINR coverage and capacity, and we provide crisp insights on the fundamental limits of wireless network densification.

## 1Introduction

With exponential increase in data traffic driven by a new generation of wireless devices, data is expected to overwhelm cellular network capacity in the near future. Heterogeneous cellular network (HetNet) deployment is a promising and effective way to provide high cellular network capacity by overlaying conventional macrocell cellular architecture with heterogeneous architectural features, such as small cellular access points (picocells and femtocells), low-power fixed relays, and distributed antennas. Ultra-dense networks (UDNs) are expected to achieve higher data rates and enhanced coverage by exploiting spatial frequency reuse, while retaining at the same time seamless connectivity and mobility. Inspired by the attractive features and potential advantages of UDNs, their development and deployment has been gaining momentum in both wireless industry and research community during the last few years. It has also attracted the attention of standardization bodies, e.g. 3GPP LTE-Advanced.

Existing cellular network analyses have focused on stochastic geometry based models, in which base stations (BSs) are located according to a Poisson point process (PPP). Most prior results have considered a standard singular power-law path loss model and Rayleigh fading, as a means to provide tractable analysis of coverage probability and other key performance metrics in downlink cellular networks. For instance, recent results using the aforementioned models with closest BS [1] and strongest BS association [?] show that the coverage probability does not depend on the network density when thermal noise is negligible. However, the singular single-slope path loss model used therein leads to unrealistic results in certain scenarios and fails to accurately capture the dependence of path loss exponent on the link distance. In [20], the authors study the impact of dual slope path loss on the performance of downlink UDNs under closest BS association and show that both coverage and capacity performance strongly depends on the network density. More precisely, it is shown that the network coverage in terms of signal-to-interference-plus-noise ratio (SINR) is maximized at some finite density and there exists a phase transition in the near-field path loss exponent with ultra densification (i.e. network density goes to infinity). In [6], the authors consider strongest BS association with bounded path loss and lognormal fading and show that the coverage attains a maximum point before going down when the network densifies. Based on system-level simulations, [13] shows that there is a fundamental limit of 1 cell/user in UDNs, although such deployments are neither cost nor energy efficient. Despite shedding light on the performance limits of UDNs and the optimal network density, previous work has mainly focused on the effect of path loss singularity [12] or boundedness [20]. The effect of fast fading has been only investigated in some particular contexts, e.g. [11], [16]. Furthermore, emerging utilization of advanced communication and signal processing techniques is expected to enhance the channel gain, which in some cases may have a regularly varying tail [19] or a diffuse power component [9].

In this work, we analyze the SINR, coverage probability and capacity of downlink UDNs under multi-slope path loss and general fading distributions, considering strongest BS association in PPP networks. In particular, using general fading distributions, we provide the ability to capture the effect of either fast fading, shadowing, or composite fast fading-shadowing. We study the tail behavior of the received signal power, as well as the SINR, coverage and capacity scaling in the ultra-dense regime. Our results provide useful insights on the fundamental limits of network densification with the following main conclusions:

Under the Poisson field assumption, the most affecting component of the path loss is its near-field exponent . Bounded path loss (obtained for ) is just a special case of with being the network dimension.

The effect of fading on the performance scaling is as significant as path loss, and regularly varying fading distributions have the same effect as path loss singularity.

In more conventional cases when the near-field path loss exponent is greater than the free-space dimension (i.e. ), or when fading is heavy tailed (i.e. , ), the coverage and capacity saturate at a limiting bound when the network density increases.

In more realistic cases with (in particular ) and fading being less heavy tailed or even truncated, both coverage and capacity exhibit an ‘inverse U’ behavior w.r.t. network density, i.e. both are maximized at finite density then vanish when the network further densifies.

Finally, all standard fading models, such as Rayleigh, lognormal, Gamma, and their composite forms, belong to the same class of fading distributions, which leads to coverage and capacity maximization at a certain finite network density.

## 2System Model

### 2.1Network Model

Consider a typical downlink user located at the origin and that the network is composed of cell sites located at positions . For convenience, cell sites are referred to as nodes, whereas the *typical user* is simply referred to as *user*. Unless otherwise stated, are assumed to be random variables independently distributed on the *network domain* according to a homogeneous Poisson point process (PPP) of intensity , denoted by . In prior work, the entire -dimensional Euclidean space , where is usually assumed as network domain. Since network domain is in reality limited, and far-away nodes are less relevant to the typical user due to path loss, we assume that the distance from the user to any node is upper bounded by some constant , which can be arbitrarily large. Each node transmits with some power that is independent to the others but is not necessarily constant.

### 2.2Propagation Model

Let represent the path loss function. The receive power is related to the transmit power by with being the transmitter-receiver distance. Physics laws require that . However, in the literature, has been usually assumed to admit a power-law model, i.e. where is the path loss exponent satisfying . This far-field propagation model has been widely used mainly due to its tractability. However, for short ranges, especially when , this model is no longer relevant and becomes singular at the origin. In the context of network densification where the inter-site distance becomes smaller, the above singular model may be unsuitable. Moreover, the dependence of the path loss exponent on the distance in emerging millimeter wave (mmWave) communications [7] advocates the use of a more generic path loss function. In this work, the path loss is modeled as follows

where is the indicator function, is a given constant characterizing the number of path loss slopes, are constants satisfying

denotes the path loss exponent satisfying

and are constants to maintain continuity of , i.e.

for . For notational simplicity, we also use the following notation

The above model (cf. (Equation 1)) captures that the path loss exponent varies with the distance while remaining unchanged within a certain range. In principle, free-space propagation in has path loss exponent equal to 2 (i.e. ), whereas in realistic scenarios, path loss models often include antenna imperfections and empirical models usually result in the general condition ( ?) for far-field propagation. Condition ( ?) models the physical property that the path loss increases faster as the distance increases. Notice, however, that this condition is not important in the subsequent analytical development. Finally, condition ( ?) is related to the near field (i.e. it is applied to the distance range ).

The multi-slope path loss function, as defined above, has the following widely used special cases:

, : , which is the standard unbounded path loss,

, : , which is the bounded path loss recommended by the 3GPP standard, in which is referred to as the minimum coupling loss.

Due to the particular importance of path loss boundedness so as to have a realistic model, we have the following definitions:

It is clear that the path loss function (Equation 1) is bounded if and only if (iff) , and is physical iff and .

Besides path loss, shadowing and fast fading are also sources of wireless link variations, which are commonly referred to as *fading* in the sequel. Let be a variable containing transmit power, fading, and any gains or attenuation other than path loss from -th node to the user. Given node location , the variables are assumed not identical to zero and independently distributed according to some distribution . To this end, the signal power, say , that the user receives from -th node is expressed as where is the Euclidean distance.

### 2.3Performance Metrics

The quality of the signal received from -th node is expressed in terms of its SINR as

where is the aggregate interference with respect to -th node’s signal, and thermal noise at the user’s receive antenna is assumed Gaussian with average power . In addition, we define to be the total interference. The main metrics used in this paper are the coverage probability and the average rate that the user experiences from its *serving cell*. Let denote the SINR that the user receives from its serving cell. The *SINR coverage probability*, denoted by , is defined as the probability that is larger than a given target , and the *capacity*, denoted by , is defined as the Shannon rate (in nats/s/Hz) assuming Gaussian codebooks, i.e.

### 2.4User Association

The above performance metrics are defined with respect to the user’s serving cell, which in turn depends on the underlying user association scheme. *Strongest cell association* is used in this work, i.e. the user is connected to the cell that provides the best/strongest signal quality (in practice, time averaging of the signal is usually performed to avoid frequent handover due to fast fading). In a longer version of this work, we also investigate the case of nearest cell association and show the effect of user association on performance scaling. Under strongest cell association, the SINR of the user is given by and can be expressed as [15]

### 2.5Notation

Quantities whose dependence on the network density is important, are denoted as , e.g. , , and . We also denote by , , and , the distance, associated fading, and received power from a random node, respectively. Let be the distribution of , and . In addition, for real functions and , we say if for , if , and if . We also use notation , , to denote the convergence in distribution, convergence in probability, and almost sure (a.s) convergence, respectively. Finally, for two random variables and defined on the same probability space, we say that is *statistically* greater than , denoted by , if . Similar definition is for , , and .

Note that if is a regularly varying function with index at , it can be represented as as for some slowly varying function .

## 3Tail Behavior of Received Signal Power

The network performance mainly depends on the received SINR and hence is a function of , which in turn depends on and (cf. (Equation 5)). The behavior of the maximum and the sum is totally determined by that of the received power . Therefore, we first study the signal power , in particular its tail behavior using tools from extreme value theory.

Under the assumption that BSs are distributed according to a homogeneous PPP, the nodes of a realization of are uniformly distributed. Thus, given and the assumption that , the distance distribution to a random node of , say , is for , and for . Then, taking , the result follows.

See Appendix Section 6.

Based on Proposition ?, we derive the following main result:

See Appendix Section 7.

Theorem ? shows that the tail behavior of the wireless link depends not only on whether path loss is bounded or not, but also on the tail behavior of the fading. More precisely, a key implication of Theorem ? is that path loss and fading have interchangeable effects on the tail behavior of the wireless link. This can also be shown using Breiman’s Theorem [5] and results from large deviation and product distributions. In particular, when the fading distribution is a regularly varying function, the wireless link is also regularly varying regardless of the path loss function’s boundedness. For lighter-tailed fading, the regular variation property of the wireless link is solely imposed by the path loss singularity.

More importantly, Theorem ? is a general result and covers all tail behaviors for the fading: case (1) covers the heaviest tails (i.e. with , e.g. Pareto distributions), as well as the moderately heavy tails (i.e. the class , e.g. exponential, normal, lognormal, Gamma distributions). Case (2) covers all remaining tails (e.g. truncated distributions). Therefore, for any realistic statistical model and distribution of fast fading and shadowing, Theorem ? enables us to characterize the tail behavior of the wireless link, which is essential to understand the behavior of the interference, the maximum received power, and their asymptotic relationship. In wireless communications, the signal distribution often involves lognormal or Gamma shadowing and Rayleigh fading, which all belong to the class , and the path loss is bounded, thus . As a result, it can be shown that in most relevant cases in wireless UDNs, belongs to the maximum domain of attraction of a Gumbel distribution [17].

Finally, the above result generalizes prior results: [11] showed that the interference is tail-equivalent with the fading if path loss is bounded and if , which is not applicable for , . [16] showed that under lognormal fading, is regularly varying for unbounded path loss, and behaves like a lognormal tail for bounded path loss.

The following result, which is a direct consequence of Theorem ?, can be provided in order to better understand the signal power scaling and the interplay between path loss function boundedness and fading.

## 4Scaling Laws

In this section, we provide the main results of this paper, namely the SINR, coverage probability and capacity scaling, when the network density is asymptotically large. Using results from Section 3, we investigate and provide insights on the performance limits of network densification.

First, we start by showing that in sparse networks, the signal quality improves for increased node density.

As , we have almost surely. Thus, for ,

where note that is intuitively evident, but a formal proof can be easily obtained using for example [2].

### 4.1SINR Scaling

We provide here the scaling of the received SINR under strongest BS association in the asymptotically large node density regime.

For ,

as since is finite and is not identical to 0.

If , due to [14]. Thus,

If with , as where has a non-degenerate distribution [3]. As a result,

where is a non-degenerate distribution.

If with or with , we have . Hence, due to [18]. Moreover,

Thus, leads to as .

Let us have a closer look at Theorem ? and on its implications in the interplay between multi-slope path loss and fading. According to Corollary ?, due to the fact that . Recall that contains the transmit power and all potential channel powers (including fading). Therefore, implies that the channel powers take large values with non negligible probability. As a result, dominates and compensates the path loss, resulting in maximum power that grows at the same rate as the aggregate interference. This provides a theoretical justification to the fact that network densification always enhances the signal quality .

When with , implies that the SINR distribution converges to a non-degenerate distribution. Moreover, from Corollary ?, this convergence is due to either large near-field exponent or heavy-tailed fading. In that case, for any SINR target , the coverage probability flattens out starting from some large value of node density. This means that further increasing the network density by deploying more BSs does not improve the network performance. In Figure ? we show the convergence of to a steady distribution for two cases: or with , where ’*Composite*’ represents the case of composite Rayleigh-lognormal fading, which belongs to the rapidly varying class , and stands for Pareto fading distribution of shape and some scale , i.e.

In practically relevant network settings, the path loss is bounded and fading is moderately heavy tailed (i.e. as in the case of lognormal shadowing and Rayleigh fading) or even truncated (i.e. where ). As a result, based on Theorem ?, we have that , hence . In other words, the SINR is proven to be asymptotically decreasing with the infrastructure density. This means that there is a fundamental limit on network densification and the network should not operate in the ultra-dense regime. In other words, deploying too many BSs would decrease the network performance due to the fact that the increased signal power cannot compensate for the faster growing aggregate interference. Figure ? confirms that with with (i.e. either or ), the tail of vanishes and converges to zero when increases (ultra-dense regime). This convergence to zero of the SINR in the ultra-dense regime further emphasizes the importance of local spatial scheduling among BSs since near-field interferers generate much stronger interference than far-field ones.

### 4.2Coverage and Capacity Scaling

To provide a complete characterization of the network performance in the ultra-dense regime, we further study the scaling of coverage probability and capacity. By the definition of the aforementioned two metrics and Theorem ?, we obtain the following result.

Case (1) directly follows Theorem ?. For case (2), by Theorem ?, we have , where has a non-degenerate distribution. Then, for constant and , as . Similarly, as .

For case (3), given the conditions and Theorem ?, we have . Hence, and by Lemma ?, s.t. . For , we first note that . Thus, by Lebesgue’s dominated convergence theorem and by the fact that . For , such that due to Lemma ?. Thus, , since is decreasing w.r.t. . In consequence, .

The following observations can be made based on Corollary ?: First, both the coverage probability (cf. Figure ?) and the capacity (cf. Figure 1) saturate (ceiling effect) when the fading is regularly varying with index in or the near-field path loss exponent is large (i.e. ). Second, when the fading is less heavy-tailed ( with ) and when the near-field path loss exponent is smaller (i.e. ), both performance metrics are maximized at a finite network density, then decrease and go to zero in the ultra-dense regime (cf. ’*inverse U*’ curves in Figure 2 and Figure 3). This suggests that there is an optimal point of network density to aim for. Furthermore, fading and path loss have equivalent and somehow interchangeable effects on the network performance. Additionally, the most affecting element of the path loss is its near-field exponent ; bounded path loss function (i.e. ) is just a special case of the class . The impact of the near-field exponent was also observed in [20] for nearest cell association.

Table 1 summarizes the behavior of network performance according to Corollary ?.

Scaling Regime |
lighter tail | |||

, |
saturation | |||

## 5Conclusion

This paper analyzes the scaling regimes of the coverage probability and average rate under multi-slope path loss model and general fading distributions, considering strongest cell association in a Poisson field network. We show that when the fading distribution is regularly varying with index in or the near-field path loss exponent is larger than the network dimension , both coverage and capacity saturate to a constant value in the ultra-dense regime. When the fading distribution is less heavy-tailed and when , both coverage and capacity are maximized at some finite density. Finally, our results show that path loss and fading have interchangeable effects on the tail behavior of SINR, coverage probability and capacity, extending previous results using dual slope path loss model and Rayleigh fading.

## 6Proof of Proposition

Using given in (Equation 1), we have

where .

(i) For : integration by parts with and for yields

where , which reduces to ( ?) after change of variable and applying the condition due to (Equation 3). On the other hand, applying again (Equation 3),

Hence, substituting (Equation 8) back to (Equation 7) yields ( ?) with .

(ii) For , we have:

Thus,

in which due to (Equation 3). Hence, ( ?) with .

## 7Proof of Theorem

Due to space limitations, we only give the proof for the case of for . The proof for the cases and are similar and would be provided in a longer version of this paper.

can be represented as for some . Then, by the monotone density theorem [4], the density function of can be written as as .

For

, as ,

If : by Karamata’s theorem [4], for :

If : similarly to the above case, we can easily obtain the same result using Karamata’s theorem. Note that if for some , we can also show that for some constant . Thus, for ,

*For *: If , then exists and as . Thus, as . Otherwise, i.e. if , then we have

where is due to [10], and is due to the representation . Thus, as ,

By substituting (Equation 9) and (Equation 10) into the expressions of given by Proposition ?, we have as ,

where , and where for the last case with and , we further have as . As a result, is regularly varying with index if , and with index if . Hence the proof for with .

### References

- J. Andrews, F. Baccelli, and R. Ganti, “A tractable approach to coverage and rate in cellular networks,”
*IEEE Trans. Commun.*, vol. 59, no. 11, pp. 3122–3134, Nov. 2011. - F. Baccelli and B. Blaszczyszyn,
*l@Anglais StochasticGeometry andWirelessNetworks,VolumeI -Theory*.1em plus 0.5em minus 0.4em NoW Publishers Inc., 2009. - N. Bingham and J. Teugels, “Conditions implying domains of attraction,” in
*Proc. 6th Conf. Probability Theory*, 1981, pp. 23–34. - N. H. Bingham, C. M. Goldie, and J. L. Teugels,
*Regular Variation*. 1em plus 0.5em minus 0.4emCambridge University Press, 1989. - L. Breiman, “On some limit theorems similar to the arc-sin law,”
*Teor. Veroyatnost. i Primenen.*, vol. 10, no. 2, pp. 351–359, 1965. - C. S. Chen, V. M. Nguyen, and L. Thomas, “On small cell network deployment: A comparative study of random and grid topologies,” in
*Proc. Veh. Technol. Conf. (VTC Fall) 2012*.1em plus 0.5em minus 0.4emIEEE, 2012, pp. 1–5. - S. Deng, M. K. Samimi, and T. S. Rappaport, “28GHz and 73GHz millimeter-wave indoor propagation measurements and path loss models,” in
*Proc. 2015 ICC Workshop*.1em plus 0.5em minus 0.4emIEEE, 2015, pp. 1244–1250. - H. Dhillon, R. Ganti, F. Baccelli, and J. Andrews, “Modeling and analysis of -tier downlink heterogeneous cellular networks,”
*IEEE J. Sel. Areas Commun.*, vol. 30, no. 3, pp. 550–560, Apr. 2012. - G. Durgin, T. Rappaport, and D. A. de Wolf, “New analytical models and probability density functions for fading in wireless communications,”
*IEEE Trans. Commun.*, vol. 50, no. 6, pp. 1005–1015, Jun. 2002. - P. Embrechts, C. Klüppelberg, and T. Mikosch,
*Modelling Extremal Events*.1em plus 0.5em minus 0.4emSpringer, 1997. - M. Haenggi and R. K. Ganti,
*Interference in large wireless networks*.1em plus 0.5em minus 0.4emNow Publishers Inc, 2009. - H. Inaltekin, M. Chiang, H. V. Poor, and S. B. Wicker, “On unbounded path-loss models: effects of singularity on wireless network performance,”
*IEEE JSAC*, vol. 27, no. 7, pp. 1078–1092, 2009. - D. Lopez-Perez, M. Ding, H. Claussen, and A. H. Jafari, “Towards 1 Gbps/UE in cellular systems: Understanding ultra-dense small cell deployments,”
*arXiv:1503.03912*, Mar. 2015. - R. A. Maller and S. I. Resnick, “Limiting behaviour of sums and the term of maximum modulus,”
*Proceedings of the London Mathematical Society*, vol. 3, no. 3, pp. 385–422, 1984. - V. M. Nguyen and F. Baccelli, “A stochastic geometry model for the best signal quality in a wireless network,” in
*Proc. 8th Int. Symp. Model. Optim. Mobile, Ad Hoc Wirel. Netw.*1em plus 0.5em minus 0.4emIEEE, May 2010, pp. 465–471. - V. M. Nguyen, “Wireless link modelling and mobility management for cellular networks,” Ph.D. dissertation, Telecom ParisTech, 2011.
- V. M. Nguyen, F. Baccelli, L. Thomas, and C. S. Chen, “Best signal quality in cellular networks: asymptotic properties and applications to mobility management in small cell networks,”
*EURASIP J. Wirel. Commun. Netw.*, vol. 2010, 2010. - G. L. O’Brien, “ l@English A limit theorem for sample maxima and heavy branches inGalton-Watson trees,”
*l@English J. Appl. Probab.*, vol. 17(2), pp. 539–545, 1980. - A. Rajan, C. Tepedelenlioglu, and R. Zeng, “A unified fading model using infinitely divisible distributions,”
*arXiv:1508.04804*, 2015. - X. Zhang and J. Andrews, “Downlink cellular network analysis with multi-slope path loss models,”
*IEEE Trans. Commun.*, vol. 63, no. 5, pp. 1881–1894, May 2015.