Joint beamforming and power splitting design for MISO downlink communication with SWIPT: a comparison between cell-free massive MIMO and small-cell deployments

Simultaneous wireless information and power transfer (SWIPT) has been advocated as a highly promising technology for enhancing the capabilities of 5G and 6G devices. However, the challenge of dealing with large propagation path loss poses a signifi‑ cant hurdle. To address this issue, massive multiple‑input multiple‑output (MIMO) is employed to enhance the efficiency of SWIPT in cellular‑based networks with mul‑ tiple small cells, and especially increase the energy for cell‑edge users. In addition, by leveraging a large set of spatially distributed base stations to collaboratively serve SWIPT‑enabled user equipment, the cell‑free massive MIMO has the potential to pro‑ vide even better performance than the conventional small‑cell systems. In this work, we extend the investigation to include the application of SWIPT technology with alter‑ nating current (AC) logic in the cell‑free networks


Introduction
Conventionally, cellular networks are constrained by inter-cell interference, wherein a user close to cell boundary would particularly suffer from strong interferences from neighbouring base stations (BSs).To alleviate this problem, the technique of network multiple-input multiple-output (MIMO), which is also referred to as distributed MIMO or coordinated multi-point transmission (CoMP), is proposed by means of a coherent cooperation among BSs to reduce such interferences [1,2].For this aim, BSs in network MIMO require advanced backhaul links to support the signal transmission in the downlink (DL) and the signal detection in the uplink (UL).In addition, it requires a central processing unit (CPU) to collect the precoded signals and channel state information (CSI) among BSs.By bearing these costs, the cooperation can yield significant improvements on spectrum efficiency (SE) and coverage probability [3,4], inspiring many related studies [5][6][7][8][9][10][11][12].
As a variant of the above, cell-free massive MIMO has been introduced and gained much interest recently, in which a large number of BSs equipped with several antennas is adopted to serve a relatively smaller number of user equipment (UEs) that share the same time-frequency resources.Herein, each UE would be surrounded by a considerable number of BSs, thereby experiencing a high degree of macro-diversity and low path losses.This implies a uniformly good level of quality-of-service (QoS) to be perceived by a UE despite its particular location in the network.Given the special merit, [13] has demonstrated that cell-free massive MIMO would pave a scalable way to implement CoMP and [14] exhibited that it has significantly superior performance than conventional small-cell networks, where each user is only served by a single BS.
Despite the differences between cell-free MIMO and small-cell, energy consumption is always an issue in networking, and energy harvesting (EH) would be one of the most important solutions to address this issue, especially for the massive amount of lowerpowered wireless devices deployed in the Internet of Things (IoT) [15,16].For efficient EH, simultaneous wireless information and power transfer (SWIPT) has been widely adopted as a special case of the wireless power transfer (WPT) system.As the name implies, SWIPT allows for the simultaneous transfer of both information and power, and it is beneficial for ensuring an acceptable level of EH without the need for costly infrastructure upgrades [17,18].
In conventional small-cell networks, SWIPT is usually studied by means of co-located massive MIMO [19][20][21][22][23][24].For example, [19] formulated achievable data rate and EH expressions for the SWIP-enabled massive MIMO systems over Rician fading channels.Similarly, [20] proposed to jointly decide transmit power and time-slot duration in order to optimize energy efficiency (EE) subject to its QoS constraints in the SWIPT-enabled system with energy beamforming.In addition, joint optimization of beamforming and power splitting has also received significant attention in SWIPT-enabled networks.Specifically, the joint optimization problem of satisfying both signal-to-inference-and-noise ratio (SINR) and harvested power constraints at each user with minimum transmission power has been addressed in [21], for instance.In that work, a feasible point pursuitsuccessive convex approximation (FPP-SCA) method was introduced to address the minimization problem for the total power of beamforming vectors (BVs).Moreover, for simultaneously optimizing data rates and harvested powers in such networks, the authors in [22] proposed a majorization-minimization (MM) approach, which could outperform the classical block-diagonalization (BD) strategy.Additionally, other optimization objectives, such as maximizing the minimum (or max-min) throughput and energy, have also been proposed [23,24].
Unlike the aforementioned progress on SWIPT, early research on cell-free massive MIMO networks primarily focused on developments without SWIPT [25][26][27][28].For example, in [25], an optimization problem was presented to jointly determine nonorthogonal pilot sequences and base stations (BSs), aiming to maximize the achievable UL and DL rates through the proposed transmit power control.In [26], a distributed conjugate beamforming method was proposed for the multi-group multicast cell-free massive MIMO network.Similarly, a minimum mean square error (MMSE)-based precoding was introduced in [27] for cell-free massive MIMO systems.Furthermore, a joint optimization model for DL beamforming and power control was proposed in [28], and a joint beamforming and power optimization framework for massive MIMO systems can be found in [29].
After the initial works, cell-free massive MIMO has made further progress with the integration of WPT.Particularly, in the scenario where UEs can utilize the energy harvested in DL to transmit UL pilot and information signal to BSs, the authors in [30,31] aimed to maximize the minimum SE of UEs while considering the transmission power constraints of both BSs and UEs.In addition, a wirelessly powered IoT system based on the cell-free massive MIMO technology has also been proposed in [32], where the UL and DL power control coefficients are jointly optimized to minimize the total transmit energy consumption while satisfying the target SINRs.Then, a long-term scheduling and power control approach to maximize the minimum time-average data rate among all sensors in WPT-aided cell-free IoT networks was further studied in [33].Additionally, the expressions for the average DL data rate for EH users based on stochastic geometry were derived in [34], including the coexistence of cell-free massive MIMO and SWIPT.
Despite the different types of networks as mentioned, most related works on SWIPT considered a linear EH model, which rectifies the received radio frequency (RF) signals into required direct current (DC) power through RF-EH circuits.However, the process of rectification can result in losses, leading to a decrease in EE.More recently, based on power splitting, a receiver architecture was proposed in [35], allowing computing circuits powered directly with alternating current (AC) voltages rather than direct current (DC) for the SWIPT receiver.Given the innovation, this work, however, studied only a simple point-to-point link configuration as an example.Here, we extend the information decoding and energy harvesting (ID-EH) receiver architecture to consider the multiple input single output (MISO) interference channel (IC) in the small-cell networks, 1 and the cell-free massive MIMO networks.Building on the extended architecture, we propose a joint beamforming and power splitting optimization framework to maximize the sum-rate while meeting energy constraints, AC logic energy supply, and total transmit power for both types of networks.
The framework, however, has intractable optimization forms mainly due to the AC computing involved.To tackle the challenge, we introduce a quadratic transformbased fractional programming (FP) model to address the non-convex optimization problem on beamforming for each type of network and provide near-optimal solutions while using a one-dimensional search algorithm to determine power splitting ratios (PSRs).To further reduce the computational complexities, we propose alternative solutions based on conventional beamforming schemes such as zero forcing (ZF), maximum ratio transmission (MRT), and signal-to-leakage-and-noise ratio (SLNR) to solve the beamforming design problem, while leaving the one-dimensional search algorithm in the above to determine PSRs.Especially unlike the previous works [14,36,37] that consider the conventional beamforming schemes adopted in non-SWIPT wireless networks for cell-free [36], small-cell [37], or both [14], we focus on developing the FP-based algorithms tailored for the SWIPT-enabled networks, along with their variants based on the conventional beamforming schemes.Furthermore, we conduct a comparative analysis of these FP algorithms and their variants, spanning both types of networks, similar to the approach in [14], which also considers both of the network scenarios.Beyond the perspective of scenarios, our study also provides insights into EH aspects, which were not addressed in the previous works.More specially, we categorized the reviewed literature based on their environment (small-cell or cell-free), energy harvesting (SWIPT or non-SWIPT), and whether they employed a comparative study (comparative or non-comparative).This categorization is now presented in Table 1.While the table reveals a considerable number of works across the considered spectrum, only a limited number of comparative studies exist apart from our work, such as that in [14].In that work ( [14]), it, however, does not consider SWIPT.In contrast, our comparative work not only encompasses beamforming methods across the spectrum but also integrates the AC-relevant SWIPT technique to tackle the non-convex optimization problem for both BV and PSR in these environments.In summary, the key characteristics and contributions of our work are as follows: • We improve EE by extending the concept of AC computing to power the computation circuits with AC power in both cell-free MIMO networks and small-cell
The rest of this paper is structured as follows: In Sect.2, we present the system and channel models for cell-free massive MIMO networks considered in this work.Based on these models, we then formulate the joint design optimization problem for cell-free networks.To address this problem, we propose a quadratic transformbased fractional programming algorithm in Sect.3, followed by a sequential beamforming and power splitting algorithm in Sect. 4. Next, we introduce the system model and problem formulation for small-cell networks in Sect. 5.In this context, a corresponding fractional programming algorithm is presented to optimize the problem in Sect.6, along with a sequential algorithm in Sect.7. Numerical experiments are then conducted in Sect.8 to compare the proposed algorithms.Finally, conclusions are drawn in Sect.9.
Notations In this work, we follow the writing convention where boldface lowercase and uppercase symbols represent vectors and matrices, respectively.The superscripts (•) −1 , (•) + , (•) † , (•) ⋆ , and (•) T denote the matrix inverse, Moore-Penrose pseudo-inverse, conjugate transpose, conjugate, and transpose, respectively.The absolute value operator is denoted by | • | , and the Euclidean norm is shown by || • || .Additionally, the operation Re (•) represents the real part of a complex number, I M stands for an M × M identity matrix, and #{K } signifies the cardinality of the set K.

