Spatial Topology Adjustment for Minimizing Multicell Network Power Consumption
Abstract
While the deployment of base stations (BSs) becomes increasingly dense in order to accommodate the growth in traffic demand, these BSs may be underutilized during most hours except peak hours. Accordingly, the deactivation of these underutilized BSs is regarded as the key to reducing network power consumption; however, the remaining active BSs should increase their transmit power in order to fill network coverage holes that result from BS switching off. This paper investigates the optimal balance between such beneficial and harmful effects of BS switching off in terms of minimizing the network power consumption, through comprehensively considering the spatial BS distribution, BS transmit power, BS power consumption behaviors, radio propagation environments, and frequency reuse. When BSs are deployed according to a homogeneous Poisson point process, the suboptimal and approximated design problems are formulated as geometric programming and the solutions lead to insightful design principles for the key design parameters including the spatial density, transmit power, and frequency reuse of remaining active BSs. The numerical results demonstrate that these solutions are very close to the optimal balances.
I Introduction
Recently, the wide dissemination of smart devices has accelerated the wireless traffic demand, and in order to accommodate this significant growth, the deployment of base stations (BSs) in wireless cellular networks continues to become more and more dense. This growth further aggravates concerns about the everincreasing energy consumption and carbon footprint [1], when considering that the power consumption of BSs in cellular networks is 60–80% of the operators’ power consumption [2] and the energy provision for the BSs is to 50% of the total operational cost [3]. Therefore, energy savings associated with the BSs are a crucial issue for decreasing not only harmful greenhouse gas emissions but also operational expenditures (OPEX). The efforts to reduce BS energy consumption are being undertaken at various levels [1, 3]: the energy source level (e.g., adopting renewable energy resources [4]), component level (e.g., improving the power amplifier efficiency [5]), link and protocol level (e.g., discontinuous transmission for long standby BSs [6]), and network topology level (e.g., heterogeneous network deployments with optimal balances of macro, micro, pico, femtocells, and relay stations [7, 8]).
In wireless cellular networks, BS deployment is typically designed to accommodate peak time traffic (e.g., evening hours) that is up to ten times higher than that of the offpeak periods (e.g., late night hours) and their traffic volume varies temporally in a comparatively regular pattern throughout the day [9, 3]. Accordingly, it is likely that BSs are underutilized during most hours except peak hours. However, a BS has poor energy efficiency particularly in low load situations; even at zero load, the direct current (DC) power consumption of a BS remains at approximately 50% of the peak power [10]. In this regard, the entire or partial deactivation of underutilized BSs could be a key to saving energy [1, 11, 12]. However, the BS switching off technology creates a coverage deficit issue due to the deactivated BSs [13, 14]. In order to fill this coverage deficit, the service area of the remaining active BSs must be expanded. Hence, the BS deactivation should be applied carefully through comprehensively considering both the energy saving effects of underutilized BSs and the increased burden on the remaining BSs that is required in order to maintain service quality. This paper investigates the minimization of the total power consumption of BSs in multicell networks through quantifying and balancing a tradeoff between the beneficial and harmful effects of BS switching off. The optimal solution of the problems will enable the derivation of the inherent design principles that depict the relationship among the BS power consumption model, radio environments, and service constraints.
Ia Related Work
Many studies have proposed various algorithms for reducing access network power consumption by adapting the number of active BSs to the dynamic traffic demands and, more specifically, for determining when and which BSs should be activated and how the remaining active BSs should expand the physical coverage and accommodate the current network load in order to maintain service quality, e.g., [15, 16, 17, 18, 19]. The suggestions in these studies have primarily focused on the design of specific algorithms and they have verified that the algorithms perform well for given BS power consumption models and radio environments regarding the minimization of the power consumption of the total BSs in a network via simulations. However, these studies have not provided comprehensive and insightful understandings about BS switching off because the performance of such specific algorithms depends heavily on the BS power consumption properties (e.g., how the transmit power contributes to the total power consumption and at which level the BSs should be deactivated) and the radio environments (e.g., the radio propagation model).
Several studies have investigated partial or entire BS switching off based on the practical BS power consumption behavior. In general, the power consumption behavior of an active BS can be expressed as the sum of the transmit powerdependent or loaddependent part due to the radio transmission and the constant part remaining due to cooling, power supply, and monitoring [20]. The work of [10] introduced the discontinuous transmission (DTX) of a BS, which was motivated by an affine BS power consumption model, and it evaluated the performance based on the third generation partnership project longterm evolution (3GPP LTE) system parameters. The work of [21] considered a practical BS power consumption model and addressed the effect of this model on the algorithm design and performance. In contrast, in [22], it was claimed that when adjusting the cell size of the remaining active BSs, the optimal cell size depended significantly on the amount of fixed power consumption, i.e., a high fixed power consumption in the active BSs resulted in a reduction in the number of active BSs and an increasing in the cell size in terms of BS energy saving. Recently, the work of [23] introduced the joint optimization of the BS transmission power reduction, DTX, and resource allocation through considering an affine BS power consumption model. However, these studies only rely on extensive simulations for the performance assessment of their algorithms.
Due to the difficulty of expressing the effect of intercell interference in wireless multicell networks, the performance of BS switching off algorithms has primarily been evaluated through simulation. The work presented in [24] mathematically analyzed the outage probability of the user equipment (UE) at the worstcase location, in order to quantify the effect of the cooperative coverage extension of the remaining active BSs. In that study, the worstcase location was considered to exclude the impact of the intercell interference, but the intercell interference is a key factor that system designers should not overlook. The studies presented in [25, 26] mathematically analyzed the effect of base station deactivation through statistically modeling the location of BSs as a homogeneous Poisson point process (PPP). The study in [25] analyzed the adverse effect of BS switching off on the uplink power consumption but only considered the intercell interference for a onedimensional cell deployment. The study in [26] examined the optimal BS density under a service outage constraint in interferencelimited homogeneous and heterogeneous cellular networks. However, neither of these studies have comprehensively considered the multicell operational parameters such as the BS switching off, BS power consumption model, BS transmit power, and frequency reuse.
IB Contributions and Organization
This paper analytically investigates the optimal balance between the power saving that results from switchedoff BSs and the load increase of the remaining active BSs in terms of minimizing the multicell network power consumption. The key contributions are highlighted in the following three aspects.
IB1 Comprehensively formulating a spatial topology adjustment problem
Optimization problems are formulated with the objective of minimizing the network power consumption under network coverage and capacity constraints. These design problems are characterized by (i) the network power consumption model incorporating the affine power consumption model of an individual BS and (ii) the network coverage and capacity constraints being expressed as functions of the spatial density (), transmit power (), and frequency reuse () of the remaining active BSs after switching off some BSs. These formulations provide baselines from which to build suboptimal or approximated design problems that enable closed form solutions and useful optimality conditions.
IB2 Deriving the design principles for green multicell planning
The suboptimal and are derived in closed forms through replacing the network coverage constraint with its lower bound (LB) that enable the design problem to be convexified. This solution leads to the design principles for greening multicell networks, where the operation for the BS switching off and transmit power adjustment can be split into four cases according to the BS sleep mode effect, i.e., the difference () between the standby power consumption of an active BS and the power consumption of a sleep mode BS. It is noteworthy that the interval of and the values of and for each case are expressed in terms of the spatial BS distribution, BS transmit power, affine BS power consumption model, required service quality, and radio propagation model.
IB3 Investigating the impact of frequency reuse on green multicell networks
The partial frequency reuse (PFR) that allows each BS to use only part of the spectrum is one of the most common methods for intercell interference mitigation. In order to investigate the effect of the PFR in terms of the network power consumption, the design problem to jointly optimize in addition to and is formulated, where denotes the number of frequency bands for the PRF. The design problem is reexpressed as geometric programming (GP) and the necessary conditions for the optimality are derived. These conditions result in insightful design principles, e.g., that indicates the spatial density of interfering BSs, i.e., the active BS density using the same subband, remains constant.
The remainder of this paper is organized as follows. Section II describes the system model for the multicell deployment and BS switching operation, and it defines the spatial performance metrics for the design problem formulation. Sections III formulates and solves the design problems for optimizing the BS switching off and transmit power adjustment in terms of minimizing the multicell network power consumption. Section IV optimizes the frequency reuse parameter in addition to the BS switching off and transmit power adjustments. Then, Section V discusses the numerical results, and Section VI concludes the paper.
Ii System Model and Performance Metric
This section describes the random multicell deployment and BS switching off models that render their performances mathematically tractable, and it defines the network performance metric used to incorporate the power consumption behaviors of an individual BS.
Iia BS Switching Off and Network Coverage Reduction
Each BS is in either an active or a sleep mode; the former means that the BS serves users within its coverage while the latter signifies that the BS does not transmit or receive signals. One simple method to select deactivated BSs is a random selection that determines their mode in a probabilistic manner. Let () denote the probability that a BS remains active. Each BS enters the sleep mode with probability while staying in the active mode with probability . In order to avoid the coverage deficits that result from sleep mode BSs, the remaining active BSs fill in these coverage deficits through expanding their cell size. There are several methods to expand cell size, e.g., cooperative beamforming and transmit power increasing [21, 24]. This paper considers the scenario where the remaining active BSs expands their service coverage through increasing their transmit power ; thus, the power saving from the deactivated BSs and the transmit power increase of active BSs should be carefully balanced, i.e., that determines the BS mode must be appropriately traded off against in order to minimize the power consumption of the entire BSs while maintaining service quality above a specific level.
In addition to the random selection of switchedoff BSs, this paper considers the random BS deployment model, where BSs are randomly distributed according to a homogeneous PPP with a density of . This random BS model provides the cumulative distribution function (cdf) of signal to interference plus noise ratio (SINR) that is not only mathematically tractable but also tracks the performance of an actual BS deployment as accurately as the grid model [27]. Accordingly, it is expected that the performance of green multicell networks that randomly select switchedoff BSs among randomly deployed BSs would be analytically tractable, which is the primary reason why this paper adopts these random models. More sophisticated methods, e.g., network adaptation accurately tracking spatial distribution of network loads [15, 16, 17, 18, 19], provide better performance, but they are difficult to use to quantitatively reveal the relationship among a variety of inherent properties of multicell networks, such as the BS power consumption behavior, radio propagation environment, and networkwide power consumption of BSs. Fig. 1 presents the coverage models of a cellular network with full density and with reduced density in random cellular networks. This figure depicts that the remaining active BSs expand their coverage due to the switching off of neighbor BSs. It is assumed that each user is served by the nearest active BS.
In order to assess the coverage served by individual BSs more quantitatively, this paper addresses the coverage probability [27], which is defined as the probability that a user experiences an SINR above a specified level, and this probability can be used as the key metric for expressing the network coverage constraint in minimizing the networkwide energy consumption of BSs. Let denote a random variable that represents the SINR at which a user receives a downlink signal from the nearest BS. The coverage probability can be expressed as where is a given value that determines the minimum link quality. In order to address the network coverage more rigorously, the uplink signal quality as well as the downlink signal quality should be considered. However, in general, users select their serving BS through measuring the downlink signal quality from the BSs, e.g., the reference signal received power (RSRP), and the BS power consumption is dominated by the transmit power radiated from the BS antennas. For this reason, this paper only focuses on the downlink coverage probability in terms of network coverage constraints, and the consideration of uplink coverage remains as future work. The study in [27] derived a simple and tractable form expression for the downlink coverage probability in multicell networks using stochastic geometry, and this result encompasses the effect of various operational factors, such as the general path loss exponent, spatial density of BSs, and BS transmit power. Multicell networks that use the same frequency among cells suffer from low performance for cell boundary users due to severe intercell interference, e.g., is below 0.6 when dB [27]. A common method of resolving this intercell interference problem is to reduce the number of interfering BSs by allocating part of the entire frequency bands to each BS. Planned and dynamic frequency reuse methods contribute to further enhancing the performance of the cell boundary users, but their performance has primarily been evaluated through simulations. The study in [27] analytically quantified the performance enhancement by simply modeling this frequency reuse as a random frequency band allocation. This random frequency reuse provides a lower performance compared with more sophisticated methods, but its analytical performance results facilitate the mathematical formulation for designing various operations of multicell networks, e.g., [28]. The following result provides the coverage probability in multicell networks where BSs are deployed according to a homogeneous PPP with and they transmit their downlink signals with transmit power using one frequency band randomly chosen among frequency bands.
Lemma II.1
If one of frequency bands is randomly allocated to each cell, the coverage probability for the random multicell networks is given by
(1) 
where denotes the received signal to noise ratio (SNR) at unit distance, which depends on and , and . In particular, when , this coverage probability has the closed form expression, as follows:
(2) 
where and denotes the function.
IiB Individual BS Power Consumption and Multicell Area Power Consumption
A BS consists of various components for power amplifying, air conditioning, signal processing and power supply; each component has an individual power consumption behavior that may be dependent on the transmit power or bandwidth [3]. Several studies have built a parameterized BS power consumption model that is expressed as an affine function of the BS transmit power [12, 20, 29], and the work of [20] verified that this model approximates its underlying complex model well. In this model, it is noteworthy that the BS power consumption in the active mode almost linearly increases with the transmit power radiated at the BS antenna and a certain level of power is consumed even in the sleep mode. In contrast, it takes time to reactivate the BSs that are in the sleep mode, and this delay is dependent on which functional components are deactivated in the sleep mode. In this sense, the power consumption in the sleep mode, which is denoted by , depends on at which level this mode is designed under the consideration of tradeoff between the amount of power saving and the reactivation delay of the components, e.g., switching off the entire BS or only core components such as the power amplifier and RF. Herein, indicates an ideal BS switching off. Let denote the instantaneous power consumption of BS located at . This paper assumes that each BS has a single sector and a single RF chain. Therefore, can be expressed as follows:
(3) 
where denotes the standby power consumption of the active BS that results from the components with transmit power independent power consumption characteristics, e.g., air conditioning and power supply, and denotes the maximum BS transmit power. Note that as a result of the sleep mode benefit.
In order to assess the power consumption of the entire access network, this paper considers the average power consumption of BSs per unit area, i.e., the area power consumption (APC), as the primary performance metric for capturing the energy efficiency of multicell networks. When BSs are distributed according to a homogeneous PPP with density and they enter sleep mode with probability , the spatial distribution of the remaining active BSs can be modeled as the independent thinning of homogeneous PPP with retention probability [30]. That is, original point process can be considered as the superposition of two independent homogeneous PPPs for active and sleep mode BSs, which are denoted by with density and with density , respectively. Therefore, the APC is given as follows:
(4) 
where denotes the expectation of random variable , (a) follows from the Campbell theorem, the stationarity of a homogeneous PPP [30], and (3), and (b) follows from the introduction of new variables, i.e., and . It is worth noting that the APC in (4) naturally connects the parameterized BS power consumption model given in (3) with the networkwide BS power consumption. The next sections quantitatively reveal how the BS power consumption property affects the minimization of this APC.
Iii Area Power Consumption Minimization
This section formulates and solves the problem for the balance between the spatial density of active BSs and their transmit power, i.e., and , in terms of minimizing the APC defined in (4). To begin with, in order to clarify the relationship between optimal and , this section only considers the universal frequency reuse (UFR) case, i.e., . The case of is addressed in Section IV.
In formulating the problem, the definition of the service quality constraint remains ambiguous. This paper addresses the service quality through discriminating between two factors: network coverage and capacity.

