Switchable Coupled Relays Aid Massive Non-Orthogonal Multiple Access Networks with Transmit Antenna Selection and Energy Harvesting

The article proposes a new switchable coupled relay model for massive MIMO-NOMA networks. The model equips a much greater number of antennas on the coupled relays to dramatically improve capacity and Energy Efficiency (EE). Each relay in a coupled relay is selected and delivered into a single transmission block to serve multiple devices. This paper also plots a new diagram of two transmission blocks which illustrates energy harvesting and signal processing. To optimize the system performance of a massive MIMO-NOMA network, i.e., Outage Probability (OP) and system throughput, this paper deploys a Transmit Antenna Selection (TAS) protocol to select the best received signals from the pre-coding channel matrices. In addition, to achieve better EE, Simultaneously Wireless Information Power Transmit (SWIPT) is implemented. Specifically, this paper derives the novel theoretical analysis in closed-form expressions, i.e., OP, system throughput and EE from a massive MIMO-NOMA network aided by switchable coupled relays. The theoretical results obtained from the closed-form expressions show that a massive MIMO-NOMA network achieves better OP and greater capacity and expends less energy than the MIMO technique. Finally, independent Monte Carlo simulations verified the theoretical results.


Introduction
Relays are proven suitable solutions for extending networking coverage. In addition, Multiple-Input-Multiple-Output (MIMO) Non-Orthogonal Multiple Access (NOMA) technologies are effective methods for enhancing the capacity of fifth-generation (5G) wireless communication networks and beyond. The NOMA technique improves spectral sharing and the capability for a greater number of connections in the same time slot/frequency [1,2]. The NOMA technique in fifth-generation (5G) wireless communication networks and beyond could therefore serve large numbers of devices. The NOMA technique's feature is the ability to broadcast a superimposed signal to devices in the network with different Power Allocation (PA) factors [3]. The farthest device with poorest Channel State Information (CSI) is allocated the largest PA factor compared to other PA factors allocated to other devices. Successive Interference Cancellation (SIC) [4] implemented at the receiver device removes interference before the device decodes its own message in the superimposed signal by treating other device information which has a stronger PA factor [5]. Major studies have significantly contributed to NOMA through the application of new techniques in 5G networks and beyond. In [6], the authors verified that the system performance of NOMA networks is dependent on power resource allocation. The main aims of NOMA networks • We present a new design for a diagram of two transmission blocks to calculate the propagation of wireless information and power (WIP). The present paper offers the potential for the practical application of wireless sensor networks (WSNs), e.g., in a water environment where relays and devices are barely powered [32]. • We maximize system throughput in a massive MIMO-NOMA network. The study deploys a TAS protocol which selects the best received signals from the pre-coding channel matrices. • The study delivers novel expressions for OP, system throughput, and EE in closedforms. We apply Monte Carlo simulation results to verify the analysis results.
This paper is structured as follows. Section 2 introduces system modeling. Section 3 presents an analysis of the system model. We describe and discuss the analysis and simulation results in Section 4. Section 5 presents a summary of the paper.

New Design for a Cooperative MIMO-NOMA Scheme
The present paper examines a new cooperative MIMO-NOMA model for emerging 5G wireless networks and networks beyond. Figure 1 depicts the system model containing a BS, coupled relay R = {R 1 , R 2 }, and multiple devices D n = {D 1 , . . . , D N }. To benefit a massive MIMO-NOMA network, the study equips a large number of antennas A 0 , A 1 , and A 2 at the network nodes and BS, coupled relays, and devices, respectively. In addition, we assume that the BS has full knowledge of the CSI [26,33]. As with the system models in the major studies [28][29][30][31], the massive MIMO underlay cooperative NOMA model shown in Figure 1 contains two time slots to complete a transmission block. It is important to illustrate the difference to previous major studies. In [16][17][18][19], the authors verified that relays are a good solution to assist in combating channel fading. In other major studies, the authors deployed a SISO relay to aid a multi-destinations [29] or a multiple SISO relay to aid a single destination [31]. The authors also proposed relay selection strategies to select the best nearest relay [17] and max-instantaneous data rate [18] (Equation (16)) to assist the destinations. The present study deploys coupled relays to assist massive MIMO-NOMA networks, however, only the relay with better power capacity among the coupled relays is selected for cooperating devices in an odd/even transmission block., i.e., while one relay is selected to forward the signal to devices, another relay has to maintain EH from the BS.
Note that the present study is designed to serve multiple devices simultaneously. The coupled relays must first have the device CSI and report this information to the BS. Based on the device CSI, the BS allocates PS factors, whereas the coupled relays feed back their own CSIs to the BS and other devices. Based on the energy capacity information, the BS selects the best powered relay to forward the signal. The devices then wait to receive the signal from the strongest powered relay. In the case of insufficient CSI, the BS may select the poorest powered relay to forward the signal, leading to a reduced lifespan of the relay. This may interrupt signal propagation because the devices are waiting for a non-cooperative relay. T r a n s m i t i n f o r m a t i o n T r a n s m i t e n e r g y