Methods
In this section, we introduce the system and channel models and formulate the joint design optimization problem for cell-free massive MIMO networks.Then, due to the extensive content of the corresponding methods, we organize them into Sects.3 and 4 and for detailed exposition.Afterwards, we transition to addressing the system model and problem formulation for small-cell networks in Sect.5, along with a detailed discussion of the methods specific to the small-cell scenario in Sects.6 and 7.

System model
As shown in Fig. 1a, a cell-free massive MIMO network is considered to consist of a set of BSs, B = BS 1 , BS 2 , . . ., BS N b , where N b = #{B} is the number of BSs in B .In this network, an ideal backhaul link is assumed to connect all BSs with a central processing unit (CPU).These BSs serve a set of UEs, U = UE 1 , UE 2 , . . ., UE N u , where N u = #{U} is the number of UEs in U .In addition, each BS can employ M BS > 1 antennas, while each UE would use only a single antenna for opera- tions.Given that, a UE-BS association is typically performed to enable each UE to be served by a subset of BSs in the network, including a set of all BSs in the conventional sense.Further, the set of BSs serving UE j is denoted by B j = BS j 1 , BS j 2 , . . ., BS j n b .Similarly, BS i can serve a subset of UEs, denoted by U i = UE i 1 , UE i 2 , . . ., UE i nu .For example, B 1 = {BS 1 , BS 2 , BS 3 , BS 4 , BSBS 5 , BS 6 , BS 7 } and U 6 = {UE 1 , UE 5 , UE 6 , UE 7 , UE 8 , UE 9 , UE 10 } as shown in Fig. 1a.

