Energy-Spectral Efficiency Trade-Offs in Full-Duplex MU-MIMO Cloud-RANs with SWIPT

The present paper investigates the trade-offs between the energy efficiency (EE) and spectral efficiency (SE) in the full-duplex (FD) multiuser multi-input multioutput (MU-MIMO) cloud radio access networks (CRANs) with simultaneous wireless information and power transfer (SWIPT). In the considered network, the central unit (CU) intends to concurrently not only transfer both energy and information toward downlink (DL) users using power splitting structures but also receive signals from uplink (UL) users. This communication is executed via FD radio units (RUs) which are distributed nearby users and connected to the CU through limited capacity fronthaul (FH) links. In order to unveil interesting trade-offs between the EE and SE metrics, we first introduce three conventional single-objective optimization problems (SOOPs) including (i) system sum rate maximization, (ii) total power minimization, and (iii) fractional energy efficiency maximization. Then, by making use of the multiobjective optimization (MOO) framework, the MOO problem (MOOP) with the objective vector of the achievable rate and power consumption is addressed. All considered problems are nonconvex with respect to designing variables comprising precoding matrices, compression matrices, and DL power splitting factors; thus, it is extremely intractable to solve these problems directly. To overcome these issues, we develop iterative algorithms by utilizing the sequential convex approximation (SCA) approach for the first two SOO problems and the SCA-based Dinkelbach method for the fractional EE problem. Regarding the MOOP, we first rewrite it as an SOOP by applying the modified weighted Tchebycheff method and, then, propose the iterative algorithmbased SCA to find its optimal Pareto set. Various numerical simulations are conducted to study the system performance and appealing EE-SE trade-offs in the considered system.


Introduction
In the context of explosive growth in mobile data traffic and in wireless connectivities, the cloud radio access networks (CRANs) have been suggested as a modern architecture for the next generation of wireless communication systems. It is promising to overcome the strict requirements for concurrently enhancing the 1000-fold data rate as well as reducing by half the power consumption [1][2][3]. Specifically, a typical CRAN composes of a central unit (CU), a set of distributed radio units (RUs), and terminal devices [1]. The CU is responsible for cooperatively processing baseband signals of all terminals and, then, communicates them to RUs through FH connections. These RUs which are inherently distributed nearby user equipment (UE) are responsible for transmitting or/and receiving data messages toward or/and from terminal devices. Based on this structure, the baseband signals are centrally processed at the CU, and thus, CRANs provide not only a cooperative data transmission but also a better interference management among RUs as compared with the conventional ones, for example, multicell models [4,5] and coordinated multipoint (CoMP) schemes [6,7]. In addition, another remarkable advantage of CRANs is the replacement of the high-cost high-consuming power base stations by the lowcost low-consuming power RUs that can significantly reduce the system power consumption.
On the other hand, with potentially increasing the degrees of freedom in the design, the multiple antenna transmission techniques can thus support to convey more data streams over a same radio resource and promisingly boost the system throughput. Hence, the multiantenna techniques have been integrated and investigated in the context of CRANs. The CRANs implemented with multiple-antenna RUs and singleantenna terminals are named as MISO CRANs and investigated in [8][9][10][11][12][13]. Meanwhile, the ones deployed with multipleantenna RUs and multiple-antenna UEs were referred to as MIMO CRANs and studied in [14][15][16][17]. Regarding the MISO CRAN deployment, [8] proposed the joint design of RU activation and beamformers to minimize power consumption in downlink (DL) MISO CRAN systems. In [9], the authors first proposed a spatial-compression-and-forward strategy for uplink (UL) MISO CRAN systems and, then, addressed the beamforming and compressing design to maximize the minimum user's quality of service (QoS). Alternatively, the work in [10] developed a nonorthogonal multiple access (NOMA) transmission in DL MISO CRANs and, then, proposed a solution to handle the maximin individual user rate subjected to transmit power budget and limited FH capacity constraints. In addition, a wireless-powered MISO CRAN system was developed in [11] where the terminal devices conduct both information decoding (ID) and energy harvesting (EH) in the first ID-EH DL phase and each user then transmits its own data in the UL data transmission phase using harvested energy. The UL rate maximization problem is tackled under minimum DL rate requirements, power budget, and FH capacity constraints. Several novel techniques including hybrid digital-analog beamforming and rate splitting multiple access were addressed in DL MISO CRANs in [12,13], respectively. As the result, the system spectral efficiency (SE) is thus considerably improved. On the other hand, by considering the MIMO CRAN scheme, the joint design of precoding and compressing matrices was investigated for the UL MIMO CRAN in [14,15] and for the DL MIMO CRAN in [16], respectively. Additionally, the authors in [17] took multicluster interference into consideration and developed an effective algorithm to handle them in the multicluster DL MIMO CRAN.
It is worth noting that the aforementioned works separately studied the UL CRAN channels or DL CRAN ones and can be referred to as half-duplex (HD) transmission that leads to inefficient spectrum usage. Recently, thanks to the advanced analog and digital self-interference (SI) cancellation strategies, the full-duplex (FD) radio is able to concurrently transmit and receive signals over the same radio resources. Then, the FD techniques potentially provides the SE approximately twice as much as that of the conventional HD counterparts [18,19]. The integration between the FD radio approach and CRAN scheme has recently drawn significant attention [20][21][22][23][24][25][26]. Reference [20] investigated a classical Wyner model based on FD-aided CRANs and illustrated that the achievable rate of the FD CRAN system is significantly improved in comparison to both HD CRANs and single-cell processing ones. Reference [21] developed a two-stage algorithm to jointly design RU beamformers, user power allocation, and user selection for the FD MISO CRAN. Furthermore, the system data rate for the multiuser scheme was analysed and optimal power allocation for the single-user case was investigated in the FD CRANs [22]. Additionally, by considering a hybrid RU set which comprises both multiantenna and singleantenna RUs, [23] aimed for joint resource allocation and interference cancellation in order to maximize the DL rate subjected to power and UL rate QoS constraints in the FD CRANs. By exploiting the difference of the convex (DC) form in the system data rate and FH capacities, the authors in [24] proposed an iterative algorithm based on the majorizationminimization method for the FD MIMO CRANs which jointly optimizes precoding and compressing matrices under transmit power and FH rate constraints. Alternatively, the work in [25] developed a low-complexity design for FD MIMO CRANs by making the use of the connection between weighted minimum mean squared error (WMMSE) problem and weighted sum rate (WSR) problem. The precoders and compressing matrices are then alternatively optimized by solving the convex problem in each iteration. Most recently, in [26], the simultaneous wireless information and power transfer (SWIPT) protocol was investigated on the FD MISO CRANs. The joint design of DL beamforming, EH-power splitting ratio, UL power allocation, and RU postcoding was addressed. Reference [27] investigated the energy-throughput trade-off by designing the beamforming for SWIPT CRANs. Alternatively, the authors in [28] developed the hybrid time-power splitting protocol for downlink user information and uplink ID-EH sensor in FD-aided massive MIMO wireless systems.
However, the aforementioned works were concentrated only on designing either the sum rate maximization (referred to as SE maximization) or the power consumption minimization in the considered CRAN. This might lead to ineffective and negative impacts on the system EE performance [29]. The EE issue was, therefore, taken into account in the transceiver design for CRAN networks [30][31][32][33][34][35][36]. Specifically, reference [30] investigated the resource allocation in order to minimize the entire system power consumption for both data sharing-and data compressing-based DL CRAN models. The authors in [31] first generalized these models to the multihop FH links MIMO CRANs and, then, proposed a successive convex quadratic programming (SCQP) algorithm to jointly design precoders, user-RU association, and RU activation to optimize the system fractional EE defined as a ratio of the system sum rate to power consumption. In order to further reduce power consumed by ineffective RUs, the issues of turning RU on/off and user-RU association were studied in various recent works. Particularly the joint design of precoding matrices, turning RUs on/off was investigated for both SE and EE optimization problems in DL MIMO CRAN-based multihop FH link networks [32]. In [33], the EE maximization problem constrained with FH capacity limitation, transmit power budget, and user's QoS requirements was tackled for a DL MISO CRAN. Alternatively, the EE beamforming design-based RU selecting and RU-UE associating designs for DL MISO CRAN were proposed either with virtual computing resource models in [34] or with FH rate-dependent power consumption models in [35]. To take advantage of FD techniques, reference [36] proposed an efficient algorithm for optimizing the system EE in FD MISO CRANbased data sharing systems. However, the FH capacity constraints and FH power consumption model were not taken into consideration.

