Macro-pico amplitude-space sharing with layered interference alignment

Inter-cell interference leads to severe performance degradation in cellular networks, and the study of multi-user interference channel is the corner stone for solving this problem. Amplitude-space layered interference alignment (IA), as an effective complementation to the vector-space IA, is a promising method to increase the data rate in static interference channels. However, recent studies of layered IA has been focused on analyzing the degrees of freedom (DoF) or the achievable rate under specific channel constraints. In this paper, we propose a layered IA scheme that can work with arbitrary channel coefficients. We develop a layer partitioning method and optimize the active layer assignment through linear programming. An implementation scheme is then introduced with multi-level nested lattice codes where the signal and interference are nested in amplitude space, and the interference from different users is nestedly aligned. The performance of the proposed scheme is finally evaluated in homogeneous and heterogeneous cellular networks with practical settings.


Introduction
Multi-tier heterogeneous networks can achieve significant areal capacity gain by intensifying spatial reuse of the spectrum [1,2]. While the high-power macro-cell base station (BS) provides basic coverage and supports high mobility, the low-power nodes like pico-cell BSs support high-capacity transmission for hotspot zones. Accordingly, the data rate requirements of macro-cell user and pico-cell user may have large difference.
Operated in the same frequency band, the interference in heterogeneous networks is complicated, including not only the co-tier interference among macro-cells and among pico-cells but also the cross-tier interference between the macro-cell and pico-cells [3]. One picouser may encounter interference from the macro-BS and the pico-BSs in the same macro-cell and also from the macro-BSs in adjacent macro-cells. The pico-cells might be deployed at any places in the macro-cell, and the transmit power of the macro-BS is much stronger than the pico-BS. Conventional methods to mitigate inter-cell interference mainly involve various orthogonalization-based processing. For example, almost blank subframe (ABS) is a time orthogonal interference coordination scheme; fractional frequency reuse (FFR) is a frequency orthogonal interference coordination scheme; coordinated beamforming (CB) is a spatial orthogonal interference coordination scheme [4][5][6][7]. The vector-space interference alignment (IA) scheme is an advanced orthogonalizationbased interference coordination scheme, where multiple interference are aligned in a subspace while leave the orthogonal signal subspace interference-free [8,9]. However, the orthogonalization-based processing schemes are only efficient when the interference has similar strength with the desired signal. If the interference is very weak or very strong, it is a waste of resource to provide orthogonal subspaces for each interference. In fact, for very weak interference channels, treating the interference as noise is optimal [10][11][12][13]. For strong interference channels, interference cancelation can achieve the capacity [14][15][16].
In multi-tier heterogeneous networks with different interference strengths in various scenarios, amplitude-space sharing of the signal and interferences is an effective way to complement the weakness of the orthogonalization-based processing. In [17], an amplitude-space sharing scheme in a two-cell network was proposed, where each BS schedules one user and the two users simultaneously scheduled in two overlapping cells consist a two-user interference channel. In different interference scenarios, the opportunities for interference cancelation are proactively created and the network sum rates are maximized through the optimized Han-Kobayashi coding [18,19]. In multi-cell networks, Han-Kobayashi coding is no longer applicable [20], we need to study amplitude-space sharing scheme that can accommodate multiple interference. In this case, the desired signal and interference should occupy different layers in the amplitude space, and multiple interference may be aligned in one layer to compress their occupied spaces.
Layered interference alignment had been studied in multi-user interference channels, but the existing results are obtained under special channel conditions and are not applicable in practical systems. In [21], the degrees of freedom (DoF) of the fully connected K-user Gaussian interference channel were proved to be arbitrarily close to K/2, but the channel coefficients should take specialized forms that the inputs are shifted by an even amount on the desired links and by an odd amount on the interference links (or vice versa). In [22], it was shown that K/2 DoF are achievable when the cross-link channel gains corresponding to the interferers are rational, whereas the direct-link channel gains corresponding to the intended signal are irrational algebraic; if the direct-link channel gains are rational as well, the DoF of the channel are strictly smaller than K/2. Later, a general formula for the DoF of real interference channels was derived in [23] by maximizing the Rényi information dimension. Recently, the authors of [24] proved that the DoF of the static interference channel are K/2 for almost all sets of channel gains through the real interference alignment. In [25], incorporating both vector and real interference alignment, the total DoF of 2KM/(K + 1) are characterized in the K × 2 and 2 × K, M antennas X channels, for almost all channel realizations. The achievability schemes in [22,24,25] rely on some results from the field of Diophantine approximation in number theory, especially the lower bounds on the approximability of irrational algebraic numbers by rationals [26][27][28][29][30][31][32]. However, these results also imply that in asymptotic high signal-to-noise ratio (SNR) region the capacity characterization of the interference channel is extremely sensitive to slight variations of the channel gains.
The generalized degrees of freedom (GDoF) of the symmetric Gaussian K-user interference channel were studied in [33], where all direct links have the same SNR and all cross-links have the same interference-tonoise ratio (INR). The authors found that the GDoF per user do not depend on the number of users, so that the characterization is identical to the two-user interference channel with the exception of a singularity when the INR equals to the SNR. The achievable rates in moderate SNR values were obtained in [34,35], where the interference alignment scheme is suitable for a class of integer-interference channel, where all cross-link channel gains are integer or rational.
In this paper, we first study the K-user layered interference alignment scheme with arbitrary channel coefficients and then apply the scheme to multi-cell heterogeneous cellular networks with practical SNR values. The main contributions of this paper are as follows. 1) With arbitrary channel coefficients, a novel layer partitioning method is designed based on the power level intersections between the signals and interference. 2) To maximize the network sum rate, a linear programming method was proposed to optimize the assignment of active transmit layers. 3) A nested lattice coding scheme is developed to implement the encoding and decoding when the power levels of each interference are not exactly aligned. 4) The performance of the proposed scheme applied in cellular networks is evaluated, and the affecting factors in different network scenarios are analyzed.
The rest of this paper is organized as follows: In Section 2, we first study the optimal active layer assignment under deterministic interference channel models. Then, based on the obtained insight, we study K-user Gaussian interference channel with arbitrary coefficients and study the layer partitioning and active layer assignment methods in Section 3. The encoding and decoding schemes are developed in Section 4, and the performance in cellular networks is evaluated in Section 5. Finally, Section 6 concludes the paper.