Channel model
Next, we consider a flat-and-block fading channel model, which is widely used for mmWave networks in 5 G and 6 G [38][39][40] to conduct the cell-free network.The model has also been investigated for massive MIMO networks in which the channel is almost flat within a resource block [41] and could be efficiently estimated [38,40,41].Furthermore, the model allows incorporating specific space-time coding schemes into massive MIMO.For example, Bell Labs layered space-time (BLAST) coding, which achieves high spectral efficiency by using multiple antennas to transmit multiple data streams simultaneously, along with dirty paper coding, has been adopted in [42] to conduct a hybrid beamforming scheme for massive MIMO under the flat-fading channel model.
Specifically, this model involves several key notations.First, the fading coefficient including the path-loss and shadowing between BS i and UE j in time t is denoted by β i,j [t] , while the complex path gain in time t is represented by α p i,j [t] .Then, the array response vector corresponding to the angle of departure and the number of channel paths in the system is given by N p i,j .Providing that BSs can receive measurements from UEs correctly and relay them to the CPU that performs the joint beamforming and power splitting (PS), the channel vector (CV) between BS i and UE j in time t can be obtained by Further, as we focus on the DL context wherein BSs are transmitting to multiple UEs, the transmitted signal from BS i in time t can be represented by Here, ω i,j [t] ∈ C MBS×1 denotes the DL beamforming vector (BV) between BS i and UE j in time t, while s j is the data symbol for UE j.Given that, the received signal at UE j in time t can be written as In the above, h i,j [t] is the CV just introduced in (1), and n A,j ∼ CN (0, ρ 2 A,j ) denotes the additive complex Gaussian noise at UE j.Given that, the total received power at this UE in time t can be obtained by

Nonlinear energy harvesting model
By means of SWIPT, each UE j can simultaneously perform ID and EH on the received signal ϒ j [t] with a PS scheme.Specifically, the PSR adopted by UE j is denoted by θ j [t] ∈ (0, 1) as that in [43].Further, a power splitter is introduced to divide the received power P r j [t] into two flows.One of them is given by (1 − θ j )P r j and directed to the decoding block while the other is given by θ j P r j and directed to the EH block, as shown in Fig. 1b.Consequently, the instantaneous SINR at UE j can be given by where ρ 2 I,j denotes the variance of random circuit noise for ID, n I,j , in addition to ρ 2 A,j just introduced.Given the SINR, the achievable data rate at UE j would be However, unlike previous works, we consider the AC computing architecture in [35] and conduct another power splitter with a PSR φ j to further divide the received power flow for EH into two parts.The first part, accounting for 1 − φ j of this flow, is used to supply the AC computational logic, while the second part, with a proportion of φ j , is dedicated to the energy management block, as shown in Fig. 1b as well.According to the power splits, the power to supply the AC computation at UE j can be given by (1) (2) Similarly, with θ j φ j of the power flow utilized to the DC energy for energy management, the EH for DC at UE j can be represented by [35] where ) , M EH is the maximum harvested energy at a UE, and a and b are constants for circuit specification [44].

Cell-free problem formulation
Based on the above system model, our objective is to optimize BVs and PSRs in order to maximize the system sum-rate while ensuring that the energy allocated for the energy management block fulfils our demand.In addition, it needs to ensure that the DC harvested energy meets a specified requirement, and the transmit power of BSs remains within a bounded range.Specifically, the optimization problem can be formulated using the following programming model, despite the time index t: where W is the set of BVs ( ω i,j , ∀i ∈ B j , ∀j ∈ U ), and and are the sets of PSRs ( θ j and φ j , ∀j ∈ U ), respectively.Given that, the first two sets of constraints, (9b) and (9c), spec- ify the requirements that the energy allocated for the AC computational logic and that for the energy management block should be higher than ǫ ac and ǫ dc , respectively.Then, the set of constraints (9d) would ensure the transmit power of BVs to be bounded within [P min , P max ] .Finally, the two sets constraints (9e) and (9f ) simply denote that PRSs, θ j and φ j , ∀j ∈ U , should lie within the interval (0, 1).

