Modeling and Analysis of Uplink NonOrthogonal Multiple Access (NOMA) in LargeScale Cellular Networks Using Poisson Cluster Processes
Abstract
Nonorthogonal multiple access (NOMA) serves multiple users by superposing their distinct message signals. The desired message signal is decoded at the receiver by applying successive interference cancellation (SIC). Using the theory of Poisson cluster process (PCP), this paper provides a framework to analyze multicell uplink NOMA systems. Specifically, we characterize the rate coverage probability of a NOMA user who is at rank (in terms of the distance from its serving BS) among all users in a cell and the mean rate coverage probability of all users in a cell. Since the signaltointerferenceplusnoise ratio (SINR) of th user relies on efficient SIC, we consider three scenarios, i.e., perfect SIC (in which the signals of interferers who are stronger than th user are decoded successfully), imperfect SIC (in which the signals of of interferers who are stronger than th user may or may not be decoded successfully), and imperfect worst case SIC (in which the decoding of the signal of th user is always unsuccessful whenever the decoding of its relative stronger users is unsuccessful). The worst case SIC assumption provides remarkable simplifications in the mathematical analysis and is found to be highly accurate for scenarios of practical interest. To analyze the rate coverage expressions, we first characterize the Laplace transforms of the intracluster interferences in closedform considering both perfect and imperfect SIC scenarios. In the sequel, we characterize the distribution of the distance of a user at rank which is shown to be the generalized Beta distribution of first kind and the conditional distribution of the distance of the intracluster interferers which is different for both perfect and imperfect SIC scenarios. The Laplace transform of the intercluster interference is then characterized by exploiting distance distributions from geometric probability. The derived expressions are customized to capture the performance of a user at rank in an equivalent orthogonal multiple access (OMA) system. Finally, numerical results are presented to validate the derived expressions. It is shown that the average rate coverage of a NOMA cluster outperforms its counterpart OMA cluster with higher number of users per cell and higher target rate requirements. A comparison of Poisson Point Process (PPP)based and PCPbased modeling is conducted which shows that the PPPbased modeling provides optimistic results for the NOMA systems.
I Introduction
Until very recently, the stateoftheart wireless communications systems have been utilizing a variety of orthogonal multiple access (OMA) technologies, in which the resources are allocated orthogonally to multiple users. These techniques include frequencydivision multiple access (FDMA), timedivision multiple access (TDMA), codedivision multiple access (CDMA), and orthogonal frequencydivision multiple access (OFDMA). Since OMA maintains the orthogonality among users in a cell, the intracell interference (i.e., interuser interference within a cell) does not exist. As a result, the information signals of users can be retrieved at a low complexity. Nonetheless, the number of served users is limited by the number of orthogonal resources.
Conversely, NOMA serves multiple users simultaneously using the same spectrum resources (i.e., radio channels), however, at the cost of increased intracell interferences. To mitigate the intracell interferences, NOMA exploits Successive Interference Cancellation (SIC) at the receivers [1]. NOMA supports low transmission latency and signaling cost compared to conventional OMA where each user is obliged to send a channel scheduling request to its serving base station (BS). With these attractive features, NOMA can be a potential access technology for 5G networks. Nevertheless, the conclusions about the performance of NOMA are largely unknown in multicell network scenarios. For instance, in uplink NOMA, a large number of transmitting users in the neighboring cochannel BSs can result in high interference at the BS of interest. Consequently, the uplink multicell interference in NOMA is directly proportional to the number of transmitting users per neighboring cochannel BS and is more severe compared to the intercell interference in OMA.
Ia Background Work
The concept of NOMA was initially proposed in [1] for downlink transmissions. Various practical aspects, such as multiuser scheduling, impact of error propagation in SIC, overall system overhead, and user mobility were discussed. System level simulations were conducted in [2] to highlight the benefits of twouser NOMA over OMA, in terms of overall system throughput as well as individual user’s throughput.
The approximate expressions for the ergodic sumrate and outage probability of a user in a given downlink NOMA cluster were derived in [3]. Later, in [4], the throughput gains of the twouser cooperative NOMA (in which the strong channel user relays the information of weak channel user) were investigated. The idea of cooperative NOMA was then applied to wirelesspowered systems in [5]. Based on users’ distances, grouping of users was performed first. Then, three user selection schemes were investigated, i.e., (i) pairing of the nearest users from each group, (ii) pairing of the nearest user from one group and the farthest user from another group, and (iii) arbitrary user pairing. The direct link was used to transfer energy from the BS. A cooperative data link was established for the lower channel gain user via the higher channel gain users. Closedform approximate expressions for the outage probability and throughput of a twouser NOMA cluster were derived. In [6], the user pairing was investigated considering fixed NOMA (FNOMA) and cognitive radio inspired (CRNOMA). In FNOMA, any two users could make a NOMA pair based on their channel gains. While in CRNOMA, a weak channel user opportunistically gets paired with the strong channel user provided that the interference caused by the strong user will not harm the rate requirement of the weak channel user. It was observed that CRNOMA pairs the strongest user with the second strongest user, whereas FNOMA pairs a strongest user with the weakest user in the system.
A general concept of uplink NOMA was discussed in [7]. An uplink power backoff policy was proposed to distinguish users in a NOMA cluster with nearly similar signal strengths (given that traditional uplink power control is applied). Closedform analysis was performed for ergodic sumrate and outage probability of a twouser NOMA cluster. Further, the problem of user scheduling, subcarrier allocation, and power control in uplink NOMA has been investigated by various researchers in [8, 9] with perfect SIC at the BS. A game theoretic algorithm for uplink power control has been designed in [10] considering a twocell NOMA system where intercell interference is assumed to be Gaussian distributed.
IB Motivations and Contributions
To date, most of the research investigations consider the throughput analysis of NOMA for singlecell downlink systems with perfect SIC at the receivers. The derived expressions generally leverage on high signaltonoise ratio (SNR) and asymptotic assumptions as well as the application of GaussianChebyshev quadrature (GCQ) technique which approximates all integrals into finite sums. Unfortunately, the conclusions about the performance gains of NOMA (compared to OMA) in singlecell/singlecluster scenarios cannot be applied directly to multicell/multicluster scenarios. The reason is the intercell^{1}^{1}1The term inter/intracell and inter/intracluster interference will be used interchangeably throughout the paper. interferences in NOMA can be quite severe as well as distinct from OMA, especially in the uplink scenarios [11]. Note that, in uplink NOMA, the intercell interference incurred at the BS of interest is directly proportional to the number of transmitting users in neighboring cochannel BSs. This is different from an equivalent uplink OMA system where only one user transmits at a time per neighboring cochannel BS.
Further, compared to downlink NOMA, we note that the performance analysis of uplink NOMA is particularly more challenging due to the mathematical structure of the intracell interferences. In the uplink NOMA, the BS receives transmissions from all users simultaneously. As such, the intracell interference to a user is a function of the channel statistics of other users within the cell. On the other hand, in downlink NOMA, the intracell interference to a user is a function of its own channel statistics [3, 4, 6]. To this end, the contributions of this paper are outlined as follows:

Using the theory of order statistics and Poisson Cluster Process^{2}^{2}2NeymanScott PCP, such as Modified Thomas Cluster process and Matern Cluster Process, are recently exploited in a set of research studies for performance evaluation of adhoc clustered networks [12], D2D systems [13, 14], and downlink multitier cellular networks [15, 16]. (PCP), we develop a framework to analyze the rate coverage probability of a user who is at rank (in terms of the distance from its serving BS) among all users in a cell and the mean rate coverage probability of all users in a cell, considering a multicell uplink NOMA system. Unlike typical stochastic geometry frameworks where the user locations are uniform over the 2D plane and are independent of the BS locations [17, 18], we exploit Matern Cluster Process (MCP) to accurately model the proximity of multiple users around a BS^{3}^{3}3The Poisson Point Process (PPP)based modeling of BSs and user devices focuses on the link between the serving BS and a typical user. Since the typical user can be located anywhere in the cell, results are averaged over all spatial positions inside the cell. Such an approach provides higher analytical flexibility. However, in practice, users are more likely clustered around a BS and are distinct due to their channel conditions. In this regard, PCP have been shown empirically to be a more accurate cellular network modeling technique [19] that allows locationspecific performance modeling of users..

NOMA systems rely on efficient SIC and the interference of a user at rank needs to be adapted according to the level of SIC. We consider three SIC scenarios that include perfect SIC (in which the signals of interferers who are stronger than th user are decoded successfully), imperfect SIC (in which the signals of interferers who are stronger than th user may or may not be decoded successfully), and imperfect worst case SIC (in which the decoding of the signal of th user is always unsuccessful whenever the decoding of its relative stronger users is unsuccessful). The worst case SIC assumption provides remarkable simplifications in the mathematical analysis and is found to be highly accurate for scenarios of practical interest.