Network coverage constraint: The network should be able to provide users with the spatial coverage with probability not less than , i.e., . In random multicell networks, the complementary cumulative distribution function (ccdf) of SINR, i.e., the coverage probability, is given by (1), which is a function of the spatial density of the active BSs and received SNR at a unit distance. The received SNR at the unit distance denoted by depends on the BS transmit power per physical resource block (PRB). Let denote the BS transmit power per PRB. This paper assumes that is the same for all PRBs. In addition, the problem concentrates on the worstcase coverage scenario where all active BSs use all available PRBs. Under this scenario, the total transmit power of a BS is expressed as , where denotes the total number of available PRBs^{1}^{1}1Practical systems require much more complex transmit power calculation due to the transmit power boosting for control channels and reference signals and the transmit power control for link adaptation. However, this paper only considers this simple transmit power model for analytical tractability.. This transmit power model pursues the reduction of the possible APC rather than an exact APC. When considering that it is almost impossible to accurately estimate the instantaneous transmit power of individual BSs, it is sensible to use this simple model that facilitates the derivation of analytically tractable results. When using this transmit power model, the received SNR at unit distance, i.e., becomes where denotes the noise power per PRB, represents the path loss at a unit distance, and . Therefore, when , the coverage probability given by (1), i.e., , can be reexpressed as follows:
(5) 
Network capacity constraint: The aggregate network throughput in a wireless cellular system almost linearly increases with the BS density [31]. That is, the spatial density of BSs is a dominant factor in determining the network capacity. Hence, it is sensible to roughly express the network capacity required for a certain period in terms of the spatial density of BSs that can accommodate a network load. Herein, there are two more justifications for the consideration of the BS density as capacity requirement during offpeak period. The first justification is the difficulty of obtaining an accurate estimate of the amount of required network capacity. The network load depends significantly on the variations in the number of users, user locations, and traffic patterns [9, 3]; thus, only a rough estimation is available. The second justification is that operators typically deploy sufficient BSs to accommodate the peak time traffic and the BS switching off for energy saving is primarily applied during low loads, i.e., off peak times. This implies that the green multicell networks focused on in this paper are coverage limited rather than capacity limited; thus, this rough expression of the network capacity is sufficient to reflect the capacity requirement in low load situations. Based on these justifications, this paper imposes a capacity constraint expressed as in addition to the coverage constraint. This network capacity constraint corresponds to in terms of the average area spectral efficiency (ASE), or in terms of the ASE associated with the percell spectral efficiency that can be spatially guaranteed with probability , when considering the above network coverage constraint^{2}^{2}2The data rate of stateofthe art communication systems, e.g., the 3GPP LTE system, can be approximated simply as an attenuated form of the Shannon bound, i.e, bps/Hz where is an attenuation factor smaller than one [32]. This paper considers an ideal link of ..
The APC minimization problem under these two constraints can be formulated as follows:
(6a)  
subject to  (6b)  
(6c)  
(6d) 
Note that problem (6) minimizes the objective function given by (6a) because in (4) is constant. The left hand side (lhs) of constraint (6b), which is given in (5), includes the integration and this term makes it difficult to solve problem (6) and express its optimal solution as a simple form. The following remark provides useful insight about the above optimization problem.
Remark III.1
Assume that problem (6) is feasible. Then,
(i) if and belong to the feasible set of problem (6) and ^{3}^{3}3Herein, denotes that and ., ;
(ii) if denotes the optimal solution of problem (6), .
(i) The increase in the coverage probability with and is obvious because the effect of the noise power deceases as and increase. For analytical proof, see Appendix A.
(ii) Assume that . From the increasing property of with respect to and , there is such that and ; thus, is also feasible. The APC for is smaller than that of , because objective function (6a) decreases with for given . This contradicts that is the optimal solution. Define and ^{4}^{4}4This section only considers the case of but the definitions of and can be extended into those of . In these definitions, in (1) is more definitely expressed in (13). as
(7)  
(8) 
Then, (resp. ) can be easily found via the bisection method because is a monotonically increasing function with respect to for fixed (resp. for fixed ). Remark III.1 implies that when one variable is given, the other variable can be readily optimized; that is, for given where (resp. for given where ), the optimal (resp. optimal ) becomes if (resp. if ).
However, the original problem in (6) should be jointly optimized for and , and it is not easy to analytically express the solution of this joint optimization problem. Note that even when , which has a closed form expression of given by (2), the analytical expression for the solution of problem (6) is not easily derived. The next subsection reformulates the nonconvex problem in (6) as GP through recasting coverage constraint (6b) and it derives the optimal active BS density and transmit power of the reformulated problem.
Iiia Formulation and Solution of Tightened Design Problem
The coverage probability in (5), which is an integral function, needs to be expressed more elegantly in order to find a good balance between and . From the property that for ^{5}^{5}5Let for . Because for and , when .,
(9) 
where denotes the Gamma function. The design problem typically aims to maintain high coverage probability, i.e., as close to one as possible, even when some BSs are switched off. In this regard, the design interest is placed on a small noise scenario where the noise effect is limited, i.e., is small for given . It is worth noting that the lower bound of in (9) gets tighter as decreases and the equality holds when . In fact, for is upper bounded by the one for ; thus, for , . In this sense, it is expected that the LB of (9) replaces well. Note that (9) is similar to equation (9) in [27]. However, therein, the coverage probability was approximated for a small but nonzero noise power by applying where denotes the higher order terms of , rather than being lower bounded. From (9), the new coverage constraint that tightens (6b) is given as follows:
(10) 
The above condition follows from (9) and therefore condition (10) is a sufficient condition for meeting constraint (6b). Because as mentioned above, representing coverage requirement needs to be less than for the sake of feasibility. This means that . It is worth noting that constraint (10) exhibits that the BS density has the equivalent effect as the transmit power to the power of in terms of the coverage probability. In twodimensional Euclidean space, scaling up the BS density by a factor of is equivalent to scaling down the distance between a UE and its serving BS averagely by a factor of ; thus, this distance reduction corresponds to increasing transmit power to in terms of the coverage probability. In this sense, constraint (10) makes sense.
Now, by using this new condition instead of constraint (6b) in problem (6), a new tightened optimization problem is yielded, as follows:
(11a)  
subject to  (11b)  
(11c)  
(11d) 
In order to guarantee the feasibility of problem (11), it is assumed that sufficient condition (10) is met when all the BSs transmit with the maximum transmit power, i.e., . Note that the optimal value of problem (11) provides the upper bound of that of problem (6).
Proposition III.1
Assume that and . When denotes the optimal solution of problem (11),
(12) 
See Appendix B. It is interesting that the operation for the spatial topology adjustment is split into four cases, according to the difference between the standby power consumption in the active mode and the power consumption in the sleep mode, i.e., that expresses the amount of power consumption saved through the BS switching off. This supports that the operation depends on the effect of the BS sleep mode. The result of Proposition III.1 can be interpreted as follows:

The case of indicates the environment where the effect of the sleep mode is not dominant and accordingly the network minimizes the power consumption through reducing the transmit power rather than through switching off BSs.

The case of denotes the instance that requires a tradeoff of the power reduction from the switchedoff BSs against the increase in the transmit power of the remaining active BSs in a coverage limited scenario. The result reveals the interesting outcome that is linearly proportional to with a slope of and is set to the solution of . This causes the setting of transmit power for balancing and to be very handy.

The case of and represents the situation where the BS switching off enables a substantial energy saving effect due to a large standby consumption in the active mode, a very small sleep mode power consumption, or both. Hence, the density of the active BSs is maintained as small as possible. denotes that the network is coverage limited when is minimized; hence, this result determines through maximizing .

The case of and designates that the network benefits from the BS switching off similar to that in (iii) but the active BSs are too sparse to accommodate the expected network traffic. Therefore, the density of the active BSs is at least equal to and the transmit power should be set to meet the coverage constraint.
In summary, the first case only reduces the BS transmit power, the third and fourth cases minimize the BS density, and the second case appropriately balances the BS density and transmit power.
The solution in (12) also reveals the effect of the radio propagation, i.e., , in addition to the BS power consumption behaviors. This affects both the interval of that determines the BS operation and the value of for the APC minimization, and its effect depends on other environmental factors, e.g., the noise power and target SINR.
IiiB Improvement of Suboptimal Solution
Recall that constraint (10) is the sufficient condition of original constraint (6b). Therefore, based on in (12) is at least larger than as long as . Because is monotonically increasing with respect to for given from Remark III.1, where defined in (7) can be readily obtained via the bisection method. Hence, is a feasible solution of original problem (6) and provides the APC less than . That is, based on in (12), a better solution can be readily yielded.
Iv Frequency Reuse and Green Multicell Networks
The partial frequency reuse (PFR) that allows each BS to use only part of the spectrum is a common and useful method for mitigating intercell interference problems. When the entire spectrum is divided into frequency bands, one of which an individual BS uses, it is known that as increases, the outercell user performance is improved while the average rate of each cell may be reduced [27]. Meanwhile, the question about the impact of this PFR on the power consumption in multicell networks is an interesting topic. This section formulates and solves the spatial topology design problem for green multicell networks through considering the number of frequency bands in addition to the BS density and transmit power as design variables.
In order to reflect the impact of ^{6}^{6}6 is a positive integer, but this paper allows to be real for analytical convenience., the coverage and capacity constraints addressed in Section III must be appropriately modified.