Quadratic transform-based fractional programming algorithm for cell-free
As shown above, the optimization problem (P0) is intractable due to the coupling of variables ω i,j , θ j , and φ j in the first two sets of constraints, resulting in a non-convex form.Additionally, the AC computing part introduces more variables and constraints compared to the linear EH model conventionally used.Given that, our objective is to efficiently address the problem (P0) in ( 9) and determine a feasible design for BVs ( W ) and PRSs ( and ) at each time t with low time complexity.
To this end, we first address the sum-logarithms present in the objective of optimization problem (P0) by employing the Lagrangian dual reformulation method in [45] to decouple the logarithms.Specifically, we have the following proposition: Proposition 1 By introducing a set of auxiliary variables, ξ = ξ 1 , ξ 2 , . . ., ξ N u , the optimization problem (P0) can be reformulated as where the objective of (P1) is In this objective, the function f i,j (W , , ) in the third term is With Proposition 1, we can then jointly optimize BVs and PSRs by solving the corresponding variables, W , , , and ξ , iteratively.For this, the optimization problem (P1) is divided into two subproblems and solved, respectively, as follows: (1) Solve W and ξ : First, by fixing and (represented by * and * here), the optimization problem (P1) can be reduced to P BV for the BV design as where ν j = 1 + ξ j .As shown readily, while (P BV ) appears to have a simpler form, f i,j in the objective is, in fact, a high-dimensional non-convex function as exhibited in (12).This function entails products of matrices and inverse matrices, and its non-convex nature is not easily mitigated by employing standard fractional programming methods, such as the classical Dinkelbach's transformation [46].
To resolve this issue, we turn to a recently proposed method, namely multidimensional complex quadratic transform (MCQT) [45], which expands scalar-form FP into a matrix-form framework, effectively tackling the challenge posed by highdimensional non-convexity.Specifically, since f i,j in ( 12) fulfils the concave-convex conditions for MCQT, the problem (P BV ) can be reformulated through the following proposition: Proposition 2 By using the method of multidimensional complex quadratic transform (MCQT) [45] and introducing a set of auxiliary variables y = [y 1 , y 2 , . . ., y N u ] , along with 1−θ j , ∀j ∈ U , the optimization problem (P BV ) can be transformed to (10 subjectto (9b), (9c), (9d), (9e), and (9f) subjectto (9b), (9c), and (9d) where With the transformed quadratic form, we can update W and ξ (whose element ξ j is shown in ν j = 1 + ξ j in the above) along with y in turn, as follows: First, by fixing W and ξ , in addition to the other variables, and , to be fixed in this subproblem, we can solve ∂g 2 ∂y j = 0, ∀j ∈ U , with the first-order optimality to obtain the optimal Similarly, if y and ξ as well as and are fixed, the optimal W * , represented by its col- umn-vector elements, ω * j = ω * i,j , ∀i ∈ B j , ∀j ∈ U , can be given by In the above, η j is a dual variable introduced for the power constraint, which limits the power of BV to be within [P min , P max ] .Here, we can efficiently find its optimal value by using a conventional bisection search method, where ω j = ω i,j , ∀i ∈ B j , ∀j ∈ U.
On the other hand, by taking partial differentiation with respect to ξ j in the objec- tive f (W , �, �, ξ ) and solving ∂f ∂ξ j = 0 , we can obtain its optimal, ξ * j , ∀j ∈ U , as Note that, although the third term of f (W , �, �, ξ ) , i.e. f i,j (W , , ) in ( 12), is similar to SINR in (5), its denominator involves desired signal, interferences, and noises, which is different from that of SINR.However, as shown in (19), the optimal auxiliary ξ * j is equal to SINR γ j in (5), verifying the fact that the optimization problems (P0) and (P1) are equivalent, as exhibited in Proposition 1.
(2) Solve and : Next, if W and ξ as well as y and η are given, the optimization problem P0 can be reduced to Recall that P r j = k∈U i∈B j |h † i,j ω i,k | 2 + ρ 2 A,j denotes the power received at UE j while ignoring the time index t, and let ǫ dc = b − 1 a ln e ab (M EH −ǫ dc ) e ab ǫ dc +M EH be the linear to nonlinear threshold transformation in [35].In terms of the above, the first two constraints, (b) and (c), can be rewritten as (15) (P2) max , j∈U R j subjectto (9b), (9c), (9e), and (9f ) In the reformulation, we can see that both requirements, θ j ≥ ǫ ac (1−φ j )P r j and θ j ≥ ǫ dc φ j P r j , should be satisfied at the same time.Thus, we require where φ j is bounded within (0, 1), and the objective function in the right-hand side is in the form of max f ( 1 (1−φ j ) ), g( 1 φ j ) .As noted in [35], the optimal φ * j in this form would satisfy the condition: Accordingly, the optimal EH can be obtained in terms of Next, we substitute φ j in (P2) with φ * to find θ * j .Specifically, if the optimal θ * j = θ * , ∀j ∈ U , and, without loss of generality, along with the given φ * j , we can transform the optimization problem (P2) into a one-dimensional search problem on θ , as follows: which can be solved with an initial feasible value.
Time complexity In summary, the quadratic transform-based FP algorithm dedicated to cell-free massive MIMO can solve problem P0 in (9) by iteratively updating the optimal y , ξ , W , , and .Its pseudocode is now shown in Algorithm 1 for easy reference.Note that, in each of L m iterations at most, an inverse operation is required to find ω * j with O((N u M BS ) 3 ) time for each column-vector element j ∈ U , where N u = #{U} is the number of UEs in U as defined.Further, the number of L η bisection-search iterations at most is also required to find η * j along with the above for ω * j as shown in line 5.In addition, the one-dimensional search for θ * can be solved in O(L θ ) time, where L θ is the length of the search interval (0, 1) for θ .Apart from the above, the other parameters can be easily obtained through the corresponding closed-form solutions.Consequently, the overall time complexity could be estimated as 3 ) in the worst-time complexity analysis.(21) Algorithm 1 The quadratic transform-based FP Algorithm for cell-free massive MIMO

Sequential beamforming and power splitting algorithm for cell-free
As can be seen readily, the non-convex BV optimization is the most time-consuming step in solving problem (P0), and the proposed FP-based algorithm can effectively reduce its time complexity.However, to further reduce the complexity, we incorporate conventional beamforming schemes such as ZF, MRT, and SLNR as alternatives to the search for optimal BV.This results in a sequential beamforming and power splitting algorithm that is sufficiently time-efficient for real-time operations.That is, despite and , the optimization problem (P0) can be reduced to focus on W , as When compared with (P0), this problem exhibits a simpler programming model.However, as noted before, it remains non-convex due to the maximization objective, which includes logarithmic functions with inputs of SINR represented as a ratio of two quadratic terms.Therefore, in addition to the introduced FP method, considering conventional beamforming schemes provides an alternative to implement a simpler yet more time-efficient sequential algorithm.

ZF beamforming
In convention, zero-forcing (ZF) is utilized to cancel the cross-interference links by using appropriate BVs.For cell-free massive MIMO, the condition for ZF BV ω i,j can be written as h † i ′ ,j ω i,j = 0, ∀i ′ ∈ B\i , while maximizing ||h † i,j ω i,j || 2 .More specifically, the optimization problem with respect to ZF can be formulated as where H i,j is the set of channels h i ′ ,j , ∀i ′ ∈ B\i for link i, j , and the transmit power ||ω i,j || 2 is normalized to the unit range [0, 1] to obtain the ZF BV as

MRT beamforming
Unlike ZF, which focuses on interference cancellation, maximum ratio transmission (MRT) is another widely used technique in wireless communication that aims to enhance the received signal quality.MRT achieves this by scaling the transmitted signal based on the ratio of the channel coefficients, resulting in an improved signal-tonoise ratio (SNR) at receiver.Specifically, without considering interference, the MRT BV can be easily realized by Note that, in the conventional sense of ID, this algorithm may result in strong crossinterference because it does not consider simultaneous transmissions.However, for SWIPT, it could exhibit different characteristics in terms of EH, in addition to those expected for ID, which will be revealed in our simulation study.

SLNR beamforming
When compared with SINR adopted in Sect.3, signal-to-leakage-and-noise ratio (SLNR) has the benefit to avoid collecting channel information from other stations.Hence, it usually serves as a convenient alternative to SINR, especially when obtaining such information poses challenges for local stations.Specifically, SLNR is defined as the ratio between the transmitted power to the desired UE and that leaked to other unintended UEs, as Here, by adding the notation Hi,j to denote the set of channels h i,j , ∀i ∈ B for link i, j , we can obtain the SLNR BV in this scenario, based on [47], as Note that the BVs obtained in the above, i.e. ω ZF i,j , ω MRT i,j and ω SLNR i,j , should be scaled back to ensure that j∈U ||ω M i,j || 2 , ∀i , where M represents the method of ZF, MRT, or SLNR, remains within the original bound [P min , P max ] .( 29) Algorithm 2 The sequential beamforming and power splitting algorithm for cell-free massive MIMO Time complexity Given these beamforming schemes, which provide W without finding ξ in Algorithm 1, the remaining PSR problem can be solved with the same way to find * and * in Sect.3. As a summary, the overall sequential beamforming and power splitting algorithm for cell-free massive MIMO is tabulated in Algorithm 2 for easy reference.As exhibited therein, this algorithm first requires O(N b N u ) to find BVs.Then, it needs O(L θ ) time to find PSRs with the one-dimensional search as that in Algorithm 1.Thus, its time complexity could be estimated as O(N b N u + L θ ).

System Model
To align with the cell-free environment, the DL wireless network is considered to comprise N b cells, with each cell containing a BS equipped with M BS antennas for transmitting to a single-antenna UE, as shown in Fig. 2. In addition, each cell can support multiple UEs by utilizing orthogonal frequency bands, thereby avoiding intracell interference.However, because there could be a frequency band adopted by, say UE i (served by BS i), to be shared by other BSs k ∈ B i \i , inter-cell interference is still possible.Therefore, when analysing a specific frequency band, we can model the channel of this system as a multi-cell multiple-input single-output interference channel (MISO-IC).Specifically, in the context of a specific orthogonal frequency band, each BS i can pre-code the transmitted signal xi with its BV, ωi , to communicate with the intended UE i.Consequently, the received signal at UE i can be formulated as Fig. 2 Small-cell network: a topology layout, and b SWIPT system with integrate ID and EH (a two-BSs example) Note that, unlike in the previous cell-free scenario wherein a CPU can assist in determining ω i,j between BS i and UE j, in the small-cell scenario, each BS i can only con- duct its ωi towards its intended UE i on a specific band, without further coordination to decide BVs towards UEs in other cells.
Given that, the SINR at UE i can be given by With this SINR, the achievable data rate at UE i in the small-cell scenario can be denoted by Similarly, with the nonlinear EH model in Sect.2.3, the power to supply the AC computation at UE i in this scenario can be given by Further, with ) , the EH for DC at UE i would be repre- sented by

