Downlink Energy Efficiency of Power Allocation and Wireless Backhaul Bandwidth Allocation in Heterogeneous Small Cell Networks
Abstract
The widespread application of wireless services and dense devices access have triggered huge energy consumption. Because of the environmental and financial considerations, energyefficient design in wireless networks becomes an inevitable trend. To the best of the authors’ knowledge, energyefficient orthogonal frequency division multiple access heterogeneous small cell optimization comprehensively considering energy efficiency maximization, power allocation, wireless backhaul bandwidth allocation, and user Quality of Service is a novel approach and research direction, and it has not been investigated. In this paper, we study the energyefficient power allocation and wireless backhaul bandwidth allocation in orthogonal frequency division multiple access heterogeneous small cell networks. Different from the existing resource allocation schemes that maximize the throughput, the studied scheme maximizes energy efficiency by allocating both transmit power of each small cell base station to users and bandwidth for backhauling, according to the channel state information and the circuit power consumption. The problem is first formulated as a nonconvex nonlinear programming problem and then it is decomposed into two convex subproblems. A near optimal iterative resource allocation algorithm is designed to solve the resource allocation problem. A suboptimal lowcomplexity approach is also developed by exploring the inherent structure and property of the energyefficient design. Simulation results demonstrate the effectiveness of the proposed algorithms by comparing with the existing schemes.
Bandwidth allocation, energy efficiency, heterogeneous network, power allocation, small cell, wireless backhaul.
I Introduction
Wireless communication networks have experienced tremendous growth in the past a few decades. It is shown that higher capacity wireless links are expected to meet the increasing quality of service (QoS) demands of multimedia applications, and these high data rate links also result in increasing device power consumption. The next generation communication systems need to provide higher data rate with limited power and bandwidth due to the rapidly increasing demands for multimedia services. Designing energyefficient wireless communication system becomes an emerging trend, due to rapidly increasing system energy costs and rising requirements of communication capacity [1, 2, 3]. According to [4] and [5], the radio access part is a major energy consumer in conventional wireless cellular networks, and it accounts for up to more than 70 percent of the total energy consumption. Therefore, increasing the energy efficiency of typical wireless networks is important to overcome the challenges raised by the rising demands of energy consumption and communication throughput.
To offload the overloaded traffics in macrocells and enhance the capacity and energy efficiency of the wireless networks, one proposed method is to shorten the distance between the base stations (BSs) and the user equipments. Small cells (e.g., picocells, femtocells and relay nodes) have been used to improve system capacity in hotspots for relieving the burden on overloaded macrocells, which is considered as a promising technique to provide an effective solution for the challenges in current macrocells [6, 7]. Therefore, there is no doubt that small cell has been paid much attention in recent years from academia and industry because it can help the system spatial reuse spectrum with low power consumption and improve the system coverage with low infrastructure cost deployment [8]. Heterogeneous small cell networks, where small cells are overlaid within a macrocell to improve coverage and increase system capacity beyond the initial deployment of macrocells, have been regarded as a promising approach to meet the increasing data traffic demand and reduce energy consumption. Although highly promising, many important problems related to heterogeneous small cell networks such as interference mitigation, resource allocation, and QoS provisioning [9, 10, 11, 12] should be addressed to fully reap the potential gains.
Resource allocation, such as power allocation and bandwidth allocation, has been widely used to maximize the energy efficiency under power limitation and QoS in heterogeneous small cell networks. Power allocation for energy efficiency has been widely studied in the literature. The distributed power control game was studied in [13] to maximize the energy efficiency of transmission for secondary users in cognitive radio networks and an optimal power control problem was formulated as a repeated game. In [14], based on the hardcore point process (HCPP), the authors investigated the maximum achievable energy efficiency of the considered multiuser multiantenna HCPP random cellular networks with the aforementioned minimum distance constraint for adjacent BSs. Different from the authors in [14], who took the minimum distance in adjacent BSs into consideration to maximize the energy efficient, we propose a suboptimal lowcomplexity approach of energyefficient backhaul bandwidth allocation by optimizing the fraction of bandwidth allocated for wireless backhauling at all small cell BSs within a macrocell range. The authors in [15] studied energyefficient power control and receiver design in cognitive radio networks, and a noncooperative power control game for maximum energy efficiency of secondary users was considered with a fairness constraint and interference threshold. The authors of [16] formulated the energyefficient spectrum sharing and power allocation in heterogeneous cognitive radio networks with femtocells as a Stackelberg game and they proposed a gradient based iteration algorithm to obtain the Stackelberg equilibrium solution to the energyefficient resource allocation problem. Some works also have been done to consider bandwidth allocation for energy efficiency. In [17], the authors studied the joint service pricing and bandwidth allocation for energy and cost efficiency at the operator level in a multitier network where an operator deploys heterogeneous small cell networks, and they formulated the problem as a Stackelberg game. The problem of joint link selection, power and bandwidth allocation for energy efficiency maximization for MultiHoming networks was investigated in [18]. A new energyefficient scheme was presented in [19] to statistically meet the demands for QoS during the bandwidth allocation for wireless networks.
In this paper, we define that wireless backhaul as the connection between macro BS and small cell BSs, and it is necessary to jointly consider the design of the radio access and backhaul network. Several related works considered the backhaul to improve energy efficiency in wireless networks. The authors of [20] studied energy efficiency of resource allocation in multicell orthogonal frequency division multiple access (OFDMA) downlink networks where the limited backhaul capacity, the circuit power consumption and the minimum required data rate are considered. The resource allocation problem for energyefficient communication with limited backhaul capacity is formulated as an optimization problem. In [21], an energy efficiency model of small cell backhaul networks with Gauss–Markov mobile models has been proposed. In [22], the authors maximized system energy efficiency in OFDMA small cell networks by optimizing backhaul data rate and emission power, and they proposed a joint forward and backhaul link optimization scheme by taking both the power consumption of forward links and the backhaul links into consideration.
To the best of the authors’ knowledge, energy efficiency of power allocation and wireless backhaul bandwidth allocation in heterogeneous small cell network has not been investigated. In this work, we study the power allocation and bandwidth allocation problem in a heterogeneous small cell network where the small cells use wireless backhauling to maximize energy efficiency of all small cell users. Similar to the paper in [23], we also use Gradient Assisted Binary Search (GABS) Algorithm to solve the energyefficient power allocation problem. Reference [24] is a conference version of this paper. Different from the conference version, we provide the detailed proof for the theorem, complexity analysis for the proposed algorithms and more simulation results in this paper. The key contributions of our work can be summarized as follows.

Design of an energyefficient OFDMA heterogeneous small cell optimization: This is a novel approach by considering energy efficiency maximization, power allocation, wireless backhaul bandwidth allocation, and user QoS into the design of OFDMA heterogeneous small cell optimization. We formulate the energyefficient wireless backhaul bandwidth allocation and power allocation problem in a heterogeneous small cell as a nonlinear programming problem, where maximum transmit power constraints of each small cell BS to each small cell user, the downlink data rate constraint of small cell BSs and the minimum data rate between each small cell BS and its corresponding user are considered to provide reliable and low energy consumed downlink transmission for small cell users. The nonconvex optimization problem is then decomposed into two convex subproblems, and an algorithm is proposed for wireless backhual bandwidth allocation and power allocation.

Support of the small cell backhauling in the context of designing power allocation schemes for heterogeneous small cell networks: We study the wireless backhaul bandwidth allocation at the small cell BS, which means a fraction of bandwidth is scheduled for backhauling and the other is assigned for communication with corresponding users. We formulate the bandwidth allocation problem as a convex problem and obtain the optimum solution.

