 Research
 Open Access
 Published:
A tractable closed form approximation of the ergodic rate in Poisson cellular networks
EURASIP Journal on Wireless Communications and Networking volume 2019, Article number: 187 (2019)
Abstract
The employment of stochastic geometry for the analysis and design of ultra dense networks (UDNs) has provided significant insights into network densification. In addition to the characterization of the network performance and behavior, these tools can also be exploited toward solving complex optimization problems that could maximize the capacity benefits arising in UDNs. However, this is preconditioned on the existence of tractable closed form expressions for the considered figures of merit. In this course, the present paper introduces an accurate approximation for the moment generating function (MGF) of the aggregate othercell interference created by base stations whose positions follow a Poisson point process of given spatial density. Given the pivotal role of the MGF of the aggregate interference in stochastic geometry and the tractability of the derived MGF, the latter can be employed to substantially simplify ensuing stochastic geometry analyses. Subsequently, the present paper employs the introduced MGF to provide closed form expressions for the downlink ergodic capacity for the interference limited case, and validates the accuracy of these expressions by the use of extensive Monte Carlo simulations. The derived expressions depend on the density of users and base stations, setting out a densification road map for network operators and designers of significant value.
Methods/Experimental
The methods used in the present paper are based on the mathematical tools of random spatial processes and stochastic geometry. The analytical framework developed for the performance evaluation of ultra dense networks is validated by extensive Monte Carlo simulations.
Introduction
The advent of multimedia interactive services and the surge in the number of interconnected devices has imposed the investigation of new approaches able to enhance wireless capacity in 5G networks. In this course, three prime axes of network flexibility have been leveraged, namely the employment of wider spectrum, the enhancement of spectral efficiency, and the employment of smaller cell sizes and, thus, of smaller transmit distances. In retrospect over the evolution of wireless networks, the efficient spatial reuse of the spectrum, through the reduction of the intersite distances, has provided, out of these three axes of flexibility, the most substantial capacity gains by a large margin [1]. Hence, the densification of networks arises as the most prominent candidate for achieving the envisaged capacity increase in the 5G era as well.
In the direction of densifying their networks, network operators employ systemlevel simulations and network measurements. However, over the last decade, the seminal work of Baccelli et al. [2] gave rise to stochastic geometry as a tractable tool for the largescale analysis and design of wireless networks. Indicatively, the formulation of tractable mathematical expressions for the expectation of the aggregate interference in wireless networks [3]—which was not analytically formulated hitherto—paved the way for the theoretical analysis of the performance of wireless networks. Thenceforth, a multitude of research works exploited these tools to provide significant insights into network densification. These insights are essential for understanding the innate features of dense networks and can be employed by network operators as densification road maps and guidelines for the use of system level simulators and of auxiliary network planning tools.
The insights provided by such theoretical analyses brought about significant changes in the understanding of wireless networks. In particular, for singleslope pathloss models and for networks comprising significantly more users than base stations (BSs), it has been demonstrated that the user signal quality is independent of the BS density [4]. Moreover, the probability of coverage (which constitutes the complementary cumulative distribution function (CCDF) of the signal to interference plus noise ratio (SINR)) is independent of the BS density and of the number of tiers [5]. Similar analysis for the uplink (UL) has demonstrated that the UL signal to interference ratio (SIR) is also invariant of the BS density [6]. As a result, it has been demonstrated that the network capacity increases linearly with the density of BSs and with the number of tiers [3].
These conclusions, however, which indeed hold for sparse wireless networks (e.g., tier of macro cells (MCs)), do not hold for extremely high BS densities. The reason for that is that after a BS densification threshold, the intersite distances become so small that the proximity of the neighboring BSs allows them to create line of sight (LOS) interference to the intended user. As a result, after this densification threshold, the probability of coverage is diminished precipitately due to the presence of LOS interference [7, 8].
As opposed to this behavior of networks comprising much more users than BSs, the system performance is not bounded by the aforementioned threshold in the case of networks with more BSs than users; which is, in fact, the case of the envisaged UDNs [9]. The reason for that is that the excess BSs that do not serve any user can be switched off, thus, reducing the system energy consumption and interference. In this setup, the probability of coverage will still decrease when LOS interference first appears in the network, however, as the BS density increases past this point, the probability of coverage will increase again with the density of BSs. That is, since for a fixed number of users, the excess BSs will remain idle, and therefore will not create any interference, regardless of the density of BSs [10].
Based on this comprehensive analysis, it has become evident that two key factors need to be taken into account by network operators in order to tap the capacity potential of UDNs. Firstly, the detrimental effect of LOS interference and, secondly, the beneficial effect of idle (i.e., nontransmitting) BSs. The incorporation of these two effects in the design of UDNs could engender extraordinary capacity gains. In particular, the investigation of BS coordination schemes to counteract LOS interference could, indeed, allow for a linear capacity increase with the density of BSs. Besides, the coordination of only LOS interferers, able to communicate and coordinate directly with the intended user, would not require intricate coordination schemes. Additionally, leveraging on the beneficial effect of the idle mode, users could be clustered dynamically under a single BS, not necessarily the one providing the best service to each user. Thus, BSs that were acting as sources of principal interference to the network could be switched off. Given the high density of BSs, the connection to a neighboring BS after the best serving BS has been switched off would entail a minimal pathloss increase that would be outweighed by the interference mitigation gain achieved by strategically switching off BSs.
The development of mathematical frameworks for complex optimization problems, like the ones mentioned above poses a great challenge. The reason for that is that the majority of the stochastic geometry approaches in the literature, including the theoretical analyses presented above, involve intractable integrations. Even though such integrals can be computed numerically, allowing the analysis of the network behavior, they cannot be employed for the investigation of complex optimization problems. In these cases, it is imperative that the considered objective functions, which evaluate the system performance, involve tractable closed form expressions. In this course, it is essential to exploit the available stochastic geometry tools to develop tractable and accurate approximations in addition to the available exact but cumbersome expressions.
In this direction, the authors have already provided approximate but accurate closed form expressions for the DL rate of coordinating BSs [11], providing a figure of merit for the investigation of coordination schemes among BSs. However, the derivation of closed form expressions for the DL ergodic rate of noncoordinating BSs had remained an open issue in the literature hitherto. Addressing this gap in the literature, the present paper provides approximate but extremely accurate closed form expressions for the DL ergodic rate for the interference limited case. The expressions account for networks comprising more users than BSs, as well as networks comprising more BSs than users. In the latter case, the DL rate is associated by a closed form expression to the density of BSs and of users allowing for the investigation of complex optimization problems like the one mentioned above. Last but not least, in the direction of obtaining the aforementioned expressions, an extremely accurate and simple approximation of the MGF of the aggregate othercell interference in the DL is provided. In the existing literature, the complexity of the exact MGF imposed inherent limitations on the extension of the stochastic geometry analysis to complex optimization problems. However, the present paper lifts these inherent limitations by introducing a simple, albeit extremely accurate expression for the MGF of the other cell interference.
To elaborate, the main contributions of the present paper can be summarized as follows:

Derivation of a tractable and extremely accurate approximation for the MGF of the aggregate othercell interference. The aggregate interference is created by BSs whose spatial distribution follows a homogeneous Poisson point process (PPP).

Extension of the MGF to account for two different scenarios. A realistic scenario where the fast fading characterizing the channels of the interferers is not known to the UE, and a scenario of academic interest, where the fast fading of the interferers can be taken into account.

Derivation of a tractable and accurate approximation for the coverage probability in the interference limited and fully loaded case (i.e., for λ_{UE}≫λ).

Derivation of a tractable and accurate approximation for the DL ergodic rate, in closed form, in the interference limited and fully loaded case.

Derivation of a tractable and accurate approximation for the coverage probability in the interference limited but not fully loaded case (i.e., the case of the envisaged UDNs). As opposed to the fully loaded case, the coverage probability depends on the density of users λ_{UE} and the density of BSs λ.

Derivation of a tractable and accurate approximation for the DL ergodic rate, in closed form, in the interference limited but not fully loaded case. The closed form of the expressions and their dependence on the density of users λ_{UE} and BSs λ, constitutes an invaluable figure of merit for network planning and network densification.

The accuracy of the proposed approximations is validated through extensive Monte Carlo simulations, and the practicality of the derived expressions for future applications is outlined.
The remainder of the paper is organized as follows. Section 3 presents the considered network architecture. Section 4 presents the approximation for the MGF of the aggregate othercell interference in the DL. Section 5 introduces for the first time in the literature closed form expressions for the ergodic rate in the DL for the interference limited case. Moreover, extensive simulation results are presented corroborating the accuracy of the derived expressions. Finally, Section 6 concludes the paper and presents perspectives.
The wireless cellular network
A wireless cellular system is considered, comprising a set of BSs, denoted by BS_{i}, whose positions \(x_{i} \in \mathbb {R}^{2}\) follow a spatial distribution given by a homogeneous PPP Ψ of density λ(BSs/m^{2}). Moreover, the positions of the overlaid user equipment (UE) follow also a spatial distribution given by a homogeneous PPP Φ of density λ_{UE}(UE/m^{2}). The reference UE, denoted by U E_{0}, is located at the origin and is served by its closest BS, denoted by BS_{0}. The UE adjoined at the origin can be singled out and the location of the other UE follows the reduced Palm distribution of Φ, which is the same as the original distribution Φ (as stated by Slivnyaks theorem [12]). Hence, adjoining U E_{0} at the origin does not change the distribution of Φ.
For the sake of simplicity in the notation, it is assumed that all UE and BSs are equipped with one antenna and all BSs transmit at the same power level. Intra cell users are assumed to be sharing orthogonal resources, as is typically the case in the literature [6], whereas all BSs use the same frequency band. If not explicitly mentioned otherwise (as will be done in the following sections), it is assumed that the network comprises significantly more users than BSs (i.e., λ_{UE}≫λ). As a result, every BS is active and transmitting, acting as an interferer in the DL. The othercell interference in the DL is mathematically defined as the interference coming from all BSs residing at a distance x_{i} from the origin greater than the distance x_{0} of BS_{0} from the origin, where · denotes the l^{2}norm. Note that the origin is where the reference UE U E_{0} is located. Last but not least, a singleslope unbounded pathloss model is assumed in the analysis.
The considered scenario is depicted in Fig. 1. U E_{0} is marked by the magenta diamond, BS_{0} by the magenta circle, and the intra cell users are depicted in red. The interfering BSs are depicted in black and reside in distances (x_{i}>x_{0}), where the circle of radius x_{0} around U E_{0} is depicted in green. Every BS residing in the region outside this circle is acting as an interferer in the DL.
MGF of the aggregate othercell interference
Having presented the considered network scenario and having defined the set of interfering BSs by means of their distance to the origin, we can proceed with the mathematical formulation of the aggregate othercell interference and the MGF of the latter.
Derivation of the MGF
The aggregate othercell interference in the DL is mathematically formulated as follows:
where P_{tx} denotes the transmit power of the BSs, β the pathloss exponent, ‖x‖ the distance of the interferer to the origin, and κ the path loss at a reference distance of 1 m.
Although a term of fast fading, denoted by h^{(x)}^{2}, is present at the propagation, it has not been introduced in the previous expression (1). The reason for that is that taking into account the fast fading of the interferers during the computation of the UE capacity would implicitly mean that the UE has perfect knowledge of the channel of all interferers. However, since this is not the case in practice, averaging over this fading would provide an upper bound for the capacity. As opposed to that, the omission of the fast fading h^{(x)}^{2} provides a lower bound for the capacity [13]. Since a worstcase scenario analysis is more sensible than an overoptimistic calculation of the achievable rate, the fast fading of the interferers is not introduced. However, if accounting for the fast fading is of interest, this can be done without increasing significantly the complexity of the presented analysis as will be demonstrated in Section 4.3.
The set of interfering BSs in (1) has been defined by means of their distances to the origin, i.e., ‖x_{i}‖>‖x_{0}‖. Equivalently, by defining the path loss as follows:
the set of interfering BSs can be defined by means of their path loss to the origin (i.e., L^{(x)}>L^{(0)}, where \({{L^{\left (0 \right)}}}={{L^{\left (x_{0} \right)}}}\) for brevity in the notation). Having defined mathematically the set of interfering BSs, the MGF of g can be obtained by employing the probability generating fuctional (PGFL) theorem according to which [3]:
where (3b) is obtained by computing the double integral of (3a) in polar coordinates for 0≤θ≤2π and for r=‖x‖, and (3c) is obtained by changing the variable of the integration to the path loss according to (2). Hence, employing (3c) the MGF of g is given by [14]:
where the second argument of M_{g}(s;L^{(0)}) denotes the dependence of the MGF on the random variable L^{(0)} and 1(·) is the indicator function. (a) holds by employing (3c), and (b) is attained by using the result of [15] according to which:
with _{1}F_{1}(a,b,z) being the Kummer confluent hypergeometric function given by:
where (·)_{k} denotes the Pochhammer function given by:
The definitions of (6) and (7) demonstrate the intractability of (4b). Since the derivation of the MGF of the aggregate interference is one of the fundamental tools of stochastic geometry, the intractability of (4b) propagates to every analysis employing the MGF, hindering the derivation of closed form figures of merit.
MGF approximation
In order to overcome the limitations imposed by (4b), the present paper introduces a simple, albeit extremely accurate, approximation of the MGF by introducing an alternative calculation of (4a). In this course, the MGF is derived as follows:
where (b) holds by employing the Taylor expansion of the exponential term, (c) holds by a simple calculation of the integral, and (d) and (e) are obtained from the definition of the generalized incomplete gamma function Γ(·,·,·).
Having defined (8e) it can be noted that the term within the integral of (8e) eventually tends to 0. When this happens, namely when \( \exp \left (\frac {s P_{\text {tx}}}{{L^{\left (0 \right)}}} \right) \approx 0\), the integral converges to a constant value. Hence, (8e) can be approximated by a piecewise function, involving a constant value when \( \exp \left (\frac {s P_{\text {tx}}}{{L^{\left (0 \right)}}} \right) \approx 0\) and a varying function when \( \exp \left (\frac {s P_{\text {tx}}}{{L^{\left (0 \right)}}} \right) \neq 0\).
In particular, by employing the constant value to which the integral converges when \( \exp \left (\frac {s P_{\text {tx}}}{{L^{\left (0 \right)}}} \right) \approx 0\), it holds that:
and when \( \exp \left (\frac {s P_{\text {tx}}}{{L^{\left (0 \right)}}} \right) \neq 0\), (8e) can be approximated by the Taylor expansion around 0 as follows:
Thus, by combining (9) and (10), (8e) can be approximated as follows:
and by employing only the first two terms of the Taylor expansion, (11) can be approximated by:
where c is a constant indicating the point after which the integral of (8e) converges to a constant value. c is the point of intersection of the two functions of (12) and can be obtained by solving the following equation:
Although (13) cannot be solved analytically, for the limited range of the pathloss exponent β (i.e., β∈[2,5]), it can be computed numerically and the value of c for any value of β∈[2,5] can be approximated by:
where log(·) in (14) and henceforth denotes the natural logarithm.
The employment of the Taylor expansion and of the exact value of the generalized incomplete function Γ(·,·,·) at the extreme cases guarantees the tightness of the approximation far from the intersection point c. However, the tightness of the proposed approximation still needs to be verified close to c. In this course, relevant figures will be provided in the following sections validating the tightness of the approximation close to c against the exact result of (4b). Given the dependence of c on β, the provided figures will demonstrate the tightness of the proposed approximation in the whole range of β∈[2,5], which is also the range of interest in wireless networks. It should also be noted, that the value of λ does not have any impact on the tightness of the approximation, since the term of (8e) that has been approximated by (9) and (10) does not involve λ. Therefore, the proposed approximation of (12) is not just far more tractable and simple than (4b), but is also tight over the whole range of values that are of interest in wireless networks.
Fast fading of interferers
The aforementioned analysis can be easily extended to account also for the fast fading of the interferers, if the latter is of interest. In this course, an independently marked point process can be employed, that is, a point process where a random variable, known as mark and denoted by M_{x}, is randomly assigned to each random point of the point process x [3]. The marks are mutually independent and the conditional distribution of mark \(M_{x} \in \mathbb {R}^{l}\) of a point x∈Ψ depends only on the location of x. For an independently marked homogeneous PPP with density λ on \( \mathbb {R}^{2}\) and marks with distribution F_{x}(dM) on \( \mathbb {R}^{l}\), the Laplace transform of a function f(x,M_{x}) is given by [16]:
Hence, if the fast fading of the interferers needs be introduced in (1), (15) can be applied directly to (4), with \(M_{x}=\left  {{h^{\left (x \right)}}} \right  \in \mathbb {R}\), in order to compute the expectation with respect to the path loss and to the fading of the interferers. That is, by setting:
and l=1 and by employing the distribution F_{H}(h) for the respective type of fading the MGF of (4a) is revised as follows:
The employment of (15) allowed for moving the expectation over the fading within the exponential term of (17), thus, simplifying the analysis to a great extent. As a result, the expectation over the fading of the interferers can also be moved within the exponential term of (12). Given the tractability of (12), the introduction of an additional integral within the exponential term has a minor effect on the complexity of the derived expressions and the analysis can be easily extended accordingly.
However, as already mentioned, the present analysis does not account for the fading of the interferers focusing on a realistic scenario where the UE does not have perfect knowledge of the channel of all interferers.
Ergodic capacity in the DL
Having defined a simple approximation of the MGF, the latter can be employed to provide closed form expressions for the DL ergodic rate for the interference limited case. In this course, the analysis will commence by employing the MGF of (4b) and (12) to derive the coverage probability for the interference limited case. The latter can be derived by both expressions, i.e., (4b) and (12), in spite of the intractability of the former, allowing for the comparison of the two results and, thus demonstrating the accuracy of the introduced approximation. Subsequently, capitalizing on the accuracy of the introduced approximation, the DL ergodic rate will be derived in closed form by employing the introduced approximate expression for the coverage probability.
The analysis which will initially account only for networks comprising much more users than BSs will then be extended to networks comprising more BSs than users.
Probability of coverage
The probability of coverage (i.e., the probability of the SINR exceeding the value γ) is defined as follows:
where \(\sigma _{N}^{2}\) denotes the noise power and h^{(0)} the fast fading of the intended user. As opposed to the fast fading of the interferers, the fast fading of the intended user is estimated and known in practice.
Assuming Rayleigh fading, then the random variable h^{(0)}^{2} follows an exponential distribution with unit mean and the P_{cov} is given by:
where (a) is obtained by the CCDF of the exponential distribution, (b) is obtained based on the definition of the MGF in (4), and (c) from computing the expectation with respect to the path loss L^{(0)} to the serving BS BS_{0}.
The probability density function (PDF) of the distance between a reference user and its closest BS for a PPP is given in [17]. By employing this PDF and the definition of (2), the PDF of the path loss between a reference user and its closest BS \(f_{L^{\left (0 \right)}}(y)\) is given by:
In the interference limited case, the exponential term of (20c) is equal to 1 and, by employing (21), the coverage probability can be calculated as follows:
where (a) is obtained by employing (4b) and (b) by employing the approximate MGF of (12).
The results of (22) verify the theoretical results presented in Section 2, since for a singleslope pathloss model and for more users than BSs, the coverage probability of (22) does not depend on the density of BSs λ, but only on the path loss exponent β and the SIR value γ (c is also a function of β given by (14)). The accuracy of (22b) and, implicitly, the accuracy of the MGF of (12), is demonstrated in Fig. 2 where the approximate coverage probability of (22b) is compared against the exact result of (22a). Apart from extremely accurate, the approximate coverage probability of (22b) is also significantly more tractable than the exact result of (22a).
Ergodic rate
Having defined the tractable and accurate approximation of (22b) for the probability of coverage, this can be employed to compute the DL ergodic rate. In particular, the probability of coverage given by (22b) constitutes the CCDF of the SIR (i.e., for SIR=w,P_{cov}=1−F_{W}(w)). Hence, the derived tractable expression of the CCDF of the SIR allows for the computation of the DL ergodic rate by averaging over the SIR as follows:
where (a) is obtained by integrating by parts and (b) is obtained by employing (22b). The employment of (22a) in (23a) would not allow the analytical computation of the above integral. However, (23b) can be computed in closed form and the rate is given by:
where the first three terms of (24) are obtained by the calculation of the first term of (23b) and the last term of (24) is obtained by the calculation of the last term of (23b). _{2}F_{1}(a,b,c,z) denotes the Gaussian hypergeometric function, \(\alpha =\sqrt {\left (\frac {2\beta 2}{\beta 2}\right)^{2}+2\beta 2}\) and c is given by (14). Again, as in the case of (22), (24) depends only on the pathloss exponent β.
The employment of (24) allows the computation of the rate for the interference limited case in closed form. The tight performance of (24) for the calculation of the ergodic rate is demonstrated in Fig. 3, where (24) is compared against the results obtained by Monte Carlo simulations. Even though the rate does not depend on the density of the BSs, the density employed for generating the simulation of Fig. 3 was λ=1.27e−06.
Several research works have focused on deriving expressions for the DL ergodic rate, since the latter constitutes the most sensible figure of merit for evaluating the performance of UDNs. Indicatively, in [18] the authors have provided expressions for the DL ergodic rate in heterogeneous cellular networks for all different types of fading. However, in all of these works, including the latter, the calculation of the DL ergodic rate involved at least one integration that had to be computed numerically. This imposed inherent limitations to the applicability of those expressions to complex optimization problems. In order to overcome this problem, the authors of [19] provided for the first time in the literature expressions for the DL ergodic rate in the interference limited case, for a fully loaded scenario (i.e., for λ_{UE}≫λ) that did not involve any numerical integration. In this direction, the authors of [19] provided lookup tables and employed the MeijerG function to provide a tight approximation of the DL ergodic rate in the fully loaded case. In particular, the approximate ergodic rate of [19] is given by:
where s^{∗} is the solution to the equation \(s^{*^{\delta }}\Gamma (\delta,s^{*})=0\), given by lookup tables that are computed a priori for all relevant values of β. Furthermore, \(\delta =\frac {2}{\beta }, D_{\delta }=\frac {s^{*}}{\log (1\text {sinc}(\delta))}, \operatorname {E_{n}}(x) = \int _{1}^{\infty } t^{n}e^{x t} dt \) is the exponential integral and
is the MeijerG function. Improving the result of (25), the present paper introduces for the first time a closed form approximation for the ergodic rate in (24), without the need for an a priori computed lookup table. Furthermore, the result of (24) is significantly more tractable than the result of (25), and Fig. 4 demonstrates the superior performance of (24) over (25) with respect to the tightness of the approximation. The accuracy and, more importantly, the tractability of the above expressions allow for the extension of the analysis to even more complex scenarios, as will be demonstrated in the following sections.
Ergodic rate over density of users and BSs
As already mentioned, the previous analysis corresponds to a scenario where the density of users is much greater than the density of BSs (i.e., λ_{user}≫λ) and, therefore, every BS is in transmission mode. However, since in the envisaged UDNs the number of BSs is expected to be higher than the number of UE [9], the analysis needs to be extended accordingly, taking into account the non transmitting mode of the excess BSs that do not have any UE in their coverage. In this course, the proposed tractable MGF is revised to account for the probability of the excess BSs to remain idle. This probability is defined by the density of UE λ_{UE} and the density of BSs λ. Thus, following a similar approach as before, we derive closed form expressions for the DL ergodic capacity (i.e., peak and divided among intracell users) which depend on the density of UE λ_{UE} and the density of BSs λ.
The probability that a randomly chosen BS does not have any UE in its Voronoi cell and, therefore, goes into idle mode is given by [20]:
and the probability of a BS being in transmission mode and, thus, acting as an interferer in the DL, is denoted by:
Since only a subset of the BSs create interference in the DL, the density of the BSs λ must be thinned out by the probability P_{active}. By introducing the thinned out density λP_{active} into (4) and (12), the MGF is revised as follows:
and, following the approach of (22) for the MGF of (29), the coverage probability is now given by:
The probability of coverage defined by the exact result of (30a) and by the approximation of (30b) is plotted in Fig. 5 for a pathloss exponent β=4 and for density ratios λ_{UE}/λ=0.1,0.4,0.7,1, demonstrating, once again, the accuracy of the derived approximation.
Following the same approach as in (23) for the results of (30b), the DL ergodic peak rate can be computed in closed form as follows:
For pathloss exponent values of β=3,4,5, the closed form expressions for the peak rate are given in Table 1. The expressions (34)–(36) provide a closed form representation of the peak DL ergodic rate over the probability P_{active} and, implicitly through (28), over the density of BSs λ and of UE λ_{UE}.
The probability that a randomly chosen UE is assigned a resource block at a given time and is served by its nearest BS is given by [20]:
Employing this probability, a more sensible figure of merit than the peak rate can be derived. This figure of merit is the actual DL ergodic rate of the reference UE, i.e., the rate of the reference UE after dividing the available resources and, thus, the peak rate among all intracell UE. The latter is given by:
The accuracy of expressions (33)–(36) is demonstrated in Figs. 6 and 7 where the peak and the actual DL ergodic rates are plotted over the ratio of the densities λ_{UE}/λ and compared against Monte Carlo simulations.
In the Monte Carlo simulations, the BSs are deployed following the PPP Ψ of density λ and the users are deployed following the PPP Φ of density λ_{UE}, whereas the reference UE resides at the origin. BSs with no UE in their coverage (i.e., within their Voronoi cell) do not create any interference. To compute the actual rate, the number of users residing within the Voronoi cell of BS_{0} are counted in each realization and the peak rate is divided among these users and the reference UE. The number of simulated BSs is fixed and the simulation area expands or contracts as the BS density λ changes, in order to accommodate the predefined number of BSs, while the density of users λ_{UE} remains fixed.
The closed form expressions derived in (33)–(36) provide, among others, a substantial computational gain when compared to the computational time of the respective Monte Carlo simulations. Especially, since the simulation of a wireless network, which comprises both users and BSs of different spatial distributions, and the calculation of their relative distances, is computationally expensive. In order to demonstrate the gain arising by the employment of the closed form expressions of (33)–(36), the computational time of the derived closed form expressions is tabulated in Table 2. Since this time is not an absolute metric, but depends on the hardware employed, Table 2 presents the computational time of the expressions as a percentage of the computational time required for simulating the respective wireless networks of Fig. 7 using the same hardware.
Table 2 demonstrates that the computational gain arising by the derived expressions is immense. The time required for the analytical computation of the rate is practically zero compared to the time required for simulating the respective scenario. At this point it should also be noted, that the variations in the values of the computational time arise due to the variations in the time required for the respective simulations. That is since the simulation area has to expand and contract as λ changes in order to accommodate a fixed number of BSs. The variations with respect to β emerge due to the different complexity of the analytical expressions of (34)–(36), with respect to β.
Apart from the computational gain that has been demonstrated by Table 2, the closed form expressions of (33)–(36) and the respective figures, i.e., Figs. 6 and 7, provide a deep understanding of the behavior of UDNs as the density of users λ_{UE} and BSs λ changes. First of all, expressions (34)–(36) and Fig. 6 verify one of the fundamental findings of stochastic geometry, already explained in Section 2, that in the fully loaded case (i.e., for λ_{UE}≫λ) the ergodic rate remains invariant while the BS density changes. However, Fig. 7 demonstrates that even if the peak rate remains invariant, the rate of the users that have to share this peak rate tends to zero as the density ratio (i.e., the expected number of users per typical cell) increases. This behavior demonstrates why the envisaged UDNs are expected to comprise more BSs than users, highlighting the importance of the nonfully loaded case. Given the importance of the nonfully loaded case, expressions (34)–(36) and Fig. 6 allow for the first time to quantify the threshold between the nonfully and the fully loaded case. In particular, in all three expressions for the different values of β, the network exhibits the behavior of a fully loaded network for λ_{UE}/λ>4, way earlier than implied by the notation λ_{UE}≫λ.
Another interesting finding can be derived from the behavior of the network in the nonfully loaded case, when the BSs that do not comprise any user in their coverage do not create any interference. In this setup, Fig. 6 demonstrates that the densification of the network can provide substantial capacity gains, since the achieved DL ergodic rate in this range is significantly higher than the rate achieved at the fully loaded case. Thus, Fig. 6 demonstrates that just by switching off the BSs that do not comprise users in their coverage, substantial capacity gains can be engendered for the network. This motivates the use of even more intricate schemes where BSs can be switched off strategically to mitigate interference. This optimization problem is only one of the problems that the derived expressions can be applied to, with additional applications being presented in the following section.
Moreover, expressions (34)–(36) and the respective figures, verify the results of the expressions (24) and (25), and of Figs. 3 and 4, according to which the DL ergodic rate increases monotonically with the pathloss exponent β. In particular, since the distance from the UE to the serving BS is always smaller than the distance to the interfering BSs, i.e., since ‖x_{i}‖>‖x_{0}‖, each term of the sum \(\sum \limits _{\substack {x \in \Psi, \\ x \neq x_{0}}}\left (\frac {\left \lVert x_{0} \right \rVert }{\left \lVert x_{i} \right \rVert }\right)^{\beta }\) decreases with β. Consequently, the SIR and the DL ergodic rate increase monotonically with β.
Applications of the derived expressions
Figures 6 and 7 corroborate the accuracy of the derived expressions, while providing a deeper understanding of the network performance as a function of the user and BS densities. Hence, these expressions, which for the first time associate the DL ergodic rate with the densities of the UE and of the BSs is a closed form, pave the way for the investigation of complex optimization problems, toward improving UDN operation and offered QoS.
In particular, given the maximum density of UE λ_{UE} during the operation of the network and the minimum rate requirement per user (imposed by the QoS constraints), Fig. 7 and (33) can be employed by network operators to define the minimum BS density λ that guarantees this rate. Hence, (33) can be employed as a rule of thumb for the lower limit of the densification of BSs that guarantees the QoS objectives, implicitly quantifying the minimum capital expenditure required by network operators.
In addition to that, given the aforementioned minimum rate requirement per user and a network of BS density λ, comprising UE whose density λ_{UE} varies during the operation of the network, Fig. 7 and (33) can be employed to dynamically define the probability of transmission P_{active} that achieves the predefined rate requirement as λ_{UE} changes. This probability can indicate the density of transmitting BSs λP_{active}, i.e., the density of the BSs comprising at least one UE in their coverage. The latter density can be input into optimization modules, to be used as a starting point for the search aiming at pinpointing the optimum set of BSs to be switched off, additionally to those that do not comprise any UE in their coverage. Given the high density of the BSs, strategically switching off the best serving BSs of a UE for the latter to be served by a neighboring BS has a only minimal impact on the path loss, while it can provide substantial capacity gains through the mitigation of the interference.
Conclusions
The present paper has demonstrated how stochastic geometry tools can be exploited to derive not just exact but cumbersome expressions, but also simple, albeit extremely accurate closed form expressions that allow for the investigation of complex optimization problems. The resolution of these problems is essential in order to reap the capacity benefits of UDNs. In this direction, the present paper presented an accurate and tractable approximation for the MGF of the aggregate othercell interference. Given the pivotal role of the MGF in stochastic geometry analyses, the derived approximation can be employed by a multitude of applications to simplify the analysis and facilitate the derivation of closed form expressions.
Building upon this result, the present paper has focused on the interference limited case, providing very tractable expressions for the coverage probability, as well as closed form expressions for the DL ergodic capacity (i.e., peak and actual) which depend on the density of UE λ_{UE} and the density of BSs λ. The derived expressions and their dependence on the user and BS densities set out a densification road map for network operators and designers of significant practical and commercial value. Moreover, such expressions can be used for the resolution of complex optimization problems in real time, improving the network performance and operation.
Availability of data and materials
Not applicable. No datasets were imported or generated while developing the presented material.
Abbreviations
 BS:

Base station
 CCDF:

Complementary cumulative distribution function
 DL:

Downlink
 LOS:

Line of sight
 MC:

Macro cell
 MGF:

Moment generating function
 PDF:

Probability density function
 PGFL:

Probability generating fuctional
 PPP:

Poisson point process
 QoS:

Quality of service
 SINR:

Signal to interference plus noise ratio
 SIR:

Signal to interference ratio
 UDN:

Ultra dense network
 UE:

User equipment
 UL:

Uplink
References
 1
V. Chandrasekhar, J. G. Andrews, A. Gatherer, Femtocell networks: a survey. IEEE Commun. Mag.46(9), 59–67 (2008).
 2
F. Baccelli, B. Blaszczyszyn, P. Muhlethaler, An Aloha protocol for multihop mobile wireless networks. IEEE Trans. Inf. Theory. 52:, 421–436 (2006).
 3
M. Haenggi, Stochastic Geometry for Wireless Networks (Cambridge University Press, Cambridge, 2013).
 4
J. G. Andrews, F. Baccelli, R. K. Ganti, A tractable approach to coverage and rate in cellular networks. IEEE Trans. Commun.59(11), 3122–3134 (2011).
 5
H. S. Dhillon, R. K. Ganti, F. Baccelli, J. G. Andrews, Modeling and analysis of ktier downlink heterogeneous cellular networks. IEEE J. Sel. Areas Commun. 30(3), 550–560 (2012).
 6