Wireless Communications and Mobile Computing
We would like to emphasize that, in the aforementioned works, only one specific system performance metric is considered as a single-objective function in design problems while the others are taken into account in the constraints. This approach can be referred to as single-objective optimization (SOO) design and can lead to an unfairness between the conflicting system measures, such as systems SE, EE and power consumption. More recently, the multiobjective optimization (MOO) framework has been introduced as an effective mathematical tool to overcome this issue as well as to provide a balance among the system performance metrics on designing the wireless communication systems [37,38]. As a result, the MOO framework was recently exploited in many studies to address various trade-offs between two or more system metrics [39][40][41][42][43][44][45][46][47][48]. Particularly, the harvested energy and information transfer trade-offs in SWIPT-enabled wireless system were studied either in cognitive FD MU-MIMO systems [39] or in MIMO broadcast channels [40,41]. In [42], the authors studied the optimal Pareto set for DL and UL transmit power fairness in the FD MISO networks. Furthermore, the balance between transmit power and achievable rate was investigated in [43,44] for secure DL MU-MISO networks and DL MISO CRANs, respectively. In addition, the EE-SE trade-offs for massive MIMO systems were addressed in [45,46]. On the other hand, the works in [47,48] have developed algorithms using MOO frameworks to study the balance among three essential and conflicting metrics. Particularly, the balance among three inconsistent design objectives of the MSE minimization, harvested energy maximization, and physical security optimization was revealed for SWIPT-aided MIMO interference channels in [47]. Furthermore, the work in [48] addressed the fairness among total transmit power minimization, EH efficiency maximization, and interference leakage-totransmit power ratio minimization for SWIPT-based MISO cognitive radio networks.
Motivated by the above viewpoints, this paper studies EE-SE trade-offs in an FD MU-MIMO CRAN-enabling SWIPT network (shown in Figure 1). The considered CRAN model includes that a set of FD RU stations, which is connected to the CU through the limited capacity FH links, concurrently communicates to the sets of EH-enabled DL users (EH-DLUs) and UL users (ULUs). Specifically, regarding DL channels, the baseband signals associating with DLUs are first centrally processed at the CU and, then, conveyed to the RUs via FH links. After radio frequency (RF) steps, the distributed RUs transmit these signals toward DLUs through wireless access links. At the DLU, the received signal is divided into two proportions using power splitting (PS) receiver structures, which are then performed for ID and EH purposes. The total amount of harvested energy can be stored at the DLUs to prolong their lifetime or to utilize for signal processing operations as discussed in previous works [26,39,47,48]. In the UL channels, the individual-received baseband signals at RUs are sent to the CU via FH links for central processing. It is noted that, in order to guarantee the rate over finite capacity FH links, the baseband signals are quantized and compressed by using FH compression strategies [14][15][16]. Hence, we aim to jointly design precoding matrices, FH compressing covariance matrices, and power splitting parameters. Moreover, to investigate the trade-offs between two conflicting EE and SE metrics, we introduce three SOO problems (SOOPs) and one MOO problem (MOOP); they are all under constraints of transmit power budget, minimum data rate, maximum FH link capacity, and minimum amount of harvested energy. Specifically, the SOOPs include (i) sum rate maximization, (ii) power consumption minimization, and (iii) fractional energy efficiency maximization whereas the MOOP is the simultaneous sum rate and power consumption optimization one. All interesting problems are the challenging nonconvex optimization ones and hard to be directly handled. This is due not only to the nonconcavity in data rate functions and nonconvexity in FH link rates caused by the coupled of UL and DL precoder variables but also to the inverse convexity in EH constraints. To overcome these issues, we exploit the efficient inner convex approximation and SCQP approaches to develop iterative algorithms with low computational complexity.
The main contributions of this work are summarized as twofold as follows: (1) This paper attempts to comprehensively study the SE-EE trade-off in the FD MIMO CRANs with SWIPT. Compared to previous related works in the literatures [24][25][26]36], we highlight the novelty of the present work as follows: (i) Although the works in [24,25]    Wireless Communications and Mobile Computing system model for enabling SWIPT at DLUs and, then, studies the additional system EE metric using two approaches, namely, the maximization of the ratio between the sum rate over the power consumption based on the SOO programming and simultaneous sum rate maximization and power consumption minimization based on MOO programming.
(ii) On the other hand, it is worth noting that the work [26] has also attempted to minimize power beamforming under the QoS of the data rate and harvested energy in SWIPT-enabled FD MISO CRANs while the authors in [36] have also proposed an algorithm to optimize the overall system EE in FD MISO CRANs. However, both works assumed single-antenna users and unlimitedcapacity FH links, therefore, those approaches cannot be directly applied for FD MIMO CRANs with SWIPT considered in this paper.
Thus, the present work is distinguished itself and can be referred to as the more comprehensive investigation (2) We also develop a computationally efficient iterative algorithm based on the SCQP approach, in which the log-determinant expressions in both objective function and FH rate constraints are recast by the quadratic convex ones within each iteration. As the result, the computational cost, which is essentially caused by log-determinant expressions, is significantly reduced as compared with the methods in [24,25].
The rest of this paper is organized as follows. Section 2 introduces the system model and transmission signaling and presents the problem formulations for SOOPs, which include sum rate maximization, power consumption minimization, and energy efficiency maximization, as well as introduces an EE-SE trade-off MOOP. Section 3 develops algorithms for considered SOOPs and MOOP. Finally, the simulation results and conclusions are provided in Section 4 and Section 5, respectively.
Operator Rfxg is the real part of the complex number x and Ef·g denotes the expectation operator.