Design of suboptimal lowcomplexity algorithm by decomposing the power allocation and bandwidth allocation: The energyefficient wireless backhaul bandwidth allocation and power allocation problem are decomposed and are optimized separately. Correspondingly, a suboptimal lowcomplexity algorithm is proposed. The effectiveness of the proposed suboptimal algorithm is demonstrated by simulations.
The rest of this paper is organized as follows. Section II describes the system model. In Section III, the energyefficient resource allocation and backhauling are presented, and in Section IV, optimization algorithms are proposed. Simulation results are discussed in Section V. Finally, Section VI concludes the paper.
II System Model
In this section, we formulate the problem of downlink energy efficiency of power allocation and unified bandwidth allocation for wireless backhauling in heterogeneous small cell networks.
We consider a heterogeneous small cell network as shown in Fig. 1 with a single macro BS, J small cells deployed within the macrocell range and K users randomly located in each small cell.
The small cells share the same spectrum with macrocell. In this work, the unified wireless backhaul bandwidth allocation is investigated. The unified bandwidth allocation factor \beta\in[0,1], which is the fraction of bandwidth allocated for wireless backhauling at all small cell BSs within a macrocell range. For simplicity, all small cells are assumed to have the same bandwidth allocation factor. We assume that the multiple antenna technology is used in the macro BS and each small cell corresponds to a beamforming group, so the interference for wireless backhaul between different small cells can be neglected. The antenna array size at macro BS is N, which is much greater than the beamforming group size B and the number of small cells, i.e., N\gg B and N\gg J. In this work, we also assume that B\geq J. Each small cell BS is equipped with single antenna. OFDMA technology is used in each small cell to support the communication between BS and users.
Let {P_{0}} be the equal transmit power of the macro BS transmit antenna targeted at corresponding small cell and {\sigma^{2}} is the additive white Gaussian noise (AWGN) power. Then the received signaltonoise ratio (SNR) in the wireless backhaul downlink of small cell j is given by
{\gamma_{j}}=\frac{{{P_{0}}{G_{j}}}}{{{\sigma^{\rm{2}}}}}.  (1) 
Let {g_{j,k}} be the channel power gain between the jth small cell BS and its corresponding kth user, where j\in\{1,2,...,J\}, k\in\{1,2,...,K\}. Let {p_{j,k}} denote the transmit power from the jth small cell BS to its corresponding kth user, and let {\bf{P}}={[{p_{j,k}}]_{J\times K}} denote the power allocation matrix.
We assume that different users in each small cell use different subchannels and cochannel interference between small cells as part of the thermal noise because of the severe wall penetration loss and low power of small cell BSs [12]. The received signaltointerferenceplusnoise ratio (SINR) of small cell user k associated with small cell j is given by
{\gamma_{j,k}}=\frac{{{p_{j,k}}{g_{j,k}}}}{{{\sigma^{\rm{2}}}{\rm{+}}{I_{j,k}}}}  (2) 
where {I_{j,k}} is the interference introduced by macro BS, {I_{j,k}}={P_{0}}{G_{j,k}}, where {G_{j,k}} is the channel power gain between macro BS and the kth user in the jth small cell. The achievable data transmission rate between the jth small cell BS and its corresponding kth user is determined by
{r_{j,k}}=\left({\frac{{1\beta}}{K}}\right){\log_{2}}\left({1+{\gamma_{j,k}}}% \right).  (3) 
Therefore, we have the relation between {r_{j,k}} and {p_{j,k}}
\displaystyle{p_{j,k}}=({2^{\frac{{K{r_{j,k}}}}{{1\beta}}}}1)\frac{{{\sigma^% {\rm{2}}}{\rm{+}}{I_{j,k}}}}{{{g_{j,k}}}}  (4)  
\displaystyle{r_{j,k}}=\left({\frac{{1\beta}}{K}}\right){\log_{2}}\left({1+% \frac{{{p_{j,k}}{g_{j,k}}}}{{{\sigma^{\rm{2}}}{\rm{+}}{I_{j,k}}}}}\right). 
Besides the transmit power during the transmission, circuit energy consumption is also incurred by device electronics in small cell BSs [25, 26]. Circuit power represents the additional device power consumption of devices during transmissions [27], such as digitaltoanalog converters, mixers and filters, and this portion of energy consumption is independent of the transmission state. If we denote the circuit power as {P_{C}} , the overall power assumption of the jth small cell BS to the kth user is {P_{C}}+{p_{j,k}}.
For energyefficient communication, it is desirable to send the maximum amount of data with a given amount of energy for small cell BSs. Hence, given any amount of energy \Delta e consumed in a duration \Delta t in each small cell BS to each user, \Delta e=\Delta t({P_{C}}+{p_{j,k}}), the small cell BSs desire to send a maximum amount of data by choosing the power allocation vector and backhaul bandwidth to maximize
\sum\limits_{j=1}^{J}{\sum\limits_{k=1}^{K}{\frac{{{r_{j,k}}(\beta,{p_{j,k}})% \Delta t}}{{\Delta e}}}}  (5) 
which is equivalent to maximizing
U(\beta,{\bf{P}})=\sum\limits_{j=1}^{J}{\sum\limits_{k=1}^{K}{{U_{j,k}}(\beta,% {p_{j,k}})}}  (6) 
where
{U_{j,k}}(\beta,{p_{j,k}})=\frac{{{r_{j,k}}(\beta,{p_{j,k}})}}{{{P_{C}}+{p_{j,% k}}}}.  (7) 
In (6), U(\beta,{\bf{P}}) is called energy efficiency for all small cells; {U_{j,k}}(\beta,{p_{j,k}}) is the energy efficiency of the kth user of the jth small cell. The unit of the energy efficiency is bits per Hertz per Joule, which has been frequently used in the literature for energyefficient communications [28, 29, 30, 31, 32].
When the downlink channel state information is estimated by the small cell BSs, the resource allocation is performed by a small cell BS under the following constraints.

Transmit power constraint of each small cell BS to each user:
0\leq{p_{j,k}}\leq{P_{\max}},\forall j,k (8) where {P_{\max}} denotes the maximal transmit power of each small cell BS to each user.

The downlink data rate constraint of each small cell BS: The throughput of the small cell is given by
{R_{j}}=\sum\limits_{k=1}^{K}{{r_{j,k}}}. (9) Due to the interuser interference within the overlapping areas of macrocell and small cell beamforming group, we use typical zeroforcing beamforming technique with equalpower allocation for each user transmission link to eliminate the interference significantly [33, 34], the capacity of the wireless backhaul downlink for small cell j is
{C_{j}}=\beta{\log_{2}}\left({1+\frac{{N{\rm{}}B{\rm{+1}}}}{B}{\gamma_{j}}}% \right). (10) The downlink wireless backhaul constraint requires
{R_{j}}\leq{C_{j}} (11) such that the downlink traffic of the jth small cell can be accommodated by its wireless backhaul.