S. Singh, X. Zhang, J. G. Andrews, Joint rate and SINR coverage analysis for decoupled uplinkdownlink biased cell associations in HetNets. IEEE Trans. Wirel. Commun.14(10), 5360–5373 (2015).
 7
M. D. Renzo, W. Lu, P. Guan, The intensity matching approach: A tractable stochastic geometry approximation to systemlevel analysis of cellular networks. IEEE Trans. Wirel. Commun.15(9), 5963–5983 (2016).
 8
M. Ding, P. Wang, D. LopezPerez, G. Mao, Z. Lin, Performance impact of los and nlos transmissions in dense cellular networks. IEEE Trans. Wirel. Commun.15(3), 2365–2380 (2016).
 9
J. G. Andrews, Seven ways that hetnets are a cellular paradigm shift. IEEE Commun. Mag.51(3), 136–144 (2013).
 10
M. Ding, D. LopezPerez, G. Mao, Z. Lin, Performance impact of idle mode capability on dense small cell networks. IEEE Trans. Veh. Technol.66(11), 10446–10460 (2017).
 11
A. I. Aravanis, O. Munoz, A. PascualIserte, M. Di Renzo, in 2019 IEEE 89th Vehicular Technology Conference (VTC2019Spring), Kuala Lumpur, Malaysia. On the Coordination of Base Stations in Ultra Dense Cellular Networks, (2019), pp. 10–6. https://doi.org/10.1109/VTCSpring.2019.8746531, http://ieeexplore.ieee.org/stamp/stamp.jsp?tp=&arnumber=8746531&isnumber=8746285.
 12
