Cooperative NOMA system with incremental relaying over Nakagami-m fading channels

In our article, we apply incremental relay to the cooperative non-orthogonal multiple access network, in which the incremental relay employs amplify and forward protocol. The main communication process of the system is that one base station transmits messages to multiple non-orthogonal multiple access users with the assistance of incremental relay. In addition, we deeply analyze the outage behavior of the incremental cooperative non-orthogonal multiple access system over Nakagami-m fading channels. Specifically, the exact and asymptotic closed-form expressions of outage probabilities are all derived for paired users u n and u f , which are nearby and distant user, respectively. On this basis, the diversity orders of the paired users u n and u f are also acquired, which are μ (n+f+1) and μ(f+1). We also investigate the conventional cooperative non-orthogonal multiple access system, which serves as a benchmark of the incremental cooperative non-orthogonal network. Through the proof, we can conclude that the outage behavior of the distant u f has been improved by adopting incremental relay. Finally, the performance of the system is simulated. The following conclusions can be drawn from the numerical simulation results: (1) the outage performance and throughput of the system have been significantly improved by using the incremental relay protocol; (2) as the multipath fading μ increases, the outage performance of the system becomes better; and (3) the greater the difference in channel gain between paired users, the incremental cooperative non-orthogonal system can achieve better outage performance.


Introduction
With the explosive growth of wireless communication devices, such as the Internet of Vehicles and multimedia equipment, massive connection becomes the main feature of wireless communication network. Subsequently, the requirements for low-latency traffic and high spectral efficiency have stimulated the latest innovations in wireless technologies. 1 Non-orthogonal multiple access (NOMA) has proven to be a revolutionary technology to enhance spectrum efficiency 2,3 in wireless communication networks. NOMA scheme includes power-domain and code-domain NOMA; we focus on power-domain NOMA and use NOMA to represent power-domain NOMA in this article, as multiple NOMA users can use the same resources (time/ frequency/code) to communicate at the same time, that is, compared with previous orthogonal multiple access 1 School of Internet of Things Technology, Wuxi Institute of Technology, Wuxi, China 2 School of Emergency Management, Henan Polytechnic University, Jiaozuo, China (OMA), NOMA improves resource utilization. When implementing the NOMA scheme, the transmitters use superposition coding for message transmission, while the receivers employ successive interference cancelation (SIC) algorithm for decoding. 4 Many studies have shown that NOMA effectively enhance the spectrum efficiency and throughput in large-scale cellular networks with massive connectivity. [5][6][7] The researchers first carried out research on the point to point NOMA communication networks, mainly including the performance analysis and optimization of the downlink and uplink point to point communication networks. Ding et al. 8 investigated a downlink NOMA system with users deployed randomly, and the research demonstrates that the outage performance of users mainly depends on the power allocation and the setting of the target rate. Under non-ideal channel state information (CSI), 9 the authors investigated the downlink NOMA system with uniformly deployed users. The above research indicates that NOMA is better than OMA in terms of interrupt performance and sum rate. Shi et al. 10 further studied the impact of power allocation and other aspects on system performance. Under different fairness constraints, the power allocation factor and decoding order were set for each user and the average outage probability has been analyzed. On the basis of the above research, the authors analyzed the outage behaviors and best diversity of downlink NOMA 11 in reference to the situation that each user feeds back one bit CSI information to the base station (BS). Liu and Petrova 12 have investigated a downlink NOMA communication network, which is based on dynamic proportional fair scheduling for users, and the researchers explored the upper limit performance of the NOMA network under ideal and actual conditions. Moreover, the authors have combined energy harvesting and NOMA technology in a network 13 with a large number of devices connected on a large scale and conducted a joint evaluation of the interference management and energy harvesting capabilities based on the actual situation of the network. These analysis results show that NOMA is superior to OMA, which also provides fruitful insights for the practical application.
To reduce the influence of wireless channel fading on the communication coverage of wireless communication systems, researchers combined cooperative communication 14 technology and NOMA in the communication system. Ding et al. 15 first proposed the concept and scheme of cooperative NOMA, where the nearby user acts as a relay to assist information transmission to the distant user. In the cooperative NOMA system of Ashraf et al., 16 the nearby user was employed as a relay with decode and forward (DF) and amplify and forward (AF) protocol, and then, the outage performance of the NOMA system was analyzed. In addition to the above-mentioned halfduplex communication, full-duplex (FD) 17 NOMA systems have also been extensively studied to improve spectral efficiency, as the information can be transmitted and received at the same time. The authors investigated the outage behaviors of a cooperative NOMA network under FD mode. 18 It has confirmed that the proposed cooperative NOMA system with FD mode is superior to that of HD mode. Furthermore, the influence of Loop Self-Interference (LI) signals in FD mode has been considered. 19 The performance analysis of the cooperative NOMA system 19 with user relays working in FD and HD modes has been carried out, respectively.
Apart from the above researched, the NOMA system with dedicated relays has also been carried out. The research of literature [20][21][22] has shown that the cooperative NOMA system using AF relay is significantly better than that of OMA including coding gain, outage behavior and throughput. In addition to the above analysis of the NOMA system under the Rayleigh fading channel, Chu and Zepernick 23 analyzed the cooperative NOMA system with DF relay under the Nakagami-m channel, including outage probability and the sum rate. Li et al. 24 studied the energy harvesting and information transmission of cooperative NOMA networks with the consideration of residual hardware damage and channel estimation errors. On this basis, the security performance of the environmental backscatter cooperative NOMA system 25 has been considered in the case of imperfect SIC; in addition, a kind of unmanned aerial vehicle-assisted cooperative NOMA network has been studied. 26 Besides, the authors have introduced a small cell network operating in a heterogeneous cellular network environment 27 by deploying FD, energy harvesting, and NOMA scheme, which have improved the spectrum efficiency and outage performance. Furthermore, the problem of relay selection in cognitive radio NOMA network 28 has been discussed under the spectrum sharing mode. Then, NOMA has been applied to the cognitive radioassisted NOMA-vehicle-to-vehicle (V2V) scene, 29 where two V2V transmissions schemes were proposed to enhance the performance of vehicles.
While the conventional cooperative NOMA (CCN) network effectively improves the spectrum utilization, there are still some shortcomings in the HD and FD CCN network. Specifically, the HD CCN network requires twice the time for information transmission compared with the non-cooperative NOMA system, which restricts the improvement of spectrum efficiency. Moreover, the loop interference (LI) signal in the FD CCN system restricts the improvement of diversity gain and spectral efficiency. Therefore, taking into account the application research of the incremental relaying (IR) protocol 14 in the cooperative communication system, since the IR protocol can effectively improve the spectral efficiency. The characteristic of incremental relay is to determine whether the system adopts the cooperative mode through one bit negative-feedback signal. When it is predicted that the direct link between the source and the receiver cannot correctly decode the received signal, relay is used to coordinate information transmission. In previous studies, most of the researches were carried out on the performance of cooperative OMA systems with the IR protocol, 30,31 which shows that the IR cooperative system can achieve higher throughput compared with the conventional cooperative systems. Driven by this, we consider introducing the IR protocol into the cooperative NOMA system for in-depth research.

