# Capacity Analysis of LTE-Advanced HetNets with Reduced Power Subframes and Range Expansion

## Abstract

The time domain inter-cell interference coordination techniques specified in LTE Rel. 10 standard improves the throughput of picocell-edge users by protecting them from macrocell interference. On the other hand, it also degrades the aggregate capacity in macrocell because the macro base station (MBS) does not transmit data during certain subframes known as almost blank subframes. The MBS data transmission using reduced power subframes was standardized in LTE Rel. 11, which can improve the capacity in macrocell while not causing high interference to the nearby picocells. In order to get maximum benefit from the reduced power subframes, setting the key system parameters, such as the amount of power reduction, carries critical importance. Using stochastic geometry, this paper lays down a theoretical foundation for the performance evaluation of heterogeneous networks with reduced power subframes and range expansion bias. The analytic expressions for average capacity and 5th percentile throughput are derived as a function of transmit powers, node densities, and interference coordination parameters in a heterogeneous network scenario, and are validated through Monte Carlo simulations. Joint optimization of range expansion bias, power reduction factor, scheduling thresholds, and duty cycle of reduced power subframes are performed to study the trade-offs between aggregate capacity of a cell and fairness among the users. To validate our analysis, we also compare the stochastic geometry based theoretical results with the real MBS deployment (in the city of London) and the hexagonal-grid model. Our analysis shows that with optimum parameter settings, the LTE Rel. 11 with reduced power subframes can provide substantially better performance than the LTE Rel. 10 with almost blank subframes, in terms of both aggregate capacity and fairness.

## 1Introduction

Cellular networks are witnessing an exponentially increasing data traffic from mobile users. Heterogeneous networks (HetNets) offer a promising way of meeting these demands. They are composed of small-size cells such as micro-, pico-, and femto-cells overlaid on the existing macrocells to increase the frequency reuse and capacity of the network. Since the base stations (BSs) of different tiers use different transmission powers and typically a frequency reuse factor of one, analyzing and mitigating the interference at an arbitrary user equipment (UE) is a challenging task.

### 1.1Related Work on Evaluation Methodology

Different approaches have been used in the literature for the performance evaluation of HetNets. The traditional simulation models with BSs placed on a hexagonal grid are highly idealized and may typically require complex and time-consuming system-level simulations. On the other hand, models based on stochastic geometry and spatial point processes provide a tractable and computationally efficient alternative for performance evaluation of HetNets [1]-[4]. *Poisson point process* (PPP) based models have been recently used extensively in the literature for performance evaluation of HetNets. However, as the macro base station (MBS) locations are carefully planned during the deployment process, PPP based models may not be viable for capturing real MBS locations, due to some points of the process being very close to each other. *Matern hardcore point process* (HCPP) provides a more accurate alternative spatial model for MBS locations. In HCPPs, the distance between any two points of the process is greater than a minimum distance predefined by hard core parameter. HCPP models are relatively more complicated due to the non existence of the probability generating functional [1]. Also, HCPP has a flaw of underestimating the intensity of the points that can coexist for a given hard core parameter [5]. Hence, HCPP models are not as tractable and simple as the PPP models.

With PPPs, using simplifying assumptions, such as Rayleigh fading channel model, and a path-loss exponent of four, we can obtain closed form expressions for aggregate interference and outage probability. Therefore, use of PPP models for performance evaluation of HetNets is appealing due to their simplicity and tractability [6]. Furthermore, the PPP based models provide reasonably close performance results when compared with the real BS deployments. In particular, results in [3] show that, when compared with real BS deployments, PPP and hexagonal grid based models for BS locations provide a lower bound and an upper bound, respectively, on the outage probabilities of UEs. Also, the PPP based models are expected to provide a better fit for analyzing denser HetNet deployments due to higher degree of randomness in small-cell deployments [2]. In this paper, due to their simplicity and reasonable accuracy, we will use PPP based models to characterize and understand the behavior of HetNets in terms of various design parameters.

### 1.2Use of PPP Based Models for LTE-Advanced HetNet Performance Evaluation

The existing literature has numerous papers based on the PPP model for analyzing HetNets. Using PPPs, the basic performance indicators such as coverage probability and average rate of a UE are analyzed in [7]-[10]. The use of range expansion bias (REB) in the picocell enables it to associate with more UEs and thereby improves the offloading of UEs to the picocells. The effect of REB on the coverage probability is studied in [11]. However, with range expansion, the offloaded UEs at the edge of picocells experience high interference from the macrocell. This necessitates a coordination mechanism between the MBSs and pico base stations (PBSs) to protect the picocell-edge UEs from the MBS interference. While [2] considers a homogeneous cellular network, [12] considers a HetNet with range expansion. The authors of [2] have obtained the information of real BS locations in an urban area from a cellular service provider. On the other hand, the authors of [13] have obtained the BS location information from an open source project [14] that provides approximate locations of the BSs around the world.

