Delay Analysis in Full-Duplex Heterogeneous Cellular Networks
Leila Marandi Mansour Naslcheraghi Seyed Ali Ghorashi Mohammad Shikh-Bahaei
cs.NI math.ST stat.AP stat.TH

	./Figures/
	




IEEE Transactions on Vehicular Technology, Vol. XX, No. XX, XXX 2019



Delay Analysis in Full-Duplex Heterogeneous Cellular Networks
     Leila Marandi, Mansour Naslcheraghi, Seyed Ali Ghorashi, Mohammad Shikh-Bahaei
	
	Copyright (c) 2015 IEEE. Personal use of this material is permitted. However, permission to use this material for any other purposes must be obtained from the IEEE by sending a request to pubs-permissions@ieee.org.
	
	L. Marandi and S. A. ghorashi are with Cognitive Telecommunication Research Group, Department of Telecommunications, Faculty of Electrical Engineering, Shahid Beheshti University, G.C., Tehran, Iran (email: le.marandi@gmail.com, aghorashi@sbu.ac.ir), M. Naslcheraghi is with Department of Electrical Engineering, Polytechnique Montreal, Montreal, Quebec, Canada (email: mansour.naslcheraghi@polymtl.ca), and M. Shikh-Bahaei is with Department of Informatics, King's College London, London, UK (email: m.sbahaei@kcl.ac.uk). 

    July 24, 2019
=============================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================



IEEEexample:BSTcontrol 


Heterogeneous networks (HetNets) as  a combination of macro cells and small cells are used to increase the cellular network's capacity, and present a perfect solution for high-speed communications. Increasing area spectrum efficiency and capacity of HetNets largely depends on the high speed of backhaul links. One effective way which is currently utilized in HetNets is the use of full-duplex (FD) technology that potentially doubles the spectral efficiency without the need for additional spectrum. On the other hand, one of the most critical network design requirements is delay, which is a key representation of the quality of service (QoS) in modern cellular networks. In this paper, by utilizing tools from the stochastic geometry, we analyze the local delay for downlink (DL) channel, which is
typically defined as the mean number of required time slots for a
successful communication. Given imperfect self-interference (SI) cancellation in practical FD communications, we utilize duplex mode (half-duplex (HD) or FD) for each user based on the distance from its serving base station (BS). Further, we aim to investigate the energy efficiency (EE) for both duplexing modes, i.e., HD and FD, by considering local delay. We conduct extensive simulations to validate system performance in terms of local delay versus different system key parameters.





	Heterogeneous Networks, Full-Duplex, Half-Duplex, Local Delay, Stochastic Geometry, Energy Efficiency.









§ INTRODUCTION


Heterogeneous networks (HetNets) provide services across both macro and small cells to optimize the mix of capabilities. As mobile broadband demand continues to grow dramatically, significant increases in capacity can be achieved by adding small cells to complement macro networks, forming HetNets. By the emergence of numerous applications in wireless networks, delay is one of the most important design measures of these networks that influences the quality of service (QoS) and accordingly the quality of experience (QoE) which directly affects the system reliability. Therefore, experts urge mobile operators to redesign their networks in more advanced ways to increase network capacity and support the growth of network traffic. Many technologies are introduced and developed in the fifth generation (5G) networks such as full-duplex (FD) communications, mmWave, Edge-Computing, Edge-Caching, etc. Meanwhile, HetNets consist of several layers of cells, generally operating at the same bandwidth. Accordingly, low power nodes (LPNs) including micro, pico, and femto BSs are distributed throughout the macro cells. In HetNets, macro cells provide widespread coverage, while LPNs use to reduce dead zones and more importantly to cover high traffic areas <cit.> to improve QoS and overall network capacity. In half-duplex (HD) communications, the mutual connection between pairs is obtained by two methods of Frequency Division Duplexing (FDD) or Time Division Duplexing (TDD). This partitioning of resources in time or frequency domain restricts the use of resources in the network and such methods do not guarantee the 5G requirements <cit.>. One of the ways recently used in HetNets is exploiting FD radios in which nodes are capable to transmit and receive simultaneously in the same frequency and time. Recent developments in dealing with the self-Interference (SI) at FD radios promises practical realization of this technology <cit.>. Recent studies on FD communication <cit.> reveal that the performance of FD outperforms its HD counterpart in terms of delay and sum throughput. Three important criteria for delivery of network information are throughput, delay and reliability<cit.>. Given delay requirements as an important design goal of the HetNets, it is not clear how system performs in terms of delay when the system is FD-capable. Also, since each tier in HetNet has distinctive characteristics, it significantly complicates the analysis as well as simulations. In the downlink, transmission to multiple users can be enabled within the same time-frequency resources <cit.> by applying superposition transmission and power control at the transmitter and successive interference cancellation (SIC) at the receiver. We are also inspired by the challenge of both FD and HetNet techniques that require high power consumption, which is in contrast with the promise of green evolution for future communication systems. Detail of these challenges is still unclear, which is the main focus of this work. Recently, stochastic geometry has been used to study delay in wireless cellular networks. This study is more crucial in two types of delay; the queueing delay and the transmission delay <cit.>. Delay of the queue mainly refers to the waiting time in the service queue while the transmission delay is the time for a successful transmission of data. Recent studies <cit.> introduced the concept of local delay that is a basic form of transmission delay. 