BS
Odd transmission block Even transmission block

Propagation and Formulations
From the model depicted in Figure 1, we designed a new propagation diagram, as shown in Figure 2. The coupled relays feed back data to the BS about their energy capacities. The BS decides which relay has more energy to forward the superimposed signal, while the remaining relay has less energy to maintain EH. Two main phases take place: EH and data transmission (DT). To illustrate, Figure 2 depicts two transmission blocks: an odd transmission block T (odd) and an even transmission block T (even) . Each transmission block T (θ) for θ = {odd, even} is separated into two equal time slots. The first time slot T (odd) 1 = T (odd) 2 in odd transmission block T (odd) is used by the BS to transmit wireless energy and superimposed information to coupled relays R = {R 1 , R 2 }. In terms of the SWIPT technique with a PS protocol, during the first time slot T (odd) 1 , a fraction λ of the power domain P S is used for EH while the remaining fraction (1 − λ) of the power domain P S is used to superimpose data from the BS. The second time slot T (odd) 2 = T (odd) 2 in odd transmission block T (odd) is used by the best powered relay to forward the superimposed signal to devices, while the worse powered relay applies EH from the BS, where the PS factor λ = 1. The present study assumes that relay R 1 is selected in the odd transmission block because relay R 1 has more energy than relay R 2 . Therefore, relay R 2 maintains EH from the BS. Similarly, the even transmission block T (even) is also separated into two equal time slots. The first time slot T (even) 1 = T (even) 2 in even transmission block T (even) is also used by the BS to simultaneously transmit wireless energy and the superimposed signal to the coupled relay R = {R 1 , R 2 }. However, in the second time slot T (even) 2 = T (even) 2 in even transmission block T (even) , the relay R 2 is selected to forward the superimposed signal to devices instead of relay R 1 because relay R 2 harvested energy from the BS, where λ = 1 during T (odd) 2 , and then relay R 2 contains more energy than relay R 1 . As a result, relay R 1 maintains EH from the BS with PS factor λ = 1 during T (even) 2 . Although the network model shown in Figure 1 is complex, it requires two time slots for signals to propagate through the network as in studies [28,29,31]. Specifically, the present study examines how multiple devices are served simultaneously. We therefore adopt the emerging NOMA technique. The BS also superimposes all the information of devices in the same signal by sharing the spectrum. The devices may be thus served simultaneously. As a result, the massive MIMO-NOMA network in the present study are low latency.