The interference power at the BS of a given cell/cluster is composed of intracluster and intercluster interferences. We derive the Laplace transform of the intracluster interference in closedform considering various SIC scenarios. In the sequel, we also derive the distribution for the distance of a user at rank which is shown to be the generalized Beta distribution of first kind and the conditional distribution of the distance of the intracluster interferers which is different for both perfect and imperfect SIC scenarios. The Laplace transform of the intercluster interference is then characterized using distance distributions from geometric probability. A lesscomplex bound is then exploited to model the Laplace transform of the intercluster interference.

The derived rate coverage expressions are customized to evaluate the performance of a user at rank in an equivalent OMA system in closedform. Numerical results are presented to validate the derived expressions. Our results indicate that the performance benefit of OMA diminishes quickly with the increase in number of users per cluster and higher user rate requirements. A comparative performance analysis of PPPbased and PCPbased modeling is conducted using simulations. It is shown that PPPbased modeling generally provides optimistic results for the NOMA systems due to the homogeneous distribution of users regardless of the BS locations which reduces the impact of intracluster interference.
The rest of the paper is structured as follows. Section II discusses the working principle of uplink NOMA along with cellular network model, channel model, and interference model. In Section III, we describe the fundamental differences between conventional SIC and SIC for NOMA. Considering perfect SIC, imperfect SIC, and imperfect worst case SIC, we model the interferences and define the desired performance metrics. In Section IV, we derive relevant distance distributions required for the characterization of the Laplace transforms of the interferences. In Section V, we derive the rate coverage expressions for both NOMA and OMA systems. Finally, Section VI discusses numerical and simulation results followed by the concluding remarks in Section VII.
Ii System Model and Assumptions
Iia Spatial Cellular Network Model
We consider a singletier cellular network composed of macrocell base stations (BSs) surrounded by user devices. The locations of transmitting user devices are modeled as a stationary and isotropic PCP.
Definition 1 (Poisson Cluster Process (PCP) [12]).
A PCP results from applying homogeneous independent clustering to a stationary Poisson Point Process (PPP). In particular, the parent points form a stationary PPP of density . The offsprings generated for a given parent are a family of independent and identically distributed (i.i.d.) finite point sets with distribution independent of the parent process. The complete PCP can thus be represented as Note that the parent points themselves will not be included in the PCP. The parents and offsprings are referred to as the cluster centers and the cluster members, respectively.
Definition 2 (NeymanScott PCP).
If the number of points per cluster follow a Poisson distribution with mean intensity , such a PCP is referred to as a NeymanScott Process.
MCP is a special case of the NeymanScott PCP where the cluster centers (BSs) are modeled by a parent homogeneous PPP in the Euclidean plane with density . Each parent point forms the center of a cluster around which daughter points (user devices) are uniformly distributed in a circle of radius as shown in Fig. 1. For ease of exposition, we consider number of users per cluster. Each daughter point located at with respect to its cluster center has a density function as:
(1) 
where is the distance of any arbitrary daughter point (user device) relative to its cluster center (serving BS) and its density function can be given as follows:
(2) 
The resulting MCP is a stationary and isotropic point process of density and can be defined as Without loss of generality, we perform analysis for a user at rank located in a randomly chosen cluster which is referred as representative cluster located at throughout the paper.
IiB Working Principle of Uplink NOMA
In an uplink NOMA cluster, each user transmits its individual signal with a transmit power such that the received signal at the BS can be defined as . Note that, to apply SIC and decode signals at the BS, it is crucial to maintain the distinctness of various signals superposed within . Since the channels of different users are different in the uplink, each message signal experiences distinct channel gain^{4}^{4}4The conventional uplink transmit power control (typically intended to equalize the received signal powers of users) may remove the channel distinctness and thus may not be feasible for uplink NOMA transmissions.. As a result, the received signal power corresponding to the strongest channel user is likely the strongest at the BS. Therefore, this signal is decoded first at the BS and experiences interference from all users in the cluster with relatively weaker channels. That is, the transmission of the highest channel gain user experiences interference from all users within its cluster, whereas the transmission of the lowest channel gain user receives zero interference from the users in its cluster.
IiC Channel and Interference Model
NOMA allows multiple users in a cluster to share the same resources by superposing their distinct message signals. All users are served on the same channel and time slot. All users/BSs are equipped with a single antenna. Within the representative cluster, any arbitrary user device located at with respect to its serving BS located at transmits its individual signal with power such that the superposed NOMA signal at the representative BS can be defined as follows:
(3) 
where is the pathloss exponent and is the exponential random variable which models Rayleigh fading associated with the channel between the node located at and the BS located at . All fading coefficients are i.i.d. and the additive noise is complex Gaussian distributed with mean zero and variance per dimension.
Consequently, the uplink NOMA transmission of a user located at within the representative cluster is vulnerable to two kind of interferences, i.e.,