In this paper, to study the delivery of information in HetNet with FD capability, the local delay is considered. Local delay can guarantee reliability and throughput of the system, by acquiring the mean number of time slots for successful transmission. For some network parameters, this local delay will is very high, which means that for a successful transmission, we need a very high number of time slots<cit.>. Therefore, evaluating local delay is significantly important to make sure that users are being served, reasonably.  The follow-up papers utilized this idea to analyze local delay in a clustered network <cit.>. The authors in <cit.> introduced several approaches to match the various traffic arrival in ultra-dense networks (UDN). They utilized tools from the stochastic geometry to study the tradeoff between delay and security of the physical layer in wireless networks. The authors in <cit.> proposed a hybrid-duplex scheme for HetNets by utilizing tools from the stochastic geometry and analyzed signal-to-interference-plus-noise ratio (SINR) distribution, area spectral efficiency (ASE), achievable rate and throughput for both DL and Uplink (UL) channels and have shown that Hybrid-duplex scheme outperforms both pure-HD and pure-FD HetNets. In another research direction <cit.>, the investigations of energy efficiency (EE) and local delay are conducted in conventional HetNets by considering HD mode for all links. However, these methods are inefficient for FD/Hybrid systems due to the SI and suffered interference from the FD users on the receiver of interest.






 §.§ Main Contributions and Outcomes

Thus far, the existing literature have not investigated the concept of local delay in cellular HetNets with hybrid HD/FD communications. Given imperfect SIC in practical FD transceivers and possible aggregate interference from the other transmitters at the receiver of interest, it is unlikely to choose FD communication as the optimal mode for all nodes. Hence, depending on the random circumstances (e.g. distance, channel, etc.) in a cellular HetNet, one duplexing mode will be optimal. In this paper, we analyze local delay by exploiting distance-based mode selection scheme to choose optimal mode for network nodes. Imperfect SI cancellation in practical FD-enabled systems and massive deployment of small cells in HetNets not only increases the expected interference on the receiver of interest, but also causes more energy consumption. Thereby, we investigate energy efficiency (EE) in the FD-enabled HetNets which is eventually more than its HD counterpart. Thus far, from the local delay point of view, to the best of our knowledge, there is no existing work to investigate EE in the FD-Enabled HetNets. Given FD capability for all nodes in the network, we aim to propose a hybrid duplex mode selection for the HetNets. Further details on the contributions are provided in the sequel.


	
  * We analyze the local delay for HetNets with Hybrid duplexing (HD/FD) and evaluate the impact of FD on the local delay by utilizing tools from the stochastic geometry. In particular, we aim to model and analyze the expected time required for a successful transmission from an arbitrary node to its respective receiver in terms of system key parameters. 
	
	
	
  * We exploit a hybrid duplex mode selection  scheme for FD-enabled HetNets in which all BSs and user devices are FD-capable. The reason behind choosing the distance-based duplex mode selection policy is that, considering imperfect self-interference cancellation (SIC), the SIR of a full-duplex user may be worse than the half-duplex one in the same place. Although the small-scale fading affects the SIR,  it could be averaged out and the path-loss would be the main factor influencing the SIR in long term. Hence, it may not be suitable for the cell edge users to establish a full-duplex communication. Based on that and similar to <cit.>, distance from the serving BS is chosen as a duplex mode selection criterion. 
	 This scheme is based on user distance from its respective serving BS. Different from <cit.>, we consider silent mode capability for all BSs to reduce the energy consumption. This capability allows BSs to be inactive and this results in reduction of the experienced interference on the receiver of interest and correspondingly, will lead to increase the chance of successful transmissions. 
	
	
	
  * Different from <cit.> in which the EE is investigated in single layer cellular networks (none-HetNet), we investigate the EE in the FD-capable HetNets by considering local delay in the analysis and utilizing tools from the stochastic geometry. In particular, we have modeled EE with respect to the behavior of local delay.

The rest of this paper is organized as follows. In section <ref>, we delineate the system model and introduce the key parameters. In section <ref>, we model and analyze the local delay. Theoretical and simulation results are provided in section <ref> and finally section <ref> concludes the paper.




§ SYSTEM MODEL
 



 §.§ Network Model
 