To mitigate the interference problems in HetNets, different *enhanced inter-cell interference coordination* (eICIC) techniques have been specified in LTE Rel. 10 of 3GPP which includes time-domain, frequency domain and power control techniques [15]. In the time domain eICIC technique, MBS transmissions are muted during certain subframes and no data is transmitted to macro UEs (MUEs). The picocell-edge users are served by PBS during these subframes (coordinated subframes) and thereby protecting the picocell-edge users from MBS interference. The eICIC technique using REB is studied well in the literature by analyzing its effects on the rate coverage [16] and on the average per-user capacity [18]. However, in the simulations of [20], the MBS transmits at reduced power (instead of muting the MBS completely) during the coordinated subframes (CSFs) to serve only its nearby UEs. Therein, the use of reduced power subframes during CSFs is shown to improve the HetNet performance considerably in terms of the trade-off between the cell-edge and average throughputs. Later on, *reduced power subframe* transmission have also been standardized under LTE Rel. 11 of 3GPP, and commonly referred therein as *further-enhanced ICIC* (FeICIC). In another study [21], simulation results show that the FeICIC is less sensitive to the duty-cycle of CSFs than the eICIC. In [22], 3GPP simulations are used to study and compare the eICIC and FeICIC techniques for different REBs and *almost blank subframe* densities. Therein, the amount of power reduction in the reduced power subframes is made equivalent to REB and its optimality is not justified.

### 1.3Contributions

In authors’ earlier work, analytic expressions using PPPs for coverage probability of an arbitrary UE is derived in [7] which has been extended to spectral efficiency (SE) derivations in [18] by considering eICIC and range expansion. Reduced power subframes, which are standardized in LTE Rel. 11 [23], are not analytically studied in literature to our best knowledge.

In the present work, generalized SE expressions are derived considering the FeICIC which includes eICIC and no eICIC as the two special cases. In this analytic framework that uses reduced power subframes and range expansion, expressions for the average SE of UEs and the 5th percentile throughput are derived. These expressions are validated through Monte Carlo simulations. Details of the simulation model are documented explicitly, and the Matlab codes can be accessed through [24] for regenerating the results. The optimization of key system parameters is analyzed with a perspective of maximizing both aggregate capacity in a cell and the proportional fairness among its users. Using these results, insights are developed on the configuration of FeICIC parameters, such as the power reduction level, range expansion bias, duty cycle of CSFs, and scheduling thresholds. The 5th and 50th percentile capacities are also analyzed to determine the trade-offs associated with FeICIC parameter adaptation. Further, we compare the 5th percentile SE results from PPP model with the real MBS deployment [25] and the hexagonal grid model.

## 2System Model

We consider a two-tier HetNet system with MBS, PBS and UE locations modeled as two-dimensional homogeneous PPPs of intensities , and , respectively. Both the MBSs and the PBSs share a common transmission bandwidth. The MBSs employ reduced power subframes, in which they transmit at reduced power levels to prevent high interference to the picocell UEs (PUEs). On the other hand, the PBSs transmit at full power during all the subframes.

The frame structure with reduced power subframes is shown in Figure 1. During uncoordinated subframes (USFs), the MBS transmits data and control signals at full power and during CSFs, it transmits at a reduced power , where is the power reduction factor. The PBS transmits the data, control signals and *cell reference symbol* with power during all the subframes. Setting corresponds to eICIC; and corresponds to no eICIC case.

Define as the duty cycle of USFs, i.e., ratio of the number of USFs to the total number of sub-frames in a frame. Then, is the duty cycle of CSF/reduced power subframes. Let and be the factors that account for geometrical parameters such as the transmitter and receiver antenna heights of the MBS and the PBS, respectively. Then, the effective transmitted powers of MBS during USFs is , MBS during CSFs is , and PBS during USF/CSF is . For an arbitrary UE, let the nearest MBS at a distance be its macrocell of interest (MOI) and the nearest PBS at a distance be its picocell of interest (POI). Then, assuming Rayleigh fading channel, the *reference symbol received power* from the MOI and the POI are given by,

respectively, where the random variables and account for Rayleigh fading. Define an interference term, , as the total interference power at a UE during USFs from all the MBSs and the PBSs, excluding the MOI and the POI. Similarly, define as the total interference power during CSFs. We assume that there is no frame synchronization across the MBSs and therefore irrespective of whether the MOI is transmitting a USF or a CSF, the interference at UE has the same distribution in both cases, and is independent of both and . Then, an arbitrary UE experiences the following four SIRs:

### 2.1UE Association