Deterministic channel model
It is very hard to directly study the interference network problem under Gaussian channels. In this paper, we first resort to the deterministic channel model to gain some insights on the layered interference alignment problem and then extend the idea to general Gaussian interference channels. The deterministic channel model is an approximation methodology developed to solve network information theory problems [36,37]. Its general principle is that the coding problem in a noisy channel is first approximated by a noiseless problem, then analyze the simplified problem, and use insights obtained from the simplified problem to find new achievable schemes and/or outer bounds of the original problem [38]. The deterministic channel model approximates the Gaussian channel as a discrete set of parallel noiseless channels.
In the single-user Gaussian channel, a real-valued input x generates a real-valued output y that is degraded by Gaussian noise z, that is where h is the channel coefficient, E[ x 2 ] = P, and the variance of z is N 0 . The SNR is thus defined as h 2 P/N 0 = γ . If x and z are normalized, then the effective channel gain is √ γ .
To transform the Gaussian channel to a deterministic channel, we first represent the normalized x in a base-2 notation as where each bit b i can be interpreted as occupying a signal level, and the most significant bit coincides with the highest level. Sincex is a normalized value,x < 1, all the bits b i are in fractional part. Given the SNR γ , the output of the deterministic channel model is where α = log 2 γ is the largest integer not exceeding log 2 γ and b α is the smallest signal level containing a transmit power larger than 1/2 α . In other words, the input bit sequence is shifted by α positions and the remained part after b α is truncated due to the degradation of noise.
In the K-user Gaussian interference channels, the inputs of K users form a vector x, and the output vector is where the entry of channel matrix H i,j stands for the channel gain from transmitter j to receiver i. The noise of different users is assumed to be independent and identically distributed (i.i.d.), and E[ zz H ] = N 0 I. The SNRs depend on the direct channel gains and are defined as γ k,k = H 2 k,k P k /N 0 , and the INRs depend on the crosschannel gains and are defined as γ i,j = H 2 i,j P j /N 0 . These SNRs and INRs compose a link quality matrix , whose Defining α i,j = log 2 γ i,j , and applying (3) to every direct and cross-links of (4), we obtain the deterministic model of interference channel. The input-output relationship is shown in the left part of Fig. 1, for the convenience of demonstration; the input levels are also shifted by α k,k positions so that the inputs and outputs of the direct links are at the same levels. For the cross-link from transmitter j to receiver i, the shifting is thus changed to α i,j − α j,j instead of α i,j . At each receiver, the direct-link outputs (signal) and the cross-link outputs (interference) are added together when they are on the same level, and the outputs below the noise level are discarded. Rather than the normal integer addition, modulo-2 addition is chosen here to avoid interaction between different levels. This is a simplified operation only applied in the deterministic model. In Gaussian channels, the addition of signals and interference will produce carry-over problem, which will be addressed in Section 3.