Motivation and contributions
Although the CCN systems has been extensively studied in improving user access density, the research on Incremental cooperative NOMA (ICN) system still needs to be in-depth to further improve the system spectrum efficiency. While the aforementioned research on cooperative NOMA systems have been laid a solid foundation for improving user access density, the ICN system is still under exploration for improving system spectrum efficiency. From previous research, it is worth pointing out that the NOMA communication networks are suitable for IoT scenarios, especially where user channel conditions are quite different. In this scenario, users with different channel conditions are paired to perform NOMA. In addition, extensive research has been carried out on cooperative NOMA to improve communication coverage. Without loss of generality, cooperative NOMA can be divided into HD and FD mode, respectively. As mentioned above, there are still some drawbacks in the traditional HD and FD cooperative NOMA system. Moreover, cooperative NOMA systems are almost researched under Rayleigh channels. Inspired by these, we conducted research on the ICN system under Nakagami-m fading channels. We have summarized the contribution of this article as follows: 1. We researched the ICN system, where the system with fixed gain AF relaying is studied over Nakagami-m fading channels. To characterize the outage behaviors of the proposed ICN system, we derive closed-form expressions of the outage probability for the paired nearby user u n and distant user u f . Furthermore, we obtain the asymptotic outage probability and diversity orders of two users paired. Through proof, it can be concluded that the outage behaviors of u f in the ICN network have been significantly improved.
2. We also do simulation on the basis of theoretical analysis. Compared with the CCN network, we confirm that the users paired of the ICN network achieve superior outage behavior and throughput. Specially, the ICN system can achieve better outage performance when the channel gain difference between the two users paired is greater.
The rest of this article is organized as follows. In the system model part, the NOMA system model with incremental relaying over Nakagami-m fading channels is presented. In the outage performance part, the outage behaviors of two users paired in the proposed system are researched. Simultaneously, we have also conducted in-depth research on the CCN system as a comparison. In the numerical results part, The theoretical analysis is verified by The theoretical analysis is verified by simulation analysis. At the end, we conclude this article in conclusion part.