Intracluster interference: is the interference received at the representative BS from all user devices located within the representative cluster (except the user located at ). However, after performing SIC, some of these interferences can be removed (details of SIC will follow in the next section).

Intercluster interference: is the interference received at the representative BS from all user devices located outside the representative cluster.
Provided that the center of the representative cluster is located at , the intracluster interference experienced by the transmission of user located at with respect to its cluster center can be modeled as follows:
(4) 
Similarly, the intercluster interference at the representative BS from the user devices located outside the representative cluster center can be modeled as follows:
(5) 
and . However, since NOMA systems rely on efficient SIC, the intracluster interference model needs to be adapted according to the level of SIC cancellations. This is elaborated in detail in the next section.
Iii Successive Interference Cancellation (SIC) and Performance Metrics
In this section, we will discuss the fundamental differences between conventional SIC and SIC for NOMA. We then describe the signaltointerferenceplusnoise ratio (SINR) modeling with perfect and imperfect SIC in uplink NOMA. Considering three possible SIC scenarios, i.e., perfect SIC, imperfect SIC, and imperfect worst case SIC, we finally define the performance metrics.
SIC is among one of the best known interference cancellation methods as (i) the SIC receiver is architecturally similar to traditional nonSIC receivers in terms of hardware complexity and cost, (ii) it uses the traditional decoder to decode the composite signal at different stages and neither complicated decoders nor multiple antennas are required, and (iii) it can achieve the Shannon capacity for both the broadcast and multiple access networks [20, 21]. Typically, SIC is used to regenerate the interfering signals and subsequently cancel them from the received composite signal to improve the SINR of the desired signal. That is, the SIC receiver first decodes the strongest signal by treating other signals as noise. Then it regenerates the analog signal from the decoded signal and cancels it from the received composite signal. The remaining signal is thus free from the strongest interfering signal. Then, the SIC receiver proceeds to decode, regenerate, and cancel the second strongest interfering signal from the remaining signal and so on, until the desired signal can be decoded.
Iiia SIC and SIC Error Propagation in Uplink NOMA
In uplink NOMA, we apply the same SIC principle at the BS, i.e., the SIC receiver first decodes the strongest signal by treating other signals as noise and so on. However, the difference is that the intracluster interfering signals are also the desired signals; therefore, it is not possible to provide the benefits of SIC (enhance the SINR) unequivocally for all users. That is, within a cluster, the user with strongest signal experiences interference from all users and the user with the weakest signal enjoys zero intracluster interference. Evidently, the intracluster interference statistics will vary for all user devices within the representative NOMA cluster.
Since the decoding of the strongest signal is performed first at the BS, its success/failure has a significant impact on the decoding of other users’ signals. Specifically, depending on the decoding result of the strongest signal, the interference used for the decoding of the second strongest signal differs, which makes the linktosystem mapping difficult. If the strongest signal is decoded correctly, its replica signal can be subtracted successfully from the superposed signal at the BS. Otherwise, the second strongest user will experience the interference from the strongest user as well as other users in the cluster. This phenomenon is referred to as SIC error propagation.
IiiB Modeling of Intracluster Interference with SIC
To model the intracluster interference with SIC, first the BS needs to rank the received powers of various users as such that with . However, note that the impact of pathloss factor is more stable and dominant compared to the instantaneous multipath channel fading effects. Therefore, the order statistics of the distance outweigh the fading effects, which vary on a much shorter time scale. As such, the ranking of users in terms of their distances from the serving BS is generally considered as a reasonable approximate of their respective ranked received signal powers [22]. This approximation provides great flexibility for the analytical purposes. Also, the SIC based on longterm channel states is more practically feasible since it requires less overheads for channel estimation. Note that the exact performance analysis of the user with strongest signal is unwieldy to solve since the distribution of is the ranked distribution of a composite uniform and exponential random variable and the joint distribution of several composite ordered random variables is required.
Approximation: As mentioned above, the impact of pathloss factor is more dominant compared to the channel fading effects. Hence, for tractability reasons, we assume that ordering of the received signal powers can be approximately achieved by ordering the distances of the users as such that with . That is, when the th strongest signal is decoded and subtracted from the composite signal, this means that the remaining interferers are located farther than the th rank user whose distance is from the representative BS.
The intracluster interference and in turn the SINR of th rank user can thus be modeled for perfect and imperfect SIC scenarios, respectively, as shown below.
IiiB1 Perfect SIC
In this case, a given user at rank receives interferences from all users with relatively weaker channel gains (or users with farther distances as per the approximation) and the BS perfectly decodes/cancels the strong interferences. The intracluster interference experienced by any user at th rank in a cluster can thus be modeled, after perfectly canceling strong interferences, as follows:
(6) 
where . Note that the ranking is applied only at the distances as per the approximation. The SINR experienced by any user at rank in the representative cluster can therefore be defined as follows:
(7) 
IiiB2 Imperfect SIC and Detection Probability
The signals from interferers (who are located closer to the BS than the rank user) may or may not be decoded perfectly; therefore, SIC may or may not be performed in a perfect fashion. In such a case, we first define the probabilities for the successful detection of the signals of the ranked users (ranked in terms of their distances), respectively, as follows:
where and , represent the probability of successful and unsuccessful detection of th ranked user’s signal, respectively. It can be seen that the BS attempts to decode the closest user without any interference cancellation. If the decoding is unsuccessful, the interference from this user remains intact. Subsequently, we can generalize the detection probability of a user at th rank as follows:
(8) 
where , is the signal detection threshold, and denotes the set of combinations in which each combination has bits. The successful detection is represented by a binary digit whereas the detection failure is given by . The intracluster interference experienced by a user at rank thus depends on whether the detections for closer users were successful or not. As such, conditioned on a given combination , the SINR experienced by a user at rank can be modeled as:
(9) 
where
Note that, even after successful detection, a given user can still experience rate outage (i.e., the achievable rate may remain below the target rate requirement).
IiiC Performance Metrics
We analyze the performance gains of uplink NOMA considering a system where each user has a target data rate requirement. For this case, the rate coverage probability of a user at rank and the mean rate coverage probability of a cluster are relevant performance metrics. These metrics are defined for different SIC scenarios in the following.
IiiC1 Rate Coverage
Rate coverage probability is the probability that a given user’s achievable rate remains above the target data rate. Mathematically, the rate coverage probability of a user at th rank can be defined as , where is the target data rate requirement of the th ranked user. Now, we define the rate coverage probability of th user in the following specific cases:

Rate Coverage with Perfect SIC: Using the definition of from (7), the rate coverage probability of a user at rank in the representative cluster can be defined as:
(10) where is the desired SINR corresponding to the rate requirement of the user at rank .

Rate Coverage with Imperfect SIC: In this case, we consider the probability of decoding/canceling interferences from closer interferers as less than one. As such, the rate coverage probability of a user at rank needs to consider all possible combinations and thus can be defined using (9) as follows:
(11) 
Rate Coverage with Imperfect SIC  Worst Case: The worstcase model assumes that the decoding of any user at rank is always unsuccessful whenever the decoding of his relative closer users is unsuccessful. Such a worstcase model is simple and allows evaluating the impact of SIC error propagation on the NOMA performance without invoking complicated NOMA specific linktosystem mapping. The worstcase detection probability of a user at rank can therefore be given as:
(12) where can be given using (7). The rate coverage probability of a user at rank can thus be given as:
(13)
As a byproduct, we can evaluate the rate coverage of a user at rank in an equivalent TDMAbased OMA system where users in a cluster are served in orthogonal time slots. The rate of a user at rank in OMA system can be defined as . For fair comparison with NOMA, we need to include the scaling factor of which shows the portion of available resources to a given user in OMA system. As such, if the target rate requirement of the user is , the desired SINR threshold can be defined as . Subsequently, substituting zero intracluster interference, for intercell interfering clusters, and replacing in the rate coverage expressions of NOMA, we can determine the rate coverage of a user at rank in an OMA system.
Further, we can also calculate the average rate of a user at rank . For perfect SIC, imperfect SIC, and worst case SIC scenarios, average rates for a user at rank can be computed, respectively, as follows:
Substituting in place of , we can integrate all coverage probability expressions over to evaluate the average rate of a user at rank numerically. In addition, substituting and , we can characterize the performance of the closest and farthest user in a NOMA cluster, respectively.
IiiC2 Mean Rate Coverage of a Cluster
Although the individual rate coverage of a user in the representative cluster is a useful metric, it does not offer a complete insight related to the cumulative performance of the users in the representative NOMA cluster. As such, we consider analyzing the collective performance of all users by defining the mean rate coverage of all users in the representative NOMA cluster as follows:
(14) 
where , and for perfect SIC, imperfect SIC, and worst case SIC, respectively.
Iv Relevant Distance Distributions for the Characterization of Interference
In this section, we characterize the relevant distributions of the distances between the BS of a representative cluster and different intracell user devices that are ordered according to their distances. These distance distributions are crucial for deriving the Laplace transforms of the intracluster interferences and the rate coverage analysis.
As mentioned in the Approximation, we consider ranking the users in terms of their distances from the serving BS, which is located at the cluster center. Consequently, we characterize the distance distributions of the closest node as well as its respective intracluster interfering nodes from the BS. Note that the distance of any arbitrary intracluster device located at with respect to is i.i.d. and follows the uniform distribution given in (1). Subsequently, the distance of any arbitrary intracluster device from the representative BS follows a sampling distribution given by in (2). Now we order devices within the representative cluster with respect to the cluster center such that . The distance distribution of the user at rank can thus be given as follows.
Lemma 1 (Distribution of the Distance of the User at rank ).
For a Matern cluster process, the distribution of the distance of a user at rank in the representative cluster, i.e., from its serving BS follows a Generalized Beta (GB) distribution of the first kind. The distance distribution can be derived as:
where denotes the Euler Beta function. GB distribution includes Beta distribution, Generalized Gamma distribution, and Pareto distribution as its special cases.
Proof.
See Appendix A. ∎
Provided that the user of interest is at rank , we now characterize the distribution of the distances of its corresponding intracluster interfering devices. Note that the possible interfering nodes for a user at rank can lie at any place (depending on the perfect and imperfect SIC) except the location of the user at rank . As such, the set of intracluster interferers can be partitioned into two subsets, i.e., and where and represent the set of interfering nodes closer and farther to the reference BS, respectively, compared to the user at rank . This setup is demonstrated in Fig. 2. Consequently, the distances from the representative BS to devices in are i.i.d as shown by Afshang et al. in [13, 14, 23, 24]. The same property holds for the users in .
Lemma 2 (Conditional Distribution of the Distances of IntraCluster Interfering Users).
Conditioned on the distance of the user of interest at rank (say ), the distribution of the distance of any device in the set from cluster center (BS) is i.i.d and can therefore be given by:
(15) 
Similarly, the distribution of the distance of any device in the set from cluster center (BS) is conditionally i.i.d. and can be given as:
(16) 
where and are conditionally independent.
Proof.
The proof can be done along the same lines as shown in [13]. For sake of completeness of the paper, we discuss it briefly in Appendix B. ∎
At this point, it is noteworthy that Lemma 2 utilizes the fact that the ordering of users in the set does not have any impact on the cumulative interference generated from all users of this set. Therefore, the interfering devices in the set can be sampled randomly without any specific ordering. Consequently, the distance distributions of the interfering users belonging to set are i.i.d. and can be given as in Lemma 2. The same statements hold for the users in set .
Nonetheless, for imperfect SIC scenarios, we need to include specific interferences from the ranked users in set whose signals go undetected. As such, for each combination , we need to include additional intracluster interferences from the specific users whose corresponding in set . Given that the user of interest is at rank , the distance distribution of a user at specific rank , can be derived by exploiting the conditional ranked distributions from the theory of order statistics as in the following.
Lemma 3 (Conditional Distribution of the Distance of an IntraCluster Interfering User at Rank ).
Given that the user of interest is at rank (say ), the conditional distribution of the distance of th rank user, such that and , is the distribution of the th order statistics of i.i.d. distance variables whose PDF is the truncated PDF of as given in (16). The conditional distance distribution of th rank user can thus be given as:
(17) 
Accordingly, the conditional distribution of the distance of th rank user, such that and , is the same as the distribution of the th order statistics of i.i.d. distance variables whose PDF is the truncated PDF of as given in (15). The conditional distance distribution of th rank user can thus be given as:
(18) 
Proof.
The proof of the first part of the lemma follows by substituting the PDF from (16) and its corresponding CDF into (A.1) and replacing with and with . Similarly, the proof of the second part of the lemma follows by substituting the PDF from (15) and its corresponding CDF into (A.1) and replacing with and with . ∎
V Rate Coverage Analysis
Using the distance distributions derived in Section IV, in this section, we derive the Laplace transforms of the intracluster interferences experienced by the transmission of th rank user considering both the perfect and imperfect SIC scenarios. We then derive the Laplace transform of the intercluster interference incurred at the representative BS by exploiting the distance distributions from geometric probability and a less complex upper bound of the Laplace Transform of the intercluster interference. Finally, we derive the rate coverage expressions for NOMA as well as OMA system.
Va Laplace Transforms of the IntraCluster Interferences
As defined in (6), it is evident that the intracluster interference with perfect SIC is due to all users that are located beyond the distance . Subsequently, the intracluster interference defined in (6) can be rewritten as follows:
(19) 
The Laplace transform of can then be given as follows.
Lemma 4 (Laplace transform of the IntraCluster Interference with Perfect SIC).
The intracluster interference experienced by the transmission of th ranked user in the cluster, with perfect SIC, can be given as follows:
(20) 
where is the generalized incomplete Beta function.
Proof.
See Appendix C. ∎
For cases of practical interest such as , the result presented in the Lemma 4 can be further simplified as follows.
Corollary 1 (Laplace transform of the IntraCluster Interference with Perfect SIC, ).
When , the Laplace transform of the intracluster interference experienced by the transmission of th ranked user in the cluster, with perfect SIC, can be further simplified as follows.
where (a) is derived by using the property . Note that, if , which represents the farthest user from the BS, there is no intracluster interference since .
Further, an accurate secondorder approximation of with a maximum absolute error of 0.0053 rad can be derived as . With this approximation, the closedform expression presented in Corollary 1 can be further simplified.
The intracluster interference incurred at a BS with imperfect SIC is defined in (9). It can be seen that the intracluster interference is composed of two parts, i.e., the interference observed from all users in set that are located beyond the distance and the interference from users in set whose signals go undetected. The Laplace transform of the intracluster interference with imperfect SIC can then be derived in closedform as follows.
Lemma 5 (Laplace transform of the IntraCluster Interference with ImPerfect SIC).
The intracluster interference experienced by the transmission of th ranked user in the cluster, with imperfect SIC, can be given as follows:
where (a) follows from the fact that, conditioned on the serving distance of th user, and are independent and the interfering distance distributions of their respective users are also independent. Note that is given in Lemma 4 and can be derived as follows:
where and is the incomplete Beta function.
Proof.
See Appendix D. ∎
VB Laplace transforms of the InterCluster Interference
Typically, the Laplace transform of the intercluster interference in PCP is characterized for a fixeddistance typical link (where the receiver is not chosen from the PCP) [12]. Recently, a more general approach is presented in [13, 14] to characterize the intercluster interference in which both the transmitter and receiver can be chosen from the PCP. Also, it was shown in [13, 14] that both the intra/intercluster interfering distances can be modeled using Rice distributions for Modified Thomas Cluster Processes.
Given the definition of the intercluster interference in (5), and following a similar approach proposed in [13], we can write the Laplace transform as follows:
where (a) follows from the definition of the Laplace transform and (b) follows from the property of the exponential function. Conditioned on the distance from the representative BS to the cluster center , the distance of each user device within the cluster (whose center is located at ) to the representative BS is i.i.d. Subsequently, (c) follows from the fact that the distances of all users in an interfering cluster to the reference BS are i.i.d., (d) follows from the PGFL of PPP since all cluster centers follow a homogeneous PPP, and (e) follows from the conversion of Cartesian to polar coordinates. The Laplace transform can thus be derived as follows.
Lemma 6 (Laplace Transform of the InterCluster Interference).
The integrals and are solvable in closedform as detailed in Appendix E and the Laplace transform as well as the rate coverage can be evaluated numerically using MAPLE and MATHEMATICA.
Note that the Laplace transform in (d) can also be attained from the PGFL of an MCP with fixed number of daughters per cluster as defined below.
Definition 3 (PGFL of the Matern Cluster Process).
The PGFL of the MCP given the number of nodes are fixed per cluster can be given as follows:
(24) 
where denotes the intensity of the parent point process which is a homogeneous PPP in case of MCP and . Since the original cluster process is stationary, the intercluster interference is independent of the position of the receiver [12].
Now observing that , we can apply Jensen inequality to derive an upper bound on the Laplace transform of the intercluster interference. That is
(25)  
(26)  
(27) 
where (a) follows from applying the Jensen inequality and (b) follows from applying the cartesiantopolar coordinate transformation and integration by parts [15].
VC Rate Coverage Probability for NOMA
In this section, we derive the rate coverage probability of a user at rank for all three cases, i.e., perfect SIC, imperfect SIC, and worstcase SIC.
VC1 Rate Coverage with Perfect SIC
Using the distance distribution of the user at rank derived in Lemma 1 and the Laplace transforms of the inter and intracluster interferences derived, respectively, in Lemma 4 and Lemma 6, the coverage of a user at rank can be derived as follows:
(28) 
where (a) follows from averaging over the distribution of and the definition . Since NOMA systems are typically interference limited, the expression in (a) can be further simplified by substituting . Also, taking and taking the upper bound of the Laplace transform of the intercluster interference in (25), we can simplify (VC1) as follows:
where
VC2 Rate Coverage with Imperfect SIC
In this case, a given user at rank is prone to the interferences from all users in the set as well as from some users in the set whose signals go undetected. As such, using the distance distribution of the user at rank derived in Lemma 1, the Laplace transforms of the intra and intercluster interferences derived, respectively, in Lemma 5 and Lemma 6, and the interference limited case, the coverage of a user at rank in the representative cluster can be derived as follows:
(29) 
VC3 Rate Coverage with Imperfect SICWorst Case
VD Rate Coverage Probability for OMA
The rate coverage performance of a user at rank in a given OMA cluster can be derived as follows.
Corollary 2 (Rate Coverage Probability of a User at Rank in OMA).
The rate coverage probability of a user at rank in TDMAbased OMA system can be given by taking zero intracluster interference, using the Laplace transform of the intercluster interference from (23) with =1, and replacing with as follows:
(30)  
(31) 
where and is the regularized confluent Hypergeometric function.
Note that (b) is derived by neglecting the noise, taking the bound of the intercluster interference in (