We consider a HetNet with K- tier of BSs that are different in transmitting powers, densities and path loss exponents. The BS locations of tier 𝒦 = {1,2, ⋯ ,K} are distributed in the Euclidean plane ℝ^2 independently based on the homogeneous Poisson point process (HPPP) denoted by Φ _k with density of λ _k. The BSs of each tier have the same transmit power, which are denoted by p_k. Due to high density of BSs in HetNets, the background noise is negligible in comparison with the amount of interference. Denoting τ_k as the system target signal-to-interference ratio (SIR) threshold for a typical user located in tier k and is operating in mode † (FD/HD), the probability of the event that this typical user satisfies system threshold τ_k is defined as 𝒫( _k^†≥τ _k). If the SIR thresholds of all tiers are the same, then we have τ_k = τ. In general, the properties of each tier can be expressed as {λ _k,p_k,τ _k}. Users are also independently distributed in the Euclidean plan based on HPPP Φ _u with density of λ _u. The transmit power of users is denoted by p_u. Both BSs and users are FD-capable. Imperfect SI cancellation is assumed with the residual SI signal power of βp_k, βp_u, for BSs and users, respectively, where 0 ≤β≤ 1 represents the residual power ratio. If β = 0, the SI cancellation is perfect and in case of β = 1, there is no SI cancellation. The network is considered to be static, i.e., the positions of BSs and users remain the same throughout time. In practical systems, the delay performance highly relies on the properties of the traffic in wireless networks. In this paper, to consider the traffic modeling, the scenario of the network is considered as backlogged <cit.>, i.e., when the transmitters are scheduled to communicate, they always have data to transmit. Furthermore, we assume the system to be open access, where the user can connect to any BS. Inspiring from the Slivniak's theorem <cit.>, we conduct the analysis on a typical user in the origin for DL channel, for this purpose, which implies that any user has the same chance of being at the origin in process. The notations of the paper are summarized in table <ref>.





 §.§ Channel Model and Duplex mode selection
 
In the proposed model, each BS and user have a single antenna for power transmission/reception. For the propagation model, we assume a combination of fading and path loss. It is assumed that all channels are Rayleigh fading with unit mean. Let us denote the power fading coefficient between BS b_k and a typical user u_o at given time slot t with h_t,b_ku_o∼exp( 1 ). We assume standard power-law path loss model l( r ) =  r ^ - α _k, where the BS of tier k is located at distance r from the origin (typical user) and α  > 2 is the path loss exponent. We consider cell association based on maximum average DL received power (DRP), where a mobile user is associated with the strongest BS in terms of long-term averaged DRP at the user. The tier index that a typical user is associated to, is given by

    k = max_i {p_ir_i^ - α _i,i ∈𝒦}.


By exploiting distance-based duplex mode selection, we denote ϑ _k as the distance threshold; if the distance between typical user and its serving BS is less than ϑ _k, the user will be treated as an FD user, otherwise, it will operate in HD mode. FD and HD user sets in tier k are denoted by  𝒰_k_ FD, 𝒰_k_ HD, respectively. The corresponding definitions can be written as

    𝒰_k_ FD = {( p_kr_k^ - α _k > p_jr_j^ - α _j,∀ j  k) ∧( r_k⩽ϑ _k)},


    𝒰_k_ HD = {( p_kr_k^ - α _k > p_jr_j^ - α _j,∀ j  k) ∧( r_k > ϑ _k)}.


Time is divided into discrete slots with the same duration δ _t and each transmission spans one time slot. BSs can transmit/receive data only if they are in the active mode. Technically, each BS of the k^th tier at any time slot with probability χ__k( 0 ≤χ__k≤ 1) is silent. The interference at typical user in time slot t consists of interference from all active BSs, except the serving BS, as well as the set of FD users, which can be expressed as

    ℐ_t = ∑_j ∈𝒦[   ∑_b_j∈Φ _ja\{b_o}p_jh_t,b_ju_or_b_ju_o^ - α _j + 
       ∑_u_i∈Φ _u^ FD\{u_o}p_uh_t,u_iu_or_u_iu_o^ - α _u],
 
where Φ _ja = Φ _ja^ FD∪Φ _ja^ HD indicate the set of active BSs that consists of active FD ( Φ _ja^ FD) and HD ( Φ _ja^ HD) BSs in tier j in time slot t. The interference caused by the active BSs can be written separately; ℐ_t,B = ∑_b_j∈Φ __ja^ FD\{b_o}p_jh_t,b_ju_or_b_ju_o^ - α _j  + ∑_b_j∈Φ __ja^ HD\{b_o}p_jh_t,b_ju_or_b_ju_o^ - α _j. Φ _u^ FD is the set of FD users. b_o represents serving BS of typical user u_o. h_t,b_ju_o indicates the fading power coefficients between BS b_j and typical user u_o, similarly h_t,u_iu_o is fading from FD user u_i and typical user u_o. The SIR of the typical user in time slot t in the k^th tier is 

    _k,t^† = p_kh_t,b_ou_or_b_ou_o^ - α _k/ℐ_t + 1_†. βp_u_o,