System model
The system model including network description and signal model will be explained in the following two subsections.

Network description
In this article, a downlink cooperative NOMA network with incremental relaying is considered over Nakagamim fading channels. As the system model shown in Figure 1, one BS communicates to two paired users with the assistance of the incremental relay R, where the two paired users are a nearby user u n and a distant user u f among all N users. All nodes in this system are equipped with single antenna. In addition, consider the scenarios where AF protocol is considered at the relay. h sr , h ru n , and h ru f . Moreover, we assume the channels of N users follow independent Nakagami-m fading distribution. In general, the order of channel gains for N users is h su 1 j j 2 ł h su 2 j j 2 ł Á Á Á ł h su N j j 2 . The communication process of ICN system is described as follows. A complete communication process consists of two time slots. BS broadcasts the superposed signal of two users to u n , R and u f in the first time slot. Once the superposed signal is received, u f compares the SINR of the decoded own signal with a predetermined decoding threshold. If the signal of u f can be decoded through the direct link, u f will feed back a positive signal to BS and R. BS employs direct NOMA transmission mode within the whole transmission process after receiving the positive feedback, that is, the second time slot repeats the communication process of the first time slot. Otherwise, BS employs cooperative NOMA transmission mode while receiving the negative feedback, that is, R will forward the received superimposed signal in the second time slot.

Signal model
In the first time slot, BS transmits superimposed signal ffiffiffiffiffiffiffiffiffi a n P s p x n + ffiffiffiffiffiffiffiffiffi a f P s p x f to R, u n and u f , where P s is the transmit power of BS, x n and x f denote the power signal for u n and u f , respectively. a n and a f are the corresponding power allocation coefficients of u n and u f . Specially, assuming a n ł a f and a n + a f = 1, 32 which is to ensure fairness and QoS among users. Hence, the signal received at u n , u f , and R can be expressed as y n = ffiffiffiffiffiffiffiffiffi a n P s p h su n x n + ffiffiffiffiffiffiffiffiffi a f P s p h su n x f + v n ð1Þ and y r = ffiffiffiffiffiffiffiffiffi a n P s p h sr x n + ffiffiffiffiffiffiffiffiffi a f P s p h sr x f + v r ð3Þ respectively; where v n , v f , and v r are the additive white Gaussian (AWGN) at u n , u f , and R, which follow the complex Gaussian distribution with zero mean and variance s 2 , respectively. According to NOMA principle, the SIC scheme is first implement at u n to decode x f ; hence, the SINR of decoding x f at x f is where r s = P s =s 2 denotes the transmit SNR of each link between BS and other node. Subsequently, u n executes the SIC scheme to eliminate the signal x f , the SNR of u n decoding its own signal is given by Finally, the SINR of u f decoding its own signal can be written as If u f estimates that it can correctly detect x f through the direct link, the communication mode is still repeated in the remaining next time slot.
Otherwise, BS employs cooperative NOMA transmission mode while receiving the negative feedback, that is, R will forward the superimposed signal in the second time slot. Consequently, assuming the transmit power of R is P R and P r = P s . The AF protocol is considered at R, and R amplifies and forwards the received superimposed signal to u n and u f . For simplifying, fixes gain factor G is employed at AF relay; then, the signals received at u n and u f are given as y AF n = Gh ru n y r + n 1 ð7Þ respectively, where G = ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi P r =P s ½Eð h sr j j 2 Þ + s 2 q and E Á f g denotes expectation operation. n 1 and n 2 are the AWGN at u n and u f and they follow the complex Gaussian distribution with zero mean and variance s 2 . Similar to the above decoding process, the SINR of x f decoded at u n can be expressed as g AF n, f = a f r s h sr j j 2 h ru n j j 2 a n r s h sr j j 2 h ru n j j 2 + h ru n j j 2 + 1=G 2 After the signal x f of u f is decoded and subtracted, the SINR of u n decoded its own signal can be written as g AF n = a n r s h sr j j 2 h ru n j j 2 Simultaneously, the SINR of x f decoded at u f can be expressed as Outage performance evaluation The following section will conduct a detailed analysis of the outage performance for the ICN system.