Problem formulation
As with the case of cell-free communication, our goal in the small-cell scenario is to optimize BVs and PSRs so that the system sum-rate can be maximized, subject to the constraints on harvested energy, AC logic energy supply, and total transmit power, as well.However, by denoting W as the set of BVs ( ωi , ∀i ), and and as the sets of PSRs ( θi and φi , ∀i ), which are specific to the small-cell scenario, the optimization problem gains its dis- tinct characteristics and can be formulated as: (33)

Quadratic transform-based fractional programming algorithm for small-cell
For this scenario, we can similarly adopt the Lagrangian dual reformulation with a set of auxiliary variables ξ = ξ1 , ξ2 , . . ., ξN b , and according to the MCQT method, rewrite the objective of (P5) to where 1−θ i , ∀i .Then, by applying the first-order optimality to f for ξ , i.e.
solving ∂ f ∂ ξi = 0, ∀i , we can obtain which is equal to the SINR for the small-cell scenario.That is, the two objectives (38a) and ( 39) are equivalent.Given that, we can similarly apply the MCQT method, similarly to that in Proposition 2, to the latter, i.e. (39), resulting in the transformed objective, where νi represents 1 + ξi , and y = [ŷ 1 , ŷ2 , . . ., ŷN b ] denotes a set of auxiliary variables for the small-cell scenario.Then, to find its optimal value, we can solve ∂ Q ∂ ŷi = 0, ∀i , with fixed W , , , and ξ , resulting in Similarly, with fixed y , , , and ξ , the optimal W * , comprising ω * i , ∀i , can be obtained by where ηi is a dual variable introduced for the power constraint in the small-cell scenario, which can be similarly obtained by using a bisection search method.
Algorithm 3 The quadratic transform-based FP Algorithm for small-cell Apart from the above for W and ξ , the steps to obtain the optimal values for and can be done like that shown in Sect.3. Finally, we present the FP-based algorithm for the smallcell scenario in Algorithm 3 for easy reference.