S. N. Chiu, D. Stoyan, W. S. Kendall, J. Mecke, Stochastic Geometry and Its Applications (John Wiley and Sons, West Sussex, 2013).
 13
G. Geordie, Devicetodevice communication and wearable networks: Harnessing spatial proximity. PhD thesis, Universitat Pompeu Fabra (2017).
 14
T. U. LamThanh, New methods for the analysis and optimization of cellular networks by using stochastic geometry. PhD thesis, Universite ParisSaclay (2018).
 15
M. D. Renzo, W. Lu., Systemlevel analysis/optimization of cellular networks with simultaneous wireless information and power transfer: Stochastic geometry modeling. IEEE Trans. Veh. Technol.66:, 2251–2275 (2017).
 16
J. Lee, F. Baccelli, in IEEE INFOCOM 2018  IEEE Conference on Computer Communications. On the effect of shadowing correlation on wireless network performance, (2018), pp. 1601–1609.
 17
S. Sadr, R. S. Adve, Partiallydistributed resource allocation in smallcell networks. IEEE Tran. Wirel. Commun.13(12), 6851–6862 (2014).
 18
M. D. Renzo, A. Guidotti, G. E. Corazza, Average rate of downlink heterogeneous cellular networks over generalized fading channels: A stochastic geometry approach. IEEE Trans. Commun.61(7), 3050–3071 (2013).
 19
