1

Optimization and Analysis of Probabilistic Caching
in N-tier Heterogeneous Networks

arXiv:1612.04030v2 [cs.IT] 28 Nov 2017

Kuikui Li, Chenchen Yang, Zhiyong Chen, Member, IEEE, and Meixia Tao, Senior Member, IEEE

Abstract—In this paper, we study the probabilistic caching for
an N-tier wireless heterogeneous network (HetNet) using stochastic geometry. A general and tractable expression of the successful
delivery probability (SDP) is first derived. We then optimize the
caching probabilities for maximizing the SDP in the high signalto-noise ratio (SNR) regime. The problem is proved to be convex
and solved efficiently. We next establish an interesting connection
between N-tier HetNets and single-tier networks. Unlike the
single-tier network where the optimal performance only depends
on the cache size, the optimal performance of N-tier HetNets
depends also on the BS densities. The performance upper bound
is, however, determined by an equivalent single-tier network. We
further show that with uniform caching probabilities regardless
of content popularities, to achieve a target SDP, the BS density
of a tier can be reduced by increasing the cache size of the tier
when the cache size is larger than a threshold; otherwise the BS
density and BS cache size can be increased simultaneously. It is
also found analytically that the BS density of a tier is inverse to
the BS cache size of the same tier and is linear to BS cache sizes
of other tiers.
Index Terms—Caching, Wireless HetNets, Density-Caching
Tradeoff, Power-Caching Tradeoff

I. I NTRODUCTION
The global mobile data traffic is estimated to increase to
30.6 exabytes per month by 2020, an eightfold growth over
2015, and the contribution by video is foreseen to increase
from 55% in 2015 to 75% in 2020 [2]. To address this mobile
data tsunami and hence meet the capacity requirement for the
future 5G network [3], an effective and promising candidate
solution is to deploy a dense network with heterogeneous base
stations (BSs), such as macro BSs, relays, femto BSs and pico
BSs [4]. The heterogeneous network (HetNet) can provide
higher throughput and spectral efficiency. In the meantime,
it also faces two challenges. One is the tremendous burden
on the backhual link due to the explosive demand for video
contents during the peak time and the other is high CAPEX
and OPEX due to the denser BSs.
Recently, caching popular contents at BSs has been introduced as a promising technique to offload mobile data
Manuscript received December 6, 2016; revised April 27, 2017 and
September 15, 2017; accepted November 15, 2017. This work is supported by
the National Natural Science Foundation of China under Grants 61571299,
61671291, 61329101, 61528103, and 61521062. This work was presented
in part at the IEEE SPAWC 2017 [1]. The associate editor coordinating the
review of this paper and approving it for publication was Prof. Chan-Byoung
Chae. (Corresponding author: Zhiyong Chen.)
The authors are with the Department of Electronic Engineering, Shanghai Jiao Tong University, Shanghai, 200240, P. R. China. Z. Chen is
also with the Cooperative Medianet Innovation Center, Shanghai, China.
M. Tao is also with Shanghai Institute for Advanced Communication
and Data Science (email: kuikuili@sjtu.edu.cn; zhanchifeixiang@sjtu.edu.cn;
zhiyongchen@sjtu.edu.cn; mxtao@sjtu.edu.cn.).

traffic in cellular networks [5], [6]. Unlike the communication
resources, the storage resources are abundant, economical,
and sustainable. By exploiting the abundance of the storage
resources in wireless networks, significant gains in network
capacity through caching can be expected [7], which enables
caching to be an essential functionality of emerging wireless
networks [8].
The aim of this work is to study how caching can address the
aforementioned challenges in a multi-tier HetNet. In specific,
we first would like to find out what is the optimal cache
placement strategy in order to alleviate the traffic burden in
backhaul links to the minimum. Second, we would like to
find out if the deployment cost of dense BSs can be traded by
BS cache storage, and if so, what are the tradeoffs and what
conditions must be met in order for it to happen.
A. Related Work
Caching has the potential to alleviate the heavy burden on
the capacity-limited backhaul link and also improves userperceived experience [9]. Utilizing the tool of stochastic geometry [10], the work [11] formulates the caching problem
in a scenario where small BSs are distributed according to
a homogeneous Poisson Point Process (HPPP). The authors
in [12] consider a two-tier HetNet and derive a closed-form
expression of the outage probability by jointly considering
spectrum allocation and storage constraints. In [13], the authors consider a 3-tier HPPP-based HetNet with caching and
theoretically elaborate the average ergodic rate, outage probability and delay. Considering an HPPP-based cache-enabled
small cell network, a closed form expression of the outage
probability and the optimal BS density to achieve a target
hit probability are derived in [14]. The work [15] proposes
a cluster-centric small cell network and designs cooperative
transmission scheme to balance transmit diversity and content
diversity. It is worth noting that these works mainly focus on
the performance analysis of cache-enabled wireless networks
for given caching strategies, such as caching the most popular
contents.
Caching strategy is an important issue for cache-enabled
wireless networks. Previous works on the optimal caching
strategy design can be classified into two trends based on
whether channel fading and interference are considered. The
early trend focuses on the connection topology only while
ignoring channel fading and interference. The authors in [16]
formulate a cache placement problem in distributed helper
stations to minimize the average download delay with both
uncoded and coded caching. It is shown in [16] that the

2

optimal caching problem in a wireless network with fixed
connection topology is an NP-hard problem (without coding).
In [17], a joint routing and caching design problem is studied
to maximize the content requests served by small BSs. By
reducing the NP-hard optimization problem to a variant of
the facility location problem, algorithms with approximation
guarantees are established. The second trend takes into account
channel fading and interference for caching optimization by
mostly utilizing the tool of stochastic geometry. The work [18]
proposes an optimal randomized caching policy to maximize
the total hit probability and overviews different coverage models to evaluate the performance. The works [19]–[21] optimize
the probabilistic caching strategy to maximize the successful
download probability in small cell networks. Further, a closedform expression for the optimal caching probabilities is obtained in the noise-limited scenario in [20]. In [22], a greedy
algorithm is proposed to find the optimal caching strategy to
minimize the average bit error rate. The work [23] studies the
problem of joint caching, routing, and channel assignment for
video delivery over coordinated multicell systems of the future
Internet.
Recently, caching strategy optimization is extended to wireless heterogeneous networks. The combination of the optimal
caching and the network heterogeneity brings more gains in
network capacity. Utilizing the tool of stochastic geometry, the
works [24], [25] investigate the optimal probabilistic caching
at helper stations while assuming deterministic caching at
macro stations to maximize the successful transmission probability in a two-tier HetNet. Based on [18], the work [26]
considers different types of BSs with different cache capacities. The cache optimization problem for the first type of
BSs is solved by assuming that the placement strategy for
other types of BSs is given. The joint probabilistic caching
optimization problem for all types of BSs is yet not considered,
and little analytical insight on the cache design and system
performance is available. In general, the joint optimization
for probabilistic cache placements in different tiers of a
HetNet is very challenging due to the different tier association
probabilities brought by the content diversity as well as the
complicated interference distribution by the nature of network
heterogeneity.
Furthermore, a tradeoff between the small BS density and
total storage size is firstly presented in [11], where each
small BS caches the most popular contents. Using the optimal
caching scheme, [25] shows that the helper density can be
traded by the cache size to achieve a target area spectral
efficiency. Note that the tradeoff studies in [11], [25] are conducted numerically only without theoretical analysis. Deriving
and analyzing the tradeoff theoretically has not been solved. In
[27], the authors address the question that how much caching
is needed to achieve the linear capacity scaling in the dense
wireless network based on scaling law method.
B. Contributions
In this work, we first investigate the optimal probabilistic
caching to maximize the successful delivery probability (SDP)
in a general N-tier (N ≥ 2) wireless cache-enabled HetNet.

We next establish an interesting connection between N-tier
HetNets and single-tier networks. We then address the tradeoffs between the BS caching capability and the BS density
analytically based on the uniform caching strategy. The main
contributions are summarized as follows:

•

•

•

•

Analyzing and optimizing the SDP for the N-tier HetNet:
We derive the tier association probability and the SDP
by modeling the BS locations in the HetNet as N-tier
independent HPPPs. The optimal probabilistic caching
problem for maximizing the SDP is then formulated. We
prove that this problem is concave in the high signal-tonoise ratio (SNR) regime. The sufficient and necessary
conditions for the optimal solution are derived.
Highlighting the connection between N-tier HetNets
and single-tier networks: We further study the optimal
caching problem in special cases, and find that the maximum SDP of single-tier networks only depends on the
cache size while that of N-tier HetNets is also determined
by the BS densities and transmit powers. Moreover, in the
high SNR regime, we prove that there exists a singletier network such that the maximum SDP of the Ntier HetNet is upper bounded by that of the single-tier
network. When all tiers of BSs have the same cache size,
the N-tier HetNet performs the same as the single-tier
network, regardless of the network heterogeneity.
Presenting insights on the impacts of the key network
parameters: We first show that the optimal performance
of single-tier networks is independent of the BS density
and transmit power. Then, under uniform caching strategy, we analytically present the impacts of the BS cache
size, density and transmit power of each tier on the system
performance. Numerical results also verify our analytical
results.
Revealing the tradeoffs between the BS density and the BS
cache size: With uniform caching strategy, our analysis
reveals that, to maintain a target SDP, the network parameters are related as follows: 1) increasing the BS caching
capability can reduce the BS density when its cache size
is larger than a threshold Q e ; 2) the BS density λi is
inversely proportional to the cache size Qi in the same
1
. Here, i denotes the index of the
tier, i.e., λi = QiK−Q
e
tier and K1 is independent of λi and Qi . For the different
tiers, we prove that the BS density λ j is a linear function
4 Qi
of the cache size Qi , i.e., λ j = K3 −K
, for i , j, where
K5
K3 , K4 and K5 are independent of λ j and Qi . Likewise,
we reveal the similar tradeoffs between the BS transmit
power and the BS cache size.

The rest of this paper is organized as follows. Section
II presents the system model. The performance metric is
analyzed in Section III. In Section IV, we formulate and
solve the optimal caching problem. Then, the impacts and
tradeoffs of the network parameters are shown in section V.
The numerical and simulation results are presented in Section
VI, and the conclusions are drawn in Section VII.

3

macro BS
pico BS

relay

2

relay
3
4

pico BS

1
2
3
4

pico BS

2
3

4

typical user
1

We adopt the probabilistic caching strategy and assume all
the BSs in a same tier use the same caching probabilities. Each
BS caches contents with the given probabilities independently
of other BSs. Define P , [pi j ] N ×M as the caching probability
matrix where pi j denotes the probability that the BSs in the
i-th tier caches the j-th content. It must satisfy
0 ≤ pi j ≤ 1, ∀i ∈ N, ∀ j ∈ M
M
Õ
pi j ≤ Qi, ∀i ∈ N .