Sequential beamforming and power splitting algorithm for small-cell
For the small-cell scenario, when and are not considered, the optimization problem (P5) can be reduced to For this problem, Ĥi denotes the set of channels, h k,i , ∀k ∈ B i \i , and B i represents the set of BSs that can interfere with UE i rather than serve it in the cell-free scenario.Additionally, the transmit power || ωi || 2 is normalized to the unit range [0, 1] as before.Given that, the ZF BV for the small-cell scenario can be obtained by Similarly, by focusing on SNR and ignoring cross-interference as before, the MRT BV can be realized in this scenario as At last, in terms of Ȟi that represents the set of channels, h k,i , ∀k ∈ B i , the SLNR BV for the small-cell scenario can be determined by Note that these BVs ( ωZF i , ωMRT i , and ωSLNR i ) would then be scaled back to make || ωM i || 2 , ∀i , where M ∈ { ZF, MRT, SLNR }, within the original [P min , P max ] .

Algorithm 4
The sequential beamforming and power splitting algorithm for small-cell (45) After addressing the BV issue, the PSR problem for the small-cell scenario can be similarly solved as described in Sect. 3 to find the optimal values for and .Finally, the overall sequential beamforming and power splitting algorithm for this scenario is summarized in Algorithm 4 for easy reference.

Results and discussion
In this section, we conduct simulation experiments to evaluate the performance of the algorithms developed for the cell-free mass MIMO networks and the small-cell networks, respectively.To this end, the simulation environment is first introduced.Then, rate-energy representations, also known as rate-energy regions, are used to illustrate the trade-off between the system sum-rate and the amount of EH realized by the different algorithms proposed.In addition, we also include the approximation solution for the max-min power control with ZF beamforming proposed in [29].This method represents a low-complexity solution with conventional beamforming in the literature.However, it does not provide a solution for the PSR problem in the SWIPT system.Thus, we randomly generate a PSR after applying this method.The resulting algorithm is called 'ZF-APX' and serves as a baseline algorithm for the sequel.
Given that, the two metrics are evaluated by varying the transmit power of BSs, and the number of cells in the system, respectively, for each of the algorithms.Through the evaluation, the performance differences of these algorithms considering both ID and EH are revealed, showcasing their characteristics and providing insights that were not previously explored in works focusing solely on ID in these networks.While showing this, we would use 'cell-free scenario' and 'small-cell scenario' to refer to 'cell-free massive MIMO network' and 'small-cell network' , respectively, for easy of notation, whenever there is no confusion.

Simulation setup
By considering both system models introduced Sects.2.1 and 5.1, we conduct the simulation environment as shown in Fig. 3 to have 19 hexagonal cells with BS 0 located at the centre, BSs 1-6 located in the first tier, and BSs 7-18 located in the second tier, similar to the environment in [48] but with different radio parameters.In the cell-free scenario, the same BSs are considered to serve the corresponding UEs as in the multi-cell environment.However, it has no cell boundaries in this scenario, which results in a similar configuration for fair comparison.In addition, it is assumed that B i = B, ∀i , and U i = U, ∀i , in the simulation study, without loss of generality.Here, the simulation codes are developed on PYTHON [49] which calls MATLAB [50] to obtain the BVs involved and are executed on a 64-bit operating system with 16 GB RAM and Intel CORE i9, 5.5 GHz.
Furthermore, to simulate the conditions of a small cell where significant energy harvesting is possible in addition to data transmission, the cell radius is limited to 20 m for SWIPT.Each UE is randomly located within each cell, and the path loss between BS i and UE j is proportional to d −k i,j , where the distance between them, d j,k , is measured in kilometres, and k = 3 is the path loss exponent.Apart from the above, the energy requirements or thresholds are given as ǫ ae = 0.04764 mW and ǫ de = 0.2 mW as that in [35], and the maximum energy harvested at UE is set as M EH = 3.9 mW according to [51], while the transmit power is limited within the range between P min = 0 and P max = 38 dBm.As a summary, the important simulation parameters are listed in Table 2.  Path loss exponent (k) 3 Cell radius for small-cell 20 m Maximum transmit power ( P max ) 38 dBm Minimum transmit power ( P min ) 0 Maximum harvested energy at UE ( M EH ) 3.9 mW [51] Number of antennas in BS ( M BS ) 3