In and , it can be noted that and are directly affected by and hence their usage will make the cell selection process dependent on . Thus, we consider and to minimize the dependence of the cell selection process on .

The cell selection process using , and the REB can be explained with reference to Figure 2. If is less than , then the UE is associated with the MOI, otherwise with the POI. After the cell selection, the UE is scheduled either in USF or in CSF based on the scheduling thresholds (for MUE) and (for PUE). In macrocell, if is less than then the UE is scheduled to USF, otherwise to CSF. Similarly, in picocell, if is greater than then the UE is scheduled to USF, otherwise to CSF (to protect it from macrocell interference). The cell selection and scheduling conditions can be combined and formulated as:

A sample layout of MBSs and PBSs with their coverage areas for the four different UE categories are illustrated in Figure 3. Note that in the related work of [16], the UE association criteria are based on the average reference symbol received power at UE, where as our model is based on the SIR at UE, it also encompasses the FeICIC mechanism. In [16], the boundary between the USF-PUEs (picocell area) and the CSF-PUEs (range expanded area) is fixed due to the fixed transmit power of PBS. On the other hand, in our approach, the boundary between USF and CSF users can be controlled using in macrocell and in picocell, the parameters which play an important role during optimization as will be shown in Section 5.3.

Using -, it can be shown that the two SIRs and could be expressed in terms of and as,

Hence, knowing the statistics of and , particularly their *joint probability density function* (JPDF), would provide a complete picture of the SIR statistics of the HetNet system. We first derive an expression for *joint complementary cumulative distribution function* (JCCDF) of and in Section 3.1. Then we differentiate the JCCDF with respect to and to get the expression for JPDF in Section 3.2, which will then be used for spectral efficiency analysis.

## 3Derivation of Joint SINR Distribution

### 3.1JCCDF of and

From , we know that and are exponentially distributed with mean and , respectively. For brevity, substitute and in and :

Using it can be easily shown that the product has a maximum value of 1.

Let, and be the random variables denoting the distances of MOI and POI from a UE. Then, the JCCDF of and conditioned on is given by,

for . Here, , , and the integration limit The integration region of is graphically represented in Figure 4. By solving the integration as shown in Appendix ?, we can obtain a closed form expression for the conditional JCCDF as

for and , where is the Laplace transform of the total interference .

Expression for can be derived as follows. We assume that the interfering MBSs of a UE are frame asynchronous and subframe synchronous. Moreover, locations of the USFs and CSFs are uniformly randomly distributed, with a USF duty cycle of for all the MBSs. Hence, each interfering MBS transmits USFs with probability and CSFs with probability and the tier of MBSs can be split into two tiers, one tier of MBSs transmitting only USFs and other transmitting only CSFs. These two tiers are independent PPPs with intensities and . Therefore, the FeICIC scenario can be modeled using three independent PPPs as illustrated in Table 1.

BS type |
PPP |
Intensity |
Tx. power |
Distance of UE to nearest BS |
---|---|---|---|---|

USF-MBSs | ||||

CSF-MBSs | ||||

PBSs | ||||

Let, , , and be the interference at UE from all interfering USF-MBSs, CSF-MBSs and PBSs. Then, the total interference is . Using [26], parameters in Table 1, and assuming , we can derive the Laplace transform of in to be,

### 3.2JPDF of and

The conditional JPDF of and ,

can be derived by differentiating the JCCDF in with respect to and . Detailed derivation of conditional probability JPDF is provided in Appendix ?. Using the theorem of conditional probability we can write

where, the PDFs of and are and , respectively. We can then express the unconditional JPDF of and as,

where, we assume that a UE is served by a BS only if it satisfies the minimum distance constraints: UE should be located at distances of at least from the MOI and from the POI.

## 4Spectral efficiency analysis

In this section, the expressions for aggregate and per-user SEs categories are derived. Considering the JPDF of an arbitrary UE in , first the expressions for the probabilities that the UE belongs to each category are derived. Then, these expressions are used to derive the mean number of UEs of each category in a cell. These are followed by the derivation of the aggregate SE. Then, per-user SE expressions are obtained by dividing the aggregate SE by the mean number of UEs.

### 4.1MUE and PUE Probabilities

Depending on the SIRs and , a UE can be one of the four types: USF-MUE, CSF-MUE, USF-PUE or CSF-PUE. Given that the UE is located at a distance from its MOI and from its POI, probabilities of the UE belonging to each type can be found by integrating the conditional JPDF over the regions whose boundaries are set by the cell selection conditions in -. Based on these conditions the integration regions for different UE categories are shown in Figure 5.

The probability that a UE is a CSF-MUE can be found by integrating the JPDF over the region R1,

To form concise equations, let us define an integral function,