(1)
(2)

j=1

Fig. 1. An example of a 3-tier wireless cache-enabled HetNet: A typical user
can obtain contents {3, 4} from relay and contents {1, 2} from macro BS.

II. S YSTEM M ODEL
A. Network and Caching Model
We consider a general wireless cache-enabled HetNet consisting of N tiers of BSs, where the BSs in different tiers
are distinguished by their transmit powers, spatial densities,
biasing factors, and cache sizes1 . The locations of BSs in each
tier are spatially distributed according to an independent HPPP,
denoted as Φi with density λi for i ∈ N , {1, 2, · · · , N }. A
three-tier HetNet including macro BSs, relays and pico BSs,
is illustrated in Fig. 1. Consider the downlink transmission.
Time is divided into discrete slots with equal duration and we
study one slot of the system. For the wireless channel, both
large-scale fading and small-scale fading are considered. The
large-scale fading is modeled by a standard distance-dependent
path loss attenuation with path loss exponent β. The Rayleigh
fading channel h is considered as the small-scale fading, i.e.,
h ∼ CN (0, 1). Each user receiver experiences an additive
noise that obeys zero-mean complex Gaussian distribution
with variance σ 2 .
Consider a database consisting of M contents denoted by
M , {1, 2, · · · , M }, and all the contents are assumed to have
equal length2 . Each user only requests one single content at
each time slot. The content popularity distribution is identical
among all users, represented by T , (ti )i ∈M , where
Í Meach user
requests the i-th content with probability ti and i=1
ti = 1.
The content popularity T is assumed to be known a prior
for cache placement. Without loss of generality, we assume
t1 ≥ t2 ≥ · · · ≥ t M . Each BS is equipped with a cache
storage. The cache capacities of N-tier BSs are denoted as
Q = {Q1, Q2, · · · , Q N }, where each BS in the i-th tier can
store at most Qi (Qi < M) contents.
When a user submits a content request, the content will be
delivered directly from the local cache of a BS that has cached
it. If the content is not cached in any BS, it will be downloaded
from the core network through backhaul links. Since the main
purpose of this paper is to optimize the caching strategy to
offload the backhaul traffic, we only consider the transmission
of the cached contents at BSs, same as [19].

Note that the conditions (1) and (2) are sufficient and
necessary for the existence of a random content placement
policy requiring no more than Qi slots of storage at each BS
in tier i for ∀i ∈ N [18]. Also note that if a BS realizes the
caching strategy by caching each file at random with the given
probability but independently of other files, the actual cache
memory in the BS can be exceeded or wasted. To strictly
meet the instantaneous cache size constraint (2) at each BS,
a novel content placement approach is proposed in Section
II-C of [18]3 . This approach brings dependency among the
caching events of different files. But such cache dependency
is irrelevant to the analysis in this work.
B. Probability of Tier Association
Without loss of generality, we carry out our analysis for a
typical user, denoted as u0 , located at the origin as in [10].
In the cache-enable HetNet, the user association policy does
not only depend on the received signal strength but also the
requested and cached contents. Specifically, when uo requests
content j, it is associated with the strongest BS among those
that have cached content j from all the tiers based on the
average received signal power. Denote the distance between
u0 and the nearest BS caching content j in the i-th tier by ri | j .
According to our tier association policy, the index of the tier
that u0 is associated with for content j is:
−β

i( j) = arg max{Bl Sl rl | j },
l ∈N

where Bl and Sl are the association bias factor and the
transmit power of BSs in the l-th tier, respectively. For notation
simplicity, we assume Bl = 1, ∀l ∈ N , in the rest of the paper.
It is essential to determine each tier’s association probability
when a user requests a content. Since each BS caches contents
independently of other BSs, the locations of the BSs caching
content j in the i-th tier can be modeled as a thinned HPPP Φi j
with density λi pi j [28]4 . Then, we have the following lemma.
Lemma 1. The probability of u0 associated with the i-th tier
for content j is given by
2

λi pi j Si β

Wi | j = Í
N

l=1

1 These model assumptions indicate that BSs in different tiers have different
traffic load to handle, and also reflect the demand heterogeneity in different
BSs.
2 Note that the extension to the general case where contents are of different
lengths is quite straightforward since the contents can be divided into chunks
with equal size.

(3)

2

λl pl j Sl β

.

(4)

3 Note that the file ordering in this content placement approach can be
varied arbitrarily at each particular realization if different file combinations
are desired.
4 Note that application of the thinning property of HPPP is based on the
assumption that each BS caches contents independently of other BSs, hence
it is not affected by the realization method of [18].

4



Proof. Please refer to Appendix A.

This lemma states that the association probability Wi | j is
determined directly by the density λi pi j and the transmit power
Si of the thinned HPPP Φi j .
III. P ERFORMANCE A NALYSIS
In this section, we analyze the SDP for a given probabilistic
caching scheme {pi j }i ∈N, j ∈M . Consider that all the BSs operate in the fully loaded state and share the common bandwidth
[10]. By using an orthogonal multiple access strategy within
a cell, the intra-cell interference is thus not considered here
and only the interference introduced by inter-tier cells and
intra-tier other cells is incorporated into analysis. Given that
u0 sends a request for content j and is associated with the ith tier, then the received instantaneous signal-to-interferenceplus-noise ratio (SINR) of u0 is given by
−β

SINRi | j =

σ2

+

ÍN

l=1

Si |hi,o | 2 Ri | j
Í

2 −β
k ∈Φl \{ni o } Sl |hl,k | dl,k

,

(5)

where Ri | j is the distance from u0 to its serving BS nio in the
i-tier tier, dl,k denotes the distance between u0 and the k-th
interfering BS in the l-th tier, |hi,o | 2 (|hl,k | 2 ) ∼ exp(1) is the
small-scale fading channel gain between u0 and the serving
BS (the k-th interfering BS). The delivery of content j from
tier i is successful when the received SINR of u0 is larger than
a threshold τ. Thus, the SDP of content j from tier i can be
expressed as5

 
(6)
Ci | j , ERi | j P SINRi | j > τ|Ri | j .
Recall that the locations of BSs caching content j in the itier can be modeled as a thinned HPPP, then the probability
density function (PDF) of Ri | j is given below.6
Lemma 2. The PDF of Ri | j is
N
2πpi j λi −π Íl=1
pl j λ l
fRi | j (r) =
e
Wi | j

Proof. Please refer to Appendix B.



Sl
Si



2
β

r2

r.

(7)


Note that different from the conventional network without
caching where each user is associated with the strongest BS
and there only exists one type of interfering BSs, in the
multi-tier cache-enabled HetNet considered in this work, the
interferences to uo when it is associated with the i-th tier for
content j can be divided into two groups. The first group of
interferences come from all the BSs (except the serving BS
nio ) in each tier that have stored content j, the locations of
which can be modeled as a thinned HPPP, denoted as Φl j with
density pl j λl for l ∈ N . The second group comes from all the
BSs in each tier that do not cache content j, the locations of
which can also be modeled as a thinned HPPP, denoted as Φl j −
5 Note that the SDP C
i | j in a cache-enabled network depends on both the
average received signal strength and the caching distribution, which is different
from the traditional coverage probability where uo is always associated with
the strongest BS.
6W
i | j and fR i | j can also be obtained by applying the thinned HPPP Φ i j
to Lemma 1 and Lemma 3 of [29], respectively.

with density (1 − pl j )λl for l ∈ N . For the first interference
group, the distance from uo to each interfering BS in Φl j is
  β1
times of the distance from uo to its serving BS
at least SSli
nio according to (30) caused by our association policy. For
the second group, the interfering BSs could be very close to
uo . By carefully handling these two groups of interferences in
N-tier HetNets, we derive an analytical expression for Ci | j in
the following proposition.
Proposition 1. The SDP Ci | j is


∫
2πpi j λi ∞
τr β σ 2
r exp −
Ci | j =
Wi | j 0
Si
!
  β2
N
Õ
 2
Sl
exp −π λl
pl j H(τ, β)+(1−pl j)D(τ, β)+pl j r dr, (8)
Si
l=1
where H(τ, β) ,
2
βτ

2
β

2τ
β−2 2 F1 (1, 1

− β2 , 2 − β2 , −τ), and D(τ, β) ,

B( β2 , 1− β2 ). Furthermore, 2 F· () denotes the Gauss hyperfunction, and B(x, y) is the Beta function defined as
∫geometric
1 x−1
µ
(1
−
µ)y−1 dµ.
0


Proof. Please refer to Appendix C.

By the law of total probability, the average SDP for u0 is
given by
M Õ
N
Õ
C,
t j Wi | j Ci | j .
(9)
j=1 i=1

Substituting (8) and (4) into (9), we obtain a tractable expression of C as follows


∫ ∞
M Õ
N
Õ
τr β σ 2
C=
2πpi j t j λi
r exp −
Si
0
j=1 i=1
!
2
 β
N
Õ
 2
Sl
pl j H(τ, β)+(1−pl j )D(τ, β)+pl j r dr. (10)
exp −π λl
Si
l=1
In the interference-limited scenario, where the noise power is
very small compared with the interference power and hence
can be neglected, the expression (10) can be simplified.
Corollary 1. In the interference-limited scenario, i.e., σ 2 →
0, the SDP C can be simplified as
0

C =

M
Õ
j=1

ÍN

i=1

2

λi Siβ pi j t j

,
2 
β T(τ, β)p
λ
S
+
D(τ,
β)
l
l
l
j
l=1

ÍN

(11)

where T(τ, β) , H(τ, β) − D(τ, β) + 1.
Proof. Substituting σ 2 = 0 to (10), we then obtain (11).



Equation (10) and (11) show a tractable expression and
a closed-form expression for the SDP in the general regime
and high SNR regime, respectively. The performance metric
depends on four main factors: the number of tiers N, the
caching probabilities {pi j }i ∈N, j ∈M , the BS densities {λi }i ∈N
and transmit powers {Si }i ∈N . In the rest of this paper, we shall
focus on the interference-limited regime with high SNR.

5

IV. C ACHING O PTIMIZATION AND A NALYSIS
In this section, we formulate and solve the optimal caching
problem for maximizing the SDP in the high SNR regime.
Further, by considering the optimal caching problem in special
cases, we establish an interesting connection between N-tier
HetNets and single-tier networks.

Í
∗
where η∗ satisfies M
j=1 p j = Q and can be found by bisection
8
method.
Based on (13), we have the following result.
Corollary 2. The optimal caching probability p∗j ( j ∈ M)
decreases with the index j, i.e, p∗j increases with the content
popularity t j . Besides, p∗j ( j ∈ M) increases with the cache
size Q.