System Model
This paper considers an FD MIMO CRAN-enabled SWIPT as shown in Figure 1, in which two sets of UL and DL users communicate with a CU through a pool of RUs. It is assumed that the limited-capacity FH links are used in CU-RU con-nections. Hence, to meet the bottleneck rate constraint, the data is compressed before being transferred over FH links. For notation simplicity, the sets of K RUs, I DLUs, and J ULUs are denoted as R ≜ fλ R : λ = 1, 2, ⋯, Kg, D ≜ fk D : k = 1, 2, ⋯, Ig, and U ≜ fℓ U : ℓ = 1, 2, ⋯, Jg, respectively. The number of transmit and receive antennas at RU λ R are N Tx λ R and N Rx λ R , respectively. Thus, the numbers of all transmit and all receive antennas at RUs are N Tx respectively. In addition, the numbers of transmit antennas of ULU ℓ U and receive antennas of DLU k D are denoted as N Tx ℓ U and N Rx k D , respectively. Moreover, the channel state information (CSI) of all links is assumed to be fully known at the CU by exploiting the joint UL and DL channel estimation algorithm [49]. It should be noticed that such an assumption of the perfect CSI is widely adopted in the recent related works [9][10][11][12][13][14][15][16][17][23][24][25][26]. In practice, the CSI can be achieved through the channel estimation, feedback, or reverse links, and thus, the imperfect CSI is available at the CU. In such cases, the results in this work can be referred to as the upper-bound performances.
At the same time and over the same frequency resource, the CU intends to transmit the message s k D ∈ ℂ d k D ×1 for DLU k D in the DL channels while the ULU also wishes to convey its own message s ℓ U ∈ ℂ d ℓ U ×1 to the CU, with d k D and d ℓ U representing the numbers of data streams of DLU k D and ULU ℓ U , respectively. All data symbols are assumed with zero mean and unit power-independent distribution, i.e., Efs i s † j g = 0 with i ≠ j and Efs i s † i g = I with an appropriate dimension for all i, j ∈ fU, Dg. In addition, to take advantages of multiantenna transmission in spatial diversity enhancement, the DL and UL signals are first beamformed by multiplying with its corresponding precoding matrices, Thus, in the UL channels, the transmit signal at ULU ℓ U is given as Then, the power constraints at ULUs, i.e., Ekx ℓ U k 2 ≤ P ℓ U , ∀ℓ U ∈ U, are determined by Concerning the DL channels, it is emphasized that the precoding matrix for DLU k D is constructed at the CU in Prior to being conveyed throughout the limited capacity FH link, the DL signal for RU λ R is quantized and compressed. Adopting a Gaussian test channel, the compressed DL signal for RU λ R is given as where the quantization noise of DL FH link λ R , q D,λ R , is independent of precoded data stream and obeys q D,λ R~C N ð0, Q D,λ R Þ as in [16,17,24,25]. This process is to guarantee the rate constraint over the FH link. Hence, the power constraint at RUs, i.e., Ekx D,λ R k 2 ≤ P λ R with λ R ∈ R, is equivalently expressed as 2.1. DL Channels with Information and Energy Transferring.
From (3), the noiseless FH rate of the link from CU to RU λ R is determined as [16,17] g D,λ R = log 2 I N Tx In addition, the received signal vector at DLU k D , , is given as represents the additive Gaussian noise vector at DLU k D with zero mean and covariance matrix σ 2 k D ,a I N Rx k D .In order to realize SWIPT via received RF signals, this paper adopts power splitting receiver structures at DLUs. Specially, we exploit perantenna power splitting (PAPS) structures in EH-DLUs. Denoting 0 ≤ w k D ,j ≤ 1 as the signal proportion factor for ID at antenna j of DLU k D with j = 1, 2, ⋯, N Rx k D , the PAPS coefficients of DLU k D are represented by a diagonal matrix Then, the received signal for ID is represents the additive circuitry Gaussian noise vector at DLU k D with zero mean and covariance matrix σ 2 k D ,c I N Rx k D . Accordingly, the residual amount of signals is used for scavenging energy. Therefore, the signal for the EH purpose is given as From (7), the DL rate at DLU k D can be expressed as the function of precoders, compression matrices, and PAPS matrices and given by where which is the covariance matrix of multiuser, cochannel, and compression interferences plus additive noise at DLU k D , where we have defined Q D ≜ diag fQ D,λ R , λ R ∈ Rg.Moreover, upon (8), the amount of harvested energy at DLU k D is expressed as follows where the DL signal covariance mapping matrix Π k D = C k D + Ω k D . In this paper, we consider a practical nonlinear model of energy harvester which is mathematically formulated as [50] where k D denotes the maximum power harvested at DLU k D at the saturation point of EH circuit and a k D and b k D are the nonlinear EH parameters by hardware limitations as illustrated in [50]. This quantity of energy, E NL k D , can be stored either in a battery or in a super capacitor and later exploited to prolong its own operation such as processing/decoding signals or transmitting control signals to the RUs [26,39,47,48,50,51].
2.2. UL Channels with Information Transferring. As the transmission strategy discussed in the early subsection, the 5 Wireless Communications and Mobile Computing received baseband signal at RU λ R in UL channels is given as where interference is referred to as SI at RU λ R ; otherwise, it is an inter-RU interference link between receive antennas at RU λ R and transmit antennas at RU θ R . The additive Gaussian noise at RU λ R is denoted as Þ is the quantization noise which is randomized and independent from the received data stream [14,15,24,25]. The rate corresponding to UL FH link λ R is thus expressed as where Define R represents the aggregative link from transmit to receive antennas between RUs. In addi- are the stacked versions of UL, DL quantization noise, and additive noise at RUs. Then, the stacked UL baseband signal vector received at all RUs, i.e., , can be expressed as Similar to [24,52], the optimal minimum mean square error and successive interference cancellation (MMSE-SIC) technique is used at the CU to decode the UL messages. Based on the MMSE-SIC decoding strategy, the user's message with the strongest channel condition is retrieved first whereas that with the worst channel link is decoded last. After detecting any user data, the CU then removes it from the received signal for interference mitigation. As the result, the later decoding ULU suppresses from the multiple-access interferences of all earlier decoding ULUs. Therefore, the data rate of ULU ℓ U is given as [14,24] where Then, the total achievable rate of UL channels, i.e., ∑ ℓ U ∈U R ℓ U , can be formulated as where is the received covariance matrix of the desired UL signal at RUs. Meanwhile, the covariance matrix of SI and compression noise plus additive noise at RUs is with Q U ≜ diag fQ U,λ R , λ R ∈ Rg being the diagonally stacking matrix of covariance compression noise of UL FH links.