where †∈{ HD, FD  }. _k,t^ FD, _k,t^ HD indicate the SIR, when typical user is FD and HD user, respectively. βp_u_o is residual SI signal at the typical user. 1_† is the indicator function and defined by

    1_† = {[   1 , † = ,;    0, † = . ].


An illustration of the system described above is delineated for a 3-tier HetNet in Fig.<ref> along with the expected interferences in DL channel. Tier 1, i.e. Macro has more transmission power than Pico (i.e. tier 2) and Femto (i.e. tier 3) cells. Shown HD/FD modes for the users are determined based on distance as described before.
 



black







§ LOCAL DELAY ANALYSIS
 

First we derive general local delay on DL channel for a typical operating in FD and HD modes, then we calculate EE based on the local delay analysis. 



 §.§ Local Delay Definition

	Let r_k be the distance between the typical user and its associated BS. Assuming the typical user at the origin, we assume that 𝒫_^† is the probability that the typical user (operating in mode †) is successfully connected to its tagged BS in a single transmission conditioned on r_k. Then, the conditional success probability is given by


    𝒫_ suc^† = ( 1 - χ__k)𝒫( _k,t^† > τ_k | r_k.).

where χ_k is the silence probability and τ_k is the given threshold defined in section <ref>. Successful transmission occurs only when the BS is active and SIR satisfies the system threshold value. The purpose of operating in silent mode is energy saving. However, if the associated user is going to transmit data and BS is in the silent mode, it should be noted that it is only for one time slot and not for the entire time duration, therefore, it can try to transmit in the next time slots. Conditioned on r_k, by considering the constant transmit power (or independent and identically distributed (i.i.d) transmit power) and i.i.d fading, the event of successful transmission in time slots and the success indicator random variables are temporally i.i.d <cit.>. Hence, the distribution of the conditional local delay is geometric with mean (𝒫_ suc)^-1. The local delay of a typical user <cit.> in tier k operating in mode † is then the expectation with respect to (w.r.t.) r_k:  

    𝒟_k^†  = 𝔼_r_k[ 1/𝒫_suc^†].

In Eqn. <ref>, 𝒟_k^† denotes the average number of slots that takes the tagged BS to successfully transmit a packet to the typical user.



Note that the typical user associates to one tier. Now, according to the law of total probability, the local delay can be written as

    𝒟_† = ∑_k ∈𝒦𝒜_k^†𝒟_k^†,


where 𝒜_k^ FD and 𝒜_k^ HD provide the per tier association probability that typical user operates in FD and HD modes, respectively. 

 
	Generally, the association probability in the k^th tier for a typical user that operates in HD or FD mode is given by

    1!𝒜_k^ FD = 2πλ _k_0^ϑ _kxexp( ∑_i  k - πλ _i( p_i/p_kx^α _k)^2/α _i - πλ _kx^2) dx.
		
	and
		
    1!𝒜_k^ HD = 2πλ _k_ϑ _k^∞xexp( ∑_i  k - πλ _i( p_i/p_kx^α _k)^2/α _i - πλ _kx^2) dx.

 The probability that the user is associated to the tier k is defined as the sum of the association probabilities in FD and HD modes, namely 
	
    1!𝒜_k = 2πλ _k_0^∞xexp( ∑_i  k - πλ _i( p_i/p_kx^α _k)^2/α _i - πλ _kx^2) dx.

	
		 
		see appendix <ref>.
	





 §.§ General Case of Local Delay
 
 
	The local delay in a K-tier HetNet with Hybrid duplex scheme and associated cell based on DRP is given by




    𝒟_ FD =    ∑_k ∈𝒦2πλ _k/( 1 - χ__k)∫_0^ϑ _kexp( sβp_u_o) F_1( r,λ_j )
       ∏_j ∈𝒦 F_2( r,λ_j ,τ_k ) F_3( r,λ_u ,τ_k ) dr,

	


    𝒟_ HD =    ∑_k ∈𝒦2πλ _k/( 1 - χ__k)∫_ϑ _k^∞ F_1( r,λ_j ) ∏_j ∈𝒦 F_2( r,λ_j ,τ_k )
        F_3( r,λ_u ,τ_k ) dr,


where s = τ _kr_k^α _k/p_k,

    F_1( r,λ_j ) = r.exp(  - π∑_j ∈𝒦λ _j( p_j/p_kr^α _k)^2/α _j),


    F_2( r,λ_j ,τ_k ) = exp(    π( 1 - χ__j)λ __j( p_j/p_kr^α _k)^2/α _j
       ∫_1^∞( τ _k/τ _k + u^α _j/2)du),


    F_3( r,λ_u ,τ_k ) = exp( 2πλ _j,u^ FD∫_e_k( j )( 1 - 1/1 + sp_uy^ - α _u)ydy).





	We provide the proof for typical user in FD mode and similar to the same trend, it is valid for the typical user in HD mode. First, we need to get the probability of successful transmission by using Eq.(<ref>). Then we put the final result of Eq.(<ref>) in Eq.(<ref>) to obtain expression for the local delay when typical user operates in FD mode.
	
	lLl
		
		
		   𝒫_suc^FD = ( 1 - χ__k )𝒫(  _k,t^FD > τ_k| r_k . )

		   = ( 1 - χ__k )𝒫( p_kh_t,b_ou_o r_b_ou_o ^ - α_k > τ_k( ℐ_t + βp_u_o )| r_k . )

		   = ( 1 - χ__k )𝒫( h_t,b_ou_o > τ_k r_b_ou_o ^α_k/p_k( ℐ_t + βp_u_o )| r_k . )

		   = ^a ( 1 - χ__k )exp(  - sβp_u_o )𝔼_ℐ_t[ exp(  - sℐ_t )| r_k . ]

		   = ( 1 - χ__k )exp(  - sβp_u_o )L_ℐ_t( s| r_k . ),
		
	

where (a) follows from the fact that h_t,b_ou_o∼exp( 1 ) and L_ℐ_t( s| r_k.) refers to the Laplace transform of interference and can be calculated as shown in Eq. (<ref>). where  ℐ_t,j = ℐ_t,B + ∑_u_i∈Φ _u^FD\{u_o}p_uh_t,u_iu_or_u_iu_o^ - α _u and for 𝔼_Φ ,h[ exp(  - s∑_j ∈𝒦ℐ_t,j)], we get Eq. (<ref>) where (a) is obtained by using probability generating functional (PGFL) of the Poisson point process (PPP). Limits of integrals in Eq.(<ref>)  are from e_k( j ) tend to ∞ and therefore, the closest interferer from tier j is located at least at the distance e_k( j ) = ( p_j/p_kr_k^α _k)^1/α _j. On the other hand, it is difficult to specify the exact range of interference from FD users to a typical user. Interfering users are distributed around their serving BSs, so it can be approximated by distance between interfering BSs and typical user, hence the limits of integrals are also e_k( j ). Since each BS specifies the duplex mode and how it connects to its user, the density of FD users is defined as λ _j,u^ FD = 𝒜_j^ FD/𝒜_j( 1 - χ__j)λ _u, where 𝒜_j^ FD is the association probability that the user operates in FD mode within the j^th tier. (b) is obtained by variables substitutions u = ( p_j/p_kr_k^α _k)^ - 2/α _jy^2 and finally, putting the Eq.(<ref>) in Eq.(<ref>) gives the final result for 𝒟_k^ FD. According to <cit.>, distance distribution of a typical user that operates in FD mode from its serving BS is given by








    f_r_k^ FD( r )= 
    
    0   r > ϑ _k,
    1/A_k^ FD2πλ _kr.exp(  - π∑_j ∈𝒦λ _jξ _k^2( r))   r ≤ϑ _k.


where ξ _k( r) = ( r^α _kp_j/p_k)^1/α _j. Hence, the local delay for a typical user associated to the k^th tier and operates in FD mode is given by  

    𝒟_k^ FD = ∫_0^ϑ _k𝒟_k^ FD( r ) .f_r_k^ FD( r )dr,




    𝒟_k^ FD = ∫_0^ϑ _k   2πλ _k/( 1 - χ__k)𝒜_k^ FDexp( sβp_u_o) F_1( r,λ_j )
       ∏_j ∈𝒦 F_2( r,λ_j ,τ_k ) F_3( r,λ_u ,τ_k ) dr.


Finally, by using Eq.(<ref>) in Eq.(<ref>), we get Eq.(<ref>). The proof of Eq.(<ref>) is similar to the same process and we avoid the repetition due to lack of space. 






  §.§.§ Special case



By considering α _k = α _j = α _u = α in Eq.(<ref>), Eq.(<ref>), Theorem <ref> can be simplified as Eq.(<ref>), Eq.(<ref>) at the page eq: special case Delay-FD.



In Eqn. (<ref>), (a) follows taking integral with respect to variable r, ρ( τ _k,α) = ∫_1^∞( τ _k/τ _k + u^α/2)du, and ρ '( τ _k,α) = ∫_1^∞( τ _kp_u/τ _kp_u + u^α/2p_j)du.







 §.§ Energy Efficiency (EE)
 

	
	For a successful communication between the typical user and its associated BS, EE is defined as the ratio of network throughput (or area spectral efficiency) to the average total power consumption <cit.>.	

 The achieved throughput in a single HD/FD link is denoted by 𝒟_†^ - 1log( 1 + τ) nats/s/Hz. The sum throughput of the hybrid HD/FD network is given by 𝒯_† = D_† ^ - 1log( 1 + τ)∑_k ∈𝒦( 1 - χ__k)λ _k^†, where λ _k^†  = 𝒜_k^†/𝒜_k. In this paper, we obtain the power consumption of each BS in each tier based on a linear approximation of the BS power modeling <cit.>. The power consumption of a BS in tier k denoted by P_k and is defined by
 
 
    P_k = 
     P_k0 + Δ _kp_k + 1_†.  P_SIC,   p_k > 0,
    P_k,sleep,   p_k = 0

 
 where the above expression indicates the power expenditure in active mode, which P_k0, Δ _k is the static power and the slope of power consumption in the k^th tier, respectively. P_SIC represents power consumption for SI cancellation. 1_† is the indicator function defined in Eq. (<ref>). P_k,sleep is the static power consumption in sleep mode <cit.>. Conclusively, we define the EE as
 
 
    η  = D_† ^ - 1log( 1 + τ)∑_k ∈𝒦( 1 - χ__k)λ _k^†/∑_k ∈𝒦λ _k^†[ ( 1 - χ__k)( P_k0 + Δ _kp_k + 1_†.P_SIC) + χ__kP_k,sleep],

 where the denominator represents power consumption of all BSs across the entire network. The unit of EE is nats/Joule/Hz. By substituting 𝒟_ FD and 𝒟_ HD from Eq.(<ref>) and (<ref>) in (<ref>), the general expression of  η can be obtained. Since the expression of EE in Eq. (<ref>) is not in closed form, we will validate it numerically in the next section. 
 
 
 

§ NUMERICAL RESULTS
 
 We consider a two-tier (macro and poico BSs) cellular HetNet (K=2). The transmit power for the BSs in first and second tiers assumed to be 46 dBm, 30 dBm, respectively. Transmit power for the user devices is 23 dBm.  ϑ _1 = 300 m and ϑ _2 = 150 m. SIC factor is β = -70 dB and power expenditures are P_10 = 139 Watts, P_20 = 9.7 Watts, Δ _1 = 5, Δ _2 = 5.7, P_SIC = 50 mWatts, P_1,sleep = 80 Watts and P_2,sleep = 6.1 Watts. Note that in these parameters, "1", "2" refer to the first and second tier, respectively. Fig. <ref> shows the local delay for a typical user FD and HD in the hybrid network and EE versus SIR threshold τ for different silent probabilities. In Fig. <ref> we can observe that the local delay for FD mode in all silence probabilities is less than the HD one. This is because the FD users are closer to their tagged BSs which leads to much stronger received power by the FD transceivers in comparison with the HD one. However, with the rising of SIR threshold, the local delay for both duplex modes is increasing. We find that in the plane of τ  < 0, by increasing the silent probability, local delay in both HD and FD communications increases. The reason is, at the low regime of τ, most of the transmissions are successful and consequently the interference from BSs on the typical user has slight impact. Conversely, in the high regime of τ (i.e., τ > 0), the impact of aggregated interference is more accented. However, by increasing the silent probability, the number of active BSs decreases, which will lead to partial reduction in the suffered interference at the typical user and this will result in local delay reduction. As can be seen from Fig.<ref>, EE is being reduced by increasing the silent probability. It is clear that by increasing the silent probability, the number of active BSs decreases, namely the number of active links between user devices and BSs will be reduced. This causes reduction in system sum throughput, while the BSs in silent mode use energy, correspondingly, the EE will be decreased according to its definition in <ref>. It can be seen that there is an optimal SIR threshold (around 0 dB) that maximizes EE. In Fig. <ref>, local delay for a typical user FD and HD in the hybrid network and EE are illustrated for a three-tier HetNet, in which the transmission power for the third tier assumed to be 23 dBm, ϑ_3 = 75 m, and power expenditures for the third tier are P_30 = 7.7 Watts, Δ _3 = 7, and P_3,sleep = 4.8 Watts. By increasing the number of tiers, due to the reduction in distance threshold, the chance of making FD communications, decreases. This is due to the fact that in the third tier, the distance threshold is much smaller than the distance thresholds in tiers one and two. Hence, the existing HD users in tire three experience lower delay due to the reduction in the interference incurred from the FD users, therefore they can receive higher power levels. Correspondingly, this will incur three-tier HetNet a lower delay, higher throughput, and consequently a better energy efficiency in comparison to two-tier HetNet. The explanation of EE in Fig. <ref> is also valid for Fig. <ref>. The effect of silent probability in three-tier HetNet is similar to the results of two-tier network demonstrated in Fig. <ref>. In Fig.<ref>, by increasing the density of tier 2, the probability that more users can operate in FD mode is raising. On the other hand, interference from FD users on a typical user increases and due to distance-based duplex selection scheme, the FD users are getting closer to each other. This phenomena causes higher aggregate interferences on the typical user and consequently will increase the chance of link failure. As a result, the local delay of the typical user that operates in FD mode will be increased. Conversely, for the HD counterpart, by increasing the density of users across the entire network, fewer users will operate in HD mode (i.e., due to distance-based duplex selection scheme), and therefore these users will be located at the far distant away. Correspondingly, the impact of aggregated interference from FD users on the typical user that operates in HD mode will be decreased. Hence, the local delay will decrease or not change with higher density. Conceding to Fig. <ref>, in the higher SIR threshold regime (τ≥0), increasing density causes more power consumption. On the other hand, due to non-satisfaction of the SIR threshold, we will have increment in local delay and reduction in EE. However, in the lower SIR thresholds (τ≤ 0), the impact of short distance association dominates against the aggregated interference, and hence the acquired SIR is improved and EE is enhanced. Fig. 5a demonstrates the local delay as a function of the silent probability of BSs of tier 2 (i.e., χ_2). By increasing the chance of being silent, the probability of operating in FD mode and consequently the number of FD transceivers decreases due to the reduction in the number of BSs in tier 2. Hence, the aggregated interference imposed on typical HD user decreases and as a result, the local delay decreases. On the other hand, reduction in the number of active BSs causes reduction in the success probability for the users, which in turn leads to higher local delay. There is an optimal silence probability to minimize local delay when a typical user operates in HD mode. In case of the impact of density, similar explanations for Fig. <ref> are also valid for Fig. <ref>. For the EE as shown in Fig.<ref>, with respect to Fig. <ref>, by increasing density of BSs, the EE decreased at 0 dB of τ. It is due to the increasing aggregated interference, high energy consumption, non-satisfaction of SIR threshold and reduction in system sum throughput. As can be seen, the curves overlap with increasing the chance of being silence in tier 2 due to reduction in the effective interference from tier 2. As a result, the throughput and consequently EE is increased. The local delay for a typical user FD in the hybrid network of DL channel with the different SIC factors is shown in Fig.<ref>. As the value of SIC increases, local delay decreases. Obviously, this is due to the better SIC. In Fig.<ref> the local delay in DL of macro and small cell users is defined as follows
  
    𝒟_k = 1/𝒜_k( 𝒜_k^ FD𝒟_k^ FD + 𝒜_k^ HD𝒟_k^ HD).

 For -10 ≤τ≤ -2, there is no significant difference in local delay between macro and small cells. However, for the higher SIR thresholds, the local delay in macro cells is less than small cells, and this is because of transmitting power of macro cell which is more than small cells (i.e. Femto and Pico cells) and eventually satisfies the SIR threshold.
 
 
 



 
 
 
 
  

  

 
 

§ CONCLUSION
  
 We analyzed the local delay of K-tier HetNet with hybrid HD/FD communications. By using tools from the stochastic geometry, we derived the local delay expression for both HD and FD modes. Analytical and simulation results gave some insights on the impact of key system parameters, including the SIR threshold, BS density and silent probability. Higher SIR threshold increases the local delay. However, from the energy efficiency (EE) point of view, we have optimal limit that increases the network EE. Higher density of the BSs in the high SIR regime causes higher local delay for FD users as well as inefficiency in the energy consumption. Rising the silent probability in the high SIR regime can improve local delay and EE. We observed that Hybrid HetNets can increase the optimal SIR range for EE. Exploring delay-constrained mode selection will be a useful extension to this work. Moreover, in this paper and all related works, the channel for the FD links is assumed to be reciprocal and exchanging data on FD nodes is assumed to be symmetric. FD communications with asymmetric data and a dissimilar channel is another worthwhile open problem to pursue.
 
 

 
 
 

 §.§ Proof for Lemma <ref>

 
The event of r_k≤ϑ _k for FD users, the association probability with the kth tier can be given as  

	
    𝒜_k^ FD  =  𝔼_r_k[∏_i  k𝒫( r_k < ( p_k/p_ir_i^α _i)^1/α _k) ∧( r_k≤ϑ _k) ] 
        = ∫_0^ϑ _k∏_i  k𝒫(r_i > ( p_i/p_kx^α _k)^1/α _i) .f_r_k( x )dx
         = ^( a ) 2πλ _k∫_0^ϑ _kx.exp(  - π∑_i  kλ _i( p_i/p_kx^α _k)^2/α _i - πλ _kx^2)dx

	 
 
 where (a) derived by using null probability and distance distribution: 𝒫(r_i > ( p_i/p_kx^α _k)^1/α _i)= exp(  - πλ _i( p_i/p_kx^α _k)^2/α _i) and f_r_k( x ) = 2πλ _kxexp(  - πλ _kx^2). By combining both equations, we obtain Eqn. (<ref>). The proof for HD mode is the same.  





IEEEtran


10
	HetNet2
	A. Ghosh, N. Mangalvedhe, R. Ratasuk and et al., “Heterogeneous cellular networks: From theory to practice," IEEE Commun. Mag., vol. 50, no. 6, pp. 54-64, Jun. 2012.
	
	HD
	A. Gohil, H. Modi and S. K. Patel, “5G technology of mobile communication: A survey," in Proc., ISSP., Mar. 2013, pp. 288-292.
	
	FD1
	S. Goyal, P. Liu, S. S. Panwar and et al., “Full duplex cellular systems: will doubling interference prevent doubling capacity?," IEEE Commun. Mag., vol. 53, no. 5, pp. 121-127, May. 2015.
	
	FD2
	G. Noh, H. Wang, C. Shin and et al., “Enabling technologies toward fully LTE-compatible full-duplex radio," IEEE Commun. Mag., vol. 55, no. 3, pp. 188-195, Mar. 2017.
	
	FD-app1
	M. Naslcheraghi, S. A. Ghorashi and M. Shikh-Bahaei, “FD device-to-device communication for wireless video distribution," IET Commun., vol. 11, no. 7, pp. 1074-1081,  2017.
	
	FD-app2
	M. Naslcheraghi, S. A. Ghorashi and M. Shikh-Bahaei, “Full-duplex device-to-device collaboration for low-latency wireless video distribution," in Proc., ICT, May. 2017, pp. 1-5.
	
	FD-app3
	M. Naslcheraghi, M. Afshang and H. S. Dhillon, “Modeling and performance analysis of full-Duplex communications in cache-enabled D2D Networks," in Proc., IEEE  ICC, May. 2018, pp. 1-6.
	
	FDSurvey
	A. H. Gazestani, S. A. Ghorashi, B. Mousavinasab and M. Shikh-Bahaei, “A survey on implementation and applications of full duplex wireless communications," Physical Commun., vol. 34, pp. 121-134, Jun. 2019.
	
	LocalDelay1
	M. Haenggi, “The local delay in Poisson networks," IEEE Trans. Inf. Theory, vol. 59, no. 3, pp. 1788-1802, Mar. 2013.
	
	EditorPaper
	L. Dai, B. Wang, Z. Ding and et al., “A survey of non-orthogonal multiple access for 5G," IEEE Commun. Surveys Tuts., vol. 20, no. 3, pp. 2294-2323, May. 2018.
	
	Delay-cause
	Y. Zhong, M. Haenggi, F. Zheng and et al., “Toward a tractable delay analysis in ultra-dense networks," IEEE Commun. Mag., vol. 55, no. 12, pp. 103-109, Dec. 2017.
	
	LocalDelay2
	M. Haenggi, “Local delay in Poisson networks with and without interference," in Proc., Allerton, Sep. 2010, pp. 1482-1487.
	
	LocalDelay3
	Z. Gong and M. Haenggi, “The local delay in Poisson networks," IEEE Trans. Wireless Commun., vol. 12, no. 9, pp. 4766-4777, Sep. 2013.
	
	LocalDelay4
	G. Alfano, R. Tresch and M. Guillaud, “Spatial diversity impact on the local delay of homogeneous and clustered wireless networks," in Proc., WSA, Feb. 2011, pp. 1-6.
	
	Reviewers
	Y. Zhong, X. Ge, H. H. Yang and et al., “Traffic matching in 5G ultra-dense networks," IEEE  Commun. Mag., vol. 56, no. 8, pp. 100-105, Aug. 2018.
	
	Reviewers2
	Y. Zhong, X. Ge, T. Han and  et al., “Tradeoff between delay and physical layer security in wireless networks," IEEE J. Sel. Areas Commun. , vol. 36, no. 7, pp. 1635-1647, Jul. 2018.
	
	Hybrid-dulex1
	J. Lee, T. Q. S. Quek, “Hybrid full-/half-duplex system analysis in heterogeneous wireless networks," IEEE Trans. Wireless Commun., vol. 14, no. 5, pp. 2883-2895, May. 2015.
	
	Hybrid-duplex2
	W. Tang and et al., “Distance-based hybrid duplex in heterogeneous networks," in Proc., IEEE GLOBECOM, Dec. 2015, pp. 1-6.
	
	Hybrid-duplex3
	W. Tang, S. Feng, Y. Liu and Y. Ding, “Hybrid duplex switching in heterogeneous networks," IEEE Trans. Wireless Commun., vol. 15, no. 11, pp. 7419-7431, Nov. 2016.
	
	LocalDelay5
	L. Liu, Y. Zhong, W. Zhang and M. Haenggi, “On the impact of coordination on local delay and energy efficiency in Poisson networks," IEEE Wireless Commun. Lett., vol. 4, no. 3, pp. 241-244, Jun. 2015.
	
	DTX-LocalDelay
	W. Nie, Y. Zhong, F. Zheng and et al., “HetNets with random DTX scheme: Local delay and energy efficiency," IEEE Trans. Veh. Technol., vol. 65, no. 8, pp. 6601-6613, Aug. 2016.
	
	EE-HD
	D. W. K. Ng and et al, “Energy-efficient resource allocation in OFDMA systems with large numbers of base station antennas," IEEE Trans. Wireless Commun., vol. 11, no. 9, pp. 3292-3304, Sep. 2012.
	
	distance-criteria
	K. T. Hemachandra and A. O. Fapojuwo, “Distance based duplex mode selection in large scale peer-to-peer wireless networks," in Proc., IEEE ICC, May. 2017, pp. 1-6.
	
	EE-Hybrid
	Z. Wei, S. Sun, X. Zhu and et al., “Energy-efficient hybrid duplexing strategy for bidirectional distributed antenna systems," IEEE Trans. Veh. Technol., vol. 67, no. 6, pp. 5096-5110, Jun. 2018.
	
	typicaluser-Palm
	M. Haenggi, Stochastic Geometry for Wireless Networks, 2013.
	
	LocalDelay6
	X. Zhang and M. Haenggi, “Delay-optimal power control policies," IEEE Trans. Wireless Commun., vol. 11, no. 10, pp. 3518-3527, Oct. 2012.
	
	Stochastic-HetNet2
	H. Jo, Y. J. Sang, P. Xia and J. G. 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.
	
	powr-consume
	G. Auer, V. Giannini, C. Desset and et al., “How much energy is needed to run a wireless network?," IEEE Wireless Commun. Mag., vol. 18, no. 5, pp. 40-49, Oct. 2011.