A. Caching Optimization for General Case

Proof. Please refer to Appendix F.

The optimal caching problem of maximizing the SDP is
formulated as

By (11) for N = 1, the maximum SDP for the single-tier
0
network, denoted by C1∗ , is

P1 :

0

max

C (P)

s.t.

(1), (2)

0

C1∗ =

P

Proposition 2. Problem P1 is a concave optimization problem.


Proof. Please refer to Appendix D.

By Proposition 2, we can use the standard interior point
method to solve P1 [30]. Let P∗ = [p∗i j ] N ×M denote the
optimal solution of P1. By the Karush-Kuhn-Tucker (KKT)
conditions, the sufficient and necessary conditions for P∗ can
be stated in the following lemma.
Lemma 3. The optimal solution P∗ of Problem P1 satisfies
the following sufficient and necessary conditions:
s


+ 
N


  1 © Vi j E Õ




ª
∗
∗
pi j = min  ­
(12)
−
G k pk j −E® , 1 ,


ηi k=1,,i
  Gi


¬
 «

2

for ∀i ∈ N, ∀ j ∈ M, where [x]+ = max{0, x}, Vi j = λi Si β t j ,
2
2
ÍN
Gi = λi Si β T(τ, β), E = D(τ, β) i=1
λ S β , and ηi is the
Í Mi i ∗
Lagrangian multiplier that satisfies j=1 pi j = Qi 7 for
∀i ∈ N .


Proof. Please refer to Appendix E.

Furthermore, according to (39), we have the following
remark to state the impact of the BS cache size of each tier.
0

Remark 1. For ∀i ∈ N , the maximum SDP C ∗ increases with
the cache size Qi (Qi < M).

M
Õ

p∗j t j

j=1

T(τ, β)p∗j + D(τ, β)

,

(14)

thus we have the following remark.
Remark 2. In the interference-limited regime, the maximum
0
SDP C1∗ of single-tier networks is independent of the BS
density, transmit power, and only depends on the cache size.
This is because the serving BS and interfering BSs have the
same caching resource, and the increase in signal power
is counter-balanced by the increase in interference power.
Similar performance independency on the BS density and
transmit power also exists for traditional networks without
cache [10].
2) Connection between 1-tier and N-tier HetNets: We
observe from (11) that the SDP in the high SNR regime
depends on the caching probabilities {pi j }i ∈N, j ∈M through
2
ÍN
β
2
ÍN
i=1 λ i Si pi j
β p . Thus, by defining x =
λ
S
, then we
ij
j
2
i=1 i i
Í
N
i=1

λ i Si

β

can formulate a new problem below,
M
Õ
xj tj
P2 : max
x
T(τ, β)x j + D(τ, β)
j=1

(15)

0 ≤ x j ≤ 1, ∀ j ∈ M
(16)
2
Í
M
N
Õ
β
i=1 λi Si Q i
.
(17)
xj = Í
2
N
β
j=1
i=1 λi Si
From (14) and (15), it is seen that P2 is identical with the
caching optimization problem in a single-tier network with x
s.t.

being the caching probability vector and Q e =

ÍN

i=1

ÍN

λ i Si

i=1

B. Caching Optimization for Special Cases
1) Optimization for N=1: When N = 1, the network
degrades to the single-tier network. Denote P1 , (p j )1×M as
the caching strategy and Q < M as the cache size for the
single-tier network, then the optimal caching probabilities in
(12) become:
s


+ 

 1


 
t
D(τ,
β)
D(τ,
β)
j
∗
 ,1 , ∀ j ∈ M (13)
p j = min 
−

∗

η
T(τ, β)  
  T(τ, β)

 

7 From

(39), it can be shown by contradiction that the maximum SDP is
achieved only when constraint (2) holds with equality. The similar proof is
given by Lemma 2 of [18]. In the rest of this paper, we use (2) with equality
as the constraint.

2
β

λ i Si

Qi
2
β

being the equivalent cache size. By further comparing P1 and
P2, we obtain the following proposition which states a general
relationship between the optimal performance of an N-tier
HetNet and that of a single-tier network.
0

Proposition 3. Let C ∗ be the optimal objective of P1 for
the considered N-tier HetNet. For the single-tier network with
cache size Q e =

ÍN

i=1

λ i Si

2
β

i=1 λ i Si

ÍN

Qi
2
β
0

, we have
0

C ∗ ≤ C1∗
8 Note that the optimal caching probability in this special case is consistent
with the prior works on single-tier networks in [19], [21]. Our work extends
the probabilistic caching strategy optimization for a single-tier network to
that for a general N -tier HetNet and contributes to presenting the impacts
and essential tradeoffs of the heterogeneous network parameters.

6

with equality if Qi = Q j for ∀i, j ∈ N .
Proof. Please refer to Appendix G.



Proposition 3 states that in the interference-limited regime,
the optimal performance of an N-tier HetNet with BS cache
sizes {Qi }i ∈N , BS densities {λi }i ∈N , and BS transmit powers
{Si }i ∈N is upper bounded by that of a single-tier network with
BS cache size Q e =

ÍN

i=1

λ i Si

2
β

i=1 λi Si

ÍN

Qi
2
β

, arbitrary BS density, and

arbitrary BS transmit power. Further, their performances are
the same when the cache sizes {Qi }i ∈N are the same for all
tiers in the N-tier HetNet.
3) Optimization for Qi = Q j , ∀i, j ∈ N : In this case,
Proposition 3 shows that P1 is equivalent to P2. Based on (13),
the optimal solution P∗ of P1 can be further given below.
Corollary 3. When Qi = Q j for ∀i, j ∈ N , there exists an
optimal solution P∗ of P1 satisfying p∗i j = x ∗j for ∀i ∈ N, where
the optimal solution x∗ of P2 follows (13) with p∗j = x ∗j and
Q = Qe .

tradeoffs of the BS density λi , transmit power Si and cache
size Qi are found.

A. The Equivalence under Uniform Cache
Consider the uniform caching strategy where P1,u = {p j =
i
and Pu = {pi j = Q
M }i ∈N, j ∈M for single-tier and Ntier HetNets, respectively. By substituting Pu and P1,u into
(11) and (15), respectively, the equivalence of the system
performance of N-tier HetNets and 1-tier Networks can be
established.
Qe
M } j ∈M

Proposition 4. In the interference-limited regime, the SDP
0
C (Pu ) of a N-tier HetNet with the caching strategy Pu equals
0
that C1 (P1,u ) of the single-tier network with the caching
strategy P1,u , and is given by
0

0

C (Pu ) = C1 (P1,u ) =

M
Õ
j=1

Qe t j
,
T(τ, β)Q e + D(τ, β)M

(18)

Proof. Problem P2 is the new constructed single-tier caching
optimization problem. Thus, its optimal solution x ∗j is the same
as (13) where let p∗j = x ∗j and Q = Q e . In the second part of
the proof of Proposition 3, we show that when Qi = Q j for
∀i, j ∈ N , P1 is equivalent to P2, and P∗ with p∗i j = x ∗j is also
an optimal solution of P1. Thus, Corollary 3 is proved.


where the cache size of the single-tier network Q e

Based on Remark 2 and Proposition 3, we have the following remark.

Based on (13) and C (Pu ) = C1 (P1,u ), we have the following
corollary.

Remark 3. In the interference-limited regime, when all the
BSs in the N-tier HetNet have the same cache size, the
maximum SDP of the HetNet is independent of the network
heterogeneity in the BS density and transmit power.
Traditionally, without caching ability at BSs, the outage
probability is independent of the number of tiers, the BS
densities and transmit powers in the interference-limited K-tier
HetNets [29]. By introducing caching resource into the system,
0
(11) states that the maximum SDP C ∗ generally depends on
the number of tiers N, the BS cache size {Qi }i ∈N , the BS
density {λi }i ∈N and transmit power {Si }i ∈N . The intuition
behind this observation is that the caching resource changes
the decision of a user to access a BS. In the multi-tier cacheenabled HetNet, the decision not only depends on the received
SINR, but also depends on the contents cached at BSs. In order
0
to further understand the impact of the cache size Qi on C (P),
we theoretically illustrate the relationship between the cache
size and the BS density/transmit power in the next section.
V. A NALYSIS ON N ETWORK PARAMETERS UNDER
U NIFORM C ACHE
In this section, with the uniform caching strategy where
each content is cached with equal probabilities regardless of
content popularities, the equivalence between the SDP of an
N-tier HetNet and that of a single-tier network is obtained.
Based on this property, we further investigate the impacts of
the key network parameters, i.e, the BS cache size, density
and transmit power on the system performance. Finally, the

2

=

β
i=1 λi Si Q i
2 .
ÍN
β
i=1 λ i Si

ÍN

Proof. Substituting Pu and P1,u into (11) and (15), respectively, we then have (18).

0

0

Corollary 4. In the interference-limited regime, consider the
scenario that the video content popularity distribution is an
1
uniform distribution, i.e., t j = M
, ∀ j ∈ M, then the uniform
caching strategies P1,u and Pu are the optimal probabilistic
caching strategies for single-tier networks and N-tier HetNets,
respectively.
1
Proof. When t j = M
, ∀ j ∈ M, we have p∗1 = p∗2 = · · · =
Í
Qe
∗
∗
p∗M from (13). Since M
j=1 p j = Q e , we thus have p j = M
0∗
0
for ∀ j ∈ M, and the maximum C1 = C1 (Pu ). Then from
0
0
Proposition 3 and Proposition 4, we have C ∗ = C (Pu ), and
hence Pu is also the optimal solution of P1.


B. Impacts of λi , Si , Qi , ∀i ∈ N
Proposition 4 further states that, with uniform caching
regardless of content popularity, the SDP performance of Ntier HetNets depends on all the system parameters through
the equivalent cache size Q e =

ÍN

i=1

λ i Si

2
β

i=1 λ i Si

ÍN

Qi
2
β

only. Then, the

impacts of the BS density λi and cache size Qi of tier i can
be obtained in the following two lemmas.
0

Lemma 4. The SDP C (Pu ) increases with Si and λi when
ÍN

Qi ≥
λi .

j=1,,i

λj Sj

2
β

j=1,,i λ j S j

ÍN

Qj
2
β

0

. Otherwise, C (Pu ) decreases with Si or

7

× (π5002)−1

25

80

S1=S2+10dBm

Tier−1 Density λ1

Tier−1 Transmit Power S1 [dBm]

Qe=20, C (Pu)=0.18

20

15

Q1>Qe>Q2=5

0≤Q1<Qe<Q2=30
10