where, is a function of and , for is the integration region as defined in Figure 5. Then, can be written as,

Similarly, the conditional probabilities that a UE is a USF-MUE, USF-PUE or CSF-PUE are respectively given as,

### 4.2Mean number of MUEs and PUEs

Since the MBS locations are generated using PPPs, the coverage areas of all the MBSs resemble a Voronoi tessellation. Consider an arbitrary Voronoi cell. Let the number of UEs in the cell be , and number of CSF-MUEs in the cell be . Then, is a random variable and the mean number of CSF-MUEs is given by,

where in we use the fact that the probability that any of the UEs in a cell being a CSF-MUE is independent of . However, it is important to note that this is itself a consequence of our assumption that there is no limit on the number of CSF-MUEs per cell. Further, the event that any one of the UEs in a cell is a CSF-MUE is independent of the event that any other UE in that cell is a CSF-MUE, and all such events have the same probability of occurrence, namely given in . Then,

Using [27], it can be shown that the mean number of UEs in a Voronoi cell is . Therefore, the mean number of CSF-MUEs in a cell are given by,

Similarly, the mean number of USF-MUEs, USF-PUEs and CSF-PUEs are respectively given by,

### 4.3Aggregate and Per-user Spectral Efficiencies

We use Shannon capacity formula, , to find the SE of each UE type. The mean aggregate SE of an arbitrarily located CSF-MUE can be found by

Similarly, the mean aggregate SEs for USF-MUEs, USF-PUEs and CSF-PUEs can be respectively derived to be

where, . Then the corresponding per-user SEs are

### 4.45th Percentile Throughput

The 5th percentile throughput reflects the throughput of cell-edge UEs. Typically the cell-edge UEs experience high interference and analyzing their throughput provides important information about the fairness among the users in a cell and the system performance.

Consider the JPDF expression in . The integration regions of the JPDF for different UE categories are shown in Figure 5. The SIR PDF of USF-MUEs can be evaluated by integrating the JPDF over in the region R2,

for . The CDF expression can be derived as

for and, the CDF of throughput of the USF-MUEs can be derived as a function of in as,

for . By using the CDF plots, the 5th percentile throughput of USF-MUEs can easily be found as the value at which the CDF is equal to 0.05. Similarly, the 5th percentile throughput of other three UE categories can also be found.

## 5Numerical and Simulation Results

The average SE and 5th percentile throughput expressions derived in the earlier sections are validated using a Monte Carlo simulation model built in Matlab. Validation of the PPP capacity results for a HetNet scenario with range expansion and reduced power subframes is a non-trivial task. In this section, details of the simulation approach used for validating the PPP analyses are explicitly documented to enable reproducibility. Matlab codes for the simulation model and the theoretical analysis can be downloaded from [24].

### 5.1Simulation Methodology for verifying PPP Model

Algorithm used in simulation to find the aggregate and per-user SEs is described below.

The X- and Y-coordinates of MBSs, PBSs and UEs are generated using uniformly distributed random variables. The number of MBS and PBS location marks are and respectively, where, is the assumed geographical area that is square in shape as illustrated in Figure 6.

The UE locations are constrained within a smaller area which is aligned at center of the main simulation area to avoid the UEs to be located at the edges. In the PPP analysis, the area is assumed to be infinite. But in simulation, this scenario is approximated by making sufficiently larger than . The number of UEs is .

The MOI (closest MBS) and POI (closest PBS) for each UE is identified. The minimum distance constraints are applied by discarding the UEs that are closer than from their respective MOIs (POIs).

The SIRs , , , are calculated for each UE using -.

The UEs are classified as USF-MUEs, CSF-MUEs, USF-PUEs and CSF-PUEs using the conditions in -.

The MUEs (PUEs) which share the same MOI (POI) are grouped together to form the macro- and pico-cells.

The SEs of all the UEs are calculated. In a cell, SE of a USF-MUE is calculated using

. The SEs of other UE types are calculated using similar formulations.The aggregate capacity of each UE type is calculated in all the cells.

Mean aggregate capacity and mean number of UEs of each type are calculated by averaging over all the cells.

The per-user SE of each UE type are calculated by (mean aggregate capacity)/(mean number of UEs).

46 dBm, 30 dBm | |
---|---|

-11 dBm | |

35 m, 10 m | |

(marks/Km) | 4.6, , 200 |

Fading model, Path-loss exponent () | Rayleigh, 4 |

0.5, 4 dB, 0dB | |

### 5.2Per-user SEs with PPPs and Monte Carlo Simulations

The system parameter settings are shown in Table 2. The per-user SE results obtained using analytic expressions of - are shown in Figure 7(a) and Figure 7(b) for macrocell and picocell, respectively. The analytic plots agree well with the simulation plots, and provide the following insights:

#### USF- and CSF-MUEs

Referring to Figure 2, USF-MUEs form the outer part and CSF-MUEs form the inner part of the macrocell. As the REB increases, some of the USF-MUEs at the macro-pico boundary which have worse SIRs are offloaded to the picocell. Consequently, the mean number of USF-MUEs decreases and their per-user SE increases as shown in Figure 7(a).

The mean number of CSF-MUEs are not affected by as long as . Considering Figure 5, it can be noted that if , the line intersects the boundary of region R1. Hence, if is increased further such that , the area of R1 decreases and thereby decreases the mean number of CSF-MUEs. Therefore, the per-user SE of CSF MUEs remains constant as long as , and increases if crosses this limit as shown in Figure 7(a).

On the other hand, as the increases, the transmit power of all the interfering MBSs increases during CSFs, hence it increases the interference power at all the UEs. This causes the SIRs of USF-MUEs (), USF-PUEs () and CSF-PUEs () to decrease, which can be noted in , , and , respectively. However, the SIRs of CSF-MUEs () would increase (despite of increased interference) because of the increase in received signal power (due to higher ) which can be noted in . Considering and , since is a constant, the degradation in causes the number of USF-MUEs to increase and CSF-MUEs to decrease. Consequently, the per-user SE of USF-MUEs decreases and that of CSF-MUEs increases for increasing , as shown in Figure 7(a).

#### USF- and CSF-PUEs

As the REB increases, the mean number of USF-PUEs remains constant if because the area of region R4 in Figure 5 is unaffected by the value of . Therefore, the per-user SE of USF-PUEs also remain constant for increasing REB as shown in Figure 7(b). With increasing REB, some MUEs are offloaded to the picocell and become CSF-PUEs. But, these UEs are located at cell-edges and have low SIRs. Hence the per-user SE of CSF-PUEs decreases as shown in Figure 7(b).

On the other hand, as the increases, the transmit power of all the interfering MBSs increases during CSFs causing , and to decrease and to increase, as explained previously. Considering and , since is a constant the degradation in causes the number of USF-PUEs to decrease and CSF-PUEs to increase. Consequently, the per-user SE of USF-PUEs increases and that of CSF-PUEs decreases for increasing , as shown in Figure 7(b).

### 5.3Optimization of System Parameters to Achieve Maximum Capacity and Proportional Fairness

The five parameters are the key system parameters that are critical to the satisfactory performance of the HetNet system. The goal of these parameter settings is to maximize the aggregate capacity in a cell while providing proportional fairness among the users.

Consider an arbitrary cell which consists of UEs. Let be the capacity of an arbitrary UE . The sum of capacities (sum-rate) and the sum of log capacities (log-rate) in a cell are respectively given by,

Maximizing the corresponds to maximizing the aggregate capacity in a cell, while maximizing the corresponds to proportional fair resource allocation to the users of a cell [28], [29]. There can be trade-offs existing between aggregate capacity and fairness in a cell. Maximizing the may reduce the , and vice versa. In this section, we try to understand these trade-offs by analyzing the characteristics of and with respect to the variation of key system parameters.

We attempt to achieve the proportional fairness by optimizing the five key system parameters to maximize the . The variation of with respect to is shown in Figure 8, for . These plots are obtained through the Monte Carlo simulations and each surface plot is the variation of with respect to and for a fixed value of and . The optimum scheduling thresholds and that maximizes the are dependent on the values of and .

Figure 9 shows the plots of and as the functions of and . The markers show the simulation results while the dotted lines show the smoother estimation obtained using the curve fitting tool in MATLAB. For small values, the optimum threshold has higher values as shown in Figure 9(a), and according to this causes very few MUEs that have to be scheduled during CSFs. This makes sense because MBS transmit power during CSFs is very low for small and hence the number of CSF-MUEs which can be covered is also less. On the other hand, for higher values, MBS transmits with higher power level during CSFs and can cover larger number of CSF-MUEs. Therefore, to improve the fairness proportionally, the optimal value decreases with increasing so that more MUEs are scheduled during CSFs.

In the picocell, with increasing the CSF-PUEs at the cell edges will experience higher interference from the MBSs. Then, more PUEs should be scheduled during USFs to improve proportional fairness. Likewise, decreasing in Figure 9(b) indicates that more PUEs are scheduled during USFs as per .

The with optimum scheduling thresholds and is plotted in Figure 10. Higher the , better is the proportional fairness. It is important to note that the range expansion bias, , has a significant effect on proportional fairness. The increases from to when is increased from 0 db to 12 dB.