Odd Transmission Block
As shown in Figure 1, it is important to note that both the BS and coupled relays (R 1 and R 2 ) are equipped with multiple antennas, where A 0 > 1 and A 1 > 1. A 0 and A 1 are the number of antennas at the BS and coupled relays, respectively. Tran et al. [23] designed a pre-coding channel matrix size of [number of transmitting antennas × number of receiving antennas]. We therefore designed the pre-coding channel matrices from the A 0 transmitting antennas at the BS to the A 1 receiving antennas on the coupled relays R 1 and R 2 , respectively, as follows: where h (a 0 ,a 1 ) 0 ∈ H 0 and g (a 0 ,a 1 ) 0 ∈ G 0 for a 0 ∈ A 0 and a 1 ∈ A 1 are channels from a transmitting antenna a 0 at the BS to a receiving antenna a 1 at the coupled relays R 1 and where d R 1 and d R 2 are the distances from the BS to coupled relays R 1 and R 2 , respectively, and the coefficient ω is the path-loss exponent factor.
By applying the PS protocol, the first time slot in the odd transmission block T (odd) 1 is used for the BS to transmit wireless energy and superimposed information simultaneously. To illustrate, two phases take place. In the first phase, the BS sends wireless energy to the coupled relay with the PS factor λ. Therefore, the EH from the best channel in the precoding channel matrices given by (1) and (2) in the first time slot T (odd) 1 in odd transmission block T (odd) at coupled relays R 1 and R 2 are expressed as follows: where η is the collection factor, λ is the PS factor, and P S is the transmission power at the BS.
In the second phase of the first time slot T (odd) 1 in odd transmission block T (odd) , in a major advantage of NOMA theories, the BS broadcasts a superimposed signal by superimposing the messages x i of devices D i for i = {1, . . . , N} to the coupled relays R 1 and R 2 . We assume that no direct down-link exists from the BS to the devices. The received signal at the coupled relays is expressed as follows: where n R ∼ CN(0, N 0 ) is the Additive White Gaussian Noise (AWGN) at the coupled relays R = {R 1 , R 2 } with zero mean and variance N 0 . The PA factors for devices D i for i = {N, . . . , 1} are, respectively, denoted by α i , constrained to α 1 < . . . < α N , α 1 + . . . + α N = 1, and given by Note that Expression (7) is derived from the feature studied [7]. However, the devices in [7] were ordered where D 1 is the farthest device. However, the system model shown Figure 1 indicates that device D N has the farthest distance from the coupled relays. In term of NOMA theory, the farthest device must be allocated the biggest PA factor. As a result, the PA factor α N for device D N is the largest value among the PA factors, whereas device D 1 is the nearest distance from the coupled relays. Therefore, device D 1 is allocated the smallest PA factor.
SIC is another feature of NOMA theories which is implemented at the user. The user therefore implements SIC to detect messages in the received signal. In [7,12], the authors investigated NOMA networks with a random number N of users. The users repeated the SIC phases until their own messages were successfully detected in the received signal. It is important to note that N devices exist in our model (Figure 1), and, therefore, N SIC phases at the relay R 1 . After selecting the received signal as given in (5), in the first SIC phase, relay R 1 detects the message x N of device D N as a result of the constraint of the PA factors α N > . . . > α 1 . In the second SIC phase, relay R 1 detects the message x N−1 of device D N−1 after removing the x N symbol from the received signals. The relay R 1 repeats SIC until it successfully detects the last symbol x 1 .
However, the present study examines a massive MIMO underlay cooperative NOMA network in contrast to the schemes presented in [29,31], where the author studied cooperative Single-Input-Single-Output (SISO)-NOMA schemes. Fortunately, the authors of [23,24] also investigated a MIMO-NOMA network with TAS and obtained Signal-to-Interference-plus-Noise Ratios (SINRs), where the devices detected information by applying SIC. To optimize system performance, our study considers the massive MIMO technique in combination with a TAS protocol, where relay R 1 selects the best received sig- for SIC. In the first SIC phase, relay R 1 decodes the x N symbol from the best received signal max Y by treating the data symbols x j = {x 1 , . . . , x N−1 } and AWGN n R 1 as interference. The SINR is therefore obtained when relay R 1 decodes the x i symbol, as follows: where ρ S is the transmission Signal-to-Noise Ratio (SNR) and ρ 0 = P S /N 0 . Maximization of the instantaneous bit rate threshold is achieved at relay R 1 when relay R 1 decodes the message x i for i = {N, . . . , 1} as follows: In the second time slot T (odd) 2 in the odd transmission block T (odd) , the relay R 1 retrieves the messages x i = {x N , . . . , x 1 } and forwards the messages to the devices in the superimposed signal while the relay R 2 continues harvesting energy from the BS. Therefore, the received signal at devices D n for n = {1, . . . , N} and the EH at relay R 2 are expressed, respectively, as follows: where P R 1 is the transmission power at relay R 1 and n n ∼ CN(0, N 0 ) is the AWGN at device D n for n = {1, . . . , N} which follows zero mean and variance N 0 . Note that the pre-coding channel matrix H n is given by where the channel h (a 1 ,a 2 ) n in the pre-coding channel matrix H n , where a 1 ∈ A 1 and a 2 ∈ A 2 , is a channel from transmitting antenna a 1 at relay R 1 to a receiving antenna a 2 at a device D n , also applying Rayleigh distribution for propagation. Each fading channel gain is given by h (.,.) where d n is the distance from relay R 1 to device D n . As a result of the combination of TAS and SIC, the SINRs are obtained at devices D n for n = {1, . . . , N} when the devices decode data symbols x i for i = {N, . . . , n} from the best received signal max{|Y 2 |} A 1 ×A 2 by treating the data symbols x j for j = {1, . . . , i − 1} and AWGN n n as interference: = max where SNR ρ R 1 = P R 1 /N 0 . The achievable bit-rate reached at device D n when it decodes the data symbol x i for i = {N, . . . , n} from the best received signal max{|Y n |} is expressed as follows:

Even Transmission Block
As with the first time slot in the odd transmission block, by applying the PS protocol, the first time slot in the even transmission block T (even) 1 is used by the BS to transmit wireless energy and superimposed information simultaneously to coupled relays. The EH from the best channel in the pre-coding channel matrix at coupled relays R 1 and R 2 is expressed as follows: where T (odd) 1 = T (odd) 2. The BS broadcasts a superimposed signal by combining the independent messages x i of devices D i for i = {N, . . . , 1}. Therefore, the received signal at coupled relays is expressed as follows: max In the second transmission block, SIC is implemented at the relay R 2 . As with the SIC in the odd transmission block, the relay R 2 has to repeat SIC until it detect all data symbols x i for i = {N, . . . , 1} in the best received signal as (20). The SINRs are obtained when relay R 2 decodes the data symbols x i = {x N , . . . , x 1 } as follows: Similar to (8) and (8), maximization of the instantaneous bit-rate threshold achieved at the relay R 2 when the relay R 2 decodes data symbols x i , where i = {N, . . . , 1} is expressed as follows: By applying a DF protocol, relay R 2 recovers the decoded data symbols x i for i = {N, . . . , 1} and forwards a beamforming superimposed signal to devices D n for n = {1, . . . , N}. The received signals in the second time slot T (even) 2 in even transmission block T (even) at devices D n while EH at relay R 1 are expressed, respectively, as follows: where P R 2 is the transmission power at relay R 2 and the pre-coding channel matrix G n is given by where the channel g (a 1 ,a 2 ) n in the pre-coding matrix channel G n , with a 1 ∈ A 1 and a 2 ∈ A 2 , is a channel from transmitting antenna a 1 at relay R 2 to a receiving antenna a 2 at a device D n , applying Rayleigh distributions for propagation. Each fading channel is represented by g (.,.) where v n is the distance from relay R 2 to device D n . The SINRs obtained at devices D n for n = {1, . . . , N} when they decode the data symbols x i for i = {N, . . . , n} from the best received signal max n are expressed as follows: where SNR ρ R 2 = P R 2 N 0 . Maximization of the achievable bit-rate thresholds achieved at devices D n when they decode the data symbols x i for i = {N, . . . , n} from the best received signal max n is expressed as follows:

System Performance Analysis
Many factors affect the system performance of wireless networks. The authors of [34] investigated the causes of OP at the BS, such as brief BS power supply variation, preventive BS activity state transition due to excessive temperature increase or decrease inside the BS rack, auto-recovery software and hardware failure, and temporal cell interference and congestion. The downlink MIMO-NOMA model with superposition transmission at BS and SIC at the terminal devices and the SIC processing is adopted in receiver side [35]. Therefore, our study considered OP at the receivers when the receivers could not successfully decode messages in the received signals. We analyze the system performance of the network model depicted in Figure 1 and delivers novel closed-forms of OP, system throughput and EE expressions at coupled relays R = {R 1 , R 2 } and devices D n for n = {1, . . . , N}.
3.1. Outage Probability at the Coupled Relays R Theorem 1. The outage event at a relay in the coupled relays R occurs when the relay cannot successfully decode at least a data symbol x i ∈ {x N , . . . , x 1 } from the best received signal for an odd transmission block or max for an even transmission block, which is the best signal after TAS. Therefore, the coupled relays R receive the best superimposed signal max by applying a TAS protocol to select the best signal from the pre-coding channel matrix H 0 for an odd transmission block or G 0 for an even transmission block to maximize the SINRs γ , as given by (8), (9), (21), and (21), and maximize the instantaneous bit-rate threshold max    , as given by (10) or (23). Maximization of the instantaneous bit rate threshold max    is less than a device's predefined bit-rate threshold R * i , an outage event will occur, i.e, the OP at coupled relays R is expressed as follows: We obtain the OP at coupled relays R = {R 1 , R 2 } in novel closed-form as follows: where β i is given by See Appendix A for proof.

Remark 1.
If the network has a large number of devices, it is challenge to apply the expression (30) in Monte Carlo simulations. Fortunately, the authors of [12] analyzed a network model with multiple relays and multiple devices. The authors presented the expressions for OP at the relays as [12] (Equation (33)). It is important to mention that the present paper extends the work of the previous study [12] by deploying massive MIMO and TAS techniques. From [12] (Equation (33)), the OP expression at the coupled relay R in (30) can be rewritten as follows: From (34), we obtain a novel expression of OP at the coupled relay R in closed-form by applying the Cumulative Density Function (CDF) as defined in [13] (Equation (71)): See Appendix B for proof.

Outage Probability at Devices
Theorem 2. The outage event in an odd or even transmission block at devices D n for n = {1, . . . , N} occurs when, on the one hand, the relays R 1 and R 2 for odd and even transmission blocks, respectively, cannot successfully decode at least data symbol x i for i = {N, . . . , n} from the best received signal max Therefore, the OP at devices in an odd or even transmission block is expressed as follows: From (36), the OP at user D n is obtained in closed-form as follows: where σ 2 n = E |H n | 2 . Note that θ = odd and R = R 1 , or θ = even and R = R 2 .
See Appendix C for proof. Similar to (30), note that Expression (36) is challenging in Monte Carlo simulations if the network has a large number of devices. Therefore, the OP at devices D n can be rewritten as follows: From (38), we obtain the expression of OP at devices D n in closed-form as follows: See Appendix D for the proof.

System Throughput
In Figures 1 and 2, the individual system model has two transmission blocks. Each odd or even transmission block is separated into two time slots. The achievable system throughput in an odd or even transmission block is the sum of the minimal device throughput at the relay and device in the same transmission block. Therefore, the system throughput is expressed as follows:

Energy Efficiency
Green wireless networks require a higher throughput yet use lower energy. To achieve this aim, the present study deploys a massive MIMO technique and SWIPT protocol. As a result, the EE performance indicates the sum of device throughput and sum of transmission power at the BS and coupled relay ratio in the same transmission block. Therefore, the EE performance of an individual network model as given by Figure 1 is expressed as follows:

Numerical Results and Discussion
To investigate the system performance of the massive MIMO-NOMA network model shown in Figure 1, we propose parameters for both a theoretical analysis and the Monte Carlo simulations, as shown in Table 1.
Let a wireless network contain coupled relays and three devices (N = 3). The coupled relays are allocated nearby. The distances from the BS to the coupled relays are d R 1 = d R 2 = 10 m, and the distances from the coupled relays to devices D 1 , D 2 and D 3 are d 1 = v 1 = 5 m, d 2 = v 2 = 7 m and d 3 = v 3 = 10 m, respectively. The path-loss exponent refers to an indoor environment with λ = 4. By modeling the more challenging indoor environment, the performance bound of the less challenging outdoor scenario is therefore also covered. In 4G Long-Term Evolution (LTE) release 8, the maximum numbers of antennas at the BS and user equipment are 4T4R and 1T2R, respectively. In 4G LTE-Advanced (LTE-A), the number of antennas is greater to allow the use of 8T8R at BSs and 2T2R at devices in LTE-A schemes. Therefore, massive MIMO networks need at least an 8T8R antenna array at the BS. However, the present study assumes a certain number of antennas equipped at the BS (A 0 = 4), coupled relays R 1 (A 1 = 4), R 2 (A 1 = 16), and devices (A 2 = 2) to prove the the benefits of massive MIMO (even transmission block) over MIMO (odd transmission block). The fading channel from the transmitting antennas at the BS to the receiving antennas at the coupled relays are distributed over Rayleigh fading channels. Based on distances and the path-loss exponent, the expected channel gains from the BS to coupled relays are σ 2 R 1 = σ 2 R 2 = 1 × 10 −4 and from the coupled relays to devices D 1 , D 2 , and D 3 are σ 2 1 = 16 × 10 −4 , σ 2 2 = 4.1649 × 10 −4 , and σ 2 3 = 1 × 10 −4 . Each fading channel randomly generates 1e6 experiments. To simplify, the devices require the same bit-rate thresholds 1bps/Hz and SNR ρ S = ρ R = {0, . . . , 30} dB. By applying (7), the PA factors for devices D 1 , D 2 , and D 3 are α 1 = 0.1667, α 2 = 0.3333, and α 3 = 0.5, respectively. The present study assumes that coupled relays may fully collect EH (η = 1). The PS factor in an odd transmission block is λ = 0.4. Therefore, 0.4P S and 0.6P S are applied to EH and DT processing, respectively, whereas the PS factor in an even transmission block is reduced (λ = 0.4) since the relay R 2 in this block is equipped with a greater number of antennas than relay R 1 in an odd transmission block.

Variables Values Units
Even transmission block A 1 16 λ 0.6 Figure 3 plots the OP performance at relay R 1 and devices in odd transmission blocks. Note the various markers and line plot analysis (Ana) and simulation (Sim) results. The analysis results of OP performance at relay R 1 are given by (30) or (34) and for devices D n by (36) or (38). The analysis results were verified with Monte Carlo simulations for relay R 1 given by (31) or (35) and for devices by (37) or (39), where R = R 1 and θ = odd. Figure 3 illustrates that device D 3 achieved the best OP results, even though device D 3 was the farthest device and therefore allocated the biggest PA factor α 3 = 0.5. When SRN ρ → ∞, the OP results of relay R 1 and devices tend toward zero. Figure 4 plots the OP performance at relay R 2 and devices in even transmission blocks. To improve the networking capacity and energy, we equipped a large number of antennas at relay R 2 and increased the PS factor λ = 0.6. By increasing the PS factor, relay R 1 was able to harvest more energy but relay R 2 could receive weak signals. However, we may observe that the OP performance at relay R 2 and devices in even transmission blocks ( Figure 4) achieved better results than OP performance at relay R 1 and devices in odd transmission blocks (Figure 3). The analysis results of OP performance for relay R 2 are given by (30) or (34) and for devices D n by (36) or (38), where R = R 2 and θ = even. The analysis results were verified with Monte Carlo simulations for relay R 2 given by (31) or (35) and for devices by (37) or (39), where R = R 2 and θ = even.
We may observe that the OP performances at relay R 2 and devices in even transmission blocks outperform those at relay R 1 and devices in odd transmission blocks at high SNRs such as ρ S = ρ R = 20 dB. The present study thus exploits the advantages of a massive MIMO-NOMA network compared to a MIMO-NOMA network. We equipped a greater number of antennas on relay R 2 (A 1 = 16) than on relay R 1 (A 1 = 4). As a result, we obtained the respective pre-coding channel matrix sizes of [4 × 16] and [16 × 2] for G 0 and G n , which, in the even transmission blocks, were much larger than the pre-coding channel matrix sizes of [4 × 4] and [4 × 2] for H 0 and H n in odd transmission blocks. From Expressions (8), (9), (14), (15), (21), (22), (27), and (28) and by applying the TAS protocol, only the best channels from the pre-coding channel matrices are selected for data decoding. Therefore, the relay R 2 and devices in even transmission blocks have better OP performance than relay R 1 and devices in odd transmission blocks under the same simulation parameters.  Figures 5 and 6 plot the system throughput performance at devices in odd and even transmission blocks, respectively. Even though device D 3 has the farthest distance from coupled relays (d 3 = v 3 = 10 m), device D 3 always achieved the best throughput performance compared to other devices. It is interesting that the system throughput results of relay R 1 and devices in odd transmission blocks ( Figure 5) are similar to the system throughput results of relay R 2 and devices in even transmission blocks ( Figure 6) at the same SNR.
To illustrate, we extract the investigated results from Matlab software. At a SNR range ρ S = ρ R = {15, . . . , 21} dB, device throughput in odd and even transmission blocks achieved TP (odd) = {0.0051, 0.0435, 0.1338, 0.2204, 0.279, 0.2983, 0.3} and TP (even) = {0.0011, 0.0261, 0.1287, 0.2169, 0.2849, 0.2999, 0.3}, respectively. We may observe that device throughput in odd transmission blocks outperformed device throughput in even transmission blocks at low SNRs, i.e., ρ S = ρ R = {15, . . . , 18} dB. However, device throughput in even transmission blocks improved and outperformed device throughput in odd transmission blocks at high SNRs, i.e., ρ S = ρ R = {19, 20} dB. Device throughput also tended toward their data rate thresholds, i.e., R * 1 = R * 2 = R * 3 = 0.1 b/s/Hz when the SNRs ρ S = ρ R → ∞. As a result, system throughput in both odd and even transmission blocks It is important to note that the PS factor in even transmission blocks (λ = 0.6) was greater than the PS factor in odd transmission blocks (λ = 0.4). Therefore, the relay R 1 in even transmission blocks harvested more energy than relay R 2 in odd transmission blocks. Certainly, even transmission blocks achieved better EE performance than odd transmission blocks. Figure 7 plot the EE performance of odd even transmission blocks. We may observe that EE performance in the even transmission block with massive MIMO had a higher peak than the odd transmission block. The massive MIMO technique therefore not only provided greater throughput but also consumed less energy.