5

λ1=λ2/5

Qe=20, C’(Pu)=0.18

75

S1=S2+20dBm

’

λ1=λ2/10
Q1>Qe>Q2=5

70

0≤Q1<Qe<Q2=30

65
60
55
50
45

0

0

5

10

15

20

25

30

35

40

40

0

5

Cache Size Q1

10

15

20

25

30

35

40

Cache Size Q

1

(a) The tradeoff between λ1 and Q1 . We consider two cases: (1). {S1 , (b) The
S1 and Q1 . We
o
n consider twoo cases: (1). {λ1 ,
n tradeoff between
1
2
S2 }={43 dBm, 33 dBm} (2). {S1 , S2 }={53 dBm, 33 dBm}. We fix λ2 }=
, 10 2 (2). {λ1 , λ2 }=
, 10 2 . We fix S2 =33
2
5
π500
π500
π5002
π500
for both cases.
λ2 =
dBm for both cases.
π5002
0

Fig. 2. Under fixed C (Pu ), the tradeoffs of the BS density λ1 , transmit power S1 and cache size Q1 in a two-tier HetNet where τ = −10 dB, t j =
0

γ = 0.8, β = 4, C (Pu ) = 0.18 (Q e =20). By Corollary 5, we have K1 < 0, K2 < 0 for Q1 ∈ [0, Q e ), and K1 > 0, K2 > 0 for Q1 ∈ (Q e , +∞).
0

Proof. From (18), we can see that C1 (P1,u ) increases with Q e .
Then we have
2
2
2
Í
ÍN
β
β
Siβ (Qi N
∂Q e
j=1 λ j S j −
j=1 λ j S j Q j )
=
2
Í
∂λi
( N λ S β )2
j=1

2
β

j=1,,i

ÍN

j=1,,i

λj Sj

ÍN

j=1,,i

2
β

λj Sj

Qj

∂Q
∂λi

, we have that

2
β

≥0

∂Q
and ∂S
≥ 0, which also means that C1 (P1,u ) increases with λi
i
0
0
and Si 9 . Due to C (Pu )=C1 (P1,u ), we thus have Lemma 4. 
0

0

Lemma 5. The SDP C (Pu ) increases with Qi (Qi < M) and
2
the increasing speed is monotonic to λi Si β , ∀i ∈ N .
∂Q e
∂Qi

2

=
0

λ i Si β
ÍN
( i=1 λi Si

2
β )2

> 0, Q e increases with Qi
0

for ∀i ∈ N. Since C1 (P1,u ) increases with Q e and C1 (P1,u ) =
0
0
C (Pu ), C (Pu ) increases with Qi . We thus have this lemma.

Remark 4. From Lemma 4 and Lemma 5, we can observe
two important features in the N-tier cache-enabled HetNet:
• Increasing the density or transmit power of the BSs
from the tier with small cache size decreases the system
performance. This somewhat surprising result is actually
intuitive since such BSs (e.g., pico or femto) with small
cache only provide little service but bring strong interferences to other BSs (e.g., macro or relay).
9 This

If the BS transmit power and density of the i-th tier are
both the largest among all the tiers, it is most effective
to increase the performance of the network by increasing
the BS cache size of the i-th tier.

ÍN

Obviously, when Qi ≥

Proof. Due to

,

j j

2
2
Í
β
λ j Sj β − N
j=1,,i λ j S j Q j )
=
,
2
Í
β 2
( N
j=1 λ j S j )
2
2
2
Í
ÍN
β −1
2
β
β
(Qi N
λ
S
i
∂Q e
j=1,,i λ j S j − j=1,,i λ j S j Q j )
β
i
.
=
2
Í
∂Si
β )2
( N
λ
S
j
j
j=1

Si (Qi

•

1/ j γ
ÍM
1/k γ
k=1

is not contradict with Remark 2 in Section IV, because the BS
densities {λi }i∈N , transmit powers {Si }i∈N and cache sizes {Qi }i∈N affect
the cache size Q in here based on Q e .

C. Tradeoffs of Qi , λi , and Si , ∀i ∈ N
In the preceding subsection, we describe the impacts of Qi ,
0
λi , and Si on C (Pu ). Now we will present the tradeoffs of
these network parameters at a target SDP.
1) Tradeoffs of one tier parameters: Lemma 4 and 5 show
that increasing the BS density, transmit power or cache size
0
influence the SDP C (Pu ) by changing Q e . This suggests
that, as long as Q e does not change, one can interchange
different types of system parameters to maintain the same
system performance, and hence can obtain the tradeoffs of
these network parameters at a target SDP. In the following,
we elaborate the tradeoff between the BS density λi (transmit
power Si ) and the cache size Qi within each tier i for ∀i ∈ N .
0
Given a target SDP C (Pu ), the communication resource (λ j
and S j ) and the caching resource (Q j ) of all the tiers without
the i-th tier, i.e., ∀ j ∈ N and j , i, we can obtain Q e based
on (18), thereby gaining the tradeoffs of the i-th tier’s cache
size Qi , BS density λi and transmit power Si , as illustrated in
the following theorem.
Theorem 1. With the uniform caching strategy, given a target
0
0
C (Pu ) determined by C1 (P1,u ), and the fixed values λ j , S j ,
Q j for ∀ j ∈ N and j , i, the network parameters Qi , λi , Si
satisfy the following tradeoffs
K1
, for given Si, ∀i ∈ N,
Qi − Qe
β
K2
Si = (
) 2 , for given λi, ∀i ∈ N,
Qi − Qe

λi =

(19)
(20)

8

where

0.8

(21)

N
Õ
2
λj
( )S j β (Q e − Q j ).
λ
j=1,,i i

(22)

λj (

Proof. Please refer to Appendix H.



Interestingly, we can observe from Theorem 1 that λi is
inversely proportional to Qi , while Si is a power function
of Qi with a negative exponent (−β/2). Accordingly, it
is natural to ask that if the BS density can be reduced by
increasing the BS caching capability. If yes, what is the
condition? Thus, we have the following corollary to answer
this question.

0.5
0.4
0.3

0.1
0

0.2

0.4

0.6

Zipf Exponent γ

0.8

1

1.2

Fig. 3. SDP with different caching strategies. Here, we set M = 1000 and
{Q1, Q2 }={200, 50}.
1
0.9

Case 1: Tier 1 cache
Case 1: Tier 2 cache
Case 2: Tier 1 cache
Case 2: Tier 2 cache
Case 3: Tier 1 cache
Case 3: Tier 2 cache
Case 4: Tier 1 cache
Case 4: Tier 2 cache
Zipf distribution

∆Q1=10

0.8

Optimal Caching Probability

•

0.6

0.2

Corollary 5. To maintain the same SDP, we have the following
results for the i-th tier
•

Theory: Optimal caching strategy
Theory: Popular caching strategy
Theory: Uniform caching strategy
Simulation: Optimal caching strategy

0.7

Sj
) (Q e − Q j ),
Si

j=1,,i

K2 =

2
β

Successful Delivery Probability

K1 =

N
Õ

The BS density λi and transmit power Si decrease with
the cache size Qi , i.e., K1 > 0 and K2 > 0,, when Qi ∈
(Q e, ∞).
The BS density λi and transmit power Si increase with the
cache size Qi , i.e., K1 < 0 and K2 < 0, when Qi ∈ [0, Q e ).

0.7
0.6
0.5
0.4
0.3

∆Q2=5

0.2
0.1

Proof. Please refer to Appendix I.



It is worth mentioning that increasing the BS caching
capability of one tier does not always reduce the BS density or
transmit power of this tier to achieve the same performance in
a multi-tier NetHet, in contrast to single-tier caching networks
[11]. Taking N = 2 for example, the tradeoff between λ1
and Q1 is shown in Fig. 2(a). The target SDP determined
0
by C1 (P1,u ) is 0.18. Note that if a tier is deployed with larger
cache size (Q1 > Q e = 20), the more the BS cache size (Q1 )
increases, the less the BS density (λ1 ) becomes, in order to
0
maintain the same C (Pu ). However, if the tier is deployed
with small cache size (Q1 < Q e = 20), increasing both the
BS cache size and density can keep the same performance,
as show in Fig. 2(a). This is because the increase in BS
density of the tier with small cache size only improves its
own tier association probability, and also causes the stronger
interference to other tiers with larger cache sizes. Therefore,
more users are associated with the tier with small cache size,
and the SDP decreases because this tier can only serve fewer
content requests. The tradeoff between S1 and Q1 shown in
Fig. 2(b) is similar to that of λ1 and Q1 .
2) Tradeoffs of different tier parameters: Similarly, we
can find the tradeoffs of the parameters in different tiers, as
specified in Theorem 2.
Theorem 2. When Qi , λ j , and S j (i, j ∈ N, i , j) satisfy the
0
following tradeoffs for any two different tiers, C (Pu ) remains
constant.
K3 − K4 Qi
,
K5
K3 − K4 Qi β
Sj = (
)2,
K6

λj =

(23)
(24)

0

20

40

60

80
100
120
Content Index

Fig. 4. The optimal caching probability
transmit powers.

pi∗ j

140

160

180

200

for different storage sizes and

for ∀i, j ∈ N, i , j, where
N
Õ

K3 = Q e

λk Sk

2
β

−

N
Õ

2

λk Sk β Q k ,

(25)

k=1,,i, j

k=1,,j
2

K4 = λi Siβ ,

(26)

2
β

K5 = S j (Q j − Q e ),

(27)

K6 = λ j (Q j − Q e ).

(28)

0

Proof. When C1 (P1,u ) keeps constant, Q e will be a fixed value.
We consider the tradeoffs between the j-th tier parameters λ j ,
S j and the i-th tier parameter Qi assuming that other parameters are constants. Then according to Q e =

ÍN

i=1

ÍN

i=1

have
Qe

λ i Si

N
Õ

2

λk Sk β −

k=1,,j

N
Õ

2

2

2
β

λ i Si

Qi
2
β

, we

2

λk Sk β Q k = λi Siβ Qi +λ j S j β (Q j −Q e ), (29)

k=1,,i, j

Substituting (25), (26), (27), (28) into (29), we thus have (23),
(24).

Different than Theorem 1, the BS density λ j is a linear
function of the BS cache size Qi , and the BS transmit
power S j is a power function of Qi with a positive exponent
β
2 , where i , j. Then, the function changing law follows the
following corollary.
Corollary 6. For different tiers, we have

9

0.5

S1=S2+10 dBm,λ1=λ2/10
Successful Delivery Probability

Successful Delivery Probability

S1=S2+10 dBm,λ1=λ2/5
S1=S2+20 dBm,λ1=λ2/10
0.45