Power Consumption Model. The whole transmit power consumed by RUs and ULUs is given by
where P Tx ℓ U and P Tx λ R are the transmitted power at ULU ℓ U and RU λ R and given in (1) and (4), respectively. In addition, 0 < ξ λ R , ξ ℓ U ≤ 1 are the power amplifier efficiency factors of RU λ R and ULU ℓ U , respectively. Otherwise, the total circuitry power is given as where P Cir λ R and P Cir ℓ U are the circuitry power consumed at RU λ R and ULU ℓ U , respectively. Additionally, the power consumption of the FH links is defined as the rate-dependent model; hence, the total power consumption over all DL and UL FH links is given as [30,31] where α FH U,λ R = P FH U,λ R ,max C −1 U,λ R and α FH D,λ R = P FH D,λ R ,max C −1 D,λ R are constant scaling factors. g D,λ R and g U,λ R are FH rates given by (5) and (14), respectively. Meanwhile, C D,λ R and C U,λ R are the maximum rate to be able to convey through the DL and UL FH link λ R , respectively. Nevertheless, P FH U,λ R ,max and

Formula of Considered Optimization
Problems. Firstly, we introduce some system requirements including minimum harvested energy, minimum user's data rate, and limited FH capacity requirements. In order to ensure the operating necessity of the energy harvester, the total amount of harvested energy at DLUs must satisfy the following condition where e req k D is the required harvested energy at DLU k D . The value of e req k D is chosen properly providing that the harvested energy is useful and can activate the EH circuitry [50]. Hence, this value should be set up to be larger than or at least equal to the minimum sensitive level of its energy harvester [26,39,48].
On the other hand, to meet the maximum limited FH capacity, the DL and UL FH rates are thus constrained as [24,25] Finally, the QoS management is predefined as the following achievable data rate constraints: where r req k D and r req R denote the minimum QoS requirement of DL and UL data rates, respectively.
For notation simplicity, we define the compact set S as In this paper, we consider three SOOPs and an MOOP which are defined below.
2.4.1. Problem SOOP-1. The sum rate maximization problem is written as a minimization one where E Tot ðSÞ ≜ R Tot ðSÞ/P Tot ðSÞ is the whole system EE.
Note that (1) and (4) present the transmit power constraints at ULUs and RUs, respectively. The nonconvex constraints in (23) guarantee the minimum harvested energy at each EH-DLU while the nonconvex constraints of DL and UL FH rate are given in (24) and (25), respectively. As a result, the three above problems are highly nonconvex due not only to nonconvexities in objective functions but also to nonconvex constraints. Thus, these problems are all extremely intractable to directly find their own optimum points. In the following sections, we develop the iterative algorithms using sequential convex approximation (SCA) methods to solve these nonconvex problems.