G. Geordie, R. K. Mungara, A. Lozano, M. Haenggi, Ergodic spectral efficiency in MIMO cellular networks. IEEE Trans. Wirel. Commun.16(5), 2835–2849 (2017). https://doi.org/10.1109/TWC.2017.2668414, http://ieeexplore.ieee.org/stamp/stamp.jsp?tp=&arnumber=7880697&isnumber=7920448.
 20
S. M. Yu, S. L. Kim, in 2013 11th International Symposium and Workshops on Modeling and Optimization in Mobile, Ad Hoc and Wireless Networks (WiOpt). Downlink capacity and base station density in cellular networks, (2013), pp. 119–124. http://ieeexplore.ieee.org.recursos.biblioteca.upc.edu/stamp/stamp.jsp?tp=&arnumber=6576422&isnumber=6576384.
Funding
The work presented in the present paper has been carried out within the framework of the project ETN5Gwireless (this project has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie SkłodowskaCurie grant agreement No. 641985). Moreover, the work has been partially funded through the grant 2017 SGR 578 (funded by the Catalan Government—Secretaria d’Universitats i Recerca, Departament d’Empresa i Coneixement, Generalitat de Catalunya, AGAUR) and the project TEC201677148C21R (AEI/FEDER, UE): 5G&BRUNNERUPC (funded by the Agencia Estatal de Investigacion, AEI, and Fondo Europeo de Desarrollo Regional, FEDER).
Author information
Affiliations
Contributions
All authors have contributed to this research work. Moreover, all authors have read and approved the final manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
About this article
Cite this article
Aravanis, A.I., Tu Lam, T., Muñoz, O. et al. A tractable closed form approximation of the ergodic rate in Poisson cellular networks. J Wireless Com Network 2019, 187 (2019). https://doi.org/10.1186/s1363801914999
Received:
Accepted:
Published:
Keywords
 Ultra dense networks
 Poisson point process
 Moment generating function
 Interference
 Ergodic capacity
 Rate
 Closed form
 Base station density
 User density