0.4

0.6
0.55
0.5
0.45

0.4
30
× (π5002)−1

25
20

0.35
20

25

30

35

40

15

Tier−2 Density λ2

Cache Size Q1

(a) Cache size Q1 . Here, we use {λ1, λ2 } = { 1 2 , 5 2 } or { 1 2 ,
π500
π500
π500
10
2 }, {S1 , S2 } = {43 dBm, 33 dBm} or {53 dBm, 33 dBm}.

6

4

2

10
8 × (π5002)−1

Tier−1 Density λ1

(b) Densities {λ1 , λ2 }.

π500

0.54

Q1=150, Q2=25

Q1=30, Q2=5

Q1=200, Q2=25

Q1=40, Q2=5

Q =150, Q =75

Q1=30, Q2=15

1

Successful Delivery Probability

Successful Delivery Probability

0.56

2

0.52
0.5
0.48
0.46
0.44

0.55

0.5

0.42
0.4
40

45

50

0.45
20

55

25

Transmit Power S1 [dBm]

30

35

Transmit Power S [dBm]
2

(c) Tier-1 transmit power S1 . Here, we set M = 1000.
(d) Tier-2 transmit power S2 .
0
Fig. 5. The impacts of the BS cache size Q1 , densities {λ1 , λ2 }, transmit powers {S1 , S2 } on the maximum SDP C ∗ .
•

•

When Q j > Q e , we have K5 > 0, K6 > 0 and K3 > 0,
K3
i.e., λ j and S j decrease with Qi where Qi ∈ [0, K
];
4
When Q j < Q e and Q e >

ÍN

k=1,,i, j

ÍN

λ k Sk

k=1,, j

2
β

λ k Sk

Qk

2
β

, we have

K5 < 0, K6 < 0 and K3 > 0, i.e., λ j and S j increase with
K3
Qi where Qi ∈ [ K
, +∞);
4
•

When Q e <

ÍN

k=1,,i, j

ÍN

λ k Sk

k=1,, j

2
β

λ k Sk

Qk

2
β

, we have K5 < 0, K6 < 0

and K3 < 0, i.e., λ j and S j increase with Qi where Qi ∈
[0, +∞).
Proof. Please refer to Appendix J.



Remark 5. To keep the target SDP, the tradeoffs of the
parameters in two different tiers are summarized as two cases:
• For tier i with cache size Q i > Q e , the BS density and
transmit power of this tier can be reduced when another
tier’s cache size increases;
• For tier i with cache size Q i < Q e , the BS density and
transmit power of this tier need to be increased when
another tier’s cache size increases.
When increasing the density of the tier with small cache

size, the SDP decreases since the tier can only serve fewer
content requests and bring the strong interference to other tiers
with larger cache sizes. Thus, the density of such a tier cannot
be reduced by increasing another tier’s cache size. Inversely,
increasing both the density of the tier with small cache size
and the cache size of another tier can keep the fixed SDP.
VI. S IMULATION AND N UMERICAL R ESULTS
In this section, simulation and numerical results are provided to evaluate the performance of the proposed caching
strategy, and present the impacts and tradeoffs of the network
parameters. We consider a two-tier cache-enabled HetNet
1
5
where the BS densities {λ1 , λ2 }={ π500
2 , π5002 }, the transmit
powers {S1 , S2 }={53, 33} dBm and the normalized cache
sizes {Q1, Q2 }={40, 10}. We consider the interference-limited
scenario and set the number of contents M = 200, the path
loss exponent β = 4, and the SINR threshold τ = −10
dB. In the simulation, the BSs are deployed according to
independent HPPPs in an area of 5000 m×5000 m, the content
popularity distribution is modeled as the Zipf distribution,
i.e., the popularity
of the j-th ranked content is given by
γ
t j = Í M1/j1/k γ , where γ ≥ 0 characterizes the skewness of
k=1

10

the popularity distribution. The parameter γ is set to 0.8 in
the simulation. The optimal caching strategy is obtained by
using interior point method. These parameters will not change
unless specified otherwise.

0.45

Successful Delivery Probability

A. Optimal Caching Strategy

0.22

0.25

0.215
0.2

0.21

0.15

0.205
1.9

2

2.1

0.1

−5

0
5
10
SIR Threshold τ [dB]

15

20

Fig. 6. Comparison between 1-tier networks and 2-tier HetNets.
Theory: Special case 1
Simulation: Special case 1
Theory: Special case 2
Simulation: Special case 2
Theory: Special case 3
Simulation: Special case 3
Theory: Special case 4
Simulation: Special case 4

0.6
0.406
0.55

0.405
0.404

0.5

0.403

6.85

6.9

3

4

6.95

0.45
0.4
0.35
0.3

1

2

5

6

Density λ (λ)
1

7

8

9

10

× (π5002)−1

Fig. 7. The impacts of tier-1 (single-tier network) BS density λ1 (λ) on the
0
maximum SDP C ∗ for four special cases.

cache size decreases the optimal performance. This is expected
because BSs with smaller cache size (e.g., pico or femto) can
only provide less service to satisfy user requests while bringing
strong interference to other BSs (e.g. macro or relay stations).
In addition, the impacts of BS transmit powers are similar to
that of BS densities, which are illustrated in Fig. 5(c) and Fig.
5(d).
C. The special cases
Connection between 1-tier networks and N-tier HetNets.
In Fig. 6, we consider four scenarios:
1) 1-tier network with Q e = 20;
2) 2-tier HetNets with {Q1, Q2 } = {20, 20};
3) 2-tier HetNets with {Q1, Q2 } = {25, 10};
4) 2-tier HetNets with {Q1, Q2 } = {28, 4}.
In scenarios 2 − 4, the two-tier BS cache sizes satisfy
Í2

λ i Si

2
β

i=1 λ i Si

Í2

Impacts of cache size. Fig. 5(a) illustrates that increasing
the cache size can greatly improve the maximum SDP, which
verifies the validity of Remark 1. It also shows that increasing
λ1 or S1 can further improve the optimal performance. That
means more user requests can be served since BSs in the
first tier have higher association probability and more cache
resource.
Impacts of BS densities and transmit powers. It can be
observed from Fig. 5(b) that the optimal SDP decreases with
the BS density λ2 (with Q2 = 10) and increases with λ1 (with
Q1 = 40), which means deploying more BSs with smaller

2−tier HetNets, {Q1, Q2}={28, 4}

0.3

0
−10

i=1

B. Impacts of network parameters with optimal cache scheme

2−tier HetNets, {Q1, Q2}={25, 10}

0.35

0.05

Successful Delivery Probability

Comparsion with other caching strategies. To investigate the performance with larger content library, we set
M = 1000 in Fig. 3. It can be seen that the optimal probabilistic caching strategy always outperforms the two considered
baselines. The first baseline is popular caching strategy where
each BS only caches the most popular contents, and the
second is uniform caching strategy where each BS caches
each content randomly with equal probabilities. The popular
caching strategy can only perform as well as the optimal
probabilistic caching with highly skewed content popularity,
e.g., γ > 1. The reason for this is that the majority of user
requests only focus on those very few popular contents and
the optimal probabilistic caching tends to cache those most
popular contents. The performance of uniform caching strategy
is not affected by γ since it caches each content with equal
probabilities. The theoretical results for the proposed optimal
caching scheme are also validated by simulation in Fig. 3.
Optimal caching probabilities for different network
parameters. In Fig. 4, we consider four cases:
10
1
1) Case 1: {Q1, Q2 }={30, 5}, {λ1 , λ2 }={ π500
2 , π5002 }, {S1 ,
S2 }={43, 33} dBm;
1
10
2) Case 2: {Q1, Q2 }={40, 10}, {λ1 , λ2 }={ π500
2 , π5002 },
{S1 , S2 }={43, 33} dBm;
1
5
3) Case 3: {Q1, Q2 }={30, 5}, {λ1 , λ2 }={ π500
2 , π5002 }, {S1 ,
S2 }={43, 33} dBm;
5
1
4) Case 4: {Q1, Q2 }={30, 5}, {λ1 , λ2 }={ π500
2 , π5002 }, {S1 ,
S2 }={53, 40} dBm.
The content popularity based on Zipf’s law is also illustrated
in Fig. 4. It can be seen that the optimal caching probabilities
highly depend on the network parameters, e.g., cache size,
content popularity, transmit power and BS density. Caching
the most popular contents is not always the best solution.
Specifically, the optimal caching probability is positively related to the content popularity. The results also show that
increasing the BS cache size has a more significant impact
on the optimal caching probabilities than increasing the BS
density and transmit power.

1−tier network, Q=20
2−tier HetNets, {Q1, Q2}=20, 20

0.4

Qi
2
β

= Q e = 20. It can be seen that the optimal

performances of the two-tier HetNets considered in scenarios
2 − 4 are always upper bounded by that of the single-tier
network considered in scenario 1. In addition, the two-tier
HetNet of scenario 2 achieves the same performance as the
single-tier network of scenario 1 when their cache sizes satisfy
Q1 = Q2 = Q e . These results verify our conclusions in
Proposition 3 and Remark 3.
Impacts of BS density for special cases. We consider
four special scenarios in Fig. 7:
1) 1-tier network with Q = 20, S = 33 dBm;
2) 2-tier HetNets with Q1 = Q2 = 15, {S1 , S2 }={53, 33}
5
dBm, λ2 = π500
2;

20

50

S1=S2+10dBm

18

1

16

Qe=20, C (Pu)=0.18

Tier−2 Density λ2

14
12
10
8 0≤Q ≤Qe<Q =30
2

6

Q ≥Qe>Q =5
1

2

4

Qe=20, C (Pu)=0.18

40
35
30

Q1≥Qe>Q2=5

0≤Q ≤Q <Q =30

25

1

e

2

20

Q1= Qe=20,
S2=−∞ dBm (0 W)

15
10

2
0

λ1=λ2/10

’

2

’

1

λ1=λ2/5

45

S =S +20dBm
Tier−1 Transmit Power S2 [dBm]

× (π5002)−1

11

0

5

10

15

20

25

30

35

40

5

0

5

10

15

Cache Size Q

20

25

30

35

40

Cache Size Q

1

1

(a) The tradeoff between λ2 and Q1 . Two cases: (1). {S1 , S2 }={43, (b) The tradeoff between S2 and Q1 . Two cases: (1). {λ1 ,
33} dBm (2). {S1 , S2 }={53, 33} dBm. We fix λ1 = 1 2 for both λ2 }={ 2 2 , 10 2 } (2). {λ1 , λ2 }={ 1 2 , 10 2 }. We fix S1 =53
π500
π500
π500
π500
π500
cases.
dBm for both cases.
0

0