Problem MOOP.
The energy-spectral efficiency optimization using the MOO framework is given as [37,38]  in which the objective function is a two-row vector whose elements are the SOOP-1 and SOOP-2 objective functions. Once again, it is very challenging to jointly find the optimal solution to these nonconvex multicriteria problems (32a) and (32b). Therefore, in the next section, we introduce an approach based on the modified weighted Tchebycheff and SCA methods to find the optimal Pareto set.

Proposed Iterative Algorithm
This section represents an approach to develop algorithms for considered nonconvex problems using SCA methods. Firstly, we provide three inner convex approximations of the data rate, FH rate, and harvested energy.

Convex Inner Approximation Analysis.
In this subsection, we provide inner approximations in order to convexify the nonconvex problems.
3.1.1. Concave Quadratic Lower Approximation for the DL and UL Rate. Following Theorem 1 of [52], [53], and Appendix A that the DLU k D rate in (9) is lower bounded by the following concave quadratic expression, 7 Wireless Communications and Mobile Computing where Φ Moreover, the concave quadratic lower bound of UL rate (18) is given as where It is shown in [52,53] Accordingly, the nonconvex QoS rate constraints (26) and (27) are innerly approximated, respectively, as the following convex constraints where g ðκÞ D,λ R = log 2 jQ ðκÞ,−1 In addition, the UL FH rate is convex upper bounded bŷ where g ðκÞ U,λ R = log 2 jI N Rx Finally, the inner convex approximation of nonconvex FH rate constraints according to the DL and UL in (24) and (25), respectively, is given aŝ

Convex Inner Approximation for DLU EH Constraints.
Firstly, the compact constraint in (23) is equivalently rewritten after some simple mathematical manipulations as where e introducing additional diagonal matrix variables X k D with k D ∈ D, which satisfies the constraints the harvested energy at DLU k D is lower, approximated as (see Appendix C for detail) with Wireless Consequently, the nonconvex constraint (43) is recast as the following convex inner approximation beloŵ

Iterative Algorithm for Considered Optimization
Problems. Upon the above analyses, the nonconvex SOOP-1 in (29a) is innerly approximated at iteration κ as the following convex one whereP Tot ðSÞ ≥ P Tot ðSÞ and given bỹ whereĝ U,λ R andĝ D,λ R are given in (39) and (41), respectively. These two recast problems are now convex ones and can be efficiently tackled with existing convex solvers, such as CVX packages [54]. The whole procedure is summarized in Algorithm 1.
On the other hand, in order to find the Pareto set, the MOO problems (32a) and (32b) are first transformed to the SOO one by utilizing the modified weighted Tchebycheff method as [37] min S max p=1,2 where F ⋆ 1 and F ⋆ 2 are the optimums of SOOP-1 and SOOP-2, respectively. The weighted factor ω p ∈ ½0, 1 with ω 1 + ω 2 = 1 represents the trade-off between two considered metrics. In addition, the small constant δ = ½0:001,0:01 should be chosen to obtain the Pareto optimality [19,37].
It is important to notice that problems (53a) and (53b) are still challenging to be directly solved due to the nonsmooth nonconvex property of the objective function and nonconvex constraints. Therefore, we exploit the analyses in the previous subsections and SCA method to develop an algorithm to iteratively handle it. Following inner convex approximations in the previous subsection, problems (53a) and (53b) are approximated at iteration κ as For given ω 1 and ω 2 , problems (54a) and (54b) are a convex one and can be handled with existing convex solvers, e.g., Initialization: Setting the small error tolerance ε > 0, the maximum number of iteration κ max . Generating an initial setS ð0Þ by solving max-min problem (58a) and (58b); Set κ = 0; repeat Find optimal setS ⋆ by solving convex problem (50a) and (50b) [or (51a) and (51b)], Algorithm 1: The iterative algorithm for the sum-SE maximization (or sum-power consumption minimization) problem. 9 Wireless Communications and Mobile Computing CVX packages [54]. The entire set of optimal Pareto points is achieved by varying value 0 ≤ ω 1 , ω 2 ≤ 1 under the condition of ω 1 + ω 2 = 1. The total procedure to find Pareto points is summarized in Algorithm 2.