Compared to , has a smaller effect on the proportional fairness. When is set to zero which corresponds to the eICIC, is at its minimum. It shows that eICIC provides minimum proportional fairness. Figure 10 moreover shows that setting which corresponds to no eICIC, also does not provide maximum . An setting between 0.125 and 0.5 maximizes the and hence the proportional fairness.

The characteristics of with optimum scheduling thresholds is shown in Figure 11. As the increases, decreases, which is the opposite effect when compared to the in Figure 10. This shows the trade-off between the aggregate capacity and the proportional fairness. Increasing the would increase the proportional fairness but decrease the aggregate capacity, and vice versa.

Comparing Figures Figure 10 and Figure 11 also explains the trade-off associated with setting . A very small value, , provides larger but smaller , which is better from an aggregate capacity point of view. Setting is better from a fairness point of view. Any value of is not recommended since it degrades the aggregate capacity as shown in Figure 11, decreases the proportional fairness as shown in Figure 10, and consumes higher transmit power by the MBSs. Setting as in the eICIC case would reduce both and drastically.

The implications in Figures Figure 10 and Figure 11 can be seen from a different perspective by using the 5th percentile SE versus 50th percentile SE graph shown in Figure 12. It shows that increasing the from 6 dB to 12 dB improves the 50th percentile SE, but degrades the 5th percentile SE. This illustrates the trade-off between the 5th and 50th percentile SEs, which is analogous to the trade-off between aggregate capacity and proportional fairness among the users, as explained in the previous paragraphs. Figure 12 also shows that has a notable effect on the 5th and 50th percentile SEs.

### 5.4Impact of the Duty Cycle of Uncoordinated Subframes

In the results of Figures Figure 9–Figure 12, was set to 0.5 and we next show the effect of varying on and . Introducing into the optimization problem makes it difficult to visualize the results due to the addition of one more dimension. Therefore, we use the optimized scheduling thresholds, and , and analyze and as the functions , and . Figures Figure 13 and Figure 14 show the versus and the versus , respectively for different values of and . The variation of with respect to is not significant, except for . Whereas, the variation of with respect to is significant.

When , the value decreases rapidly for . Nevertheless, is shown to have poor performance in the previous paragraphs and hence it is not recommended. For other values of , variation in does not affect the significantly, which shows that by using a fixed value of , proportional fairness can be achieved by optimizing (to maximize ) the scheduling thresholds. Figure 14 shows that fixing approximately to 0.43 maximizes the irrespective of and , provided the scheduling thresholds are optimized to maximize .

In [16], the boundary of CSF-PUEs that form the inner region of picocell (excluding the range expansion region) is fixed due to the fixed transmit power of PBS. The *association bias* and *resource partitioning fraction* parameters are used as the variables to be optimized. It is analogous for us to have a fixed and optimize and . But in contrast, we fix the for simplicity and optimize the other four parameters, since coordinating among the cells through the X2 interface is complex and adds to communication overhead in the backhaul.

### 5.55th Percentile Throughput

Using the expressions derived in Section 4.4, the 5th percentile throughput versus for different is shown in Figure 15(a) for MUEs, and in Figure 15(b) for PUEs. As the increases, MBSs transmit at higher power level during CSFs and the UEs of all types experience a higher interference power. However, the received signal power at CSF-MUEs increases with and results in improved 5th percentile throughput as shown in Figure 15(a). But, the SIRs of USF-MUEs and USF/CSF-PUEs degrade due to higher interference and therefore their 5th percentile throughput decreases with increase in as shown in Figures Figure 15(a) and Figure 15(b).

Increasing the REB, , causes the USF-MUEs with poor SIR, located at the edge of macrocell, to be offloaded to the picocell and thereby increasing the 5th percentile throughput of USF-MUEs as shown in Figure 15(a). The offloaded UEs in picocell are scheduled during CSFs and due to their poor SIR the 5th percentile throughput of CSF-PUEs decreases as shown in Figure 15(b).

### 5.6Comparison with Real BS Deployment

We obtained the data of real BS locations in United Kingdom from an organization [25] where the mobile network operators have voluntarily provided the information of location and operating characteristics of individual BSs. The data set in [25] was last updated in May 2012, and it provides exact locations of the BSs. Also, the BSs of different operators can be distinguished.

In this section, we compare the 5th percentile SE results from the PPP model with that of the real BS deployment and hexagonal grid model. The real MBS locations of two different operators in a km area of London city were obtained from [25] as shown in Figure 16. In this area, the average BS densities of the two operators were found to be 1.53 MBSs/km and 2.04 MBSs/km. To have a fair comparison, the MBS locations for hexagonal grid and PPP models were also generated with the same densities. The PBS locations were generated randomly using another PPP model. The parameters dB, , , dB, dB, and dBm were fixed while the PBS density was varied to analyze its effect on the 5th percentile SE.