Fig. 8. Under fixed C (Pu ), the tradeoffs of the BS density λ2 , transmit power S2 and cache size Q1 in a two-tier HetNet where C (Pu ) = 0.18 (Q e =20).
By Corollary 6, we have K3 > 0, K5 > 0, Q1 ∈ [0, Q e ] for Q2 = 30 > Q e , and K3 > 0, K5 < 0, Q1 ∈ (Q e , +∞) for Q2 = 5 < Q e .

3) 2-tier HetNets with {Q1, Q2 } = {15, 5}, {S1 , S2 }={73,
5
33} dBm (S1  S2 ), λ2 = π500
2;
4) 2-tier HetNets with {Q1, Q2 } = {15, 5}, S1 = S2 = 33
dBm, λ2 = λ1 .
We can see that the optimal performances are independent of
the BS density in the considered four special cases, which
also verifies our conclusions in Remark 2 and Remark 3. The
optimal performance of scenario 2 is almost the same as that of
scenario 3 because both them are almost entirely determined
by the tier with cache size Q = 15. Besides, from scenario 4, it
can be seen that when two tiers of BSs are deployed with the
same density and transmit power, however, still have different
cache sizes, performance independency on the BS density still
exists.
D. Tradeoffs of network parameters
In Fig. 8, we present the tradeoffs of the network parameters
0
in different tiers. To keep the fixed C (Pu ), increasing the BS
cache size Q1 of the first tier cannot always bring the reduction
of the BS density λ2 or transmit power S2 of the second tier.
The crucial condition is that the cache size of the second tier
should satisfy Q2 > Q e . It can be further seen in Fig. 8(a)
that the slope that λ2 varies with Q1 is positively related to
S2 . Similarly, the slope that S2 varies with Q1 is also positively
related to λ1 , as illustrated in Fig. 8(b).
VII. C ONCLUSION
In this paper, we contributes to developing a general Ntier wireless cache-enabled HetNet framework and proposing
an optimal probabilistic caching scheme for maximizing the
SDP. Our analysis results show that the optimal performance
of such a network not only depends on the cache size, but
also the BS densities and transmit powers. Numerical results
show that the proposed optimal caching strategy achieves a
significant performance gain.
Toward analyzing the optimal performance of N-tier HetNets, we establish an interesting connection between N-tier
HetNets and single-tier networks. Furthermore, we quantify

the tradeoff between the BS cache size and the BS density with
the uniform caching strategy. We show that: i) increasing the
BS caching capability does not always reduce the BS density;
ii) the BS density decreases with the BS caching capability
when the BS cache size is larger than a threshold; iii) the BS
density is inversely proportional to the BS cache size within
the same tier, and becomes a linear function of the BS cache
sizes of different tiers.

A PPENDIX
A. Proof of Lemma 1
By definition, we have




Wi | j = P Sr,i | j > max Sr,l | j
l,l,i
" N
#
i
Ö h
−β
−β
P Si r > Sl rl | j |ri | j = r
= Eri | j
l=1,l,i

=

∫

∞

0

"

N
Ö

P rl | j

l=1,l,i



Sl
>r
Si

 β1 #
fri | j (r)dr,

(30)

where fri | j (r) is the probability density function (PDF) of ri | j ,
  β1
means that it does not exists any BS of the
and rl | j > r SSli
   1 2
β
l-th tier caching the content j in the area A = π r SSli
.
Using the fact that the null probability of a 2-D Poisson point
process with density λ in an area A is e−λA, we can deduce
"
P rl | j



Sl
>r
Si

 β1 #

=e

−π pl j λl



Sl
Si



2
β

r2

,

(31)

Thus, fri | j (r) can be given by


d 1 − P[ri | j > r]
2
fri | j (r) =
= 2πpi j λi re−πλi pi j r .
dr

(32)

12

Substituting (31) and (32) into (30), we have
Wi | j =

∞

∫

2πpi j λi re

−π

ÍN

−π

ÍN

l=1, l,i

pl j λ l







Sl
Si



2
β

r 2 −πλi pi j r 2

e

dr

0
(a)

= 2πpi j λi

∫

∞

e

l=1

pl j λ l

Sl
Si

2
β

r2

rdr

0
2

λi pi j Si β

=Í
N

l=1

2

λl pl j Sl β

,

(33)

  β2
where (a) follows that pl j λl SSli
= λi pi j for l = i. We
complete the proof of Lemma 1.
B. Proof of Lemma 2
According to the definition, we have


d 1 − P[Ri | j > r]
.
fRi | j (r) =
dr

(34)

Let k j denote the tier associated with the typical user requesting content j. Therefore, we have
P[Ri | j > r] = P[ri | j > r |k j = i]


P ri | j > r, Sr,i | j > maxl,l,i Sr,l | j
=
Wi | j
"
  β1 #
∫ ∞ Ö
N
Sl
1
P rl | j >r
=
fri | j (r)dr
Wi | j r l=1,l,i
Si
 2
∫ ∞
ÍN
S
β 2
(b) 2πpi j λi
−π l=1
r
pl j λ l S l
i
e
rdr, (35)
=
Wi | j
r
where (b) is obtained by substituting (31) and (32). Substituting (35) into (34), we thus have Lemma 2.
C. Proof of Proposition 1
Í
−β
According to (6) and let Il = k ∈Φl \{ni o } Sl |hl,k | 2 dl,k , we
then have
∫ ∞


Ci | j =
P SINRi | j (r) > τ|r fRi | j (r)dr
0
#
∫ ∞ "
N
Õ
=
P |hi,o | 2 > τSi −1 r β (σ 2 +
Il )|r fRi | j (r)dr
0

(c)

=

=

∫

l=1
∞

0
∫ ∞

h

EIl e

ÍN
−τSi −1 r β (σ 2 + l=1
Il )



exp −τSi −1 r β σ 2

0

N
Ö

i
|r, Il fRi | j (r)dr



 Ö


1


= EΦl j 
−β 
k ∈Φl j \{ni o } 1 + τSi −1 r β Sl dl,k 


!


∫ ∞
1
(d)
= exp −2πpl j λl  S  1 1 −
xdx
β
1 + τSi −1 r β Sl x −β
r Sl
i
! !
2
  β2
∫ ∞
β −1
2
µ
2
S
(e)
l
2
dµ
= exp −πpl j λl τ β
r
Si
β τ −1 1 + µ


Sl β2 2
= exp −πpl j λl ( ) r H(τ, β) ,
(37)
Si
where (d) follows from the probability generating functional (PGFL) of the PPP and (e) is obtained using µ =
2
2
2τ
τ −1 Si r −β Sl−1 x β , H(τ, β) denotes β−2
2 F1 (1, 1 − β , 2 − β , −τ).
2

Similarly, let D(τ, β) be β2 τ β B( β2 , 1− β2 ), the Laplace transform
L Il j − is calculated as


L Il j − τSi −1 r β



∫ 
1
= exp −(1 − pl j )λl
xdx
1−
1 + τSi −1 r β Sl x −β
R2


Sl 2
(38)
= exp −π(1 − pl j )λl ( ) β r 2 D(τ, β) .
Si
Substituting (7), (37), (38) into (36), we have Proposition 1.

D. Proof of Proposition 2
It is easy to see that the feasible set {P|0 ≤ pi j ≤ 1,
ÍM
∀i ∈ N, ∀ j ∈ M} is a convex set.
j=1 pi j ≤ Q i ,
0
The objective function C (P) is a summation of M polynomials.Í For ∀ j ∈ M, the j-th polynomial is given by
N
2
2
V pi j
Fj = Í Ni=1G ipj +E
, where Vi j = λi Si β t j , Gi = λi Si β T(τ, β),
i=1 i i j
2
ÍN
and E = D(τ, β) i=1
λi Si β . The first derivative of Fj is
∂Fj
Vi j E
∂C
=
= 
 2 > 0.
Í
∂pi j
∂pi j
N
G
p
+
E
i
i
j
i=1
0

We then have



L Il τSi −1 r β fRi | j (r)dr,

l=1

(36)
|2

  β1
respectively. By the condition dl,k ≥ rl | j ≥ r SSli
given in
(30), the Laplace transform L Il j can be calculated as


h
i
−1 β
L Il j τSi −1 r β = EIl j e−τSi r Il j

where (c) follows |hi,o ∼ exp(1). As discussed in Section III,
there exist two groups of interferences to u0 from the l-th tier:
the interference Il j from the BSs storing the requested content
j and the interference Il j − from the BSs without storing content
j. According to [28], the locations of these two groups of interfering BSs can be modeled as two thinned and independent
HPPPs Φl j with density pl j λl and Φl j − with density (1−pl j )λl ,

2T (τ,β)t j E
ÍN
, ai
( i=1
Gi pi j +E)3
matrix H N ×N of

∂2 F j
∂pi j ∂pk j

2
β

= −mai ak

(39)

< 0, where m =
2

= λi Si , and ak = λk Sk β . Then, the Hessian
Fj is given by

H N ×N , −m[ai ak ] N ×N = −m(aaT ),

(40)

where we define the column vector a , (a1, a2, ..., a N )T .
Obviously, the Hessian matrix H is a real symmetric negative
semi-definite matrix, i.e., Fj is concave of pi j . According to
the property that the summation of concave functions is also
0
a concave function, C is also a concave function of P. Thus,
Problem P1 is a concave optimization problem.

13

E. Proof of Lemma 3
By constructing the Lagrangian function of P1, we have
0

L(P, µ, λ, η) =C (P)+

N Õ
M
Õ

µi j p i j +

i=1 j=1

+

N
Õ

ηi (Qi −

N Õ
M
Õ

λi j (1 − pi j )

i=1 j=1
M
Õ

pi j ),

(41)

j=i

i=1

0

where µi j , λi j are the Lagrange multipliers associated with
(1), ηi is the Lagrange multiplier associated with (2), and we
use µ , (µi j ) N ×M , λ , (λi j ) N ×M , η , (ηi )i ∈N . Thus, we have
Vi j E
∂L(P, µ, λ, η)
+ µi j − λi j − ηi .
= ÍN
∂pi j
( k=1 G k pk j + E)2
∂L(P∗, µ, λ, η)
= 0, ∀i ∈ N, ∀ j ∈ M,
∂p∗i j

(43)

µi j p∗i j = 0, λi j (1 − p∗i j ) = 0, ∀i ∈ N, ∀ j ∈ M,

(44)

ηi (Qi −

p∗i j )

= 0, ∀i ∈ N,

(45)

j=1
M
Õ

p∗i j = Qi, 0 ≤ p∗i j ≤ 1, ∀i ∈ N, ∀ j ∈ M

(46)

j=1