Outage probability analysis
First according to the description of the system model, we can express the probability density function (PDF) of y = h j j as in which m is the multipath fading parameter, v 0 is the control spread and G Á ð Þ is the Gamma function. So we can get the PDF and cumulative distribution function (CDF) of t = h j j 2 as follows where h j j 2 is unsorted channel gain. Combined with the order statistics, 33 we can obtain the PDF and CDF of h m j j 2 as follows and in which h m j j 2 is the mth user's channel gain after channel sorting.
Next, we analyze the outage behaviors of the two users paired in the system.
Outage probability of u n . Considering the communication scheme of the ICN network, the outage event for u n is divided into two categories. One is when the u f can successfully detect its own information by direct link, that is, the system is in direct link communication mode. On the contrary, the system works in cooperative communication mode. Hence, the outage probability of u n can be written as P ICN n, AF = P r g f ø g thf À Á 1 À P r g n, f ø g thf , g n ø g thn À Á Â Ã where the decoding threshold of u n and u f under direct link transmission mode are g thn = 2 R n À 1 and g thf = 2 R f À 1. Moreover, R n and R f are the target rate of u n and u f , respectively. Additionally, the outage probability of u n for CCN network is P CCN n, AF , which can be given as follows P CCN n, AF = ½1 À P r ðg n, f ø g 0 thf , g n ø g 0 thn Þ 3 ½1 À P r ðg AF n, f ø g 0 thf , g AF n ø g 0 thn Þ Similarly, the decoding threshold of u n and u f under cooperative transmission mode is g 0 thn 2 2R n À 1 and g 0 thf = 2 2R f À 1.
The closed-form expression of P ICN n, AF is presented in Theorem 1 as follows.
Theorem 1. In the ICN network, the closed-form expression of outage probability for nearby user u n is obtained as In addition e = g thf r s a f À a n g thf À Á b = h su n j j 2 ø g thn a n r s O = max e, b ð Þ e 0 = g 0 thf r s a f À a n g 0 denotes the modified Bessel function of the second kind with order x. Note that equation (19) is derived on the condition of a f .a n g th f and a f .a n g 0 th f .
In the derivation process of Theory 1, the closedform expression of P CCN n, AF for CCN system is given in the following corollary, which can be used as a comparison for ICN network performance. Corollary 1. In the CCN system, the closed-form expression of the outage probability for nearby user u n can be given by 3 K iÀt + 1 2m The condition for the above equation (20) to be true is a f .a n g 0 th f .
Outage probability of u f . Similarly, considering the communication scheme of ICN network, the outage event of u f occurs when u f cannot detect x f correctly through direct link. Therefore, the outage probability of u f is given as As g thf \g 0 thf , equation (21) can be further rewritten as As the derivation process of Theory 1, the closedform expression of P ICN f , AF can be given in the following Theorem 2.
Theorem 2. In the proposed ICN network, the closedform expression of outage probability for distant user u f is derived as Proof. See Appendix 2.
It is worth noting that the conditions for the establishment of equation (23) are a f .a n g thf and a f .a n g 0 thf , otherwise P ICN f , AF = 1. In addition, the outage probability of u f for CCN network is P CCN f , AF , which can be written as In the derivation process of Theory 2, the closedform expression of P CCN f , AF for CCN system is given as Corollary 2.
Corollary 2. In CCN system, the closed-form expression of the P CCN f , AF for distant user u f is given as Where and c 0 f = me 0 =v su f . Note that equation (25) is also derived on the condition of a 2 =a 1 . g 0 thf , otherwise P CCN f , AF = 1. (22) and (24), it is obvious that P ICN f , AF \P CCN f , AF , which shows that the performance of u f in the ICN network exceeds that of the CCN network.

