I Introduction
In the fifth generation (5G) communication systems and internet of things (IoT) networks, massive connectivity is required to support large number of devices in various scenarios with limited spectrum resources [2, 3]. Unmanned aerial vehicles (UAVs) communication is an effective approach to provide connectivity during temporary events and after disasters. Some initial research contributions in the field of UAVs communication have been made by researchers in [4, 5, 6, 7, 8]. UAVs can provide connectivity to multiple users as wireless relays for improving wireless coverage. On the other hand, ground users can have access to UAVs, which can act as aerial base stations, for reliable downlink and uplink communications [4, 5, 6]. Additionally, UAV networks are capable of enhancing spectrum efficiency by having lineofsight (LoS) connections towards the users, which provide higher received power for the users [7, 8].
Nonorthogonal multiple access (NOMA) has been considered as a promising technique in 5G mobile networks because of its superior spectrum efficiency [9, 10]. The key idea is that multiple users are served within the same frequency, time and code block [11, 12], which relies on the employment of superposition coding (SC) and successive interference cancellation (SIC) at the transmitter and receiver, respectively [13]. Moreover, SIC technique at receivers allows that users who have better channel conditions to remove the intrachannel interference. Ding et al. [14]estimated the performance of NOMA with fixed power allocation (FNOMA) and cognitive radio inspired NOMA, which demonstrated that only the user with higher channel gain influences the system performance. Timotheou and Krikidis optimized the userpower allocation problem to improve user fairness of a NOMA system [15]. To further exploit NOMA networks, the outage performance and capacity in the downlink and uplink transmission scenario with dynamic power allocation factors was estimated by Yang et al. [16], which can guarantee the quality of service in dynamic power allocation NOMA (DNOMA).
Ia Motivations and Related Works
1) Studies on MIMONOMA systems:
Current multipleinput multipleoutput (MIMO) NOMA applications can be primarily classified into a pair of categories, namely, beamformer based structure and cluster based structure. In beamformer based structure, one centralized beamformer serves only one user in a beam
[17]. Choi proposed a coordinated multipoint transmission scheme in [18], where two base stations (BS) serve paired NOMA users simultaneously by beamformer based structure. The multiple antenna scenario of beamformer based structure was proposed by Shin et al. in [19], where a joint centralized and coordinated beamforming design was developed for suppressing the intercell interference. Some related works on cluster based MIMONOMA have been investigated in [20, 21, 22, 22, 23]. More specifically, Ding et al. [20] proposed a MIMONOMA model with transmit power control and detection scheme. By adopting this design, the MIMONOMA model can be separated to multiple independent SISONOMA arrangements. Then, outage probabilities and the sum rate gap of multiple SISONOMA arrangements were estimated, which the model assumes that global channel state information (CSI) is unknown at the BS. The sum rate and multiuser capacity were investigated by Zeng et al. [24], which show that multiuser MIMONOMA is not a preferable solution of NOMA due to high computational complexity. In order to establish a more general framework of MIMONOMA, and circumvent the restrictive assumption of receiver antennas, the signal alignment technique was proposed for both downlink and uplink transmission scenarios in [21].2) Studies on NOMA in stochastic geometry systems: Stochastic geometry is an effective mathematical tool for capturing the topological randomness of networks. As such, stochastic geometry tools were invoked to model the impact of the locations for NOMA users [25]. Some research contributions with utilizing stochastic geometry approaches have been studied in [22, 23]. More particularly, Liu et al. [22] proposed an innovative model of cooperative NOMA with simultaneous wireless information and power transfer (SWIPT). In this model, the wireless power transfer technique was employed at users, where near users acted as energy harvesting relays for supporting far users. Ding et al. [25] evaluated the performance of NOMA with randomly deployed users. The analytical results show that it is more preferable to group users whose channel gains are more distinctive to improve the diversity order in NOMA system. With the goal of enhancing the physical layer security of NOMA networks, Liu et al. [23] proposed a NOMA assisted physical layer security framework in largescale networks. The secrecy performance of both single antenna and multiple antenna aided BS scenarios have been investigated.
3) Studies on UAV: In UAV aided wireless networks, the probability that each device has a LoS link is dependent on the environment, location of the device, carrier frequency and the elevation angle [4, 8, 26]. Sharma and Kim estimated the outage performance of single antenna assisted NOMA in UAV networks [27], where UAVground channels are characterized by LoS transmission. Recently, the general form of both LoS and NLoS transmission scenarios, namely, Nakagamim fading channels, have been proposed in the literature. Hou et al. [28] estimated the outage performance of FNOMA downlink transmission scenario in both LoS and NLoS scenarios.
While the aforementioned research contributions have laid a solid foundation with providing a good understanding of single input single output (SISO) NOMA terrestrial networks, how MIMONOMA technique is capable of assisting UAV networks is still unknown. To the best of our knowledge, there has been no existing work intelligently investigating the effect of the network performance of MIMONOMA assisted UAV networks, which motivates us to develop this treatise.
IB Contributions
The novel structure design in this work–by introducing the MIMONOMA assisted UAV networks–can be a new highly rewarding candidate, which will contribute the following key advantages:

We propose a general MIMONOMA aided UAV framework with interference, where stochastic geometry approaches are invoked to model the locations of users and interference sources. Utilizing this framework, LoS and NLoS links are considered to illustrate the general case of NOMA assisted UAV networks.

We derive closedform expressions for outage probability of paired NOMA users in the proposed framework. We provide tractable analytical upper bounds for both LoS and NLoS scenarios. Diversity orders are obtained for the paired NOMA users based on the developed outage probability. The obtained results confirm that for the case that interference power is proportional to the transmit power, diversity orders of paired NOMA users are zero.

We derive exact analytical expressions for ergodic rate in both LoS and NLoS scenarios. We provide tractable analytical lower bounds for the general case. We also provide closedform expressions for the special case when the path loss is three. We obtain high SNR slopes for the paired NOMA users based on the developed ergodic rate. The obtained results confirm that the high SNR slopes of far users are zero in both LoS and NLoS scenarios.

We demonstrate that 1) the outage performance and ergodic rate can be enhanced by the LoS propagation; 2) the ergodic rates of near users have the same rate ceiling in both LoS and NLoS scenarios; 3) diversity orders and high SNR slopes of paired NOMA users are not affected by the LoS transmission of the proposed framework.
IC Organization and Notations
The rest of the paper is organized as follows. In Section II, a model of UAVaided transmission scenario is investigated in wireless networks, where NOMA users are uniformly allocated on the ground. Then precoding and detection strategies have proposed for UAV networks. Analytical results are presented in Section III to show the performance of UAVaided MIMONOMA networks. Our numerical results are demonstrated in Section IV for verifying our analysis, which is followed by the conclusion in Section V. Table I lists all notations used in this article.
and  Power allocation factors 