Alignment in amplitude space
From the left part of Fig. 1, we can see that, if all the transmit levels are active, there must be collisions between signals and interference. To avoid collision, only part of the transmit levels can be active. Therefore, we need to study an optimal active level assignment algorithm to avoid collision and simultaneously maximize the network sum rate.
In the K-user deterministic interference channel, each active level conveys one bit. To avoid collision, each level at the receiver side should be occupied either by a signal bit or by an interference bit. Actually, this interference bit can be a modulo-2 addition of several interference bits, since it is not necessary to decode each of these interference bits separately. In order to leave as many receive levels as possible "free" for reception of signal bits, it will be advantageous if multiple interference bits fall on one receive level. This is the principle of interference alignment in amplitude space: the alignment of interference saves amplitude levels for desired signals and thus the throughput can be increased.
The optimal active level assignment problem can be formulated as where b k,l indicates whether the l-th transmit level of user k is active or not, and l = l + (α k,j − α j,j ). When the level is active, b k,l = 1, otherwise b k,l = 0. The second constraint means that, if the l -th level of user j interferes with the l-th level of user k through the cross-link, b k,l and b j,l can not be "1" simultaneously; otherwise, the received signal is not decodable. Since the transmit levels of user j are shifted by α k,j − α j,j positions when they arrive at receiver k, the l -th transmit level of user j will fall on level l +(α k,j −α j,j ) at receiver k. Thus,when l = l +(α k,j −α j,j ), b j,l is an interference to b k,l , and these two bits are not allowed to be "1" simultaneously. If there is another user i and its l -th level b i,l also interfere with b k,l , similarly there will be a constraint that b k,l +b i,l ≤ 1. However, we do not restrict the summation of the interference bits, i.e., b j,l +b i,l can be any value even when they fall on the same level at receiver k. If both b i,l and b j,l are active, where l +(α k,i −α i,i ) = l +(α k,j −α j,j ) = l, the constraint can still be satisfied as long as b k,l = 0.
The optimization problem in (5) is a binary integer programming problem, which can be solved efficiently by a linear programming (LP)-based branch-and-bound algorithm [39]. The algorithm searches for an optimal solution to the binary integer programming problem by solving a series of LP relaxation problems, in which the binary integer requirement on the variables is replaced by the weaker constraint 0 ≤ b k,l ≤ 1.
One example of the optimized results is shown in the right part of Fig. 1. We can see the features of an amplitude-aligned interference network. From the transmitter side, each user has different active level assignment, and the active levels may not be contiguous. At the receiver side, several interference links might be aligned on one level, and the interference levels might be above, below, or interlaced with the signal levels.