Diversity analysis
Based on the analysis of the outage probability of u n and u f in the above part, we further analyze the asymptotic outage probability in high SNR region and obtain diversity order of two paired users. The diversity order 32 can be defined as Diversity order of u n and u f . When r s ! ', the approximate expression of P ICN n, AF for ICN system is obtained as Corollary 3.
Corollary 3. In the ICN system, the outage probability of u n at high SNR can be approximated as Proof. See Appendix 3.
Substituting equation (27) into equation (26), we can obtain d ICN u n , AF as where d ICN u n , AF is the diversity order of u n for ICN network.
Similarly, the approximate expression of P ICN f , AF for ICN system is obtained as Corollary 4.
Corollary 4. In the ICN system, the outage probability of u f at high SNR can be approximated as Substituting equation (29) into equation (26), we obtain d ICN u f , AF as where d ICN u f , AF is the diversity order of u f for ICN network.
According to Yue et al., 22 the diversity order of u n and u f in the CCN mode are d CCN u n , AF = m n + 1 ð Þ and d CCN u f , AF = m f + 1 ð Þ, respectively.
Remark 2. The diversity order of the nearby user u n in the ICN system is significantly higher than that of the CCN system.

Throughput analysis
For the proposed ICN network, as the information transmission rates are fixed value, the throughput of the network in delay-limited transmission mode mainly depends on the outage behaviors. Therefore, under delay-limited transmission mode, the throughput of ICN system can be expressed as Similarly, the throughput of the system under CCN mode is given as follows, which can be used as a benchmark of the ICN system throughput Numerical results Simulation evaluate and verify the performance analysis of the above ICN system, where the transmission is performed under Nakagami-m fading channels.
Assuming the total number of users is N = 5. The power distribution coefficients for u f and u n are a f = 0:8 and a n = 0:2, 32 respectively. Moreover, we assume that R f = 1 BPCU and R n = 1:5 BPCU are the target rate of u f and u n , where BPCU refers to bit per channel use. To simplify, the gain of the AF relay is fixed at C = 0:9. 22 Besides, v sr = 1=d a sr and v ru = 1=(1 À d sr ) a are the average power, which can be obtained by setting the pathloss exponent a = 3 and normalized distance d sr = 0:5. 21 In the following simulation analysis, Monte Carlo simulations are performed to validate the derived exact analytical results. The simulation results for the outage probability and system throughput are averaged over 10 6 channel realizations. It is emphasized that the CCN and OMA networks have also been investigated and simulated, which are used as a comparison for the ICN network performance. Figure 2 plots the outage probability curves of two users paired in the ICN network versus transmit SNR. Meanwhile, the outage probability curves under CCN and OMA networks are also simulated for comparison. The black, blue, and red curves represent the outage behavior of u n and u f for ICN, CCN, and OMA networks. According to equations (19) and (23), the exact theoretical curves of u n and u f in the ICN system are presented. It is obvious that the Monte Carlo curves coincide with the exact theoretical curves, which shows that our theoretical derivation is accurate. In addition, the asymptotic outage behavior curves of u n and u f in ICN network are presented as the black dashed line, which are plots according to equations (27) and (29). As shown in the figure that the approximate and precise curves of outage probability tend to be consistent under high SNR. Moreover, the exact theoretical curves of u n and u f in the CCN system are presented as blue curves, which are plots according to equations (20) and (25). By comparison, it is shown that the outage performance of the ICN system surpasses that of CCN and OMA system. Simultaneously, Remark 1 has also been verified. Figure 3 plots the exact outage probability curves of u n and u f in ICN network versus transmit SNR, where m = 2, 3, 4, respectively. From the curve in the figure, it can be concluded that an increase in the value of m will bring about an improvement in outage behaviors. Figure 4 plots the outage probability versus transmit SNR curves of two paired users with different path loss exponent. In this simulation, the path loss exponent a takes 3, 4, and 5 respectively; other parameters remain unchanged. When path loss exponent a = 3, 4, 5, the curves of the outage probability for the two paired users are black, red, and blue, respectively. From the simulation curve, it can be concluded that the outage behaviors of u n and u f is greatly affected by a. Specifically, the system performance becomes better as a becomes larger, mainly because the average power becomes larger as a increases. However, the diversity order of two paired users u n and u f does not change with the change of a, which also verifies the conclusions shown in equations (28) and (30). Figure 5 plots the outage probability versus transmit SNR curves of two paired users for different targeted data rates. In this simulation, the distant and nearby users have taken three pairs of different rates, which are R n = 1:5 BPCU, R f = 0:8 BPCU; R n = 1:5 BPCU, R f = 1:0 BPCU; and R n = 2:0 BPCU, R f = 1:0 BPCU, respectively. It can be seen from Figure 5 that the  outage behaviors of u n and u f deteriorate with the increase of R n and R f , respectively. These are mainly because the decoding threshold enhances as the rate increases. From the simulation results in Figure 5, another conclusion can be drawn that the performance of u f has nothing to do with R n , because u f does not decode the signal of u n when decoding x f , only takes the signal x f of u f as interference noise. These also verify the conclusions as equation (22). Moreover, the outage performance of u n is related to the change of R f ; because SIC is performed at u n , the signal x f is decoded and subtracted before the self-signal x f is decoded. Figure 6 plots the outage probability versus transmit SNR curves for different paired users in ICN network. In this simulation, a different total number of users are used and three different pairing methods for the distant and nearby users have adopted, which are N = 5, f = 1, n = 5; N = 10, f = 1, n = 5 and N = 10, f = 1, n = 10; respectively. From Figure 6, we can conclude that the outage behavior of u f is mainly affected by the total number of users N, when the order of u f in all users remains unchanged. Specifically, the outage behaviors of u f deteriorates as the N increases. As the order of u f among N users is fixed, the channel gain is mainly affected by the total number of users according to equations (15) and (16). Therefore, the decoding SINR of u f will change with the number of users, which in turn affects the outage performance of u f . Furthermore, it is shown from Figure 6 that the outage behaviors of u n is affected by its self-ordering and the value of N . When the value of N is fixed, as the sort value of the nearby user u n and distant user u f increases, the greater difference between the channel gains of them, the better outage performance of the nearby user u n . Similarly, if the ranking of u n and u f is fixed and the value of N is smaller, it means that the increase in the channel difference between u n and u f brings the performance improvement of u f . In general, it proves that NOMA is suitable for scenarios with large differences in user channel conditions and can significantly improve user outage behaviors in this scenario. Figure 7 plots the outage probability versus transmit SNR curves in ICN network for different power allocation coefficients. In this simulation, three different power allocation coefficients for the distant and nearby users have adopted, which are a n = 0:2, a f = 0:8;   a n = 0:15, a f = 0:85 and a n = 0:1, a f = 0:9; respectively. When a f = 0:8, a f = 0:85 and a f = 0:9, the outage probability curves of u n and u f are black, red and blue respectively. From Figure 7, it can be shown that the outage performance of u n deteriorates as the value of a f increases (a n decreases). Similarly, the outage behavior of u f improves as the value of a f increases. As a result, the optimization of user power allocation can also improve system performance. Figure 8 plots the ICN and CCN system throughput curves versus transmit SNR, which are plotted in delay-limited transmission mode according to equations (31) and (32), respectively. It can be concluded from Figure 8 that the system throughput of the ICN network exceed that of CCN system in limited delay mode. This is because that the system throughput is mainly affected by the probability of outage, and the outage behaviors of the ICN network has been improved. Similarly, the system throughput of ICN and CCN networks increases with the increase of m.

Conclusion
In this article, the IR protocol has been applied in the cooperative NOMA network, where the transmission was performed over Nakagami-m fading channels. Specially, the outage behavior of the ICN network were investigated, which included analysis of outage probability and diversity gain of two paired users. We obtained the closed form solutions of the outage probability and diversity order for the two paired u n and u f . In addition, the performance of CCN system has been a benchmark for ICN system. Finally, the performance of the ICN network has been simulated, the simulation results has shown that the outage performance and the throughput in the delay limited mode of the ICN network exceed the CCN network.

Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: Thank you for the support of the Natural Science Foundation of Jiangsu Higher Education Institutions. This work was support by Basic Science (Natural Science) Research Project of Higher Education Institutions in Jiangsu Province under grant 21KJB510035.
Q 1 = 1 À P r h su f 2 \ g th f r s a f À a n g th f