Transmit power of the UAV  
Transmit power of interference sources  
and  Detection vectors 
Path loss exponent  
and  The radius of small disc and large disc 
and  The target rate of the th user and of the th user 
Smallscale fading parameter 
Ii MIMONOMA Assisted UAV Network Model
Consider a MIMONOMA downlink communication scenario in which a UAV equipped with antennas is communicating with multiple users equipped with antennas each. Fig. 1 illustrates the wireless communication model with a single UAV, which is supported by MIMO beamforming, namely the clusterbased MIMONOMA. Multiple users are grouped into one cluster, and the UAV serves multiple clusters simultaneously in the clusterbased MIMONOMA, which can perfectly improve the system performance [17, 29].
Iia System Description
The UAV cell coverage is a disc area, denoted by , which has a coverage radius of
. It is assumed that the users are uniformly distributed according to HPPPs distribution, which is denoted by
and associated with the density , within large disc and small disc with radius and , respectively. For simplicity, we only focus our attention on investigating a typical user pairing in this treatise, where two users are grouped to deploy NOMA transmission protocol.The downlink users also detect signals sent by interferers which are distributed in according to HPPPs distribution with density [30]. It is assumed that the interference sources are equipped with one antenna each and use identical transmission power, which is denoted by .
IiB Channel Model
Consider the use of a composite channel model with two parts, largescale fading and smallscale fading. denotes the channel coefficient between the th antenna of the UAV and the th antenna of the th user, and the channel gain is modelled as
(1) 
where and denotes the largescale fading and smallscale fading, respectively. It is assumed that and with , are independent and identically distributed (i.i.d.). In this paper, largescale fading represents the path loss and shadowing between the UAV and users. Generally speaking, the largescale fading between the th user and the UAV changes slightly for different antenna pairs, which can be expressed as
(2) 
For simplicity, the smallscale fading matrix from the UAV to user is defined as
(3) 
where is a matrix whose elements represent Nakagamim fading channel gains. The density function of the elements is , where denotes the fading parameter. In this paper, the UAV can be projected to the coverage disc by projection theorem. Thus, the distance between the UAV to user can be written as
(4) 
where denotes the height of the UAV, and is the horizontal distance between the th user and projective point of the UAV. Therefore, the largescale fading can be expressed as
(5) 
where denotes the path loss exponent. For simplicity, we still use to represent the largescale fading between the UAV and the th user in Section II. Thus, the received power for the th user from the UAV is given by
(6) 
where denotes the transmit power of the UAV. Besides, it is assumed that the CSI of users is perfectly known at the UAV. The proposed MIMONOMA assisted UAV network model for downlink transmission is described in the following subsection.
IiC Downlink transmission
In this subsection, we estimate the downlink quality of MIMONOMA assisted single UAV networks, where two users, i.e., the th user and the th user, are grouped to perform NOMA. We consider a more practicable method that the number of antenna equipped on each user is assumed to greater than , i.e., . In the beginning, it is assumed that the maximum clusters is . Therefore, the UAV sends the following informationbearing vector to users as follows:
(7) 
where is the signal intended for the th user, is the power allocation coefficient, and . The cochannel interference can be further expressed as follows:
(8) 
where denotes an all one matrix, denotes the distance from the th user to the th interference source. and are the HPPPs distribution and the single antenna transmission power of interference sources, respectively. The interference model omits the small scale fading, since the effect of path loss is the dominant part for the long distance interferences.
Similar to [31, 21], the detection vectors can be designed as follows:
(9) 
where denotes a vector, which is normalized to 2, i.e., . denotes a matrix containing the right singular vectors of
corresponding to its zero singular values.
Define as the effective channel vector shared by users in the cluster, , and denotes the th element on the main diagonal of the matrix. Thus, with the design of precoding matrix and detection vector in [31, 21], the signaltointerferenceplusnoise ratio (SINR) of two users can be derived as
(10) 
and
(11) 
where is the transmit SNR, and . Again, it is assumed that the th user is the user with higher channel gain. Thus, based on NOMA protocol, the th user treats the signal of user as noise, and the th user can neglect the signal from th user by SIC technique. The use of signal alignment strategy reforms smallscale fading gains of two users into one matrix. Hence, it can be considered that the th user and the th user are sharing the same smallscale fading matrix.
Iii Performance Evaluations
In this section, we discuss the performance of downlink MIMONOMA assisted UAV networks. In this paper, fixed power allocation is employed at the UAV. New channel statistics, outage probabilities and ergodic rates are illustrated in the following three subsections.
Iiia New Channel Statistics
In this subsection, we derive new channel statistics for MIMONOMA assisted UAV networks, which will be used for estimating the outage probabilities and ergodic rates in the following subsections.
Lemma 1.
Assuming 2 randomly located NOMA users in the disc of Fig. 1, the detection vector of signal alignment assisted UAV networks over Nakagamim fading channels is
(12) 
where is the fading parameter of Nakagamim fading channels between the NOMA user and the UAV.
Proof.
Please refer to Appendix A. ∎
Remark 1.
The result in (12) indicates that the channel coefficients for the considered user pair are determined by the fading parameter of Nakagamim fading channels.
Remark 1
provides insightful guidelines for estimating the outage performance of the proposed frameworks. The paired users share the same small scale fading channel coefficient, which follows exponential distribution with variance
, the performance of NOMA user pair can be calculated by their stochastic geometric approach.IiiB Outage Probabilities
In this subsection, we first focus on the outage behavior of the far user , who is the user with poorer channel gain. The fixed power allocation strategy is deployed at the UAV, which the power allocation factors and are constant during transmission. It is assumed that the target rate of user is . Therefore, the outage probability of the th user is given by
(13) 
The detection vector and modified channel factor are difficult to evaluate, which renders the evaluation of the outage probability very challenging. Therefore, we assume that is 2, which is the greatest value under the constraint of (9). In addition, . Hence, the outage probability of the th user can be transformed into
(14) 
where . Therefore, is the upper bounds of the outage probability for and .
Constraint 1 (Target rate). The target rate of the th user satisfies , where . Otherwise, the outage probabilities of paired NOMA users are one.
Constraint 1 provides the basic guideline of parameter setting for MIMONOMA assisted UAV networks with fixed power allocation. In NOMA, the far user treats the signal of the th user as noise, and decodes the signal directly by SIC technique. Thus, the outage probability of the th user is one unless the target rate of the far user must under the constraint . On the other hand, the near user needs to decode the signal intended for the far user before decoding its own signal. Thus, Constraint 1 is needed to satisfy for both NOMA users.
Then we turn our attention to calculating the outage probability of the far user with interference source, which is given in the following Theorem.
Theorem 1.
Given any random vector , the upper bound of the outage probability for the far user can be expressed as
(15) 
where , , , , and represents the lower incomplete Gamma function.
Proof.
Please refer to Appendix B. ∎
In the proposed framework, the lower bound of the outage probability, which is meaningless, is only affected by power allocation factors due to that the minimum detection vector approaches zero. Thus, in this article, the lower bound of outage probability will not be estimated.
It is challenging to solve the integral in (15) directly due to the lower incomplete Gamma function. Thus, in order to derive the diversity order and gain further insights into the system’s operation in the high SNR regime, the asymptotic behavior is analyzed, usually when the SNR of the channels between the UAV and users is sufficiently high, i.e., when the transmit SNR obeys .
Corollary 1.
Assuming that , the asymptotic outage probability of the th user is given by
(16) 
where .
Proof.
Please refer to Appendix C. ∎
Remark 2.
The derived results in (16) demonstrate that the outage probability of the far user can be decreased by increasing the fading parameter or decreasing the target rate of the far user itself.
Proposition 1.
From Corollary 1, one can yield the diversity order by using the high SNR approximation, and the diversity order of the th user in the proposed MIMONOMA assisted UAV networks is given by
(17) 
The exact expression in the case of is also worth calculating, which is given in the following corollary.
Corollary 2.
For the special case that the interference sources does not exist, the exact expression can be written as
(18)  
In NOMA, the th user needs to decode the signal intended for user before decoding its own signal via SIC technique. Therefore, the outage probability of the th user is given by
(21)  
where denotes the target rate of the near user . Again, we focus on the upper bound of the outage probability for user . By some algebraic handling, the upper bound of the th user’s outage probability can be rewritten as
(22) 
where , , and . The expression of outage probability for the near user is provided in the following theorem.
Theorem 2.
Under Constraint 1, the upper bound of outage probability for the near user can be written as
(23) 
if
Proof.
Please refer to Appendix D. ∎
Based on Theorem 2, we can calculate the asymptotic result of the near user in the following corollary.
Corollary 3.
In the case where approaches infinity, it is straightforward to attain that goes to zero. Therefore, the asymptotic result can be obtained as
(24) 
where .
Proof.
Similar to Appendix C, the asymptotic outage probability of the near user can be proved. ∎
Remark 3.
Due to the fact that two users share the same small scale fading matrix by signal alignment technique, one can obtain that the diversity order of the th user is also one.
Proposition 2.
Remark 4.
The diversity order in (25) demonstrates that the power of interference sources, which is proportional to transmit SNR, influences the outage performance dramatically.
Corollary 4.
For the special case that the interference is limited, the exact expression for the near user can be rewritten as
(26)  
IiiC Ergodic Rates
The ergodic rate is a critical metric for performance evaluation. Therefore, in this subsection, we focus on analyzing the ergodic rates of individual users, which are determined by their channel conditions and geometry parameters in the proposed framework. The ergodic rates of the far user is given by
(27) 
Thus, the ergodic rate for the far user can be obtained in the following theorem.
Theorem 3.
The closed expression of ergodic rate for the th user is given by
(28) 
where , and .
Proof.
Please refer to Appendix E. ∎
Remark 5.
One can observe from the result in (28) that the ergodic rate of the far user is entirely dependent on the power allocation factors in the high SNR regime.
The integration of lower incomplete gamma function is quit challenging to estimate. Thus, we only derive the ergodic rate of the near user in two special cases, i.e., high SNR approximation and exact expression in the case of fixed path loss exponent.
Theorem 4.
When transmit SNR obeys , the approximation ergodic rate for the th user is given by
(29)  
where .
Proof.
Please refer to Appendix F. ∎
The exact expression of the ergodic rate for the near user is difficult to calculate. Therefore, the special case, which path loss exponent is fixed, is shown in the following corollary.
Corollary 5.
In the special case in which path loss exponent , the exact expression of the ergodic rate for the near user can be attained to
(30)  
where represents MeijerG function.
Proof.
Please refer to Appendix G. ∎
Remark 6.
To gather deep insights for the system performance, the high SNR slope, as the key parameter determining the ergodic rate in high SNR regime, is worth estimating. Therefore, we first express the high SNR slope as
(31) 
Proposition 3.
Remark 7.
In the proposed framework, the ergodic rate for the far user can be consider as a constant in the high SNR regime from Proposition 3.
Iv Numerical Studies
In this section, numerical results are provided to facilitate the performance evaluation of MIMONOMA assisted UAV networks. Mont Carol simulations are conducted to illustrate the correctness of analytical results. In the considered network, it is assumed that the power allocation factors are and . Some simulation parameters are summarized in Table II.
Mont Carlo Simulations repeated  times 
The radius of the large plane  20 m 
The radius of the small plane  10 m 
The height of the UAV  10 m 
Path loss exponent  
The interference sources density  
The radius of the interference sources  1000 m 
The noise power  dBm 
The antenna number of the UAV  4 
The antenna number of users  3 
Iva Outage Probabilities
In this subsection, the outage probability achieved by paired NOMA users with different choices of interference power in both LoS and NLoS scenarios is demonstrated in Fig. 2(a), Fig. 2(b), Fig. 3(a), and Fig. 3(b).
In Fig. 2(a), we evaluate the performance of the proposed approach for a network with different interference power levels in NLoS scenario. In Fig. 2(a), the black and blue curves are the outage probability for far users, and the red and green curves are the outage probability for near users. We can see that, as interference power increases, the outage probability of both the NOMA users increases. This is due to the fact that, as higher power level of interference sources are deployed, the received SINR decreases. It is also confirmed the close agreement between the simulation and analytical results in high SNR regime for the fixed interference power scenario. Besides, it is worth noting that all curves in fixed interference power scenario have the same slopes, which indicates that the diversity orders of the schemes are all one. This phenomenon validates the insights from Proposition 1 and Remark 3. Note that error floors appear when is proportional to the UAV power, which meet the expectation due to the strong cochannel interference. In this case, the dynamic interference power, which is a function of transmit SNR, is much greater than noise power in the high SNR regime. Thus, the error floor occurs in the case of high SNR situation. It is also worth noting that the asymptotic result does not exist due to the fact that lower incomplete gamma function cannot be approximated if the interference power is too large.
Fig. 2(b) shows the outage probability achieved by NOMA users in LoS scenario. In order to better illustrate the performance affected by LoS transmission, the NLoS case is also shown in the figure as a benchmark for comparison. In Fig. 2(b), we can see that as the fading parameter of small scale fading channels, , increases, the outage probability decreases for difference interference power levels. This is because that the LoS link between the UAV and users provides higher received power level. It is also shown that the diversity order for the LoS links is one in the case of fixed interference power. This is due to the fact that NOMA users share the same small fading coefficient matrix generated by signal alignment technique. The light blue curves are the outage probability of near users with interference in LoS and in NLoS scenario, respectively. Besides, several observations are drawn as follows: 1) When the interference power level is fixed, the diversity order of the transmission is one in both LoS or NLoS scenarios. 2) For the case that interference power is proportional to the UAV power, diversity orders are zero in both LoS and NLoS scenarios for NOMA users.
Fig. 3(a) plots the outage probability of NOMA users in LoS scenario versus the height of the UAV and its target rate. In this case, the transmit SNR is fixed to 60dB. One can observe that an outage ceiling exists even if the height goes to zero. This is because the fact that the target rate of far users is higher than the threshold, which is . It is also shown that when the disc radius is fixed, as the distance between the UAV and users increases, the outage probability of near users increases much faster than far users. This is a byproduct of the fact that when the height of the UAV is much larger than disc radius, the large scale fading between the UAV and NOMA users can be considered as the same. Hence, far users yield better outage probabilities because of higher received power from the UAV.
Fig. 3(b) plots the outage probability of paired NOMA users versus both the height of the UAV and the disc radius, with fixed target rate BPCU and fixed interference power dBm. Fig. 3(b) clearly illustrates the impact of disc radius to the proposed framework. We can see that the outage probability of near users increases as disc radius increases. It is also worth to mention that when the height goes to zero, the proposed UAV framework is degenerated to the traditional BS framework. Thus, the outage of far users occurs more frequently than near users, which is similar to previous literature. Besides, it is also demonstrated that the outage probabilities of NOMA users are decreased dramatically when decreasing target rates of paired users. The simulation results also confirms the insights from Fig. 3(a) that the outage probability of near users is higher than far users in the case of remote UAV.
IvB Ergodic rates
Fig. 4(a) compares the ergodic rates of individual users versus transmit SNR with different power levels of interference sources. Several observations can be drawn as follows: 1) An ergodic rate ceiling for far users exists even if the transmit SNR goes to infinity. This is due to the fact that the ergodic rate of far users is entirely affected by power allocation factors of paired NOMA users, which is a constant in FNOMA. 2) The black line and the red line are the ergodic rate of near users and far users with interference sources, respectively. As transmit SNR increases, the ergodic rate of near users increases because the received signal is increased, which improves the performance of ergodic rate. 3) The asymptotic simulations and analytical simulations are provided to confirm the close agreement between analytical results and the simulations. The dotted lines and dashed lines show the precise agreement between the approximation results and exact results. Thus, for the case of , the approximation results and exact results can be considered as the same. 4) As transmit SNR increases, the ergodic rates with and without interference sources for far users are approaching. This is due to the fact that the fixed interference power level is getting relatively smaller when SNR increased, which leads the received signal power relatively higher.
Fig. 4(b) depicts the ergodic rates of individual users versus transmit SNR in LoS scenario. The performance of NLoS scenario is also shown in the figure as a benchmark for comparison. As can be seen from the figure, far users in both NLoS and LoS scenarios have the same ergodic rate ceiling. This phenomenon validates the insight from Remark 5, where the ergodic rate of far users is only dependent on power allocation coefficients. Besides, the ergodic rate of near users with LoS link is higher than the NLoS situation. This is because the LoS propagation increases the received power level, which increases the ergodic rate of near users.
Fig. 5 plots the high SNR slope of paired NOMA users versus transmit SNR. It is observed that the high SNR slope of far users and near users goes to zero and one, respectively. This behavior can be explained as follows. The ergodic rate of far users, which changes slightly in the high SNR regime, is entirely determined by the power allocation factors. Thus, in the low SNR regime, as increases, the high SNR slope of far users is increased. On the other hand, the high SNR slope is zero for far users in the high SNR regime. For near users, the SNR slope increases monotonously, which shows the high SNR slope for near users is one in the high SNR regime. Another insight is that the LoS link accelerates the increasing rate and the decreasing rate of the paired NOMA users. On the other hand, the slopes of two scenarios are the same, which confirms Proposition 3.
V Conclusions
In this paper, the application of MIMONOMA assisted UAV networks was designed. Specifically, stochastic geometry based techniques were used for modeling both the locations of NOMA users and of the interference sources. Additionally, new analytical expressions for outage probability and ergodic rate were derived for characterizing the performance in MIMONOMA enhanced UAV frameworks. Diversity orders and high SNR slopes were also obtained to evaluate the system performance. It was analytically demonstrated that the diversity order was one due to the property of signal alignment technique in both LoS and NLoS scenarios. It was also shown that the ergodic rate of far users was only dependent on power allocation factors.
Appendix A: Proof of Lemma 1
In order to estimate the outage probability and ergodic rate of NOMA users, the channel factor can be rewritten as , where . Recall that the effective vector is generated by (9), i.e.,
(A.1) 
where all elements of
are independent and identically complex Gaussian distributed, and
is a random vector. Based on two random matrices, we can define two matrix and . It is quite challenging to estimate the general case. Thus, for the special case that , and the realvalue channel coefficients, two matrices can be written as(A.2) 
and
(A.3) 
where and denote the th and th channel coefficient between the UAV and users in the case of . Therefore, the base vector can be expressed as
(A.4) 
Thus, the variable can be transformed into
(A.5) 
Accordingly, we can rewrite the above matrix by BoxMuller transform [31, 32] to derive two variable as
(A.6) 
and
(A.7) 
where and
are independent random variables, and uniformly distributed in the interval
. Substituting (A.6) and (A.7) into (A.5), one can obtain(A.8) 
Replacing with will not change the density function because is also uniformly distributed in the interval . Therefore, we can consider the following distribution holds
(A.9) 
Hence, is a complex Gaussian vector with zero mean and variance, and the lemma is proved.
Appendix B: Proof of Theorem 1
First, the upper bound of outage probability can be written as follows:
(B.1) 
where . After some algebraic manipulations, the above function can be transformed into
(B.2) 
if . Denote , which is a constant in FNOMA. Recall that is exponentially distributed. Therefore, the upper bound of outage probability can be expressed as
(B.3) 
where .
In order to deduct the outage probability, we first focus on the interference sources, which are allocated under HPPP distribution with density . The HPPP distribution can be considered as stationary, and hence the summation of interference sources can be equivalently calculated by focusing on the interference reception seen at the node located at the th user. Accordingly, the expectation of interference sources can be obtain as by Laplace functionals [30]:
(B.4)  
where denotes the interferers allocated in the disc . For notation simplicity, denotes . It is assumed that the interference sources does not distributed in the disc to prevent infinite received signal power. After changing to polar coordinates, the above expectation can be converted as