Network coverage constraint: The PFR does not only decrease the density of the downlink interferers by a factor of , but it also enables the boosting of the SNR on one PRB by concentrating the BS transmit power on an available frequency band. This impact is described by the results in Lemma II.1, and, herein, is set to , in pursuit of the possible APC minimization justified in Section III. That is,
(13) for is equal to in (5); thus, encompasses the coverage probability for the UFR as well as the PFR.

Network capacity constraint: For the ease of formulation, the capacity constraint is expressed in terms of the BS density, which is similar to that in Section III. Because the network capacity is proportional to the bandwidth as well as the BS density for a given SINR distribution, this constraint can be simply cast into , where in the denominator reflects the reduction of the available bandwidth due to the PFR. However, this model tends to underestimate the average ASE given by , where the distribution of is also improved as increases. In contrast, it abstracts the performance well in terms of the ASE associated with the percell spectral efficiency that can be spatially guaranteed with probability , i.e, for specific constant . Considering the difficulty of its accurate estimation and the bias to the coverage limited scenario, as described in Section III, this simple capacity constraint is imposed on the design problem for green cellular networks.
Iva Formulation and Solution of Tightened or Approximated Design Problems
In order to solve the design problem, the coverage constraint based on coverage probability (13) needs to be recast into a simpler form. Similar to Section III, the coverage probability in (13) can be lower bounded by using , as follows:
(14) 
However, (14) does not provide the relationship among , , and that facilitates the convexification of the design problem. For this reason, another LB or approximation of the coverage probability for PFR is considered based on the LB of (14), as follows:
(15) 
(16) 
In contrast, only when , the LB that enables to convexify the design problem can also be obtained from , as follows:
(17) 
Note that from , and this upper bound is equal to for . In addition, the LB in (17) becomes tighter and eventually approaches as increases. In this regard, in (17) approximates the original coverage probability given by very well for close to one. Moreover, it can be shown that when , by comparing (17) with (14).
In summary, the four methods for recasting the original coverage probability were introduced, and it is worth noting that the latter three equations, i.e., , , and , enables to convexify the design problem. In addition, and are lower bounds for while is not. Their relations are summarized in Table I.
Scenario  Relationship 