Heterogeneous QoS guarantee: The QoS requirement {R_{t}} should be guaranteed for each user in each small cell to maintain the performance of the communication system
{r_{j,k}}\geq{R_{t}}. (12)
Our target is to maximize the energy efficiency of power allocation and unified bandwidth allocation for wireless backhauling in heterogeneous small cell networks under power constraint and data rate requirements. Thus, the corresponding problem for the downlink can be formulated as the following nonlinear programming problem
\mathop{\max}\limits_{\beta,{\bf{P}}}U(\beta,{\bf{P}})=\mathop{\max}\limits_{% \beta,{p_{j,k}}}\sum\limits_{j=1}^{J}{\sum\limits_{k=1}^{K}{{U_{j,k}}(\beta,{p% _{j,k}})}}  (13) 
s.t.  \displaystyle{C1:0\leq{p_{j,k}}\leq{P_{\max}}}  (14)  
\displaystyle{C2:{R_{j}}\leq{C_{j}}}  
\displaystyle{C3:{r_{j,k}}\geq{R_{t}}}  
\displaystyle{C4:0\leq\beta\leq 1}. 
III Energyefficient Resource Allocation and Backhauling
Since the optimization problem formulated in (13) and (14) is nonconvex and we notice that the continuous variables \beta and {p_{j,k}} are separable in (13). Therefore, we consider a decomposition approach to solve the energyefficient resource allocation problem. We decompose the nonconvex optimization problem into two convex subproblems: one for energyefficient power allocation and one for energyefficient wireless backhaul bandwidth allocation. Then, we solve the subproblems of energyefficient power allocation and energyefficient backhaul bandwidth allocation individually.
IIIA EnergyEfficient Power Allocation
The concept of quasiconcavity will be used in our discussion and is defined in [35].
Definition 1.
A function f, which maps from a convex set of real ndimensional vectors, D, to a real number, is called strictly quasiconcave if for any {x_{1}},{x_{2}}\in D and {x_{1}}\neq{x_{2}},
f(\lambda{x_{1}}+(1\lambda){x_{2}})>\min\{f({x_{1}}),f({x_{2}})\}  (15) 
for any 0<\lambda<1.
Given a value \beta for unified wireless backhaul bandwidth allocation, the optimization algorithm begins with the power allocation subproblem P1.1 that is formulated as
{\rm{P1}}{\rm{.1}}:\mathop{\max}\limits_{\bf{P}}U({\bf{P}})=\mathop{\max}% \limits_{{p_{j,k}}}\sum\limits_{j=1}^{J}{\sum\limits_{k=1}^{K}{{U_{j,k}}({p_{j% ,k}})}}  (16) 
s.t.  \displaystyle{C1:0\leq{p_{j,k}}\leq{P_{\max}}}  (17)  
\displaystyle{C2:{R_{j}}\leq{C_{j}}}  
\displaystyle{C3:{r_{j,k}}({p_{j,k}})\geq{R_{t}}} 
where
{r_{j,k}}({p_{j,k}})=\left({\frac{{1\beta}}{K}}\right){\log_{2}}\left({1+% \frac{{{p_{j,k}}{g_{j,k}}}}{{{\sigma^{\rm{2}}}{\rm{+}}{I_{j,k}}}}}\right)  (18) 
is strictly concave and monotonically increasing in {p_{j,k}} with {r_{j,k}}(0)=0, when {p_{j,k}}=0.
The optimal energyefficient power allocation achieves the maximum energy efficiency, i.e.
{{\bf{P}}^{*}}=\mathop{\arg\max}\limits_{\bf{P}}U({\bf{P}}).  (19) 
It is proved in Appendix A that U({\bf{P}}) has the following properties.
Lemma 1.
If {r_{j,k}}({p_{j,k}}) is strictly concave in {p_{j,k}}, {U_{j,k}}({p_{j,k}})\in U({\bf{P}}) is strictly quasiconcave. Furthermore, {U_{j,k}}({p_{j,k}}) is first strictly increasing and then strictly decreasing in any {p_{j,k}}, i.e. the local maximum of U({\bf{P}}) for each {p_{j,k}} exists at a positive finite value.
For strictly quasiconcave functions, if a local maximum exists, it is also globally optimal [35]. Hence, a unique globally optimal transmission rate vector always exists and its characteristics are summarized in Theorem 1 according to the proofs in Appendix A.
Theorem 1.
If {r_{j,k}}({p_{j,k}}) is strictly concave, there exists a unique globally optimal transmission power vector {{\bf{P}}^{*}}=\{{p^{*}}_{j,k};(j,k)\in J\times K\} for {{\bf{P}}^{*}}{\rm{=}}\mathop{\arg\max}\limits_{\bf{P}}U({\bf{P}}), for each element in {{\bf{P}}^{*}}, p_{j,k}^{*}{\rm{=}}\mathop{\arg\max}\limits_{{p_{j,k}}}{U_{j,k}}({p_{j,k}}) where p_{j,k}^{*} is given by
\begin{array}[]{*{20}{l}}{{{\left.{\frac{{\partial{U_{j,k}}({p_{j,k}})}}{{% \partial{p_{j,k}}}}}\right}_{{p_{j,k}}=p_{j,k}^{*}}}=0,\;f({p_{j,k}})=0,}\\ {i.e.,\;{U_{j,k}}(p_{j,k}^{*})=\frac{{{r_{j,k}}(p_{j,k}^{*})}}{{{P_{C}}+p_{j,k% }^{*}}}={{\left.{\frac{{\partial{r_{j,k}}({p_{j,k}})}}{{\partial{p_{j,k}}}}}% \right}_{{p_{j,k}}=p_{j,k}^{*}}}.}\end{array} 
In order to solve the problem P1.1 for power allocation, we rewrite the objective function in (16) as
\mathop{\max}\limits_{{p_{j,k}}}{U_{j,k}}({p_{j,k}})=\mathop{\max}\limits_{{p_% {j,k}}}\frac{{{r_{j,k}}({p_{j,k}})}}{{{P_{C}}+{p_{j,k}}}}.  (20) 
If each small cell user could reach the maximum energy efficiency, all small cell users could reach the maximum energy efficiency. The total data rate in each small cell could not exceed the capacity of the wireless backhaul downlink for small cell j, {R_{j}}\leq{C_{j}}, we can approximate that the data rate for each user to be less than \frac{{{C_{j}}}}{K}, {r_{j,k}}({p_{j,k}})\leq\frac{{{C_{j}}}}{K}, and the maximum of power for user k in small cell j is \frac{{{P_{S}}}}{K}. Thus, P1.1 is equivalent to
{\rm{P1}}{\rm{.2}}:\mathop{\max}\limits_{{p_{j,k}}}{U_{j,k}}({p_{j,k}})  (21) 
s.t.  \displaystyle{C1:0\leq{p_{j,k}}\leq\frac{{{P_{S}}}}{K}}  (22)  
\displaystyle{C2:{r_{j,k}}({p_{j,k}})\leq\frac{{{C_{j}}}}{K}}  
\displaystyle{C3:{r_{j,k}}({p_{j,k}})\geq{R_{t}}}. 
We can rewrite C2 in (22) according to (4) as
{p_{j,k}}\leq\left({\frac{{{\sigma^{\rm{2}}}{\rm{+}}{I_{j,k}}}}{{{g_{j,k}}}}}% \right)\left({{{\rm{2}}^{\left({\frac{\beta}{{1\beta}}}\right){{\log}_{2}}% \left({1+\frac{{N{\rm{}}B{\rm{+1}}}}{B}\frac{{{P_{0}}{G_{j}}}}{{{\sigma^{\rm{% 2}}}}}}\right)}}{\rm{1}}}\right).  (23) 
We can rewrite C3 in (22) according to (4) as
{p_{j,k}}\geq\left({\frac{{{\sigma^{\rm{2}}}{\rm{+}}{I_{j,k}}}}{{{g_{j,k}}}}}% \right)\left({{{\rm{2}}^{\frac{{K{R_{t}}}}{{1\beta}}}}{\rm{1}}}\right).  (24) 
Therefore, we have
{L_{j,k}}\leq{p_{j,k}}\leq{H_{j,k}}  (25) 
where
{L_{j,k}}=\left({\frac{{{\sigma^{\rm{2}}}{\rm{+}}{I_{j,k}}}}{{{g_{j,k}}}}}% \right)\left({{{\rm{2}}^{\frac{{K{R_{t}}}}{{1\beta}}}}{\rm{1}}}\right)  (26) 
{H_{j,k}}=\min\left\{{\left({\frac{{{\sigma^{\rm{2}}}{\rm{+}}{I_{j,k}}}}{{{g_{% j,k}}}}}\right)\left({{{\rm{2}}^{\left({\frac{\beta}{{1\beta}}}\right){{\log}% _{2}}\left({1+\frac{{N{\rm{}}B{\rm{+1}}}}{B}\frac{{{P_{0}}{G_{j}}}}{{{\sigma^% {\rm{2}}}}}}\right)}}{\rm{1}}}\right),{P_{\max}}}\right\}  (27) 
only if the following inequality is satisfied
{L_{j,k}}\leq{H_{j,k}}.  (28) 
The energyefficient power allocation is given by
{{\hat{p}}^{*}}_{j,k}=\mathop{\arg\max}\limits_{{p_{j,k}}}\frac{{{r_{j,k}}({p_% {j,k}})}}{{{P_{C}}+{p_{j,k}}}}  (29) 
subject to
{L_{j,k}}\leq{p_{j,k}}\leq{H_{j,k}}.  (30) 
We can solve (20) by using Theorem 1 to find the optimal power allocation solution. We can also use the lowcomplexity iterative algorithms based on the GABS algorithm proposed in [23] to realize the energyefficient power allocation for the kth user in the jth small cell BS as follows.
If the output {\hat{p}^{*}}_{j,k} satisfies the power constraint, i.e. \hat{p}_{j,k}^{*}{\rm{=}}p_{j,k}^{*}; otherwise, we can obtain the maximum {U_{j,k}}({p_{j,k}}) by
p_{j,k}^{*}={L_{j,k}}  (31) 
if {{\hat{p}}^{*}}_{j,k}<{L_{j,k}}, or we can get the maximum {U_{j,k}}({p_{j,k}}) by
p_{j,k}^{*}={H_{j,k}}  (32) 
if {{\hat{p}}^{*}}_{j,k}>{H_{j,k}}, since {U_{j,k}}({p_{j,k}}) is first strictly increasing and then strictly decreasing in any positive finite {p_{j,k}}.
IIIB EnergyEfficient Wireless Backhaul Bandwidth Allocation
Once the optimal solution {{\bf{P}}^{*}}=\{p_{j,k}^{*};(j,k)\in J\times K\} is obtained for the convex subproblem P1.2 parameterized by \beta, it can be used in the following subproblem P1.3 for the unified wireless backhaul bandwidth allocation
{\rm{P1}}.3:\mathop{\max}\limits_{\beta}U(\beta,{{\bf{P}}^{*}})=\mathop{\max}% \limits_{\beta}\sum\limits_{j=1}^{J}{\sum\limits_{k=1}^{K}{{U_{j,k}}(\beta,p_{% j,k}^{*})}}  (33) 
s.t.  \displaystyle{C1:0\leq\beta\leq 1}  (34)  
\displaystyle{C2:{R_{j}}(\beta,{{\bf{P}}^{*}})\leq{C_{j}}(\beta,{{\bf{P}}^{*}})}  
\displaystyle{C3:{r_{j,k}}(\beta,p_{j,k}^{*})\geq{R_{t}}} 
where {R_{j}}(\beta,{{\bf{P}}^{*}}) is the function value of {R_{j}} evaluated at {{\bf{P}}^{*}}, {C_{j}}(\beta,{{\bf{P}}^{*}}) is the function value of {C_{j}} evaluated at {{\bf{P}}^{*}}. In order to obtain the solution to the original problem in (13) and (14), the two subproblems P1.2 and P1.3 are solved iteratively until convergence.
Maximizing the objective function of P1.3 with respect to \beta is equivalent to maximizing \left({1\beta}\right) only, since (33) is a monotonically decreasing function of \beta. Problem P1.3 reduces to a feasibility problem whose solution is the smallest feasible value of \beta given constraints (34).
According to C2, {R_{j}}(\beta,{{\bf{P}}^{*}})\leq{C_{j}}(\beta,{{\bf{P}}^{*}}), we have
\beta\geq\frac{{\sum\limits_{k=1}^{K}{{{\log}_{2}}\left({1+\frac{{p_{j,k}^{*}{% g_{j,k}}}}{{{\sigma^{\rm{2}}}+{I_{j,k}}}}}\right)}}}{{K{{\log}_{2}}\left({1+% \frac{{N{\rm{}}B{\rm{+1}}}}{B}\frac{{{P_{0}}{G_{j}}}}{{{\sigma^{\rm{2}}}}}}% \right)+\sum\limits_{k=1}^{K}{{{\log}_{2}}\left({1+\frac{{p_{j,k}^{*}{g_{j,k}}% }}{{{\sigma^{\rm{2}}}+{I_{j,k}}}}}\right)}}}.  (35) 
According to C3, {r_{j,k}}(\beta,p_{j,k}^{*})\geq{R_{t}}, we have
\beta\leq 1\frac{{K{R_{t}}}}{{{{\log}_{2}}\left({1+\frac{{p_{j,k}^{*}{g_{j,k}% }}}{{{\sigma^{\rm{2}}}+{I_{j,k}}}}}\right)}}.  (36) 
So we have
\beta=\max\left\{{{\phi_{j}}}\right\}  (37) 
where
{\phi_{j}}=\frac{{\sum\limits_{k=1}^{K}{{{\log}_{2}}\left({1+\frac{{p_{j,k}^{*% }{g_{j,k}}}}{{{\sigma^{\rm{2}}}+{I_{j,k}}}}}\right)}}}{{K{{\log}_{2}}\left({1+% \frac{{N{\rm{}}B{\rm{+1}}}}{B}\frac{{{P_{0}}{G_{j}}}}{{{\sigma^{\rm{2}}}}}}% \right)+\sum\limits_{k=1}^{K}{{{\log}_{2}}\left({1+\frac{{p_{j,k}^{*}{g_{j,k}}% }}{{{\sigma^{\rm{2}}}+{I_{j,k}}}}}\right)}}}  (38) 
only if {R_{t}} satisfies the following condition
{R_{t}}\leq\min\left\{{{\varphi_{j}}}\right\}  (39) 
where
{\varphi_{j}}=\frac{{{{\log}_{2}}\left({1+\frac{{N{\rm{}}B{\rm{+1}}}}{B}\frac% {{{P_{0}}{G_{j}}}}{{{\sigma^{\rm{2}}}}}}\right){{\log}_{2}}\left({1+\frac{{p_{% j,k}^{*}{g_{j,k}}}}{{{\sigma^{\rm{2}}}+{I_{j,k}}}}}\right)}}{{K{{\log}_{2}}% \left({1+\frac{{N{\rm{}}B{\rm{+1}}}}{B}\frac{{{P_{0}}{G_{j}}}}{{{\sigma^{\rm{% 2}}}}}}\right)+\sum\limits_{k=1}^{K}{{{\log}_{2}}\left({1+\frac{{p_{j,k}^{*}{g% _{j,k}}}}{{{\sigma^{\rm{2}}}+{I_{j,k}}}}}\right)}}}.  (40) 
IV Algorithm Design
According to the analysis of power allocation and wireless backhaul bandwidth allocation discussed above, we propose an iterative optimization algorithm and a suboptimal lowcomplexity algorithm.
IVA Iterative Resource Allocation Algorithm
The proposed iterative resource allocation algorithm is shown in Algorithm 1.
In Algorithm 1, each small cell BS calculates {\phi_{j}} according to (38) and then sends {\phi_{j}} to macro BS. The macro BS chooses the maximum {\phi_{j}} to be the optimal bandwidth allocation factor \beta and broadcasts \beta to all small cell BSs.
IVB LowComplexity Optimization Algorithm
To reduce the complexity of Algorithm 1, we propose a lowcomplexity optimization algorithm where bandwidth allocation factor is calculated from the equal power allocation and we fix \beta to calculate the power allocation according to the scheme proposed in Section IV. This lowcomplexity optimization algorithm is shown in Algorithm 2.
IVC Complexity Analysis
Since the problem formulated in (13) and (14) is not convex, the only way to obtain the optimal solution is to use the method of exhaustion. If we assume that it costs P operations to calculate {r_{j,k}} and it costs Q operations to calculate {C_{j}}, the complexity of checking C2 and C3 in (14) entails KP+K+Q operations and P+1 operations, respectively. If we assume that it costs S operations to calculate {U_{j,k}}, the complexity of obtaining the total energy efficiency of all small cell users entails JKS+\left({J1}\right)\left({K1}\right) operations. The total complexity of getting the value of objective function in (13) under the constraints in (14) entails KP+K+P+1+JKS+\left({J1}\right)\left({K1}\right) operations for specific {p_{j,k}} and \beta values. If we assume the value step size for {p_{j,k}} is a and the value step size for \beta is b, there are \frac{1}{b}{\left({\frac{{{P_{\max}}}}{a}}\right)^{JK}} choices for the values of {p_{j,k}} and \beta. Therefore, the complexity for the method of exhaustion is O\left({\frac{{JKS}}{b}{{\left({\frac{{{P_{\max}}}}{a}}\right)}^{JK}}}\right).
In Algorithm 1, the worstcase complexity of calculating bandwidth allocation factor \beta from (37) entails J operations in each iteration. If we assume that it costs \Omega operations in each GABS to search the optimum power allocation without power constraint, then the worstcase complexity of finding the power allocation for every user in each small cell entails JK\left({\Omega{\rm{+4}}}\right) operations in each iteration. Suppose the Algorithm 1 needs \Delta iterations to converge, so the total complexity of Algorithm 1 is O\left({JK\Omega\Delta}\right). Since iteration is not applied in Algorithm 2, the total complexity of Algorithm 2 is O\left({JK\Omega}\right), which is less than that of Algorithm 1. In the simulation, the typical value for \Delta is around 16, the typical value for \Omega is less than 500, and the typical values for \frac{1}{b} and \frac{{{P_{\max}}}}{b} are both 100. So the complexities of Algorithm 1 and Algorithm 2 are always less than that of the method of exhaustion. When the number of small cells J and the number of users in each small cell K increase, the complexity of the method of exhaustion increases exponentially, so the complexity of the method of exhaustion is much larger than the complexities of two proposed algorithms.
V Simulation Results
Simulation results are given in this section to evaluate the performance of the proposed power allocation and backhaul bandwidth allocation algorithms. In the simulations, it is assumed that small cells are uniformly distributed in the macrocell coverage area, and small cell users are uniformly distributed in the coverage area of their serving small cell. AWGN power {\sigma^{\rm{2}}}{\rm{=}}3.9811\times{10^{14}} W. The coverage radius of the macrocell is 500 m, and that of a small cell is 10 m. Small cell has a minimal distance of 50 m from the macro BS. The minimal distance between small cell BSs is 40 m. We assume that the channel fading is composed of path loss, shadowing fading, and Rayleigh fading. The pathloss model for small cell users is based on [36]. The lognormal shadowing between small cell BS and small cell users is 10 dB. At the macro BS, we assume that the transmit power is 33 dBm, the antenna array size N=100 and beamforming group size is B=20. We consider that all the small cell users have the same QoS requirement.
Figure 2 shows the convergence in terms of the energy efficiency of all small cell users for the proposed Algorithm 1 versus the number of iterations, where J=5, {R_{t}}=0.01 bps/Hz, {P_{\max}}=20 dBm. It can be observed that the proposed resource allocation algorithm takes nearly 16 iterations to converge to stable solutions. This result, together with the previous analysis, ensures that the proposed Algorithm 1 is applicable in heterogeneous small cell networks.
Figure 3 shows the total energy efficiency of all small cell users when the number of users per small cell is increased from 2 to 10, for Algorithm 2 under {P_{\max}}=7 dBm, {P_{\max}}=10 dBm and {P_{\max}}=20 dBm compared with Algorithm 1 under {P_{\max}}=20 dBm. The simulation parameters are set as J=5, {R_{t}}=0.01 bps/Hz. Fig. 3 shows that the energy efficiency performance of Algorithm 1 is 20% more superior to that of Algorithm 2. It also can be seen from Fig. 3 that the more number of users in small cell is, the better performance is obtained because of the multiuser diversity.
Figure 4 shows the total energy efficiency of all small cell users when the number of small cells is increased from 3 to 15, for Algorithm 2 under {P_{\max}}=7 dBm, {P_{\max}}=10 dBm and {P_{\max}}=20 dBm when compared with Algorithm 1 under {P_{\max}}=20 dBm. The simulation parameters are set as K=5, {R_{t}}=0.01 bps/Hz. Fig. 4 indicates that more number of small cell is, the better performance is obtained. It can also be seen from Fig. 4 that the energy efficiency performance of Algorithm 1 is always better than that of Algorithm 2 and the gap between them becomes larger when the number of small cells increases. The energy efficiency performance of Algorithm 1 is 30% superior to that of Algorithm 2 when the number of small cells is 10.
Figure 5 shows the total downlink capacity of all small cell users when the number of users per small cell is increased from 2 to 10, for Algorithm 2 under {P_{\max}}=7 dBm, {P_{\max}}=10 dBm and {P_{\max}}=20 dBm compared with Algorithm 1 under {P_{\max}}=20 dBm. The simulation parameters are set as J=5, {R_{t}}=0.01 bps/Hz. Fig. 5 shows that the total downlink capacity of Algorithm 1 is more than 3 bps/Hz higher than that of Algorithm 2. It can also be seen from Fig. 5 that the more number of users in small cell is, the better performance is obtained due to the multiuser diversity. The total downlink capacity of Algorithm 1 is 21% higher than that of Algorithm 2 when the number of users in each small cell is over 10.
Figure 6 shows the total downlink capacity of all small cell users when the number of small cells is increased from 3 to 15, for Algorithm 2 under {P_{\max}}=7 dBm, {P_{\max}}=10 dBm and {P_{\max}}=20 dBm compared with Algorithm 1 under {P_{\max}}=20 dBm. The simulation parameters are set as K=5, {R_{t}}=0.01 bps/Hz. Fig. 6 illustrates that Algorithm 1 is superior to Algorithm 2 in terms of the total downlink capacity and the gap between them becomes larger when the number of small cells increases. The total downlink capacity of Algorithm 1 is 29% larger than that of Algorithm 2 when there are 14 small cells in the heterogeneous network.
Figure 7 shows the total energy efficiency of all the small cell users when using Algorithm 2 for power constraint {P_{\max}} ranging from 0 dBm to 12.79 dBm where the number of users in each small cell is 3, 4, 5. The simulation parameters are set as J=5, {R_{t}}=0.01 bps/Hz. Fig. 7 presents that the more users in each small cell are, the higher total energy efficiency can be obtained, which has already been shown in Fig. 3. It can also be seen from Fig. 7 that the larger power constraint is, the better performance is obtained. This is because the larger power constraint leads to the larger region of the optimizing variable.
Figure 8 shows the total energy efficiency of all the small cell users when the number of users per small cell is increased from 2 to 10, for different algorithms. Algorithm 1 and Algorithm 2 are the iterative optimization algorithm and the lowcomplexity optimization algorithm, respectively, which we have proposed in Section IV. Algorithm 3 is an existing energy efficiency optimization algorithm with equal power allocation and Algorithm 4 is an algorithm that uses the optimum power allocation we proposed given a random \beta to optimize energy efficiency. All the algorithms are under the constraint {P_{\max}}=20 dBm. Fig. 8 indicates that more users in each small cell are, the better performance can be obtained, which has already been shown in Fig. 3. It also can be seen from Fig. 8 that Algorithm 1 has the best performance, and then it follows by Algorithm 2, Algorithm 3 and Algorithm 4. The energy efficiency performance of Algorithm 1 is 30.5% and 56.6% higher than that of Algorithm 3 and Algorithm 4, respectively.
Figure 9 shows the total energy efficiency of all the small cell users when the number of small cells is increased from 2 to 5, for the optimal solution and the two proposed algorithms. Since the complexity of the method of exhaustion is high, we only consider the situation with small dimension where there are two users located in each small cell, K=2. All the algorithms are under the setting of {P_{\max}}=20 dBm and {R_{t}}=0.01 bps/Hz. From Fig. 9, we can observe that the difference between the optimal solution and Algorithm 1 in terms of energy efficiency is small, which ensures the effectiveness of the proposed algorithms. The energy efficiency performance of the optimal solution is only about 7% and 24% higher than that of Algorithm 1 and Algorithm 2 respectively when the number of small cells is 3. The difference between the optimal solution and the proposed Algorithm 1 is mainly caused by the approximation of C2 in (17). We can also observe that the performance of Algorithm 1 is slightly better than that of the existing algorithm, which is the backhaul bandwidth allocation in conjunction with the resource allocation algorithm in [14]. This phenomenon can be explained as follows. As the QoS requirement of small cell users increases, the more power is required to meet the higher QoS requirement.
VI Conclusion
In this paper, we investigated the energyefficient wireless backhaul bandwidth allocation and power allocation in a heterogeneous small cell network. We demonstrated the existence of a unique globally optimal energy efficiency solution and provided an iterative algorithm to obtain this optimum. For the downlink scenario, we first found the near optimal energyefficient resource allocation approach, and then developed a lowcomplexity suboptimal algorithm by exploring the inherent structure of the objective function and the feature of energyefficient design. From the simulation results, we observed that energy efficiency is improved by increasing the number of small cells and the number of users per small cell, and the capacity is also improved by increasing the number of small cells and the number of users per small cell. Simulation results showed great energy efficiency improvement of the proposed iterative optimization algorithm than that of the lowcomplexity optimization algorithm and the existing schemes. The proposed lowcomplexity algorithms can achieve a promising tradeoff between performance and complexity. If the future, we will investigate the nonunified backhaul bandwidth allocation and intersmallcell interference in heterogeneous small cell networks.
Appendix A Proof of Lemma 1
We first focus on {U_{j,k}}({p_{j,k}}) and then we can obtain the properties of U({\bf{P}}). If every user in each small cell can reach the maximum energy efficiency, all small cell users can reach the maximum energy efficiency.
Denote the \alpha–superlevel sets of {U_{j,k}}({p_{j,k}}) as
{S_{\alpha}}=\{\left.{{p_{j,k}}\geq 0}\right{U_{j,k}}({p_{j,k}})\geq\alpha\}  (41) 
where {{p_{j,k}}} is nonnegative. Based on the propositions in [35], {U_{j,k}}({p_{j,k}}) is strictly quasiconcave if and only if {S_{\alpha}} is strictly convex for any real number \alpha. In this case, when \alpha<0, no points exist on the contour {U_{j,k}}({p_{j,k}})=\alpha. When \alpha=0, only {p_{j,k}}=0 is on the contour {U_{j,k}}(0)=\alpha. Hence, {S_{\alpha}} is strictly convex when \alpha\leq 0. Now, we investigate the case when \alpha>0. We can rewrite the {S_{\alpha}} as {S_{\alpha}}=\{\left.{{p_{j,k}}\geq 0}\right\alpha{P_{C}}+\alpha{p_{j,k}}{r_% {j,k}}({p_{j,k}})\leq 0\}. Since {r_{j,k}}({p_{j,k}}) is strictly concave in {{p_{j,k}}}, {r_{j,k}}({p_{j,k}}) is strictly convex in {{p_{j,k}}}; therefore, {S_{\alpha}} is strictly convex. Hence, we have the strict quasiconcavity of {U_{j,k}}({p_{j,k}}).
Next, we can obtain the partial derivative of {U_{j,k}}({p_{j,k}}) with {p_{j,k}} as
\frac{{\partial{U_{j,k}}({p_{j,k}})}}{{\partial{p_{j,k}}}}=\frac{{({P_{C}}+{p_% {j,k}}){{r^{\prime}}_{j,k}}({p_{j,k}}){r_{j,k}}({p_{j,k}})}}{{{{\left({{P_{C}% }+{p_{j,k}}}\right)}^{2}}}}{\rm{=}}\frac{{f({p_{j,k}})}}{{{{\left({{P_{C}}+{p_% {j,k}}}\right)}^{2}}}}  (42) 
where f({p_{j,k}})=({P_{C}}+{p_{j,k}}){{r^{\prime}}_{j,k}}({p_{j,k}}){r_{j,k}}({p_{% j,k}}), {{r^{\prime}}_{j,k}}({p_{j,k}}) is the first partial derivative of {r_{j,k}}({p_{j,k}}) with respect to {p_{j,k}}. If p_{j,k}^{*} exists such that {\left.{\frac{{\partial{U_{j,k}}({p_{j,k}})}}{{\partial{p_{j,k}}}}}\right_{{p% _{j,k}}=p_{j,k}^{*}}}=0, it is unique, i.e. if there is a p_{j,k}^{*} such that f(p_{j,k}^{*})=0. In the following, we investigate the conditions when p_{j,k}^{*} exists.
The derivative of f({p_{j,k}}) is
f^{\prime}({p_{j,k}})=({P_{C}}+{p_{j,k}}){{r^{\prime\prime}}_{j,k}}({p_{j,k}})  (43) 
where {{r^{\prime\prime}}_{j,k}}({p_{j,k}}) is the second partial derivative of {r_{j,k}}({p_{j,k}}) with respect to {p_{j,k}}. Since {r_{j,k}}({p_{j,k}}) is strictly concave in {{p_{j,k}}}, so {{r^{\prime\prime}}_{j,k}}({p_{j,k}})<0, f^{\prime}({p_{j,k}})<0. Hence, f({p_{j,k}}) is strictly decreasing.
\begin{split}\displaystyle\mathop{\lim}\limits_{{p_{j,k}}\to\infty}f({p_{j,k}}% )&\displaystyle=\mathop{\lim}\limits_{{p_{j,k}}\to\infty}\left({({P_{C}}+{p_{j% ,k}}){{r^{\prime}}_{j,k}}({p_{j,k}}){r_{j,k}}({p_{j,k}})}\right)\\ &\displaystyle=\mathop{\lim}\limits_{{p_{j,k}}\to\infty}\left({{P_{C}}{{r^{% \prime}}_{j,k}}({p_{j,k}})+{p_{j,k}}{{r^{\prime}}_{j,k}}({p_{j,k}}){r_{j,k}}(% {p_{j,k}})}\right)\end{split}  (44) 
where
{{r^{\prime}}_{j,k}}({p_{j,k}})=\left({\frac{{1\beta}}{K}}\right)\left({\frac% {{{g_{j,k}}}}{{{\sigma^{\rm{2}}}+{I_{j,k}}}}}\right)\left({\frac{1}{{\ln 2}}}% \right)\left({\frac{1}{{1+\frac{{{p_{j,k}}{g_{j,k}}}}{{{\sigma^{\rm{2}}}{\rm{+% }}{I_{j,k}}}}}}}\right)  (45) 
and
\mathop{\lim}\limits_{{p_{j,k}}\to\infty}{{r^{\prime}}_{j,k}}({p_{j,k}})=0.  (46) 
So we have
\mathop{\lim}\limits_{{p_{j,k}}\to\infty}{P_{C}}{{r^{\prime}}_{j,k}}({p_{j,k}}% )=0.  (47) 
According to the L’Hopital’s rule, it is easy to show that
\begin{split}\displaystyle\mathop{\lim}\limits_{{p_{j,k}}\to\infty}{p_{j,k}}{{% r^{\prime}}_{j,k}}({p_{j,k}})&\displaystyle=\mathop{\lim}\limits_{{p_{j,k}}\to% \infty}\left({\frac{{1\beta}}{K}}\right)\left({\frac{{{g_{j,k}}}}{{{\sigma^{% \rm{2}}}+{I_{j,k}}}}}\right)\left({\frac{1}{{\ln 2}}}\right)\left({\frac{{{p_{% j,k}}}}{{1+\frac{{{p_{j,k}}{g_{j,k}}}}{{{\sigma^{\rm{2}}}{\rm{+}}{I_{j,k}}}}}}% }\right)\\ &\displaystyle=\mathop{\lim}\limits_{{p_{j,k}}\to\infty}\left({\frac{{1\beta}% }{K}}\right)\left({\frac{{{g_{j,k}}}}{{{\sigma^{\rm{2}}}+{I_{j,k}}}}}\right)% \left({\frac{1}{{\ln 2}}}\right)\left({\frac{1}{{\frac{{{g_{j,k}}}}{{{\sigma^{% \rm{2}}}+{I_{j,k}}}}}}}\right)\\ &\displaystyle=\mathop{\lim}\limits_{{p_{j,k}}\to\infty}\left({\frac{{1\beta}% }{K}}\right)\left({\frac{1}{{\ln 2}}}\right)\end{split}  (48) 
\mathop{\lim}\limits_{{p_{j,k}}\to\infty}\left({{r_{j,k}}({p_{j,k}})}\right)=% \mathop{\lim}\limits_{{p_{j,k}}\to\infty}\left[{\left({\frac{{1\beta}}{K}}% \right){{\log}_{2}}\left({1+\frac{{{p_{j,k}}{g_{j,k}}}}{{{\sigma^{\rm{2}}}{\rm% {+}}{I_{j,k}}}}}\right)}\right]=\infty.  (49) 
So we have
\mathop{\lim}\limits_{{p_{j,k}}\to\infty}f({p_{j,k}})<0.  (50) 
Besides,
\begin{split}\displaystyle\mathop{\lim}\limits_{{p_{j,k}}\to 0}f({p_{j,k}})&% \displaystyle=\mathop{\lim}\limits_{{p_{j,k}}\to 0}\left({({P_{C}}+{p_{j,k}}){% {r^{\prime}}_{j,k}}({p_{j,k}}){r_{j,k}}({p_{j,k}})}\right)\\ &\displaystyle={P_{C}}{{r^{\prime}}_{j,k}}(p_{j,k}^{(0)}){r_{j,k}}(p_{j,k}^{(% 0)})\\ \end{split}  (51) 
where p_{j,k}^{(0)} denotes {p_{j,k}}=0
\begin{split}\displaystyle{{r^{\prime}}_{j,k}}(p_{j,k}^{(0)})&\displaystyle={% \left.{\left({\frac{{1\beta}}{K}}\right)\left({\frac{{{g_{j,k}}}}{{{\sigma^{% \rm{2}}}+{I_{j,k}}}}}\right)\left({\frac{1}{{\ln 2}}}\right)\left({\frac{1}{{1% +\frac{{{p_{j,k}}{g_{j,k}}}}{{{\sigma^{\rm{2}}}{\rm{+}}{I_{j,k}}}}}}}\right)}% \right_{{p_{j,k}}=0}}\\ &\displaystyle=\left({\frac{{1\beta}}{K}}\right)\left({\frac{{{g_{j,k}}}}{{{% \sigma^{\rm{2}}}+{I_{j,k}}}}}\right)\left({\frac{1}{{\ln 2}}}\right)\end{split}  (52) 
{r_{j,k}}(p_{j,k}^{(0)})=0  (53) 
\mathop{\lim}\limits_{{p_{j,k}}\to 0}f({p_{j,k}})=\left({\frac{{1\beta}}{K}}% \right)\left({\frac{{{P_{C}}{g_{j,k}}}}{{{\sigma^{\rm{2}}}+{I_{j,k}}}}}\right)% \left({\frac{1}{{\ln 2}}}\right)>0.  (54) 
So, together with \mathop{\lim}\limits_{{p_{j,k}}\to\infty}f({p_{j,k}})<0, we see that {p^{*}}_{j,k} exists and {U_{j,k}}({p_{j,k}}) is first strictly increasing and then strictly decreasing in {p_{j,k}}.
Lemma 1 is readily obtained. \hfill\blacksquare
References
 [1] C. Jiang, H. Zhang, Y. Ren, and H.H. Chen, “Energyefficient noncooperative cognitive radio networks: Micro, meso and macro views,” IEEE Commun. Mag., vol. 52, no. 7, pp. 14–20, July 2014.
 [2] F. R. Yu, X. Zhang, and V. C.M. Leung, Green Communications and Networking., CRC Press, 2012.
 [3] C. Xu, M. Sheng, C. Yang, X. Wang, and L. Wang, “Pricingbased multiresource allocation in OFDMA cognitive radio networks: an energy efficiency perspective,” IEEE Trans. Veh. Technol., vol. 63, no. 5, pp. 2336–2348, June 2014.
 [4] T. Edler and S. Lundberg, “Energy efficiency enhancements in radio access networks,” Ericsson Review, 2004, http://www.ericsson.com/ericsson/corpinfo/publications/review/2004 01/files/2004015.pdf
 [5] Y. Chen, S. Zhang, S. Xu, and G. Y. Li, “Fundamental tradeoffs on green wireless networks,” IEEE Commun. Mag., vol. 49, no. 6, pp. 30–37, June 2011.
 [6] H. Zhang, X. Chu, W. Guo, and S. Wang, “Coexistence of WiFi and heterogeneous small cell networks sharing unlicensed spectrum,” IEEE Commun. Mag., vol. 53, no. 3, pp. 158–164, Mar. 2015.
 [7] D. LopezPerez, X. Chu, A. V. Vasilakos, and H. Claussen, “Power minimization based resource allocation for interference mitigation in OFDMA femtocell networks,” IEEE J. Sel. Areas Commun., vol. 32, no. 2, pp. 333–344, Feb. 2014.
 [8] H. Zhang, C. Jiang, J. Cheng, and V. C. M. Leung, “Cooperative interference mitigation and handover management for heterogeneous cloud small cell networks,” IEEE Wireless Commun., vol. 22, no. 3, pp. 92–99, June 2015.
 [9] H. Zhang, C. Jiang, N. C. Beaulieu, X. Chu, X. Wang, and T. Quek, “Resource allocation for cognitive small cell networks: A cooperative bargaining game theoretic approach,” IEEE Trans. Wireless Commun., vol. 14, no. 6, pp. 3481–3493, June 2015.
 [10] W. Cheng, X. Zhang, and H. Zhang, “StatisticalQoS driven energyefficiency optimization over green 5G mobile wireless networks,” IEEE J. Sel. Areas Commun., vol. 34, no. 12, pp. 30923107, Dec. 2016.
 [11] H. Zhang, C. Jiang, X. Mao, and H.H. Chen,“Interferencelimited resource optimization in cognitive femtocells with fairness and imperfect spectrum sensing,” IEEE Trans. Veh. Technol., vol. 65, no. 3, pp. 1761–1771, Mar. 2016.
 [12] H. Zhang, C. Jiang, N. C. Beaulieu, X. Chu, X. Wen, and M. Tao, “Resource allocation in spectrumsharing OFDMA femtocells with heterogeneous services,” IEEE Trans. Commun., vol. 62, no. 7, pp. 2366–2377, July 2014.
 [13] M. Le Treust and S. Lasaulce, “A repeated game formulation of energyefficient decentralized power control,” IEEE Trans. Wireless Commun., vol. 9, no. 9, pp. 2860–2869, Sep. 2010.
 [14] X Ge, B Du, Q Li, D. S. Michalopoulos, “Energy efficiency of multiuser multiantenna random cellular networks with minimum distance constraints,” IEEE Trans. Veh. Technol., vol. 66, No. 2, pp. 1696–1708, Feb. 2017.
 [15] S. Buzzi and D. Saturnino, “A gametheoretic approach to energyefficient power control and receiver design in cognitive CDMA wireless networks,” IEEE J. Sel. Topics Signal Proc., vol. 5, no. 1, pp. 137–150, Feb. 2011.
 [16] R. Xie, F. R. Yu, H. Ji, and Y. Li, “ Energyefficient resource allocation for heterogeneous cognitive radio networks with femtocells,” IEEE Trans. Wireless Commun., vol. 11, no. 11, pp. 3910–3920, Nov. 2012.
 [17] C. M. G. Gussen, E. V. Belmega, and M. Debbah, “Pricing and bandwidth allocation problems in wireless multitier networks,” Signals, Systems and Computers (ASILOMAR), 2011 Conference Record of the Forty Fifth Asilomar Conference on, pp. 1633–1637, Nov. 2011.
 [18] Q. Vu, L. Tran, M. Juntti, and E. Hong, “Energyefficient bandwidth and power allocation for multihoming networks,” IEEE Trans. Signal Process., vol. 63, no. 7, pp. 1684–1699, Apr. 2015.
 [19] W. Wang, X. Wang, and A. A. Nilsson, “Energyefficient bandwidth allocation in wireless networks: Algorithms, analysis, and simulations,” IEEE Trans. Wireless Commun., vol. 5, no. 5, pp. 1103–1114, May 2006.
 [20] D. W. K. Ng, E. S. Lo, and R. Schober, “Energyefficient resource allocation in multicell OFDMA systems with limited backhaul capacity,” IEEE Trans. Wireless Commun., vol. 11, no. 10, pp. 3618–3631, Oct. 2012.
 [21] X. Ge, S. Tu, T. Han, Q. Li, and G. Mao, “Energy efficiency of small cell backhaul networks based on GaussMarkov mobile models,” IET Networks, vol. 4, no. 2, pp. 158–167, Mar. 2015.
 [22] G. Nie, H. Tian, and J. Ren, “Energy efficient forward and backhaul link optimization in OFDMA small cell networks,” IEEE Commun. Lett., vol. 19, no. 11, pp. 1989–1992, Nov. 2015.
 [23] G. Miao, N. Himayat, and G. Y. Li, “Energyefficient link adaptation in frequencyselective channels,” IEEE Trans. Commun., vol. 58, no. 2, pp. 545–554, Feb. 2010.
 [24] H. Liu, H. Zhang, J. Cheng, and V. C. M. Leung, “Energy efficient power allocation and backhaul esign in heterogeneous small cell networks,” Proceedings of IEEE International Conference on Communications (ICC 2016), Kuala Lampur, Malaysia, May 2327, 2016.
 [25] S. Cui, A. Goldsmith, and A. Bahai, “Energyconstrained modulation optimization,” IEEE Trans. Wireless Commun., vol. 4, no. 5, pp. 2349–2360, Sep. 2005.
 [26] A. Y. Wang, S. Chao, C. G. Sodini, and A. P. Chandrakasan, “Energy efficient modulation and MAC for asymmetric RF microsensor system,” Proc. Int. Symp. Low Power Electronics Design, Huntington Beach, CA, pp. 106–111, Aug. 2001.
 [27] S. Cui, A. Goldsmith, and A. Bahai, “Energyefficiency of MIMO and cooperative MIMO techniques in sensor networks,” IEEE J. Sel. Areas Commun., vol. 22, no. 6, pp. 1089–1098, Aug. 2004.
 [28] S. Verdu, “Spectral efficiency in the wideband regime,” IEEE Trans. Inf. Theory., vol. 48, no. 6, pp. 1319–1343, June 2002.
 [29] F. Meshkati, H. V. Poor, S. C. Schwartz, and N. B. Mandayam, “An energyefficient approach to power control and receiver design in wireless networks,” IEEE Trans. Commun., vol. 53, no. 11, pp. 1885–1894, Nov. 2005.
 [30] R. G. Gallager, “Power limited channels: Coding, multiaccess, and spread spectrum,” in Proc. Conf. Inf. Sci. Syst., vol. 1, Mar. 1988.
 [31] D. Goodman and N. Mandayam, “Power control for wireless data,” IEEE Wireless Commun., vol. 7, no. 2, pp. 48–54, Apr. 2000.
 [32] N. Feng, S. C. Mau, and N. B. Mandayam, “Pricing and power control for joint networkcentric and usercentric radio resource management,” IEEE Trans. Commun., vol. 52, no. 9, pp. 1547–1557, Sep. 2004.
 [33] N. Wang, E. Hossain, and V. K. Bhargava, “Joint downlink cell association and bandwidth allocation for wireless backhauling in twotier HetNets with largescale antenna arrays,” IEEE Trans. Wireless Commun., vol. 15, no. 5, pp. 3251–3268, May 2016.
 [34] D. Bethanabhotla, O. Y. Bursalioglu, H. C. Papadopoulos, and G. Caire, “User association and load balancing for cellular massive MIMO,” Inf. Theory Appl. Workshop (ITA), Feb. 2014, pp. 1–10.
 [35] E. Wolfstetter, Topics in Microeconomics: Industrial Organization, Auctions, and Incentives., Cambridge University Press, 1999.
 [36] Further Advancements for EUTRA, Physical Layer Aspects.,3GPP Std.TR 36.814 v9.0.0, 2010.