Interference alignment in Gaussian channels
The deterministic channel model described above assumes that all SNRs and INRs are integer on the log scale, i.e., the signals fall exactly on the evenly spaced "levels" at the receiver. In Gaussian channels, however, this assumption is not valid. In the following, we will address the interference alignment problem taking into account arbitrary SNR and INR values.
Given the link quality matrix , at each receiver, there is one collision pattern in amplitude space, i.e., part of the desired signal may collide with part of the interference. The collision pattern changes at different receivers because they experience different SNR and INRs. The layer partitioning requires to find all the possible collision areas in the amplitude space between the signals and interference.
The layer partitioning procedure is shown in Fig. 2. For each user k, the SNR of the direct channel is γ k,k , which is represented by a bar with a height log γ k,k at the transmitter side. At receiver k, there are K bars representing the K received signals. Their relative positions in the amplitude space are affected by the corresponding SNRs and INRs. For user i's bar, the upper and lower boundaries are log γ k,i and log γ k,i − log γ i,i , respectively. If the bar of user j, j = k, overlaps with the bar of user k, both bars are split (separated into different layers) at the intersecting boundary positions. Repeating these processes at every receiver, we will obtain at most 2K(K − 1) layers for each transmit signal. The results shown in Fig. 2   After layer partitioning, user k has L k layers in the transmitter side, and the upper and lower boundaries of the l-th layer are respectively P k,l and P k,l−1 . Thus, the l-th layer has an amplitude space represented by ρ k,l = P k,l /P k,l−1 . Assume now that we have a layered transmission scheme where the bits in the highest layers are decoded first. Thus, the decisions of the current layer can at most be disturbed by as-yet undecoded layers, i.e., the layers with a lower power than the current one. Thus, ρ k,l can be viewed as a signal-to-interference ratio (SIR) of this layer. If Shannon capacity achieved transmission were used, R k,l = 1/2 log(1 + ρ k,l ) bits could be transmitted in this layer. Considering the encoding and decoding methods to be introduced in Section 4, which have a rate loss of at most 0.5 bit, we use R k,l = 1/2 log(ρ k,l ).
As in the deterministic model of Section 2, some layers must be inactive to avoid collision between the layered signals of different users. We can still apply the binary integer programming algorithm to search the optimal active layer assignment, so that the total throughput can be maximized. In Gaussian channels, since each layer has different amplitude space, the objective function should be weighted by the transmission capability of each layer, i.e., R k,l . The optimization problem is therefore formulated as where b k,l denotes the transmit state of layer l of user k. When this layer is active, b k,l = 1, otherwise b k,l = 0. The second constraint is to avoid collision between the desired signal and interferences. If the l -th layer of user j interferes with the l-th layer of user k through the cross-link, b k,l and b j,l can not be "1" simultaneously. Unlike in the deterministic model, the relation between l and l does not have an explicit expression here. This optimization problem can also be solved by the LPbased branch-and-bound algorithm. The search result is shown in Fig. 3, where the contiguous active layers are combined. We can see that some interference layers overlap and occupy the same part of the amplitude space, but they may not be strictly aligned with respect to the upper and lower boundaries. Since we need to decode the superimposed interference, this kind of alignment complicates the encoding designs. We will discuss in detail the encoding and decoding schemes in Section 4. Different from the deterministic channels, the addition of two layers of interference has a higher signal level than any of the two layers, i.e., the carry-over problem must be considered. The active layers assigned through (7) might be tightly connected, which means that at one receiver the lower boundary of an upper signal layer might be the upper boundary of a lower interference layer. If the interference layer is superimposed from two or more users, the carry-over part might collide with the upper signal layer. Therefore, at the bottom of the signal layer, we need to reserve some amplitude space for the carry-over interference. In practice, this can be done by retaining the transmit power and reducing the data rate of the upper signal layer.
For example, at receiver k, suppose there are two interference layers below the signal layer l. One is from user j, and the upper boundary is h 2 k,j P j,l ; the other is from user i, and the upper boundary is h 2 k,i P i,l . Thus, the lower boundary of the signal layer l at the transmitter side is changed intõ P k,l−1 = max P k,l−1 , where the first term inside the maximum operation is the upper boundary of the (l − 1)-th layer of user k, and the second term is the sum power of two layers of interference divided by the square of the direct-link channel gain. If the second term is larger, at receiver k, the lower boundary of the l-th layer of user k is changed from h 2 k,k P k,l−1 to h 2 k,j P j,l + h 2 k,i P i,l . The reserved space is at most 3 dB and thus the data rate loss is at most 0.5 bit. When the lower layer is occupied by the superimposed interference from K users, the reserved space is at most 10 log 10 (K) dB and the data rate loss is at most 1/2 log(K) bits. This data rate loss can be neglected compared to the increasing sum rate when SNR goes to infinity.