µi j ≥ 0, λi j ≥ 0, ∀i ∈ N, ∀ j ∈ M.

(47)

According to (42) and (43), we have that
Vi j E
ηi = Í N
+ µi j − λi j
( k=1 G k p∗k j + E)2

ηi ≥

(

Vi j E
G k pk∗ j +E)2

ÍN

k=1,,i

0

1 ≤ j ≤ n, p∗j > p∗j for n + 1 ≤ i ≤ m − 1, and p∗j = p∗j = 0 or
0
0
p∗j > p∗j = 0 for m ≤ j ≤ M. As a result, we have p∗j ≥ p∗j ,
i.e., p∗j increases with Q.
G. Proof of Proposition 3
By letting x j =

2
β

λ i Si pi j

ÍN

i=1

2
β

ÍN

i=1

, it is clear that the objective

λ i Si

function of P1 is the same as that of P2. Hence, to prove
the first part of the proposition, it suffices to show that the
feasible set of P1 is a subset
Í Mof that of P2. Denote P = {P =
[pi j ] N ×M | 0 ≤ pi j ≤ 1,
j=1 pi j = Q i , ∀i ∈ N, ∀ j ∈ M}
as the feasible set of P1 and X = {x = (x j )1×M | 0 ≤ x j ≤ 1,
2
ÍN
β
ÍM
i=1 λi Si Qi
2 , ∀ j ∈ M} as the feasible set of P2.
j=1 x j =
β
i=1 λ i Si
2
β
λ i Si
2 1×N
ÍN
β
i=1 λ i Si

ÍN

∀k ∈ N, ∀ j ∈ M. (48)

As a result, we consider the following three cases to analyze
the optimal solution:
•

0

that when Q increases to Q , ηq∗ will decrease to η∗ to
Í
0
t j D(τ,β)
1
∗0
satisfy M
− D(τ,β)
j=1 p j = Q , and T (τ,β)
η∗
T (τ,β) increases
r
t j D(τ,β)
1
∗0
∗
− D(τ,β)
to T (τ,β)
0
T (τ,β) . Thus, we have p j = p j = 1 for
η∗
0

(42)

The KKT conditions are then written as follows

M
Õ

p∗i > p∗j , p∗i = p∗j = 0, i.e., p∗i ≥ p∗j . Thus, p∗j decreases with
j, i.e., p∗j increases with t j .
Then we prove that p∗j increases with the cache size Q.
Based on (13), without loss of generality, we assume p∗i = 1
for 1 ≤ i ≤ n and p∗j = 0 for m ≤ j ≤ M with
Í
Í
∗
1 ≤ n < m < M. Then, we have M
p∗j = n + m−1
j=1
j=n+1 p j =


q
Í
t j D(τ,β)
1
n + m−1
= Q. It can be seen
− D(τ,β)
j=n+1 T (τ,β)
η∗
T (τ,β)

: If 0 < p∗i j ≤ 1, then we have

0 ≤ λi j < µi j , which is in conflict with µi j = 0 due to
(44). If p∗i j = 0, we have λi j = 0 and µi j ≥ 0. Thus, we
have p∗i j = 0.

Let a = (

)

, then for any given element P ∈ P,

we can define x = aP. Since

∂x j
∂pi j

=

2
β

λ i Si
ÍN

l=1

2
β

> 0 and 0 ≤

λ i Si

pi j ≤ 1 for ∀i ∈Í N, ∀ j ∈ M, we have 0 ≤ x jÍ≤ 1 for
M
∀ j ∈ M. Due to M
j=1 pi j = Q i for i ∈ N , we have
j=1 x j =
ÍN

2
β

λ i Si

ÍM

j=1
2
ÍN
β
λ
S
l=1 i i

i=1

pi j

=

ÍN

i=1

ÍN

2
β

λ i Si Q i

i=1

2
β

. Thus, we have x ∈ X. This

λ i Si

first case,
and µi j = 0, λi j ≥ 0.
Summarizing above results,
have (12). Based on (46),
Í we thus
∗
we have that ηi satisfies M
j=1 pi j = Q i for ∀i ∈ N . Thus, we
complete the proof.

means that for any P ∈ P, it always holds that x = aP ∈ X.
Hence, the feasible set X of P2 includes the transformation of
the feasible set P of P1 as a subset. As a result, the optimal
value of P2 provides an upper bound of that of P1.
We next prove the second part of the proposition. If Qi = Q e
for ∀i ∈ N , then for any given element x ∈ X, we can always
construct an P0 where pi j = x j for
Í i ∈ N . Due
Í M to 0 ≤ pi j =
x j ≤ 1 for ∀i ∈ N, ∀ j ∈ M and M
p
=
i
j
j=1
j=1 x j = Q e for
∀i ∈ N , we thus have P0 ∈ P. This suggests that if x∗ ∈ X
is an optimal solution of P2, then P0 with pi j = x ∗j is also an
optimal solution of P1. Hence, the optimal values of P1 and
P2 are equal when Qi = Q e , ∀i ∈ N . The proposition is thus
proved.

F. Proof of Corollary 2

H. Proof of Theorem 1

•

(

ÍN

k=1,,i

Vi j E
G k pk∗ j +Gi +E)2

< ηi <

(

ÍN

Vi j E
G k pk∗ j +E)2

k=1,,i

< 1, we have µi j = λi j = 0 due to (44); If
have 0 ≤ µi j < λi j , which is in conflict with
to (44); If p∗i j = 1, we have 0 ≤ λi j < µi j and
which are contradictory. Thus, we have 0 < p∗i j
q

ÍN
Vi j E
∗ −E .
p∗i j = G1i
−
G
p
k
k=1,,i
ηi
kj
p∗i j

•

ηi ≤

(

ÍN

Vi j E
:
G k pk∗ j +Gi +E)2
∗
we have pi j = 1

: If 0 <

p∗i j
λi j

= 0, we
= 0 due
µi j = 0,
< 1 and

Similar to the analysis for the

k=1,,i

p∗j

First, we prove that
decreases with the index j. For any
two contents i, j ∈ M and i < j, their popularities satisfy ti ≥
t j . By (13), there are only three possible cases: p∗i = p∗j = 1,

0

According to (18), when C1 (P1,u ) is maintained constant
and λ j , S j and Q j are fixed values for ∀ j ∈ N, j , i, the
equivalent cache size Q e also remains unchanged. We only

14

need to consider the interactions of the i-th tier’s parameters,
thus we can derive that
2
2
ÍN
Í
β
β
2
Qe N
j=1,,i λ j S j −
j=1,,i λ j S j Q j
β
λi Si =
.
(49)
Qi − Qe
Substituting (21) into (49), (22) into (49) respectively, we have
(19) , (20).
I. Proof of Corollary 5
Equation (49) can be rewritten as following
!
2
ÍN
β
2
ÍN
j=1,,i λ j S j Q j
β
Qe − Í
2
j=1,,i λ j S j
N
j=1,,i

2

λi Siβ =

λj Sj

β

Qi − Qe

Obviously, when Q e >

ÍN

j=1,,i

ÍN

λj Sj

j=1,,i

2
β

λj Sj

Qj
2
β

.

(50)

, we have K1 > 0.

2
2
2
Í
ÍN
β
β
λi (Si ) β Qi + N
j=1,,i λ j S j Q j
j=1,,i λ j S j Q j
,
Q
>
e
2 Í
2
2
ÍN
β
β
λi (Si ) β + N
j=1,,i λ j S j
j=1,,i λ j S j
2
ÍN
λj Sj β Qj
to Qi > Íj=1,,i
. Meanwhile, because of
2
N
β
j=1,,i λ j S j

Based on Q e =
is equivalent
1

λi (Si ) β > 0, we have Qi > Q e . As a result, we have K1 > 0
and Qi − Q e > 0 when Qi ∈ (Q e, +∞), then K1 < 0 and
Qi − Q e < 0 when Qi ∈ [0, Q e ). Because the relationship of
λi and Qi is inversely proportional, thus λi decreases as Qi
increases when K1 > 0 and Qi − Q e > 0, while λi increases as
Qi increases when K1 < 0 and Qi −Q e < 0. Similarly, we have
Si decreases as Qi increases when K2 > 0 and Qi − Q e > 0,
while Si increases as Qi increases when K2 < 0 and Qi −Q e <
0. Finally, we have the corollary.
J. Proof of Corollary 6
According to (25)-(28), we consider the following cases:
1) When K3 > 0, from (25), we have Q e >
ÍN

k=1,,i, j

ÍN

λ k Sk

k=1,, j

2
β

λ k Sk

Qk

2
β

.

2) When K5 > 0 (K6 > 0), i.e., Q j > Q e , then based
on Q e =

ÍN

k=1,, j

λ k Sk

2
β

Q k +λ j S j

k=1,, j λ k Sk
2
2
ÍN
β Q +λ S β Q
λ
S
j j
e
k
k=1,, j k k
2
2
ÍN
β +λ S β
λ
S
j j
k=1,, j k k
2
ÍN
β
k=1,,i, j λ k Sk Q k
2
ÍN
β
k=1,, j λ k Sk

ÍN

2
β

+λ j S j

2
β

Qj

2
β

, we have Q e >
ÍN

, i.e., Q e >

k=1,, j

λ k Sk

ÍN

k=1,, j

2
β

λ k Sk

Qk
2
β

>

. Further, we have K3 > 0. Since λ j ≥ 0

3
and S j ≥ 0, we thus have Qi ∈ [0, K
K4 ].
3) When K5 < 0 (K6 < 0) and K3 > 0, i.e., Q j < Q e and

Qe >

ÍN

k=1,,i, j

ÍN

λ k Sk

k=1,, j

2
β

λ k Sk

Qk

2
β

3
, we thus have Qi ∈ [ K
K4 , +∞) due

to λ j ≥ 0 and S j ≥ 0.
4) When K3 < 0, i.e., Q e <
have Q e <

ÍN

k=1,, j

ÍN

λ k Sk

k=1,, j

2
β

λ k Sk

ÍN

k=1,,i, j

ÍN

λ k Sk

k=1,, j

Qk
2
β

2
β

λ k Sk

Qk

2
β

, then we

. Further, according to Q e =

ÍN

k=1,, j

ÍN

λ k Sk

k=1,, j

2
β

Q k +λ j S j

λ k Sk

2
β

+λ j S j

2
β

Qj

2
β

, we have Q j < Q e , i.e., K5 < 0