Conclusions
This paper proposes a design for a switchable coupled relay model to assist a massive MIMO-NOMA wireless network in serving multiple devices and extending a network's lifespan. A diagram of two transmission blocks illustrates signal propagation and EH processing. Propagation and formulations are analyzed. We derive the closed-form novel expressions for OP at the coupled relays and devices. The theoretical results show that the massive MIMO technique in combination with TAS and SWIPT protocols in an underlay cooperative NOMA network provides higher throughput and consumes lower energy. The obtained results verify the massive MIMO technique as effective for 5G wireless networks. The present paper offers the potential for the practical application of a massive MIMO-NOMA network model assisted by switchable coupled relays, for example in a water environment where relays and devices are barely powered. Our massive MIMO-NOMA assisted by switchable coupled relays in combination with TAS and EH protocols can not only improve OP and system throughput performance but also extend the network's lifespan. Specifically, EH at the relay may be used to forward signals without consuming the relay's own energy. This is promising as a potential solution in extending a network's lifespan.
Author Contributions: These authors contributed equally to this work. All authors have read and agreed to the published version of the manuscript. Data Availability Statement: Data was randomly generated when we were running program. There is no data to share.

Acknowledgments:
We are especially grateful for our reviewers' comments and all the constructive suggestions that have helped us improve the quality and presentation of our manuscript.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the study's design, the collection, analysis and interpretation of data, writing of the manuscript, nor in the decision to publish the results.

Abbreviations
The following abbreviations are used in this manuscript:

Notations
The following notations are used in this manuscript: Distances from BS to relays R R = {R 1 , R 2 } Set of relays D n n ∈ N Devices d n n ∈ N Distances from relay R 1 to devices v n n ∈ N Distances from relay R 2 to devices Data of devices n R , n n , n ∈ N AWGN at relays and devices After some algebraic manipulation, we obtain From PDF expressions [23] (Equations (59) and (60)), we obtain By substituting (A3) into (30), we obtain the expression of OP at coupled relays R, as shown in (31).

Appendix C. Proof of Theorem 2
From (36), we obtain Expression (A8) can be solved and obtained in closed-form, as shown in (37).