Implementation by multi-level nested lattice codes
In the layered interference alignment scheme, since the superimposed interference layers need to be decoded, random coding is no longer applicable. We thus present a coding scheme with multi-level nested lattice codes in this section. Lattice coding is a structural coding; if there are two codewords that are selected from a lattice, their sum and difference are also within the same lattice. Thus, we can directly decode the superposition of the aligned interferences, instead of decoding the interference signals one by one. A lattice is an n-dimensional discrete subgroup of the Euclidean space R n under vector addition. Thus, if λ 1 and λ 2 are in , their sum and difference are also in . A lattice 2 is said to be nested in a lattice 1 if 2 ⊆ 1 . The lattice 1 is often referred to as a fine lattice and 2 as a coarse lattice. A nested lattice code L is the set of all points of the fine lattice that are within the fundamental Voronoi region V 2 of the coarse lattice, i.e., L = 1 ∩ V 2 .
For any two powers P a ≥ P b ≥ 0, [40] shows that there exist nested n-dimensional lattices 2 ⊆ 1 , as n → ∞, their second moment σ 2 ( 2 ) = P a , σ 2 ( 1 ) = P b , and the coding rate of L satisfies where |L| is the cardinality of set L and Vol(V i ) is the volume of Voronoi region V i . In [41], it is shown that nested lattice codes can achieve the capacity of point-to-point AWGN channels. In [42], a doubly nested lattice coding scheme was provided which can approach the capacity region of a two-way relay channel within 0.5 bit. Reference [43] provides a practical implementation scheme for nested lattice coding, where turbo coding and trellis shaping (multidimensional quantization) are involved. Due to space limitation, we refer the interested reader to [44,45] for the detailed definitions and general construction methods of lattice codes.
The codebook used by level l is a nested lattice code C k,l = k,l−1 ∩ V k,l , and the codeword c k,l ∈ C k,l . The rate of this code is The transmitted signal of user k is the summation of all the codewords where the corresponding layers are active, i.e., The received signal at user k is the superposition of the signals from all transmitted users, which is However, it has decoding problem if we just use L k layers of nested lattice code for each user k. As shown in Fig. 3, the received signals and interference are in layers and might be interlaced. If an interference layer is above a signal layer, we need first to decode and cancel the interference layer before decoding the signal layer. If the interference from two or more users exist above a signal layer, we need to decode the superimposed interference irrespective of whether the boundaries of interference layers are aligned or not.
For example, at receiver k, above the l-th signal layer there is a superimposed interference layer that is added by the l -th layer of user j and the l -th layer of user i. The upper boundary of the signal layer is h 2 k,k P k,l , and the lower boundaries of two interference layers are respectively h 2 k,j P j,l −1 and h 2 k,i P i,l −1 , where we have h 2 k,j P j,l −1 ≥ h 2 k,k P k,l and h 2 k,i P i,l −1 ≥ h 2 k,k P k,l . The received signal codeword is h k,k c k,l ∈ h k,k k,l−1 , and the received interference codewords are h k,j c j,l ∈ h k,j j,l −1 and h k,i c i,l ∈ h k,i i,l −1 . To decode the superimposed interference h k,j c j,l + h k,i c i,l , we require that the lattices h k,j j,l −1 and h k,i i,l −1 are aligned or nested. Through adjusting the transmit powers of related layers, we can make h k,j j,l −1 and h k,i i,l −1 aligned at receiver k, but it is impossible to simultaneously make h k ,j j,l −1 and h k ,i i,l −1 aligned at another receiver k .
To solve this problem, we just require a relationship that h k,j j,l −1 and h k,i i,l −1 are nested at receiver k. Since h 2 k,j P j,l −1 ≥ h 2 k,k P k,l and h 2 k,i P i,l −1 ≥ h 2 k,k P k,l , we can add a level of nested lattice j,m with second moment h 2 k,k /h 2 k,j P k,l at transmitter j, and add a level of nested lattice i,m with second moment h 2 k,k /h 2 k,i P k,l at transmitter i. The lattice nesting relation becomes j,l −1 ⊆ j,m ⊆ j,l −2 and i,l −1 ⊆ i,m ⊆ i,l −2 . Thus, both h k,j c j,l and h k,i c i,l are nested in h k,k k,l . Obviously, their sum is also nested in h k,k k,l . The superimposed interferences can then be decoded. Similarly, at receiver k , if there are interference layers above the signal layers and the interference layer is superimposed from multiple users, we can add more levels of nested lattice at the transmitter side to make the received interference codewords nested.
After these level insertion operations, the nested lattice levels increase to L + k for user k, but the number of active signal layers does not change. The decoded codeword of the l-th signal layer iŝ i.e., the received signal r k is first quantized to lattice k,l−1 and then taken a modulus operation relative to lattice k,l . The quantization operation is to mitigate the interference of lower layers, and the modulus operation can mitigate the interference of the upper layers by decoding and cancelation.
The developed implementation scheme is based on the idea of nested alignment and can guarantee the superimposed interference decodable at every receiver with any kind of power relationships.