(K6 < 0). Thus, Qi ∈ [0, + ∞).
We thus have this corollary.
R EFERENCES
[1] K. Li, Z. Chen, and M. Tao, “Probabilistic caching in multi-tier heterogeneous networks: Optimization and tradeoff,” in Proc. IEEE 18th Int.
Workshop Signal Process. Adv. Wireless Commun. (SPAWC), Jul. 2017.
[2] Cisco, “Cisco visual networking index: Global mobile data traffic
forecast update 2015-2020,” white paper, Feb. 2016.
[3] J. G. Andrews, S. Buzzi, W. Choi, S. V. Hanly, A. Lozano, A. C. K.
Soong, and J. C. Zhang, “What will 5G be?” IEEE J. Sel. Areas
Commun., vol. 32, no. 6, pp. 1065–1082, Jun. 2014.
[4] H. ElSawy, E. Hossain, and M. Haenggi, “Stochastic geometry for
modeling, analysis, and design of multi-Tier and cognitive cellular
wireless networks: A survey,” IEEE Commun. Surveys and Tutorials,
vol. 15, no. 3, pp. 996–1019, Jun. 2013.
[5] X. Wang, M. Chen, T. Taleb, A. Ksentini, and V. Leung, “Cache in the
air: Exploiting content caching and delivery techniques for 5G systems,”
IEEE Commun. Mag., vol. 52, no. 2, pp. 131–139, Feb. 2014.
[6] E. Baştuğ, M. Bennis, and M. Debbah, “Living on the edge: The role
of proactive caching in 5G wireless networks,” IEEE Commun. Mag.,
vol. 52, no. 8, pp. 82–89, Aug. 2014.
[7] M. A. Maddah-Ali and U. Niesen, “Fundamental limits of caching,”
IEEE Trans. Inf. Theory, vol. 60, no. 5, pp. 2856–2867, May 2014.
[8] H. Liu, Z. Chen, and L. Qian, “The three primary colors of mobile
systems,” IEEE Commun. Mag., vol. 54, no. 9, pp. 15–21, Sep. 2016.
[9] H. Liu, Z. Chen, X. Tian, X. Wang, and M. Tao, “On content-centric
wireless delivery networks,” IEEE Wireless Commun., vol. 21, no. 6, pp.
118–125, Dec. 2014.
[10] J. G. Andrews, F. Baccelli, and R. K. Ganti, “A tractable approach to
coverage and rate in cellular networks,” IEEE Trans. Commun., vol. 59,
no. 11, pp. 3122–3134, Nov. 2011.
[11] E. Baştuğ, M. Bennis, M. Kountouris, and M. Debbah, “Cache-enabled
small cell networks: Modeling and tradeoffs,” EURASIP J. Wireless
Commun. Netw., vol. 2015, no. 1, p. 41, Feb. 2015.
[12] S. Tamoor-ul Hassan, M. Bennis, P. H. J. Nardelli, and M. Latva-Aho,
“Caching in wireless small cell networks: A storage-bandwidth tradeoff,”
IEEE Commun. Lett., vol. 20, no. 6, pp. 1175–1178, Jun. 2016.
[13] C. Yang, Y. Yao, Z. Chen, and B. Xia, “Analysis on cache-enabled
wireless heterogeneous networks,” IEEE Trans. Wireless Commun.,
vol. 15, no. 1, pp. 131–145, Jan. 2016.
[14] S. Tamoor-ul-Hassan, M. Bennis, P. H. J. Nardelli, and M. LatvaAho, “Modeling and analysis of content caching in wireless small cell
networks,” Proc. 2015 Int. Symp. Wireless Commun. Systems (ISWCS),
pp. 765–769, Aug. 2015.
[15] Z. Chen, J. Lee, T. Q. S. Quek, and M. Kountouris, “Cooperative caching
and transmission design in cluster-centric small cell networks,” IEEE
Trans. Wireless Commun., vol. 16, no. 5, pp. 3401–3415, May 2017.
[16] K. Shanmugam, N. Golrezaei, A. Dimakis, A. Molisch, and G. Caire,
“Femtocaching: Wireless content delivery through distributed caching
helpers,” IEEE Trans. Inf. Theory, vol. 59, no. 12, pp. 8402–8413, Dec.
2013.
[17] K. Poularakis, G. Iosifidis, and L. Tassiulas, “Approximation algorithms
for mobile data caching in small cell networks,” IEEE Trans. Commun.,
vol. 62, no. 10, pp. 3665–3677, Oct. 2014.
[18] B. Blaszczyszyn and A. Giovanidis, “Optimal geographic caching in
cellular networks,” in Proc. IEEE Int. Conf. Commun. (ICC), Jun. 2015,
pp. 3358–3363.
[19] Y. Chen, M. Ding, J. Li, Z. Lin, G. Mao, and L. Hanzo, “Probabilistic
small-cell caching: Performance analysis and optimization,” IEEE Trans.
Veh. Technol., vol. 66, no. 5, pp. 4341–4354, May 2017.
[20] S. H. Chae and W. Choi, “Caching placement in stochastic wireless
caching helper networks: Channel selection diversity via caching,” IEEE
Trans. Wireless Commun., vol. 15, no. 10, pp. 6626–6637, Oct. 2016.
[21] Y. Cui, D. Jiang, and Y. Wu, “Analysis and optimization of caching
and multicasting in large-scale cache-enabled wireless networks,” IEEE
Trans. Wireless Commun., vol. 15, no. 7, pp. 5101–5112, Jul. 2016.
[22] J. Song, H. Song, and W. Choi, “Optimal content placement for wireless
femto-caching network,” IEEE Trans. Wireless Commun., vol. 16, no. 7,
pp. 4433–4444, Jul. 2017.

15

[23] A. Khreishah, J. Chakareski, and A. Gharaibeh, “Joint caching, routing,
and channel assignment for collaborative small-cell cellular networks,”
IEEE J. Sel. Areas Commun., vol. 34, no. 8, pp. 2275–2284, Aug. 2016.
[24] Y. Cui and D. Jiang, “Analysis and optimization of caching and multicasting in large-scale cache-enabled heterogeneous wireless networks,”
IEEE Trans. Wireless Commun., vol. 16, no. 1, pp. 250–264, Jan. 2017.
[25] D. Liu and C. Yang, “Caching policy toward maximal success probability and area spectral efficiency of cache-enabled hetnets,” IEEE Trans.
Commun., vol. 65, no. 6, pp. 2699–2714, Jun. 2017.
[26] B. Serbetci and J. Goseling, “On optimal geographical caching in
heterogeneous cellular networks,” in Proc. 2017 IEEE Wireless Commun.
Netw. Conf. (WCNC), Mar. 2017, pp. 1–6.
[27] A. Liu and V. K. N. Lau, “How much cache is needed to achieve
linear capacity scaling in backhaul-limited dense wireless networks?”
IEEE/ACM Trans. Netw., vol. 25, no. 1, pp. 179–188, Feb. 2017.
[28] M. Haenggi and R. K. Ganti, “Interference in large wireless networks,”
Now: Foundations and Trends in Networking, vol. 3, no. 2, pp. 127–248,
2009.
[29] H.-S. Jo, Y. J. Sang, P. Xia, and J. Andrews, “Heterogeneous cellular
networks with flexible cell association: A comprehensive downlink SINR
analysis,” IEEE Trans. Wireless Commun., vol. 11, no. 10, pp. 3484–
3495, Oct. 2012.
[30] S. Boyd and L. Vandenberghe, Convex optimization.
Cambridge
university press, 2004.

Kuikui Li received the B.E. degree in communications engineering from the School of Telecommunications Engineering, Xidian University, Xi’an,
China, in 2016. He is currently pursuing the Ph.D.
degree with the Department of Electric Engineering,
Shanghai Jiao Tong University, China. His research
interests include cache-enabled heterogeneous networks, mobile edge computing, and cooperative
communications.

Chenchen Yang received the B. Eng. degree in
the School of Electronics and Information in 2013
from Northwestern Polytechnical University (NPU),
Xi’an, China. He is currently a Ph.D student in
the Institute of Wireless Communications Technology (IWCT), Department of Electronic Engineering,
Shanghai Jiao Tong University (SJTU), Shanghai,
China. From 2016 to 2017, he visited the Department of Electrical Engineering, Columbia University, USA, as a Research Scholar. His research
expertise and interests include caching networks,
cooperative communications and mobile edge computing.

Zhiyong Chen received the B.S. degree in electrical
engineering from Fuzhou University, Fuzhou, China,
and the Ph.D. degree from the School of Information
and Communication Engineering, Beijing University
of Posts and Telecommunications (BUPT), Beijing,
China, in 2011. From 2009 to 2011, he was a visiting
Ph.D. Student at the Department of Electronic Engineering, University of Washington, Seattle, USA.
He is currently an Associate Professor with the
Cooperative Medianet Innovation Center, Shanghai
Jiao Tong University (SJTU), Shanghai, China. His
research interests include mobile communications-computing-caching (3C)
networks, mobile VR/AR delivery and mobile AI systems. He currently serves
an Associate Editor of IEEE Access, and served as the Publicity Chair for the
IEEE/CIC ICCC 2014 and a TPC member for major international conferences.

Meixia Tao (S’00-M’04-SM’10) received the B.S.
degree in electronic engineering from Fudan University, Shanghai, China, in 1999, and the Ph.D.
degree in electrical and electronic engineering from
Hong Kong University of Science and Technology
in 2003. She is currently a Professor with the Department of Electronic Engineering, Shanghai Jiao
Tong University, China. Prior to that, she was a
Member of Professional Staff at Hong Kong Applied
Science and Technology Research Institute during
2003-2004, and a Teaching Fellow then an Assistant
Professor at the Department of Electrical and Computer Engineering, National
University of Singapore from 2004 to 2007. Her current research interests
include wireless caching, physical layer multicasting, resource allocation, and
interference management.
Dr. Tao currently serves as a member of the Executive Editorial Committee
of the IEEE T RANSACTIONS ON W IRELESS C OMMUNICATIONS and an
Editor for the IEEE T RANSACTIONS ON C OMMUNICATIONS. She was on the
Editorial Board of the IEEE T RANSACTIONS ON W IRELESS C OMMUNICA TIONS (2007-20110), the IEEE C OMMUNICATIONS L ETTERS (2009-2012),
and the IEEE W IRELESS C OMMUNICATIONS L ETTERS (2011-2015). She
also served as the TPC chair of IEEE/CIC ICCC 2014 and as Symposium
Co-Chair of IEEE ICC 2015.
Dr. Tao is the recipient of the IEEE Heinrich Hertz Award for Best
Communications Letters in 2013, the IEEE/CIC International Conference
on Communications in China (ICCC) Best Paper Award in 2015, and the
International Conference on Wireless Communications and Signal Processing
(WCSP) Best Paper Award in 2012. She also receives the IEEE ComSoc
Asia-Pacific Outstanding Young Researcher award in 2009.