Dinkelbach-Based
ðSÞ being already defined below as (50a), (50b), (51a), and (51b), respectively. Problems (55a) and (55b) are now a concave-convex fractional program with convex constraints. In specific, the objective is a pseudoconcave function that inspires us to exploit the nonlinear fractional programming approach [55]. Based on the Dinkelbach algorithm, we first reformulate (55a) and (55b) into a parameterized polynomial subtractive form with respect to υ as follows Following [56], it is emphasized on parametrized problem (56) that GðυÞ is a monotonically decreasing and continuous function with respect to υ. This implies that GðυÞ = 0 is fulfilled with a unique solution which is given as On the other hand, the optimal solution of (56) at υ ⋆ is also an optimal solution of (55a) and (55b) and υ ⋆ is the optimal objective value. The detailed explanation can be refereed in [56]. It implies that one can first fix υ to find the optimal solution of convex problems (55a) and (55b). Afterward, we tackle identity GðυÞ = 0 to get the unique υ ⋆ . Finally, the whole procedure to achieve the optimal solution is summarized in Algorithm 3. and ω 2 = 1 − ω 1 ; Generating initial set ofS ð0Þ by solving max-min problem (58a) and (58b); Find the optimum F ⋆ 1 and F ⋆ 2 by solving (29a), (29b), (30a), and (30b), respectively, with Algorithm 1; repeat Update ρ ← ρ + 1; Set κ = 0; repeat Given ω ðρÞ 1 , ω ðρÞ 2 , find the optimal setS ⋆ by solving convex problem (54a) and (54b); Update κ ← κ + 1 andS ðκÞ ←S ⋆ ; Compute Δη = 1 − η ðκÞ /η ðκ−1Þ with η is objective value in (54a); until reaches jΔηj ≤ ε or κ = κ max ; Update the ρ-th component of the Pareto set: S opt ðρÞ ← S ⋆ ; until reaches ρ = ρ max ; Output: The optimal Pareto set fS opt ðρÞ : ρ = 1, 2, ⋯, ρ max g.   At the beginning of the proposed algorithms, it is hard to find one feasible set. Thus, we consider the following max-min problem which helps to find an initial feasible set. The initial feasible point is generated by solving the following max-min problem The solution can be iteratively found by solving a inner convex approximation problem at iteration κ as whereR k D ,R R , andÊ k D are given in (33), (35), and (45), respectively. The algorithm terminates as the objective value is larger than or equal to 1.

Numerical Results
In this section, various numerical simulation results are conducted to illustrate EE-SE trade-offs for the interested FD MU-MIMO CRANs with SWIPT.

Simulation Parameters.
We adopt the CRAN model in [24,25], in which a square cell with a length of 100 meters is divided into four equally small square cells with a length of 50 meters. One RU is located at the center of each small

13
Wireless Communications and Mobile Computing cell and a pair of a DLU and a ULU is in each cell. Specifically, the DLU is randomly distributed within a circle with a radius of r 0 = 10 meters from RU while the ULU is randomized within its square cell area. In addition, the minimum distance between RU and any user is set as 5 meters. This cell setup implies that 4 DLUs, 4 ULUs, and 4 RUs are simulated. Furthermore, the transmit and receive antennas at each user is set as N Tx ℓ U = N Rx k D = 2 antennas for all ℓ U and k D while each RU possesses N Tx λ R = N Rx λ R = 2 antennas for all λ R ∈ R. The numbers of data streams are d ℓ U = N Tx ℓ U and d k D = N Rx k D .For the sake of convenience, an identical power budget is set for all ULUs and all RUs, i.e., P ℓ U ,max = P U,max , ∀ℓ U ∈ U and P λ R ,max = P R,max , ∀λ R ∈ R, respectively. In addition, the power budgets are set as P R,max = 26 dBm and P U,max = 23 dBm [24,25]. The power amplifier factors are ξ −1 λ R = ξ −1 ℓ U = 2:8 and P Cir λ R = 84 (W) and P Cir [31]. The noise floor is set as σ 2 k D ,a = σ 2 k D ,c = σ 2 R = − 90 dBW [51]. Regarding energy harvesters, the nonlinear EH parameters are set P DC k D = 24 (mW) and ða k D , b k D Þ = ð 150,0:014Þ for all k D ∈ D [57]. The required harvested energy is e req k D = e req for all k D and we also set e req = −20 (dBm). In addition, we also set the required data rate of DL and UL channels as r req k D = r req for all ∀k D ∈ D and r R = J × r req . The maximum capable FH capacities and FH power scale factor are set as C D,λ R = C U,λ R = C FH and α FH U,λ R = α FH D,λ R = α FH R , respectively, for all ∀λ R ∈ R. The required minimum rate is set to r req = 0:5 (bps/Hz). Without other statements, we also set P FH U,λ R ,max = P FH D,λ R ,max = 50 (W) and C U,λ R = C D,λ R = 100 (Mbps) with λ R ∈ R over the channel bandwidth of 10 (MHz), which implies α FH R = 5 (W/(bps/Hz)) [30][31][32].
The channel matrix from node A to node B with distance d A,B meters are simulated as M A,B = d −a/2 A,BMA,B with α being the path loss exponent. Meanwhile, the SI channel matrix at RU is conducted as G λ R ,λ R = 10 −C SI /20G λ R ,λ R where C SI (in decibel) denotes the SI suppression level. We also set α = 4 and C SI = 90 dB [52]. The nominal small-scale fading matrix follows the Rician fading model asM = ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ð1 + KÞ/KM LoS p + ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi 1/ð1 + KÞM NLoS p , where K = 10 dB is the Rician factor, the line-of-sight (LoS) term M LoS is a deterministic matrix, and each element of M NLoS is an independent and identically distributed complex Gaussian variable with zero mean and unit variance, i.e., M NLoS~C N ð0, 1Þ. In this paper, we utilize the far-field uniform linear antenna array model [51] M LoS = ½1, e jθ r , e j2θ r , ⋯, e jðN r −1Þθ r † × ½1, e jθ t , e j2θ t , ⋯, e jðN t −1Þθ t , where N r and N t refer to the receiving and transmitting antennas, θ r = 2πdλ −1 sin ðϕ r Þ and θ t = 2πdλ −1 sin ðϕ t Þ with antenna spacing d = 0:5λ, carrier wavelength λ, the angle-of-arrival ϕ r , and the angle-of-departure ϕ t . The values of ϕ r , ϕ t are randomly created in a range of ½0, 2π. Regarding the iteration algorithm, we set the error tolerance to be ε = 10 −4 and the maximum iteration to be κ max = 150. In the MOO algorithm, we set δ = 0:005.