Rate-energy tradeoff
In the first set of experiments, we evaluate the proposed algorithms on the tradeoff between the system sum-rate and the energy harvesting with the given energy requirements, in the two network scenarios.To show their names concisely, the quadratic transform-based FP algorithm is denoted by FP , the sequential beamforming and power splitting algorithms with ZF, MRT, and SLNR are denoted by ZF , MRT , and SLNR , respectively, which are also the legend names for the small-cell scenario in the following figures.However, to distinguish the applied scenarios, we use ZFcf , MRT-cf , and SLNR-cf , with the suffix cf , as the legend names for the algorithms applied to the cell-free scenario.In addition, we use 'data rate' and 'system sum-rate' interchangeably in the following, whenever there is no potential for confusion.
Aided by these legends, it can be seen clearly in Fig. 4a that FP can achieve the highest data rate while trading off the amount of EH in the small-cell scenario.On the other hand, the MRT and SLNR algorithms achieve relatively lower data rates but offer trade-offs similar to FP until they reach their rate limits.In contrast, ZF limits its power towards other UEs, resulting in a reduction in the amount of EH.Moreover, as each ZF BV, ω k , ∀k ∈ B i \i , is only conducted to transmit towards its own UE k, it could not completely eliminate the interference to a specific UE i with other BVs not defined.As a consequence, the data rate of ZF is also reduced, resulting in the poorest rate-energy trade-off among the algorithms in this scenario.
However, the performance trend of ZF changes in the cell-free scenario, as shown in Fig. 4b.In this scenario, all BSs k will coordinate their transmissions towards a specific UE i, with their other BVs, ω k,i ′ , ∀i ′ ∈ U k \i , that would eliminate the interferences towards UE i.Given that, although ZF still suffers from a reduction in EH due to its zero-forcing nature, it can achieve a rate-energy trade-off that is higher than the other algorithms but lower than FP in the high rate region.
Apart from the above, it can be also seen that by controlling the transmit power leaked to other unintended UEs, SLNR could achieve higher data rate than MRT which only scales the transmitted signal to its intended UE.The performance difference in terms of data rate is further amplified in the simulated cell-free scenario where all BSs would coordinate their transmissions to minimize such leakages, resulting in reduced inferences and increased SINRs as expected.Finally, the results for both scenarios are summarized in Fig. 4c, allowing for simultaneous comparison of these algorithms in the different network scenarios.

Impact of transmit power
In this set of experiments, we aim to show the impact of the transmit power of BSs on the proposed algorithms in the different scenarios.That is, instead of the trade-off between ID (or data rate) and EH, we consider them individually with respect to the transmit power of BSs.Specifically, by varying the transmit power from 30 to 46 dBm with a step size of 2 dBm, we show the performance of these algorithms on the two metrics (data rate and EH) in Figs 5 and 6, respectively.Here, ZF-APX and ZF-APX-cf denote the baseline algorithm applied to the small cell scenario and the cell-free scenario, respectively, in addition to the others noted before.Specifically, the performance trend observed in Fig. 5 demonstrates that increasing the transmit power leads to a higher data rate for all the compared algorithms, as expected.On one hand, the data rate performance of ZF, which is the lowest among the algorithms in the small-cell scenario (shown in Fig. 5a), increases and approaches that of FP in the cell-free scenario (shown in Fig. 5b), reflecting the same trend as discussed in the previous subsection.On the other hand, it can be seen in Fig. 5b that, unlike the significant increasing trend observed for both FP and ZF, increasing the transmit power provides only marginal benefits to MRT and SLNR.This confirms the result shown in Ch. 3 of [47] that MRT is good at very low SNR and ZF is good at high SNR when global interference coordination is considered.The above is particularly noteworthy as it also reveals the fact that, without mechanisms to reflect SINR, increasing the transmit power may not effectively improve the data rate, especially for MRT, in the cell-free scenario.
In further numerical terms, the data rates of FP, ZF, MRT, and SLNR in the cell-free scenario are 3.25, 5.26, 0.85, and 1.83 times their values in the small-cell scenario, respectively, when comparing Fig. 5a and b at 46 dBm.This confirms the performance trend reported in [52] that cell-free massive MIMO can provide five-to-tenfold improvement  in throughput over small-cell counterpart.Apart from the expected trend, our results also demonstrate that not all conventional beamforming schemes are suitable for the cell-free scenario, as the numerical examples suggest.Further, our results also show that the amounts of EH obtained by FP, ZF, MRT, and SLNR in the cell-free scenario can be 1.15, 1.70, 1.42, and 1.05 times their values in the small-cell scenario, respectively, at maximum, when comparing Fig. 6a and b at 30 dBm.Clearly, the benefits provided by EH are relatively marginal compared to those for data rate.Despite the difference, these results can compensate for the previous works by showing the distinct trends to be observed in the different scenarios when taking SWIPT into account.In addition, as observed in [29], the rate performance of ZF-APX in the cell-free scenario surpasses that in the small-cell scenario, as expected.However, due to potentially lower path loss in the small-cell scenario, the energy harvested by UE from its AP with ZF-APX using random PSR in the small-cell scenario could exceed that in the cell-free scenario, especially when the transmit power is high.However, maximizing EH is not the objective of this work.In comparison, with a special concern to allocate PSR, our algorithms are able to approach the objective of maximizing the data rate while satisfying the EH requirements, yielding superior performance to the baseline in a majority of the cases in both scenarios.Finally, through the overall comparisons shown in Figs.5c and 6c, it could be summarized that augmenting the transmit power leads to a substantial enhancement in the data rate of both FP and ZF, transitioning from the small-cell to the cell-free scenario, as exhibited in Fig. 5c.However, for the other algorithms, increasing the transmit power may not have a substantial impact during the same transition, as evident from the performance metrics (data rate and EH) shown in Figs.5c and 6c.