The plots of 5th percentile SE versus PBS density are shown in Figure 17 for the two operators. The 5th percentile SE of operator-2 is better than that of operator-1 since the former has higher MBS density. As expected, the 5th percentile SE improves with the increase in PBS density. It can also be observed that increasing the PBS transmit power from 10 dBm to 30 dBm will result into almost twice the 5th percentile SE. Since hexagonal grid model is an ideal case, it has the best 5th percentile SE and forms an upper bound. The PPP model has the worse 5th percentile SE and forms a lower bound. The real MBS deployment is usually planned and hence it is not completely random in nature. On the other hand, it is also not equivalent to the idealized hexagonal grid model due to the practical constraints involved during the deployment. Hence, the 5th percentile SE of real MBS deployment lies in between the two bounds of hexagonal grid and random deployments.

## 6Conclusion

In this paper, spectral efficiency and 5th percentile throughput expressions are derived for HetNets with reduced power subframes and range expansion. These expressions are validated using the Monte Carlo simulations. Joint optimization of the key system parameters, such as range expansion bias, power reduction factor, scheduling thresholds, and duty cycle of reduced power subframes, is performed to achieve maximum aggregate capacity and proportional fairness among users. Our analysis shows that under optimum parameter settings, the HetNet with reduced power subframes yields better performance than that with almost blank subframes (eICIC) in terms of both aggregate capacity and proportional fairness. However, transmitting the reduced power subframes with greater than half the maximum power proved to be inefficient because it degrades both the aggregate capacity and the proportional fairness. Increasing the range expansion bias improves the proportional fairness but degrades the aggregate capacity. In case of eICIC, the duty cycle of almost blank subframes has a significant effect on the fairness, but with reduced power subframes and optimized scheduling thresholds, duty cycle has a limited effect on fairness. Hence, fixing the duty cycle and optimizing the scheduling thresholds is preferable since it avoids the overhead of coordinating the duty cycle among the cells through the X2 interface. We also compared the 5th percentile SE results from PPP model with that of real BS deployment and hexagonal grid model. We observed that the hex grid model forms the upper bound while the PPP model forms the lower bound. Increasing the PBS density or the PBS transmit power would improve the 5th percentile SE.

### References