System
Performance Results. Firstly, the system SE and EE performances versus FH capacity with several design methods are illustrated in Figures 2 and 3, respectively. As can be seen in Figure 2, in the phenomenon of increasing FH capacity, the achievable sum rate first rapidly grows and then saturates for all design methods. From Figure 3, the system EE performance of the sum rate maximization algorithm is significantly degraded as the FH capacity is higher than a specific value, i.e., 12 bps/Hz for such system parameters. It is because that almost all of the quantity of the power budget is utilized to improve the system sum rate in the SE design (see also Figure 4), as well as the amount of consumed power is exponentially climbing with respect to the value of the FH rate as the rate-dependent FH power consumption model (see also Figure 5). Consequently, the improved gain of the data rate substantially drops compared to the steep acceleration in the amount of consumed power, and thus, system EE should be dramatically decreased. On the contrary, the EE optimization and EE-SE MOO approaches tend to stabilize system EE metric in the high regime of FH capacity even though a lower throughput should be obtained.
On the other hand, the average systems SE and EE as the function of the minimum required harvested energy are shown in Figures 6 and 7, respectively. It is noted that at the high value of required EH e req , the problem might not be feasible. In this scenario, the systems SE and EE are set to equal zeros. It can be seen from these figures that both system performances degrade as the required EH grows. This is due to the fact that as the EH constraints are more restricted, the proportion signal for ID is reduced or/and more numbers of infeasibility of the related optimization problem can likely incur for random channel realizations. In order to obtain more insights into how algorithms worked on allocating power, various power consumption results are conducted and shown in this part. In Figure 4, the power consumption ratio is computed as the ratio of consumed power P Tot to the power budget P budget Tot From this figure, the power consumption ratio of SE optimization design approximately reaches at one for all FH configurations; this argument implies that the total power consumed by the SE maximization design and power budget is almost the same value. This can be concluded that the SE algorithm tends to completely utilize the overall available power to increase the system sum rate while the others reduce power consumption against ineffective power usage and to remain the system EE. Figures 5 and 8 plot individual power expenditures consumed by the data transmission in (20) and the FH connection in (22), respectively, with different values of FH capacity. In general, there are upward trends in consuming power for both power consumption categories with the growing in FH capacity, and the separate design approaches have reached different power expenditures. It is because more data rate can be conveyed through the FH link possessing a higher capacity level and, thus, the maximum FH rate constraints are significantly slacker that results in more power being allocated. Accordingly, we can conclude from these two figures that the power consumption designed with the SE maximization algorithm reaches the highest value in all FH capacity setups whereas that of the power minimization design is the lowest one. Otherwise, the power figures corresponding either to the EE maximization or to EE-SE trade-off approaches constitute middle values between powers consumed by two former approaches since they are facilitated to achieve a balance between the rate and power utilities.
On the other hand, the power expending for FH links is remarkably greater than that for data transmission. For example, at C FH = 10 (bps/Hz), the maximum power for data transmission only makes up at nearly 3 (W) while that for FH

15
Wireless Communications and Mobile Computing links accounts for approximate 400 (W). These different levels are more significant at the high FH capacity. Therefore, it is necessary to take rate-dependent FH power consumption into consideration in investigating the CRAN systems.

Average Pareto for SE-EE Trade-Off Results.
In this subsection, the optimal Pareto sets of various trade-offs between two among three considered metrics with several parameters are given. Specifically, the EE-power consumption, the energyspectral efficiency, and the spectrum-power consumption trade-offs with C FH = 20 bps/Hz or C FH = 30 bps/Hz are depicted in Figures 9 and 10, respectively. In these figures, the Pareto curves are obtained by changing the MOOweighted parameter, i.e., 0 < ω 1 < 1, and the optimal solutions of SE, EE, and power consumption optimization problems are also depicted. As can be seen in the figures, as expected, a better Pareto region can be achieved either with a smaller value of α FH or with a larger FH capacity C FH . Furthermore, the energy-spectral efficiency and SE-power consumption Pareto sets are concave curves and have a maximum point; this opti-mum is almost the same as the optimal solution with the EE maximization design. Moreover, we can see in these figures that the SE maximization approach consumes the largest amount of power and reaches the greatest rate while the opposite is true for the power consumption minimization design method. Meanwhile, the EE maximization algorithm tends to obtain the maximum EE with an appropriate sum rate and affordable consumed power. Alternatively, the EE-SE MOO manner permits to adjust the priority in designing either with the SE maximization or with the power consumption minimization.
Finally, Figure 11 shows the EE-SE surface which is the trade-off between systems EE and SE with respect to various FH capacities C FH . As we can see, the more the capacity of the FH link, the larger EE-SE region we can obtain. This is because, in this case, the FH rate constraints are more relax and, thus, more power can be allocated. On the other hand, the figure also shows that there are three different optimum behaviors corresponding to SE maximization, EE maximization, and power consumption minimization design methods.   Figure 10: The optimal Pareto sets with α FH = f3, 5, 8g (W/bps) and C FH = 30 (bps/Hz). 16 Wireless Communications and Mobile Computing Particularly, the system SE performance obtained by the SE maximization algorithm has a tendency to reach the highest sum rate even if its EE metric is considerably degraded in high FH capacity values. Alternatively, the EE approach tends to maximize the system EE with a suitable sum rate, whereas the power minimization method provides the lowest performance in terms of both SE and EE metrics for all FH capacity configurations since the only enough amount of power is used to achieve at the minimal system requirements.
It is worth noting that the MOO algorithm provides an optimal Pareto curving surface in which the optimal solutions of the other approaches are considered as its particular points. Hence, the EE-SE trade-off-based MOO approach is actually a more general framework as compared to conventional SOO-based design methods which are sum rate maximization, EE maximization, and power consumption minimization approaches.

