1. Introduction
To cope with the mobile traffic explosion, upcoming 5G focuses on the gigabit-scale peak data rate with higher capacity. Deploying hyper-dense small cell networks is a strong candidate to enhance the capacity up to 100-times [
1]. While dense small cells can meet the capacity requirement, total power consumption on the network increases in proportion to the number of small cell base stations (BS) [
2]. In addition, we should consider that network infrastructures take 3% of the world’s electrical energy consumption, and BSs consume about 80% of the energy consumption of mobile network operators (MNOs) [
3,
4]. In light of these illustrations, it is indispensable to enhance energy efficiency. Thus, many efforts are continued in the standards of 3GPP and the research projects of EARTH, Green Touch, 5GrEEn, etc. [
5]. In the perspective of 5G, the authors of [
6] recently suggested a design framework considering both energy efficiency (EE) and spectral efficiency (SE), called EE-SE co-design.
In the literature, many on-going research works investigate green BS operations to enhance energy efficiency. In [
7], the authors proposed a BS turning-off algorithm using the optimal cell coverage that minimizes BS area power consumption. The authors in [
8] proposed a cell zooming technique, such that the BS turning-off operation reduces the power consumption, and other turned-on BSs increase their coverage. In [
9,
10,
11], the on/off methods based on flow-level dynamics are presented. The authors in [
9] solved the problem encompassing dynamic BS operation and the related problem of user association. Related to [
9], the authors in [
10] presented the distributed threshold-based BS off algorithm based on an overlay network using Delaunay triangulation. The authors in [
11] proposed the energy-efficient user association method by using a population game approach. In [
12,
13,
14], stochastic geometry is used. The authors in [
12] proposed macro BS turning-off strategies in homogeneous and heterogeneous networks and showed that the deployment of small cells can increase energy efficiency, but the energy efficiency gain is saturated when small cell density is increased. In a two-tier heterogeneous network, the authors in [
13] proposed the repulsive cell activation scheme that activates small cells according to user density and offloading macro BS traffic. A deployment strategy considering the density and transmission power of BSs is proposed in [
14]. The authors in [
15] considered a separation architecture substituting the conventional macro BS into a coverage-providing BS and multiple small cells in order to increase energy efficiency. The authors derived the optimal intensity of small cells for the optimal deployment. In [
16], the energy-saving problem is investigated while considering multiple sleep base station modes. Furthermore, a survey in [
17] presents various existing green cellular techniques to show that energy saving can be achieved by turning on or off BSs.
In this paper, we study a small cell BS on/off operation to enhance energy efficiency. Our approach is applicable to various types of small cells, e.g., pico BS, femto BS or even Wi-Fi BS. Hereafter, we use the terminology access point (AP) instead of BS to emphasize the small cell BS. Since we assume that users can directly deploy small cell APs, intractable randomness is engendered. In the densely-deployed small cell networks, the AP on/off operation brings changes to the network environment. While flow-level analysis considers the time-averaged information [
9,
10,
11], we mainly focus on the state information of APs because a more frequent and dynamic AP on/off operation is expected. Thus, to know the state, devices’ channel measurement is crucial. These assumptions motivate us to take a new approach for the dynamic AP on/off operation using user equipment (UE) feedbacks. Our approach basically exploits belief propagation (BP) to solve an optimization problem on a factor graph. This paper extends the results of its conference version [
18]. The main difference is that we further consider the case of ultra-dense networks and investigate the impact of the proposed algorithms under the extremely dense scenario. For example, we consider the case of 100 APs with 50 UEs for dense network, as well as 10,000 APs with 1000 UEs for ultra-dense networks. In doing this, we have modified the conventional BP algorithm to the approximated BP and also verify that the proposed online algorithms of DANCE perform well even under extreme densification. Furthermore, in
Section 4.3, we provide an operational procedure and a framework to practically implement the online algorithms of DANCE.
The main contributions of this paper are summarized as follows. We apply the BP optimization framework to the AP on/off operation. To target the small cell environment, a model exploiting UE feedbacks is proposed. Then, by constructing a factor graph representing APs and UEs, the optimization problem is solved by the BP optimization framework. Since BP is an offline algorithm, however, we also propose algorithms that can run in an online manner where these algorithms are named DANCE (Device-Assisted Networking for Cellular grEening) that is a collection of low complexity algorithms inspired by BP. As the complexity of the BP approach can increase with the number of BSs and UEs, DANCE algorithms can be used in dense networks. Then, we evaluate our algorithms by performing numerical simulations. Our extensive simulation results show that BP and DANCE increase energy efficiency significantly. For instance, in a small-sized network, BP increases energy efficiency 129% compared to the baseline where all APs are always on. Furthermore, in ultra-dense networks, DANCE algorithms increase energy efficiency from 25-times up to 44-times compared to the baseline in the case of ultra-dense networks.
The remainder of this paper is organized as follows. In
Section 2, the assumptions of the system model and the problem formulation are presented.
Section 3 is devoted to the graphical modeling of our framework using the BP algorithm. In
Section 4, we propose a family of online algorithms. In
Section 5, the performances of the proposed algorithms are demonstrated with extensive simulations. Finally, we conclude the paper in
Section 6.
3. Belief Propagation
We adapt BP to solve our optimization problem expecting that a simple, but practical algorithm can provide an approximated solution of our problem. BP originally proposed by Pearl [
21] is a message-passing algorithm widely used in inference problems. BP has successfully been demonstrated in many applications, such as error-correcting codes or artificial intelligence. By using BP, the marginals of a joint probability distribution can be calculated with iterations of message passing on graphical models including Bayesian networks and Markov random fields. We specifically focus on running BP on a factor graph, which also can be readily converted into other forms of graphical models, like Bayesian networks. For details, please refer to [
22], which shows a general overview of BP, and [
23,
24,
25], which provide explanations about BP on a factor graph. The optimization method using BP was used by the authors of [
26]. The authors applied BP to inter-cell interference coordination in a femtocell network to solve the problems approximately.
By using the elements of the AP on/off state indication vector, we define a set of discrete random variables . Then, can be seen as a Bernoulli random variable denoting the turned-on or turned-off state. The realizations of the random variable is denoted by as in (1). Let us use the abbreviated notation denoting . By using these definitions above, the joint probability distribution is expressed in a shorter form .
To define the joint probability distribution
, the concept of Boltzmann’s law in thermal equilibrium from statistical mechanics is used. Our framework can be seen as a system that is assumed to have
J particles, and a state of each particle is denoted by
. By substituting the negative of free energy part of Boltzmann’s law into our global function,
is defined as:
where
z is a normalization constant and
τ denotes the temperature. Also, by changing the exponential of sum to product of exponentials,
can be shown as
when BP is applied to the non-physical system, the temperature can be seen as a constant parameter appropriately chosen.
3.1. Factor Graph
Observing (7) and (8), we see that (7) including the global function
is factorized into the product of the local functions
. Thus, (8) can be shown in a factor graph. A factor graph
is a bipartite graph with vertices
V and edges
E. Vertices, i.e., nodes of the graph, correspond to APs and UEs. There are two kinds of nodes: the variable nodes and the factor nodes. The variable nodes representing APs, shown as circle in the graph, are referred to
. The factor nodes denoting UEs, shown as the square in the graph of
Figure 2, are referred to
, where the local function is included. Edges mean the channels between APs and UEs. To solve our optimization problem, we exploit the constructed factor graph in BP message passing.
3.2. Belief Propagation Optimization Algorithm
Following the results of [
27], when
,
concentrates around the maxima of
, and then, we have:
where
, and
means the marginal expectation of the random variable
. In principal, standard BP can be seen as a process to calculate the estimated marginal probability distribution with respect to the random variable
. Thus, when the marginal probability distribution of
is estimated, an approximated solution of the optimization problem can be found.
Definition 2 (Belief).
The belief is the estimated marginal probability distribution of with respect to .
Thus, we leverage the BP algorithm to find
. On the constructed factor graph, belief messages are passed through the edges between APs and UEs. After the belief message is repeatedly passed between APs and UEs, the set of turned-on APs is determined according to the estimated marginal probability distribution. It is well known that when a factor graph is acyclic,
converges to the true marginal probability distribution of
with respect to
[
26]. However, as shown in
Figure 2, the constructed factor graph usually has cycles, e.g., AP 2, UE 3, AP
j, UE
i, and AP 2. In this case, BP yields an approximated result; thus, we compare our simulation results with the optimal value obtained from the exhaustive search in
Section 5.
Now, we present the algorithm below.
Initialization. At time
, messages
for
and
are set to have an arbitrary distribution. For instance, we use the uniform distribution, such as:
UE update. After the initialization, UE
generates and sends the belief messages to its neighboring APs
where
denotes the neighbor APs of UE
i. At
, the message update is defined by:
Therefore, when the expectation (11) is calculated over independent
,
, the message update can be presented by
where for given
. Furthermore, in (12)
is calculated for given
. Since (12) has both a sum and a product, BP is called the sum-product algorithm. Equation (12) is the product of the factor with belief messages from AP
,
, to UE
i, marginalized over all random variables, except
. To calculate
, at
, we assume that
,
. At
,
,
, is initialized to one and incremented by the number of associated UEs where each UE is associated with AP
.
AP update. AP
generates and sends the belief messages to all neighboring UE
, where
denotes the neighbor UEs of AP
j. The message is defined by:
where
is the normalizing constant making
. The belief message (13) is the product of the belief messages from UE
, all neighboring factor nodes, except
i. Then, the steps ‘UE update’ and ‘AP update’ are repeated until the number of iteration is satisfied.
Decision-making. After finishing the iteration procedure, the estimated marginal probability distribution with respect to
is finalized as:
where
is the normalizing constant to make
as a probabilistic distribution.
Then, as a final step, to describe our association rule
S, we define an integer-valued function
that determines the position of the
j-th element when the elements of the input vector are sorted in increasing order. Then, we deliberately use the following rank-based association rule of UE
i given by:
where
,
, and
η is a weighting parameter considering the tradeoff between the channel condition and the estimated marginal probability from BP. The intuition of (15) is as follows. If each UE is associated with AP
j solely considering energy efficiency, then most UEs are associated with a few APs having a higher value of
. High
may indicate that there is a global consensus that energy efficiency can be increased when the BS
j is turned on. However, connecting too many UEs to one AP can degrade the data rate, which may decrease energy efficiency. In this sense, user association cannot be determined solely relying on
. Thus, we exploit
to determine the user association, so each UE can decide its association along with the global perspective of
and also personal perspective
to preserve its service quality.
Our user association rule can be seen as an extension of the convention user association rule that is based on SINR. For example, if , each UE is associated with the AP that provides the highest SINR among APs in , which is the conventional user association rule. On the other hand, if , the UE is associated with the AP j whose is the greatest. It could be possible that SCC determines η to maximize energy efficiency by using adaptive learning methods. After user association is done for all UEs, the APs with no UEs are turned off.
While the BP framework is readily well defined on many graphical models, the drawback of BP framework is its complexity, which is given by from the view point of the UE. Hence, the complexity of the BP could be a barrier to large-scale ultra-dense networks of APs and UEs. Thus, we next propose online algorithms for dense cellular networks in the next section.
4. Online Algorithm: DANCE
We propose a family of online algorithms that are inspired by the BP framework. In BP, the node having a high certainty of belief can influence the final estimated marginal probability distribution. To mimic this behavior, DANCE considers priority to the AP or UE whose impact can be great. Since the ultra-dense network (UDN) environment is considered where APs are densely deployed, DANCE is designed to operate with low computational complexity.
4.1. Feedback from UE
We assume that UEs can send feedback messages as shown in
Figure 1. In the SCC, the collected feedback messages constitute the feedback matrix
. The feedback message can contain merely the connectivity information when using only one bit or can have more detailed information of connectivity, such as the maximum achievable data rate or appropriate adaptive modulation and coding (AMC) when using multiple bits. Thus,
has two variations: a matrix
is when one-bit feedback is used, and
is when
N-bit feedback is used. An element of a matrix,
or
, means that the feedback message from UE
i to AP
j is either one-bit or
N-bit. If a one-bit message is used,
represents only the connectivity. It is assumed that connectivity exists if
is greater than some threshold. If an
N-bit message is used,
represents the achievable data rate or the quantized level of AMC. Examples of the UE feedbacks are illustrated in
Table 1. We assume that the achievable data rate can be represented by using multiple
N bits. In our study,
N-bit message implies the achievable data rate. A study about the number of bits considering the tradeoff between performance and feedback overhead is beyond the scope of this paper and remains as a future work. However, when comparing to the BP algorithm where the repeated message exchanging is required, the DANCE algorithms only require one-time feedback messages, which significantly suppresses overall signaling overhead.
4.2. DANCE Algorithms
DANCE solves the same problem given in (3) to (6). By using matrix-based algorithms, DANCE determines the turned-on/off APs and the user association together. Once the target AP is chosen, then the association is automatically decided on the matrix , so that the result of user association is shown on the matrix . DANCE includes three algorithms: AP-first, UE-first and the proximity-ON (Prox-ON) algorithm.
In DANCE algorithms, two main algorithms, AP-first and UE-first, are designed to choose the most influential AP or UE, respectively. Depending on the type of feedback, both AP-first and UE-first have several variations: AP-N/AP-1 and UE-N/UE-1, where the number implies the length of the feedback bits. Following the categorization, AP-first and UE-first are described in Algorithms 1 and 2, respectively.
Algorithm 1 AP-first |
- 1:
Initialize , , . - 2:
while - 3:
. - 4:
for all - 5:
if , then , - 6:
else . - 7:
end for - 8:
Eliminate the column of . - 9:
Eliminate rows of where . - 10:
end while
|
Algorithm 2 UE-first |
- 1:
Initialize , , . - 2:
while - 3:
. - 4:
where . - 5:
for all - 6:
if , then , - 7:
else . - 8:
end for - 9:
Eliminate the column of . - 10:
Eliminate rows of where . - 11:
end while
|
The algorithm works as follows. AP-first algorithms, AP-1 or AP-N, firstly choose an AP on the matrix as previously mentioned. AP-1 algorithm chooses the AP that has the greatest number of association-possible UEs. The AP-N algorithm otherwise chooses the AP that has the greatest sum of the maximum achievable data rate or the AMC of each UE depending on whether stands for the data rate or AMC levels. Then, the chosen AP is turned on, and all UEs having connectivity with the chosen AP are associated with the chosen AP. This is repeated until all UEs are associated.
As a counterpart, UE-first algorithms, UE-1 and UE-N, consider UEs foremost. For example, in the UE-1 algorithm, we firstly find the UE having the smallest number of association-possible APs, i.e., potentially the worst case UE is considered first. Then, among the APs that can connect to the chosen UE, we select the AP
where
in order to minimize the number of turned-on APs. Note that UE-first is a generalization of [
28] where only a single bit feedback was assumed.
The other is the proximity-ON algorithm. In this algorithm, each UE associates with the AP that provides the highest SINR, while non-selected APs are turned off. It can be seen as a basic heuristic algorithm in order to determine the set of the turned-off APs. Therefore, among DANCE algorithms, the results of the proximity-ON algorithm can provide baselines to measure performance gains of AP-first and UE-first, respectively.
4.3. Practical Implementation
The message passing procedures for both the turning off and on operation in DANCE are illustrated in
Figure 3. As mentioned in
Section 2, each AP periodically wakes up and broadcasts a beacon signal, and each UE senses the APs’ beacon signals. Information gathered at each UE is sent to SCC. Then, SCC runs one of the DANCE algorithms in order to make a decision for each AP’s on/off status by using the feedback information from UEs. Then, SCC broadcasts the final decision to all APs, and APs follow the direction of SCC.
6. Conclusions
In this paper, we have studied the mechanisms to enhance energy efficiency in the ultra-dense small cell environment. We have formulated the energy efficiency optimization problem and solved the problem by BP. To apply BP, we have constructed a factor graph, and BP has been used to derive the estimated marginal probability. Then, according to our association rule, the set of turned-on APs has been determined. Whereas BP is an offline algorithm, we also have proposed its online version called DANCE, which is inspired by BP. Numerical simulations have confirmed that the proposed BP algorithm significantly enhances energy efficiency, and DANCE requiring low complexity achieves a close performance to BP. We also have applied BP and DANCE in UDN and showed that the proposed algorithms can significantly increase energy efficiency. In UDN, BP and DANCE commonly have increased energy efficiency, e.g., up to 44-times greater than that of the baseline, but also showed different properties about the power consumption and data rate. With BP and DANCE, MNOs can achieve high energy efficiency and/or data rate by selecting one of the algorithms considering the AP and UE density and operational purposes.