- H. Elsawy, E. Hossain, and M. Haenggi, “Stochastic geometry for modeling, analysis, and design of multi-tier and cognitive cellular wireless networks: A survey,”
*Communications Surveys Tutorials, IEEE*, vol. 15, no. 3, pp. 996–1019, Third 2013. - J. Andrews, F. Baccelli, and R. Ganti, “A tractable approach to coverage and rate in cellular networks,”
*IEEE Transactions on Communications*, vol. 59, no. 11, pp. 3122–3134, Nov. 2011. - R. Ganti, F. Baccelli, and J. Andrews, “A new way of computing rate in cellular networks,” in
*Proceedings of the IEEE Int. Conf. Commun. (ICC)*, Kyoto, Japan, June 2011, pp. 1–5. - J. Andrews, R. Ganti, M. Haenggi, N. Jindal, and S. Weber, “A primer on spatial modeling and analysis in wireless networks,”
*Communications Magazine, IEEE*, vol. 48, no. 11, pp. 156–163, November 2010. - M. Haenggi, “Mean interference in hard-core wireless networks,”
*Communications Letters, IEEE*, vol. 15, no. 8, pp. 792–794, August 2011. - M. Haenggi, J. Andrews, F. Baccelli, O. Dousse, and M. Franceschetti, “Stochastic geometry and random graphs for the analysis and design of wireless networks,”
*Selected Areas in Communications, IEEE Journal on*, vol. 27, no. 7, pp. 1029–1046, September 2009. - S. Mukherjee, “Distribution of downlink SINR in heterogeneous cellular networks,”
*IEEE J. Select. Areas Commun. (JSAC), Special Issue on Femtocell Networks*, vol. 30, no. 3, pp. 575–585, Apr. 2012. - R. Heath, M. Kountouris, and T. Bai, “Modeling heterogeneous network interference using poisson point processes,”
*IEEE Trans. Signal Process.*, vol. 61, no. 16, pp. 4114–4126, Aug. 2013. - T. Novlan, R. Ganti, and J. Andrews, “Coverage in two-tier cellular networks with fractional frequency reuse,” in
*Proceedings of the IEEE Global Telecommun. Conf. (GLOBECOM)*, Houston, TX, Dec. 2011, pp. 1–5. - H. S. Dhillon, R. K. Ganti, F. Baccelli, and J. G. Andrews, “Modeling and analysis of K-tier downlink heterogeneous cellular networks,”
*IEEE J. Select. Areas Commun. (JSAC), Special Issue on Femtocell Networks*, vol. 30, no. 3, pp. 550–560, Apr. 2012. - S. Mukherjee, “Downlink SINR distribution in a heterogeneous cellular wireless network with biased cell association,” in
*Proceedings of the IEEE Int. Conf. Commun. (ICC)*, Ottawa, Canada, June 2012, pp. 6780 –6786. - H.-S. Jo, Y. J. Sang, P. Xia, and J. Andrews, “Outage probability for heterogeneous cellular networks with biased cell association,” in
*Proceedings of the IEEE Global Telecommun. Conf. (GLOBECOM)*, Houston, TX, Dec. 2011, pp. 1–5. - =2plus 43minus 4 C.-H. Lee, C.-Y. Shih, and Y.-S. Chen, “ l@English Stochastic geometry based models for modeling cellular networks in urban areas,”
*l@English Springer Wireless Networks*, vol. 19, no. 6, pp. 1063–1072, 2013. [Online]. Available: http://dx.doi.org/10.1007/s11276-012-0518-0 =0pt - =2plus 43minus 4 website. [Online]. Available: www.opencellid.org =0pt
- D. López-Pérez, I. Güvenç, G. de la Roche, M. Kountouris, T. Q. Quek, and J. Zhang, “Enhanced inter-cell interference coordination challenges in heterogeneous networks,”
*IEEE Wireless Commun. Mag.*, vol. 18, no. 3, pp. 22–31, June 2011. - S. Singh and J. G. Andrews, “Joint resource partitioning and offloading in heterogeneous cellular networks,”
*CoRR*, vol. abs/1303.7039, 2013. - S. Singh, H. S. Dhillon, and J. G. Andrews, “Offloading in heterogeneous networks: Modeling, analysis, and design insights,”
*IEEE Trans. Wireless Commun.*, vol. 12, no. 5, pp. 2484–2497, May 2013. - S. Mukherjee and I. Guvenc, “Effects of range expansion and interference coordination on capacity and fairness in heterogeneous networks,” in
*Proceedings of the IEEE Asilomar Conf. Sig., Syst., Computers*, vol. 1, Monterey, CA, Nov. 2011, pp. 1855–1859. - A. Merwaday, S. Mukherjee, and I. Guvenc, “On the capacity analysis of hetnets with range expansion and eICIC,” in
*Proceedings of the IEEE Global Commun. Conf. (GLOBECOM)*, Atlanta, GA, Dec 2013. - Panasonic, “Performance study on ABS with reduced macro power,” 3GPP TSG-RAN WG1, Tech. Rep. R1-113806, Nov. 2011.
- A. Morimoto, N. Miki, and Y. Okumura, “Investigation of inter-cell interference coordination applying transmission power reduction in heterogeneous networks for LTE-advanced downlink,”
*IEICE Trans. on Commun.*, vol. E96-B, no. 6, pp. 1327–1337, June 2013. - M. Al-Rawi, J. Huschke, and M. Sedra, “Dynamic protected-subframe density configuration in LTE heterogeneous networks,” in
*Proceedings of the IEEE Int. Conf. Computer Commun. Net. (ICCCN)*, Munich, Germany, 2012, pp. 1–6. - =2plus 43minus 4 “Overview of 3gpp release 11 v0.1.7,” Dec. 2013. [Online]. Available: http://www.3gpp.org/ftp/Information/WORK_PLAN/Description_Releases/Rel-11_description_20131224.zip =0pt
- =2plus 43minus 4 Mpact lab data management. [Online]. Available: http://www.mpact.fiu.edu/data-management/ =0pt
- =2plus 43minus 4 Sitefinder website. [Online]. Available: http://www.sitefinder.ofcom.org.uk =0pt
- S. Mukherjee, “UE coverage in LTE macro network with mixed CSG and open access femto overlay,” in
*Proceedings of the IEEE Int. Conf. Commun. (ICC) Workshops*, Kyoto, Japan, June 2011, pp. 1–6. - S. M. Yu and S.-L. Kim, “Downlink capacity and base station density in cellular networks,” in
*Proceedings of the IEEE SpaSWiN workshop (in conjunction with WiOpt)*, Tsukuba Science City, Japan, May 2013, pp. 119–124. - P. Viswanath, D. Tse, and R. Laroia, “Opportunistic beamforming using dumb antennas,”
*IEEE Trans. Inf. Theory*, vol. 48, no. 6, pp. 1277–1294, Aug. 2002. - M. R. Jeong and N. Miki, “A comparative study on scheduling restriction schemes for lte-advanced networks,” in
*Proceedings of the IEEE 23rd Int. Symp. Personal Indoor and Mobile Radio Communications (PIMRC)*, Sydney, Australia, Sept. 2012, pp. 488–495.