Conclusion
In this paper, we have considered a problem of EE-SE tradeoffs in FD MU-MIMO CRANs with SWIPT. The problem of interest is to design the transmit precoders, FH compression matrices, and EH power splitting factors to investigate the trade-offs of the system EE-SE. We have first formulated three common optimization problems by using SOO programming including the system sum rate maximization, power consumption minimization, and fractional sum rate to power consumption ratio maximization. Furthermore, we also have made use of the efficient MOO framework to simultaneously optimize both system sum-rate and overall power consumption. Each problem is constrained by the minimum QoS of the user data rate, transmit power budget, limited FH rate, and minimum harvested energy require-ment. All design problems are highly nonconvex because of the nonconcavity in user data rate expressions, nonconvex FH rate conditions, and nonconvex EH constraints and, thus, extremely intractable to tackle directly. For the nonconvex SOO problems, the inner convex quadratic approximations are conducted in the data rates, FH rates, and harvested energy to exploit the convexity in the optimization problems.
Regarding the MOO EE-SE trade-off problem, it has first been recast to an SOO problem by applying the modified weighted Tchebycheff method, and then, the inner convex quadratic approximations are applied. As a result, SCQPbased iterative algorithms guaranteeing convergence are developed. The simulation results have revealed the design effectiveness in terms of the system EE in cases of sufficient large FH capacity as well as unveiled the interesting tradeoffs between the system energy and spectral efficiencies. Furthermore, the results have also illustrated that the conventional problems relying on SOO programming can be referred to as a special case of the considered MOO-based algorithm.

Appendix
A. Concave Lower Bound of the System Rate  (9) is thus lower bounded by a quadratic function in (33).
Regarding UL rate (18), upon defining H R,U ≜ ½H R,1 U , H R,2 U , ⋯H R,K U and F U ≜ diag fF ℓ U , ℓ U ∈ Ug, this implies that Then (18) can be rewritten, with d U = ∑ ℓ U ∈U d ℓ U , as ðA:4Þ Next, by utilizing inequality (A.1) for (A.4), the UL rate (18) is globally lower bounded by the concave quadratic function in (35). This is end of proof.

B. Convex Upper Bound of the FH Link Rate
Firstly, equation (5) can be rewritten as the difference of two concavities g D,λ R = log 2 Θ λ R − log 2 Q D,λ R , ðB:1Þ where Θ λ R = ∑ k D ∈D E † λ R F k D F † k D E λ R + Q D,λ R . Accordingly, its upper bound can be achieved by taking upper and lower bounds of first and second terms in the right-hand side (RHS) of (B.1), respectively. To this end, the upper bound of concave function log 2 jΘ λ R j around fF ðκÞ k D , Q ðκÞ D,λ R : k D ∈ Dg is obtained by taking the first-order Taylor approximation as Hence, by plugging these bounds into (B.1), we obtain the convex expression in (39) which is the upper bound of the DL FH rate. On the other hand, with a similar procedure, the upper bound of UL FH rate is given in (41). This is end of proof.

C. Linear Lower Bound of Harvested Energy at DLU
Starting from (11), we first introduce auxiliary diagonal matrix variables X k D ⪰ I N Rx with k D ∈ D; the harvested energy of DLU k D is then lower projected by using the epigraph form as ðC:1Þ It is noted that the first three terms in the right-hand side of (C.1) are convex that makes it nonconcave. Therefore, in the following, we represent in detail the method to linearise these terms. Firstly, by using the matrix inequality ( [59], Appendix B and [58], Appendix A), for all Y 0 , Y, M 0 ≻ 0 and M ≻ 0, we get with φðH k D F i D , X k D Þ being already defined in (46). That yields Similarly, with the same steps, the second trace part in (C.1) can be lower bounded as with φðG k D ,ℓ U F ℓ U , X k D Þ being given in (47). Turning our attention into the third nonconcave nonconvex trace component in (C.1), by [60], (Eq. B.3), for all Y 0 ≻ 0, Y ≻ 0, M 0 ≻ 0, and M ≻ 0, we get Wireless Communications and Mobile Computing Finally, by plugging the bounded expressions in (C.4), (C.5), and (C.7) into (C.1), the linear lower bound formula for harvested energy is completed in (45). This is end of proof.

Data Availability
The data used to support the findings of this study are included within the article.

Conflicts of Interest
The authors declare that they have no conflicts of interest.