Simulation results
We first study the performance of the layered IA scheme in deterministic interference channels and compare it with the orthogonal transmission scheme, and a virtual IA scheme where each user can achieve 1/2 DoF of the channel.
The virtual IA scheme demonstrated here is not a real transmission scheme. Since currently there is only DoF result for the layered IA transmission with arbitrary channel coefficients, we use this result to provide a performance benchmark of the achievable sum rate. In the virtual IA scheme, we compute the achievable rate of each user based on its SNR and without taking into account the impact of interference, and the sum rate of K users is obtained by the summation of each user's achievable rate and divided by 2. This method reflects the fact that each user can exploit 1/2 DoF of the channel. If there is only one user, this user will use the full DoF of the channel. The sum rate obtained in this way is labeled as "K/2 DoF" in Figs. 4 and 5. When SNR approaches infinity, the DoF result can reflect the upper bound of channel capacity. However, in moderate SNR values the achieved rate obtained from the DoF may have large difference with the capacity.
In deterministic models, the logarithmic channel gains are integer. We set the logarithmic SNRs and INRs as random integer selected from 1∼6. In K-user orthogonal transmission scheme, each user occupies 1/K of the time or frequency resources no matter what the SNRs and INRs are. The corresponding results are shown in Fig. 4. We can see that the sum rate of the virtual IA scheme linearly increases with the number of users when K ≥ 2, but its performance is not the best for all cases. In twouser or three-user interference channels, the proposed layered IA scheme can achieve higher sum rate. According to the optimization process in Section 2, the active level assignment is not an equal distribution among users. By contrast, the user with better channel conditions (e.g., large SNR and small INR) might be assigned more active levels to maximize the sum rate. Hence, when the number of users is not too much, the optimized level assignment scheme outperforms the scheme that each user exploits one half of the channel resource. In interference channels with more than three users, the virtual IA scheme will gradually dominate the performance due to the advantage that each user can exploit 1/2 DoF of the channel. No matter how many users exist, the achieved sum rate of the orthogonal transmission scheme always keeps unchanged.
The comparisons of the layered IA scheme and other two schemes in Gaussian channels are shown in Fig. 5, where the SNRs and INRs are all randomly selected from 0∼40 dB. Because of the carry-over effect, in moderate SNR levels, the achieved sum rate of the layered IA scheme grows not as fast as in deterministic channels. Yet with two and three users, the layered IA scheme still outperforms the virtual IA scheme. Similarly, with more users the virtual IA scheme will dominate the performance again since its achieved sum-rate keeps linearly increased.
We then investigate the performance of the layered IA scheme in cellular systems. Two network topologies, i.e., homogeneous networks and heterogeneous networks, are respectively studied. In homogeneous networks, every cell has similar coverage and the user is probably interfered by all other adjacent BSs. In heterogeneous networks, the macro-BS may interfere with a lot of pico-cell users, but the pico-BS only interferes with a few macro-cell users. Both kinds of cellular deployment are shown in Fig. 6.  where D is the distance between the BS and user. The celledge SNR is set as 5 dB, and small-scale Rayleigh fading is also considered.
We first observe the performance of layered IA scheme in three-cell homogeneous networks. Although the proposed scheme can be applied with any number of cells, for the homogeneous deployment, the interferences to one cell are mainly from the adjacent two cells. The outer cells with further distance only contribute noisy interference and will not affect the layer participation and alignment scheme. In this configuration, two kinds of user distribution are simulated. The first is symmetric distribution, where the macro-users are moved along the line connecting the macro-BS and the central vertex of the three cells, as the asterisks ("*") shown in Fig. 6, and their distances to macro-BSs are the same. The second is random distribution, where the macro-users are randomly distributed in the cell-edge areas, as the dots ("·") shown in Fig. 6, and the cell-edge boundaries are changed symmetrically. Conventionally, a FFR scheme is used to avoid inter-cell interference in the cell edge areas, where each cell uses 1/3 of the available bandwidth. In the layered IA scheme, we use frequency reuse factor of one, i.e., every cell uses all the bandwidth, and the interference are coordinated in the amplitude space. The sum rate performances are shown in Fig. 7. We can see that the layered IA scheme has approximately two times of the data rate over the FFR scheme both for the symmetric and random user distribution scenarios. The performance gain increases as the users approach the BSs. Although with small-scale fading, when the users are close to the BSs, the average SNRs are higher and the average INRs are lower; thus, there are more opportunities to assign active layers. The random distribution scenario is more practical; when the cell-edge boundary moves far away from the BS, the average sum Compare the average sum rates of the layered IA scheme and the FFR scheme in three-cell homogeneous networks rate achieved in this scenario can be greater than that in symmetric distribution scenario. The reason is that random user distribution diverges the SNR and INR values, which also creates more opportunities to assign active layers.
The performance of the layered IA scheme in heterogeneous networks is evaluated in Fig. 8, where 1∼3 picocells coexisting with one macro-cell are considered. In this kind of deployment, the macro-user is randomly distributed in one sector of the macro-cell, as shown in Fig. 6, and there is also one pico-user randomly distributed in each pico-cell. The layered IA scheme demonstrates great potential in this scenario, and the performance gain keeps increasing along with more coexisting pico-cells. When the distance between the pico-BS and macro-BS changes,  Fig. 8 Compare the average sum rates of the layered IA scheme in heterogeneous networks with different number of pico-cells the sum rate varies in different patterns with different number of pico-cells. With only one pico-cell coexists with the macro-cell, it constitutes a two-user interference channel. Although the position of the macro-user is random, the interference scenarios that the pico-user experienced have certain rules to follow. When the picocell moves from cell center to cell edge, the interference caused by the macro-BS to the pico-user will change from strong to weak. In this kind of two-user interference channel, the sum rate is higher both in strong interference or weak interference scenarios and is lower when the interference has a similar strength with the signal. From the explanation in Section 3, we also know that the strong interference can share the upper amplitude space above the signal layer, and the weak interference can share the lower amplitude space below the signal layer. If the interference has a similar strength as the signal, both the signal and the interference will reduce their occupied amplitude space for coexisting. In the simulation results of Fig. 8, the average sum rate for one pico-cell configuration is a concave curve, which subsidizes our analysis.
When the number of pico-cells increases to two and three, the sum rate curves are monotonically increasing. The reason behind this phenomenon is the interference among the pico-cells. In the simulations, the SNRs and INRs of all links are computed from the transmit power and the simulated channel gains, which include the largescale path loss and small-scale Rayleigh fading. When the pico-cells are close to the macro-BS, as can be seen from Fig. 6, the distances between these pico-cells are shorter. In this scenario, not only the interference from macro-BS to pico-users is stronger but also the interference from pico-BSs to pico-users. When the pico-cells move away from the macro-BS, the interference from macro-BS and pico-BS are all weaker. As observed in Section 3, the layer partitioning and assignment are complicated in the multiuser case. Although we do not provide an explicit analysis, the simulation results show that the average sum rate will increase when we move the pico-cells towards the macro-cell edge.

Conclusions
In this paper, we proposed a layered interference alignment scheme for Gaussian interference networks with arbitrary channel coefficients and applied the idea of amplitude-space sharing to homogeneous and heterogeneous cellular networks. We introduced a layer partitioning method and optimized the active layer assignment based on the insights obtained in deterministic channels. The transmission is implemented by multiple-level nested lattice codes, where the encoding method is judiciously designed to guarantee the superimposed interference layers keep aligned at all receivers. Simulation results show that in Gaussian interference channels the achieved sum rate grows with the number of users, and with two or three users, the proposed scheme outperforms the virtual scheme that each user occupies half of the channel resource. In practical cellular systems, the layered IA scheme provided evident rate gain over the orthogonalbased transmission schemes and showed great potential to mitigate the complicated co-tier and cross-tier interferences in heterogeneous networks.