,  
,  
,  
, 
New coverage constraints based on , , and for replacing original coverage probability (13) can be reexpressed as the following common form:
(18) 
where , , and are nonnegative constants and given as follows:
(19)  
(20)  
(21) 
It is sensible to assume that , because close to one is considered in this design problem. This assumption guarantees that all of , , and in (19), (20), and (21) are positive.
In problem (6) for , i.e., the UFR, by imposing new coverage constraint (18) and capacity constraint instead of constraints (6b) and (6c), the design problem for the PFR, i.e., , can be formulated as follows:
(22a)  
subject to  (22b)  
(22c)  
(22d)  
(22e)  
(22f) 
Note that objective function (22a) and constraint (22b) are posynomial functions. That is, problem (22) is GP in posynomial form [33]. GP can be transformed to a convex problem, which can be solved by the interiorpoint method [33] or available solver [34]. However, the solution based on these iterative methods does not reveal the analytical relationship among the design parameters such as , , and . In this regard, the necessary conditions for the optimal solution of problem (22) are derived.
Proposition IV.1
Assume that , and let denote the optimal solution of problem (22). These optimal variables are related as follows:
(23) 
where . To be more specific, satisfies one of the following conditions:

, , .

, , , where .

, , , where .

, , .
In (ii), is given by with coefficients (19) for in (15) (resp. (20) for in (16)) when , whereas it can also be given by with coefficients (21) for in (17) only when , where
(24)  
(25) 
See Appendix C. The interpretation of the necessary conditions (i)–(iv) in Proposition IV.1 corresponds to that of the four cases in (12) for the UFR, even though the results for the PFR do not explicitly offer the conditions that are expressed in terms of the interval of . One notable difference is regarding the additional variable , and this variable is conditioned on , like (23). Note that the density of the interfering BSs is expressed as in the PFR. Therefore, condition signifies that the density of the interfering BSs is always maintained as a fixed value of in terms of minimizing the APC in problem (22). Conditions (i), (iii), and (iv) in Proposition IV.1 correspond to the extreme cases where one of constraints (22d), (22e), and (22f) is active. In contrast, condition (ii) is relevant to the balance , , and , and the selection of these design variables begins with solving in (24) or in (25).
IvB Derivation of Suboptimal Solutions
The candidate solution for problem (22) can be obtained through comparing the APC for values derived from conditions (i)–(iv) in Proposition IV.1 and then choosing the feasible one with the least APC among them. With a slight misuse of notation, this candidate solution is denoted as . Recall that problem (22) uses the LB or approximation instead of the original coverage probability for its coverage constraint. Furthermore, it can be readily proven that the two properties in Remark III.1 also hold for in (13) by replacing with given . In this sense, similar to the method presented in Section IIIB, in order to return to the original coverage constraint, new solution is derived from candidate solution , where is defined in (7). For the design problem based on the LBs of , i.e., and , this new solution reduces the APC compared with that from . In contrast, for the design problem based on the approximation of , i.e., , this method renders the solution feasible if while decreasing the APC if .
IvC Design of and for Given
Unlike the previous subsections that considered as a design variable in addition to and , one can consider the optimization of and for given . Given , the design problem can be formulated by replacing and with and , respectively, in problem (6). From (14), the coverage constraint in integral form, i.e., for given , can be tightened by the following constraint in monomial form:
(26) 
As a result, the optimal solution of the tightened design problem for given , i.e., , can be obtained from Proposition III.1 by replacing and with and , respectively. In addition, because in (13) monotonically increases with respect to for given and as addressed in Section IVB, brings the APC less than that of .
Variable  Description  Value 