Impact of number of cells
In this set of experiments, we vary N b ∈ {6, 9, 12, 15, 19} while setting the transmit power to 38 dBm and keeping the other parameters fixed to investigate the impacts of these algorithms on the network size in the different network scenarios.It is also assumed that B i = B, ∀i , and U i = U, ∀i , as the previous, regardless of the network size.In addition, given a network size composed of N b cells, the cell-free scenario is similarly obtained by placing the BSs at the same positions as in the small-cell scenario, which eliminates the cell boundaries and results in a similar network setup.Now, by taking a broader perspective on examining and comparing the results shown in Figs.7 and 8, we can see that increasing the number of cells, or network size, has a more pronounced impact on the performance on EH than on the data rate.This could reflect the fact that the performance on EH is more related to the total received power in (4) which could more proportionally increase as the network size increases, as shown in Fig. 8.In contrast, the performance on data rate is related to SINR in (5) in the cellfree scenario or (34) in the small-cell scenario.In either one, SINR is a ratio between the desired power and the sum of interference power and noise, and its numerator and enumerator could increase with varying degrees as the network size increases, depending on the different algorithms applied.Thus, the increased network size is not necessary to proportionally increase the data rate, as shown in Fig. 7. by taking a closer look to compare these results, we can see that while the relative performance of each algorithm remains consistent within each network scenario, there are discernible performance trends across the two metrics as the network size increases.For example, when compared with itself in the small-cell scenario, MRT becomes the worst to provide the data rate in the cell-free scenario, with very limited enhancements by increasing the network size.This result aligns with the findings shown previously and is further confirmed here by varying the network size.Similarly, the other algorithms also exhibit a trend of increasing the data rate as the network size increases, but yield more improvements with the same increase in the size.Moreover, with the middle of the transmit power range, i.e. 38 dBm, the proposed algorithms consistently outperform the baseline algorithm, ZF-APX, in terms of both data rate and energy harvesting, regardless of the network size.
According to the above including those in the previous sets of experiments, we could conclude that the sequential beamforming and power splitting algorithms would be significantly influenced by the adopted beamforming schemes while reducing the computational complexity.In this context, ZF is a good option that provides a reasonable trade-off in the cell-free scenario, and especially it can approach the near-optimal data rate realized by FP in this scenario.In contrast, MRT and SLNR can serve as good strategies in the small-cell scenario, offering a high rate-energy trade-off similar to FP until reaching their maximum data rates.Among them, SLNR would be a more versatile strategy as it combines the respective benefits of MRT and ZF, as noted in [47].However, if feasible, FP is the preferred choice for near-optimal solutions, even if it may take more time to iteratively approach the optimum.

Conclusion
In this we have focused on addressing joint optimization problems related to downlink (DL) beamforming vectors (BVs) and power splitting ratios (PSRs) in simultaneous wireless information and power transfer (SWIPT)-enabled cell-free massive multiple-input multiple-output (MIMO) networks and small-cell networks.Our objective has been made to maximize the sum-rate while considering the constraints on harvested energy, AC logic energy supply, and total transmit power in these networks.To tackle the non-convex nature of these optimization problems, we have proposed quadratic transform-based fractional programming (FP) algorithms that can iteratively provide near-optimal solutions.To further reduce the time complexity, we have also applied conventional schemes such as zero forcing (ZF), maximum ratio transmission (MRT), and signal-to-leakage-and-noise ratio (SLNR) for beamforming vector (BV) design, while incorporating the one-dimensional search algorithm for power splitting ratio (PSR) design within the FP-based framework.
Through simulation studies, we have examined the performance differences among these algorithms in terms of system sum-rate, energy harvesting (EH), and rate-energy region as the trade-off between the two metrics.More specifically, we have shown that under the condition of this work that data rate is the objective while energy harvesting is considered as constraint, the proposed algorithms enable cell-free massive MIMO to achieve data rates up to nearly five times higher than small-cell MIMO, with limited energy harvested through SWIPT.They also outperform the baseline in both scenarios in general.Among them, FP emerges as the preferred option for near-optimal solutions, while the remaining algorithms offer versatility for diverse situations.
As a conclusion, our findings shed light on the aspects that have not been extensively explored in previous works.However, due to space constraints, there are still further research directions that warrant investigation.For instance, exploring algorithms for long-term optimization, based on the FP algorithms and their variants in the different network scenarios, would be an interesting avenue for future research.

Fig. 1
Fig. 1 Cell-free mass MIMO network: a system model, and b SWIPT architecture with integrate ID and EH (a two-BSs example)

Fig. 4
Fig. 4 Rate-energy region: a small-cell, b cell-free, and c all

Fig. 6
Fig. 6 Energy harvesting with varying transmit power: a small-cell, b cell-free, and c all

Fig. 7 Fig. 8
Fig. 7 Data rate with varying number of cells: a small-cell, b cell-free, and c all j || Liu et al.J Wireless Com Network (2024) 2024:19

Table 2
Simulation parameters