Total density of deployed BSs  
Minimum required BS density  
Maximum BS transmit power  
Standby power consumption of an active BS  
Slope of affine BS power consumption model  
Path loss at unit distance  
Total noise power over bandwidth  
Normalized noise power over total bandwidth ()  
Path loss exponent  4 or 5  
Minimum required SINR 
V Numerical Results and Discussion
This section evaluates the APC performance of the green multicell networks designed in Sections III and IV, and it examines the effect of the affine BS power consumption model given by (3). The system parameters for the performance evaluation are presented in Table II, and these parameter values are used to obtain results in this section, unless stated otherwise. The parameter values for the BS power consumption and path loss models are taken from [23, 35].
Figs. 2 and 3 examine the coverage probability for the design problems proposed in Sections III and IV, prior to evaluating the APC performance of multicell networks.
Fig. 2 compares the coverage probability and its alternatives when 20% and 50% of the total BSs are functioning. In this evaluation, each BS sets its transmit power to be 10% of . For the exact coverage probability, i.e., , it is observed that the analysis results coincide precisely with the empirical results; thus, the correctness of all analytical results is validated. This figure also demonstrates that the coverage probability is quite well approximated using its alternatives, which become closer for higher BS densities. This results from the noise power becoming increasingly overwhelmed by the aggregate interference as the BS density increases, and , , , and provide better lower bounds and a better approximation for smaller values. More specifically, for , , , and tightly approximate for the entire range of coverage probability, while only approximates well for high coverage probabilities. However, when considering that systems typically support high coverage probabilities, e.g., above 0.8 or 0.9, also remains a tight LB of .
Fig. 3 presents the contour graphs of the coverage constraints based on , , , , and , for high coverage probabilities. Unlike Fig. 2 for fixed and , these contour graphs enable to more precisely examine the relationship between and and compare and its alternatives by focusing on a specific . As summarized in Table I, it is observed that is lower bounded more and more tightly in the order of , , and , whereas may tighten or relax depending on when . Recall that when , does not only enable the convexification of the coverage constraint but also provides tighter LB than that of , but it is only available for . The increase in causes a quicker decay in both the desired and interference signals, and accordingly the impact of the noise power becomes more significant. For this reason, the LBs of in Fig. 3(b) tend to have larger deviations than those in Fig. 3(a).
Fig. 4 demonstrates that the design of and proposed in Section III functions well in terms of minimizing the APC. In the subfigures, the ordinates on the left and right denote the APC values and the ratios of the APC to the maximum APC with , respectively. The curves in this figure depict the APC at the optimal transmit power for given and defined in (7) and they indicate that this optimization decreases the APC to 15–34% of that for . In addition, these results support that the BS power consumption behavior parameterized by has a dominant impact on the optimal topology adjustment. That is, the BS switching off is less effective for small , while it is the key to reducing the APC as increases. In this figure, is normalized by with the value in Table II; thus, where close to zero indicates the quite small effect of the BS sleep mode while indicates an ideal sleep mode with zero BS power consumption. This figure also demonstrates that derived in Proposition III.1 provides an excellent suboptimal solution for the joint optimization of and . For example, when , nearly minimizes the APC and is almost the same as , which further improves as proposed in Section IIIB. In contrast, when , as a result of the harsh propagation loss, the noise power impact increases relatively and therefore has a slight deviation from as shown in Fig. 3(b). Thus, through further optimizing for given , the APC performance can be improved. Note that in Fig. 4, when , the curves according to meet at one point because there is no BS in sleep mode.