Enabling Wireless Power Transfer in Cellular Networks: Architecture, Modeling and Deployment
Abstract
Microwave power transfer (MPT) delivers energy wirelessly from stations called power beacons (PBs) to mobile devices by microwave radiation. This provides mobiles practically infinite battery lives and eliminates the need of power cords and chargers. To enable MPT for mobile charging, this paper proposes a new network architecture that overlays an uplink cellular network with randomly deployed PBs for powering mobiles, called a hybrid network. The deployment of the hybrid network under an outage constraint on data links is investigated based on a stochasticgeometry model where singleantenna base stations (BSs) and PBs form independent homogeneous Poisson point processes (PPPs) with densities and , respectively, and singleantenna mobiles are uniformly distributed in Voronoi cells generated by BSs. In this model, mobiles and PBs fix their transmission power at and , respectively; a PB either radiates isotropically, called isotropic MPT, or directs energy towards target mobiles by beamforming, called directed MPT. The model is applied to derive the tradeoffs between the network parameters under the outage constraint. First, consider the deployment of the cellular network. It is proved that the outage constraint is satisfied so long as the product is above a given threshold where is the pathloss exponent. Next, consider the deployment of the hybrid network assuming infinite energy storage at mobiles. It is shown that for isotropic MPT, the product has to be above a given threshold so that PBs are sufficiently dense; for directed MPT, with denoting the array gain should exceed a different threshold to ensure short distances between PBs and their target mobiles. Furthermore, for directed MPT, has to be sufficiently large as otherwise PBs fail to deliver sufficient power to target mobiles regardless of powertransfer distances. Last, similar results are derived for the case of mobiles having small energy storage.
Power transmission, cellular networks, energy harvesting, stochastic processes, adaptive arrays, mobile communication
I Introduction
One of the most desirable new features for mobile devices is wireless charging that eliminates the need of power cords and chargers. To realize this feature, the paper proposes a novel network architecture (see Fig. 1) where stations called power beacons (PBs) are deployed in an existing cellular network for charging mobiles via microwave radiation known as microwave power transfer (MPT). Charging a typical smartphone requires average received microwave power (MP) of tens to hundreds of milliwatts [1]. For instance, a conventional smartphone powered by a battery supplying mAh at V and a battery life of hr can be selfsustaining by WPT if the average transferred power is above mW assuming lossless RFDC conversion. The challenges in realizing wireless charging by MPT are threefold, namely creating lineofsight (LOS) links from PBs to mobiles to enable closetofreespace power transfer, forming sharp energy beams at PBs to counteract propagation loss, and reducing power consumption of mobiles. These challenges can be tackled by the advancements of three technologies, respectively. First, PBs require neither backhaul links nor complex computation, allowing lowcost and dense deployment. Specifically, PBs can be installed wherever there are connections to the electrical grid or even in vehicles and airplanes. Deploying many PBs shortens transfer distances, resulting in LOS powertransfer links. Second, massive antenna arrays with tens to hundreds of elements are a technology currently under active development [2] and can be used at PBs for forming sharp energy beams towards mobiles to achieve the WPT efficiency close to one. As discussed in the sequel, gigantic antenna arrays have been used to achieve an efficiency of for MPT over km (from a solarpower satellite to the earth). Third, cellular networks are evolving towards the deployment of dense base stations that will substantially reduce the required transmission power of mobiles [3]. Consequently, PBs’ transmission power can be reduced or the MPT distances increased. Furthermore, the advancements in the areas of lower power electronics and displays (e.g., reflective displays) will also decrease the power consumption of mobile devices and increase the practicality of MPT.
In this paper, an uplink cellular network overlaid with PBs, called the hybrid network, is studied using a stochasticgeometry model where base stations (BSs) and PBs form independent homogeneous Poisson point processes (PPPs) and mobiles are uniformly distributed in corresponding Voronoi cells with respect to BSs. In the model, PBs are deployed for powering uplink transmissions under an outage constraint by either isotropic radiation or beamforming towards target mobiles, called isotropic or directed MPT, respectively. The model is used to derive the tradeoffs between the network parameters including the PB/mobile transmission power and PB/BS densities for different network configurations accounting for isotropic/directed MPT and mobiles having large/small energy storage. The results provide insight into the hybridnetwork deployment and throughput.
Ia Microwave Power Transfer
Wireless power transfer is usually implemented in practice based on one of the three different technologies, namely inductive coupling [4], magnetic resonance coupling [5], and MPT [6] corresponding to short range (tens of centimeters), mid range (several meters), and long range (up to tens of kilometers), respectively. Besides limited transfer distances, both inductive and magneticresonance coupling require alignment and calibration of coils at transmitters and receivers, making these technologies unsuitable for mobile charging. MPT has no such limitations but suffers from potentially severe propagation loss over a long transfer distance.
Such loss can be overcome by freespace beaming and intercepting microwave energy using antennas with large apertures. Based on this principle, efficient pointtopoint MPT over long distances has been demonstrated in largescale systems. In 1964, W. C. Brown first showed the practicality of MPT by demonstrating a microwavepowered small helicopter hovering at an altitude of feet, where a rectenna was used to intercept power beamed from the ground using the frequency of GHz and supply DC power of W [7]. This experiment motivated researchers in Canada and Japan to develop microwavepowered airplanes in s [8] and s [9], respectively. An extreme MPT system was designed for the National Aeronautics and Space Administration (NASA) of the United States for transferring power of about GW from a solarpower satellite to a station on the surface of the earth [10]. The powertransfer efficiency is predicted to be over a transfer distance as long as km by using gigantic transmit and receive antenna arrays with diameters of m and km, respectively. Large antennas with the dimensions in the aforementioned experiments are impractical and unnecessary for MPT in the hybrid network where efficient MPT relies on deploying dense PBs to shorten powertransfer distances.
Recently, simultaneous transfer of information and power has emerged to be a new research area, which traditionally was treated as two separate problems [11, 12, 13, 14]. In [11], the author investigated the information capacity of a pointtopoint wireless channel under the constraint that the average received power is above a given threshold, establishing the relation between the transferred power and information rate over the same channel. For the channel in an inductivecoupling circuit, the efficiency of power/information transfer increases/decreases with the inputsignal bandwidth. Note that the powertransfer efficiency is the highest at the resonance frequency of the circuit and decreases as the transfer frequency deviates from the resonance point. This results in a tradeoff between the transferred power and information capacity that is characterized in [12]. A twouser broadcast channel is considered in [14] where a multiantenna base station transmits information to one multiantenna terminal and energy to the other. A tradeoff between the power and information rate delivered over the broadcast channel is achieved by precoding at the base station. Different from prior work, this paper considers coupled wireless power and information transfer in a largescale network and the resultant relation between the hybridnetwork parameters can be interpreted as the network counterparts of the existing capacityandpower tradeoffs for smallscale systems.
IB Modeling the Access and PowerBeacon Networks
Cellular networks traditionally are modeled using hexagonal grids but such models lack tractability and fail to account for ad hoc BS deployment that is a trend of network evolution. These drawbacks can be overcome by using stochastic geometry for modeling the cellularnetwork architecture [15, 16, 17, 18]. In [15], a cellular network was modeled as a bivariate PPP that comprises two independent PPPs representing BSs and mobiles, respectively. As a result, the network architecture is distributed as a Poisson Voronoi tessellation. This random network model was applied in [16] to study the economics of cellular networks, and combined with a channelfading model to study the downlink coverage in [17] and uplink coverage in [18], yielding analytical results consistent with the performance of practical networks. An uplink network is considered in this paper where singleantenna mobiles transmit data to singleantenna BSs under an outageprobability constraint for a target received signaltointerferenceandnoise ratio (SINR). Aligned with existing approaches, BSs in the uplink network is modeled as a homogeneous PPP with density and each cell of the resultant Voronoi tessellation comprises a uniformly distributed active mobile. Data links have path loss with the exponent but no fading.
Mobiles in the current cellularnetwork model are powered by MPT rather than reliable power supplies assumed in prior work. Randomly deployed PBs are modeled as a homogeneous PPP similar to the existing models of wireless ad hoc networks (see e.g., [19]). MPT uses frequencies outside the data bandwidth (e.g., in the ISM band) and hence causes no interference to uplink transmissions. Given fixed PB transmission power, the MPT bandwidth is required to be sufficiently large so that the powerspectrum density meets regulations on microwave radiation. Like data links, MPT links have no fading and their pathloss exponent need not be equal to due to the difference in transmission ranges and frequencies. Each PB radiates energy either isotropically or along the directions of targeted mobiles by beamforming, referred to as isotropic and directed MPT, respectively. Mobiles intercept energy transferred from PBs continuously and store it for powering subsequent uplink transmission.
IC Contributions and Organization
The main contributions of the paper are summarized as follows. Let denote a set of constants to be derived in the sequel.

First, consider the deployment of the cellular network. Define the feasibility region as all feasible combinations of the network parameters under the outage constraint. It is proved that
(1) where denotes the set of nonnegative numbers. It is observed that the minimum increases with decreasing and vice versa. Consider the special case of an interferencelimited network with noise omitted and under an additional constraint that requires received signal power at each BS to be above a threshold for given probability. It is shown that outageprobability is independent with and , making the outage constraint irrelevant. The corresponding feasibility region is found to have a similar form as in (1):

Next, consider the deployment of the hybrid network. Define the feasibility region for the hybrid network as all combinations of that satisfy the outage constraint constraint for the cellular network. Note that is determined by . Assume that mobiles have large energy storage. For isotropic MPT, it is proved that
(2) For directed MPT, an inner bound on is obtained that reduces to the following expression if is sufficiently large:
(3) where denotes the array gain, the first condition for the inner bound ensures that PBs are sufficiently dense, and the second condition corresponds to that the maximum power transferred from a PB to an intended mobile is sufficiently high. The results quantify the tradeoffs between and show that the gain of MPT beamforming is equivalent to increasing by a factor of .

Assume that mobiles have small energy storage. The potentially random mobiletransmission power can be stabilized by enforcing the constraint that the received raw power at each mobile exceeds a given threshold with high probability. Under this and the outage constraint, inner bounds on are obtained for isotropic MPT as
(4) and for directed MPT as
(5) Similar remarks as on the results in (3) also apply to those in (4) and (5). Moreover, compared with their counterparts for mobiles with large energy storage, the current tradeoffs between and depend on the pathloss exponents for both data and MPT links. In particular, with other parameters fixed, decreases inversely with increasing if .
The remainder of the paper is organized as follows. The models, metrics and constraints are described in Section II. The feasibility regions for the cellular and hybrid networks are analyzed in Section III and IV, respectively. Numerical and simulation results are presented in Section V followed by concluding remarks in Section VI.
Notation: The notation is summarized in Table I.
Symbol  Meaning 

,  PPP of base stations, density of 
, 
PPP of PBs, density of 

Typical base station, mobile and PB 

Mobile served by base station 
, 
Pathloss exponent for a data channel, a MPT channel 
, 
Outage probability and its maximum for the typical base station 

Transmission power for mobiles, PBs 

Outage threshold 

Set of base stations serving mobiles whose interference to is canceled 

Target received signal power at a BS in an interferencelimited network 

Maximum probability that received signal power at a BS is below in an interferencelimited network 
,  Required received raw power at a mobile and the maximum poweroutage probability for the case of mobiles having small energy storage 
Set of nonnegative numbers  
Number of interfering mobiles each BS cancels 
Ii Models, Metrics and Constraints
Iia AccessNetwork Model
In the uplink network, BSs are modeled as a homogeneous PPP with density , where represents the coordinates of the corresponding BS. Given that mobiles are associated with their nearest BSs, the Euclidean plane is partition into Voronoi cells with points in being the nuclei as illustrated in Fig. 2. Mobiles in the same cell are independent and uniformly distributed, and time share the corresponding BS. All active mobiles transmit signals with fixed power . By applying Slyvnyak’s Theorem [20], a typical base station, denoted as , is assumed to be located at the origin without loss of generality. The active mobile served by is called the typical active mobile and denoted as . It is assumed that uses a multiuser detector [21] to perfectly cancel interference from nearest interfering mobiles. For ease of notation, let denote the active mobile served by BS and represent the set of BSs serving nearest interfering mobiles for . Channels are characterized by path loss based on the longrange propagation model such that signals transmitted by a mobile with power is received by BS with power where is the pathloss exponent. It follows that the signal power received at is and the interference power is given as
(6) 
Last, time is slotted and all channels and point processes including BSs, mobiles and PBs (discussed in the sequel) are assumed fixed within one slot and independent over different slots.
IiB PowerBeacon Network Model
PBs are modeled as a homogeneous PPP, denoted as , with density and independent with the BS and mobile processes. Each mobile deploys a MP receiver with a dedicated antenna as shown in Fig. 3 to intercept microwave energy transmitted by PBs. Inspired by the design in [22], the power receiver in Fig. 3 deploys two energystorage units to enable continuous MPT. When a mobile is active, one unit is used for MPT and the other for powering the transmitter; their roles are switched when the latter unit is fully discharged. When the transmitter is inactive, the MP receiver remains active till all units are fully charged.
Consider isotropic MPT. A shortrange propagation model [23] is used that avoids singularity caused by proximity between PBs and mobiles and thereby ensures finite average raw (unmodulated) power received by a mobile. To be specific, the raw power received at the typical active mobile is given as
(7) 
where is a constant, the pathloss exponent need not be equal to for data links, and the received signal power from other mobiles is neglected as it is much smaller than raw power from PBs.
Next, consider directed MPT. Each mobile is constrained to be charged by only the nearest PB using beamforming. To simplify analysis, the beamforming responses are assumed to have only two fixed levels, namely the mainlobe and sidelobe levels. A PB can serve multiple mobiles simultaneously where the total transmitted raw power is multiplied to ensure the reliability of MPT. Specifically, the transmission power of is where is the random number of mobiles for which is the nearest PB. Based on these assumptions, the raw power received by can be written as
where the constants represent the beam mainlobe and sidelobe responses, respectively. In contrast with the isotropic counterpart, directed MPT requires handshaking between PBs and mobiles e.g., using the following protocol. By estimating the relative distances and identifications of nearby PBs from their continuously broadcast signals, a mobile sends a charging request and pilot symbols to the nearest beacon that then estimate the mobile direction and beams raw power in this direction.
IiC Metrics and Constraints
The cellularnetwork performance is measured by the outage probability that also gives the fraction of mobiles outside network coverage [17]. Define the outage probability for the typical BS as
where is given in (6), denotes the channelnoise variance, and is the target SINR.
Constraint 1 (Cellular network).
The cellularnetwork deployment satisfies the outage constraint: with .
The special case of an interferencelimited network is also considered. To this end, besides the outage constraint, the following one is applied to ensure sufficiently high received signal power at BSs for the cellular network to operate in the interferencelimited regime.
Constraint 2 (Interferencelimited cellular network).
The received signal power at a BS is required to exceed a threshold with probability larger than with and .
Mobiles store energy in either large storage unites (e.g., rechargeable batteries) or small ones (e.g., super capacitors). Large storage at a mobile removes the randomness of received raw power and provides fixed output power (see Lemma 2). However, given small storage, the transmission power of a mobile may fluctuate due to the random locations of PBs. The fluctuation is suppressed by ensuring sufficiently large and stable instantaneous power transferred to each mobile as follows.
Constraint 3 (Power beacons).
Given mobiles having small energy storage, the PB deployment satisfies the constraint that the received raw power at each mobile is below a given threshold (power outage) with probability no larger than a constant with .
The probability is chosen to be sufficiently small such that a powershortage event at an active mobile can be coped with by drawing energy buffered in the small storage units to ensure uninterrupted transmission with power .
Iii AccessNetwork Deployment
In this section, the requirements for the cellularnetwork parameters are investigated by deriving the feasibility region . It is found that for the case of interferencelimited networks is independent with the outage constraint and is determined only by the constraint on the BS received signal power (see Constraint ).
Iiia AccessNetwork Deployment with Nonzero Noise
To characterize the feasibility region, a useful result is obtained as shown in Lemma 1 that is proved in Appendix A using Mapping Theorem [24].
Lemma 1.
Consider the BS process as a function of the density denoted as . The two random variables and follow an identical distribution, where and depends only on .
Furthermore, define such that
(8) 
The main result of this section is shown in the following proposition.
Proposition 1.
Under Constraint , the feasibility region for the cellular network is given as
(9) 
with defined in (8).
Proof: Consider the outage probability as a function of . To facilitate analysis, the outage probability can be rewritten as
(10) 
Using (10) and Lemma 1, it can be obtained that
(11) 
Since the summation in (11) is a continuous random variable, there exists as defined in (8) such that for all . The desired result follows.
Two remarks are in order.

The result in Proposition 2 is consistent with the intuition that deploying denser BSs shortens transmission distances and hence requires lower transmission power at mobiles. In particular, doubling decreases by a factor of , which is more significant for more severe propagation attention (larger ).

It is challenging to derive a closedform expression for defined in (8) and its relation with is evaluated via simulation in the sequel. However, some simple properties of the relation can be inferred including that a) is a strictly monotone decreasing function of , b) corresponds to interference limited networks and c) .
IiiB AccessNetwork Deployment with Zero Noise
Consider an interferencelimited cellular network where noise is negligible.
Proposition 2.
For an interferencelimited cellular network with zero noise, is independent with the BS density and mobiletransmission power .
Proof: By substituting into (10), the outage probability can be written as
(12) 
Let and the BS process be functions of , dented as and , respectively. Using (12) and given , can be written as
(13)  
(14) 
where (13) applies Lemma 1 and (14) follows from (12). The proposition statement is a direct result of (14), completing the proof.
A few remarks are made.

The dual of the result in Proposition 2 for downlink networks with fading and without interference cancelation is shown in [17] using the method of Laplace transform. These results are the consequence of the fact that in both uplink and downlink interferencelimited networks, transmission power has no effect on the signaltointerference ratio (SIR), and increasing the BS density shortens the distances of data and interference links by the same factor and hence also has no effect on the SIR. Therefore, such results hold regardless of if fading is present and interference cancelation is applied.

For interferencelimited networks, the outage probability is largely affected by the outage threshold and the passloss exponent that determines the level of spatial separation.

Though outage probability is independent with in an interferencelimited network, the network throughput grows linearly with increasing since denser active mobiles can be supported.
The result in Proposition 2 implies that Constraint has no effect on the feasibility region for an interferencelimited network. To derive the region under Constraint , the distribution function of the received signal power at the typical BS has no closedform but can be upper bounded as
(15) 
that results from the assignment of mobiles to the nearest BSs. Note that follows the same distribution as based on Slyvnyak’s Theorem [20]. The distribution of the distance between two nearest points in a PPP is well known [25], yielding that
(16) 
Combining (15) and (16) and applying Constraint lead to the following proposition.
Proposition 3.
Consider the deployment of an interferencelimited cellular network under Constraint and . The feasibility region is independent with Constraint and inner bounded as
(17) 
with .
By comparing Proposition 1 and 3, the conditions on and for both the cases of interferencelimited and nonzeronoise networks have the same form, namely with being a constant, and differ only in . The similarity is a consequence of the fact that increasing and affects only the received signal power but not the SIR (see Remark on Proposition 2).
Iv HybridNetwork Deployment
In this section, the hybridnetwork deployment is analyzed by deriving the feasibility region for different PBnetwork configurations combining isotropic/directed MPT and large/small energy storage at mobiles.
Iva HybridNetwork Deployment: Mobiles with Large Energy Storage
In this section, it is assumed that energystorage units (see Fig. 3) at mobiles have infinite capacity. Large storage provides active mobiles reliable transmission power as specified in the following lemma directly following from [26, Theorem 1] that studies energy harvesting in wireless ad hoc networks.
Lemma 2 ([26]).
For mobiles with infinite energy storage and powered by MPT, the probability that a mobile can transmit signals with power is given as
(18) 
where is the duty cycle defined as the probability that a mobile transmits a packet in a slot and is the received raw power in (7) and (31) for isotropic and directed MPT, respectively.
Lemma 2 suggests that an active mobile can transmit signals continuously with power up to . Closedform expressions for are derived by applying Campbell’s Theorem [24]. The results are shown in Lemma 3.
Lemma 3.
The average received raw power for a mobile is given as follows:

for isotropic MPT,
(19) 
for directed MPT,
(20) where
The proof of Lemma 3 is provided in Appendix B. For a sanity check, substituting into (20) gives the result for isotropic MPT in (19). Two remarks are given as follows.

For isotropic MPT, a mobile is powered by all nearby PBs and thus the expected received raw power given in (19) is proportional to the PB density . The value of is large if the pathloss exponent for MPT links is closed to , which corresponds to the scenario where there are LOS links between a mobile and many PBs. In addition, is also proportional to the PBtransmission power and duty cycle .

For directed MPT, there is only one PB transferring significant raw power to each mobile. Consequently, the first term at the righthand side of (20) is dominant given that and hence where quantifies the propagation loss. Furthermore, if , and hence that is equal to the counterpart for isotropic MPT in (19) with the factor replaced with . It is important to note that these two factors represent respectively the advantages of isotropic and directed MPT, namely that strong raw power at a mobile comes form a large number of PBs for the case of isotropic MPT with low propagation loss ( is close to ) and from sharp beamforming for the case of directed MPT (large ).
The main result of this section is shown in Proposition 4, which follows from combining Proposition 1, Lemma 2 and 3, and the fact that the last term in (20) is positive (see the proof of Lemma 3).
Proposition 4.
Given mobiles with infinite energy storage, the feasibility region for the hybrid network deployed under the outage constraint (Constraint ) is given as follows.

For isotropic MPT,
(21) 
For directed MPT, is inner bounded as
(22) with .
Two remarks are in order.

For isotropic MPT, the required PB density decreases with the increasing PB transmission power and BS density as shown in (21). The result also shows the effects of the pathloss exponents and for data and MPT links, respectively. Specifically, the required PB density diminishes as approaches corresponding to feespace propagation and increases with growing for which datalink path loss is more severe and higher transmission power at mobiles is needed. In addition, for both isotropic and directed MPT, the required PB density is proportional to the mobile duty cycle .

Consider the inner bound on in (22). The first condition ensures that PBs and their served mobiles are sufficiently close (within a distance of ) with high probability. Under the second condition, a single PB is able to provide sufficient power to an intended and nearby mobile by either sharp beamforming (large ) or/and high transmission power (large ).
If the ratio is large, the inner bound on can be simplified as shown in the following corollary.
Corollary 1.
Given mobiles with infinite energy storage and powered by directed MPT, as , the feasibility region for the hybrid network deployed under the outage constraint (Constraint ) can be inner bounded as
(23) 
IvB HybridNetwork Deployment: Mobiles with Small Energy Storage
For the current case, the deployment of the hybrid network has to satisfy both Constraint and . Th derivation of the resultant feasibility region requires analyzing the cumulativedistribution function (CDF) of the received raw power at a mobile that is a summation over a PPP. Though the function has no closed form [27], it can be upper bounded as shown in the following lemma by considering only the nearest PB for a mobile.
Lemma 4.
Given mobiles with small energy storage, the poweroutage probability is upper bounded as follows:

for isotropic MPT, if ,
(24) 
for directional MPT, if ,
(25)
The proof of Lemma 4 is provided in Appendix C. Combining the results in Proposition 1, Lemma 4 and applying Constraint give the main result of this section as shown below.
Proposition 5.
Given mobiles with small energy storage, the feasibility region for the hybrid network deployed under Constraint and can be inner bounded as follows:

for isotropic MPT,
(26) 
for directional MPT,
(27)
It can be observed from the results that the product of the PB parameters, namely , is proportional to the logarithm of the maximum powershortage probability and hence insensitive to changes on . Moreover, the results in Proposition 5 show that beamforming reduces the required by the factor of . Like the case of mobiles having large energy storage, the required value of for the current case also decreases with increasing BS density but at a different rate.
Last, the results in Proposition 5 have their counterparts for the interferencelimited network as shown in the following corollary obtained similarly as Corollary 5.
Corollary 3.
By replacing with , the results in Proposition 5 also apply to the case of an interferencelimited cellular network with zero noise and mobiles having small energy storage.
V Numerical and Simulation Results
The results in the section are obtained based on the following values of parameters unless specified otherwise. The pathloss exponents are and and the outage threshold is . Each BS cancels interference from nearest unintended mobiles. Considering Constraint , the maximum outage probability for each BS is . The parameters in Constraint and are set as dB and while is a variable. Last, the PB transmission power is fixed as dB.^{1}^{1}1In the pathloss models, the mobile/PB transmission power is normalized by a reference value measured at a reference distance that is used to normalize all propagation distances [28].
Va AccessNetwork Deployment
The curve of versus defined in (8) is plotted in Fig. 4 and observed to be consistent with Remark on Proposition 1. In particular, approaches as increases and its minimum value at gives the outage probability () for an interferencelimited network. In addition, it is observed that corresponds to considered for subsequent simulations.
The feasibility regions for the cellular network as given in Proposition 1 and 2 are plotted in Fig. 5 for both the cases of nonzero and zero noise (an interferencelimited network). The boundaries of two feasibility regions are parallel as they have the same form. The minimum in dB is observed to decrease linearly with increasing .
VB HybridNetwork Deployment
Fig. 6 shows the curves of mobiletransmission power versus BS density for different combinations between large/small energy storage at mobiles and isotropic/directional MPT. The curves for the case of large energy storage are computed numerically using Lemma 3 and those for the case of small storage are obtained by simulation. On one hand, the values of for directional MPT are observed to be insensitive to changes on if it is sufficiently large. The reason is that is dominated by the maximum raw power transferred to a mobile from the nearest PB that is a constant . On the other hand, the values of (in dB) for isotropic MPT do not saturate and increase approximately logarithmically with growing . The reason is that for this case depends on power transferred by many PBs instead of a dominant one and hence are sensitive to changes on .
In Fig. 7, the feasibility regions for the hybrid network are plotted for different cases accounting for large/small energy storage at mobiles and isotropic/directional MPT. The feasibility region for the case of isotropic MPT with mobiles having large energy storage is obtained using Proposition 1 and those for other cases by simulation. Also shown are inner bounds on the feasibility regions based on Proposition 4 and 5 with their boundaries plotted with dashed lines. The bounds are observed to be tight except for the case of isotropic MPT and mobiles with small energy storage due to the lack of a dominant PB for each mobile in terms of received raw power. By inspecting the feasibility regions for directed MPT, there exists a threshold on the BS density below which the required PB density grows extremely rapidly. The BS density above the threshold corresponds to the regime where the power transferred to a mobile from the nearest PB is larger than the minimum required for satisfying the outage constraint on the cellular network. Otherwise, besides receiving power from the nearest PB, mobiles rely on accumulating power from unintended PBs via their sidelobe transmissions so as to satisfy the outage constraint, resulting in a high required PB density. This is the reason that isotropic MPT enlarges the feasibility region for relatively small BS densities and directed MPT is advantageous for relatively large BS densities as observed from Fig. 7
Vi Concluding Remarks
The deployment of PBs for powering a cellular network via MPT has been investigated based on a stochasticgeometry network model, resulting in simple tradeoffs between network parameters. First, considering the cellular network under an outage constraint, the minimum mobiletransmission power has been shown to increase superlinearly with the decreasing basestation density . However, the absence of noise renders outage probability independent with . Next, building on these results, the requirements on the PB deployment have been studied by analyzing the tradeoff between and the PB transmission power and density . Given isotropic MPT and mobiles having large capacity for storing transferred energy, the product is required to be above a given threshold so as to satisfy the outage constraint. For the case of mobiles having small energy storage, both the products and should be sufficient large. It has been found that compared with isotropic MPT, directed MPT by beamforming effectively increases by the array gain.
This work relies on a stochasticgeometry network model and various simplifying assumptions to derive tradeoffs between network parameters. The results provide firstorder guidelines for realizing MPT in cellular networks. To derive more elaborate insight into MPT implementation, further investigations in practical settings are necessary by considering the dependance of the directedMPT efficiency on the array configurations, using realistic channel models for smallcell networks, modeling the hotspot deployment of PBs using a clustered point process, characterizing the effect of mobility on MPT, and taking into account the dynamics of energylevels at mobiles.
References
 [1] D. GrahamRowe, “Wireless power harvesting for cell phones,” Tech. Review, Jun. 2009.
 [2] T. L. Marzetta, “Noncooperative cellular wireless with unlimited numbers of base station antennas,” IEEE Trans. on Wireless Comm., vol. 9, pp. 3590–3600, Sep. 2010.
 [3] V. Chandrasekhar, J. G. Andrews, and A. Gatherer, “Femtocell networks: A survey,” IEEE Comm. Magazine, vol. 46, pp. 59–67, Sep. 2008.
 [4] R. Want, “An introduction to RFID technology,” IEEE Pervasive Computing, vol. 5, pp. 25–33, Jan. 2006.
 [5] A. Kurs, A. Karalis, R. Moffatt, J. D. Joannopoulos, P. Fisher, and M. Soljačić, “Wireless power transfer via strongly coupled magnetic resonances,” Science, vol. 317, pp. 83–86, Jul. 2007.
 [6] W. C. Brown, “The history of power transmission by radio waves,” IEEE Trans. on Microwave Theory and Techniques, vol. 32, pp. 1230–1242, Sep. 1984.
 [7] W. C. Brown, “Experiments involving a microwave beam to power and position a helicopter,” IEEE Trans. on Aerospace and Electronic Sys., vol. AES5, pp. 692–702, Sep. 1969.
 [8] J. J. Schlesak, A. Alden, and T. Ohno, “A microwave powered high altitude platform,” IEEE MTTS Digest, pp. 283–286, 1988.
 [9] H. Matsumoto, “Research on solar power satellites and microwave power transmission in Japan,” IEEE Microwave Magazine, vol. 3, pp. 36–45, Dec. 2002.
 [10] J. O. Mcspadden and J. C. Mankins, “Space solar power programs and microwave wireless power transmission technology,” IEEE Microwave Magazine, vol. 3, pp. 46–57, Apr. 2002.
 [11] L. R. Varshney, “Transporting information and energy simultaneously,” in Proc., IEEE Intl. Symposium on Information Theory, pp. 1612–1616, Jul. 2008.
 [12] P. Grover and A. Sahai, “Shannon meets Tesla: Wireless information and power transfer,” in Proc., IEEE Intl. Symposium on Information Theory, pp. 2363–2367, Jun. 2010.
 [13] X. Zhou, R. Zhang, and C. Ho, “Wireless information and power transfer: Architecture design and rateenergy tradeoff,” in submitted to IEEE Globecom (Available: http://arxiv.org/abs/1205.0618).
 [14] R. Zhang and C. Ho, “MIMO broadcasting for simultaneous wireless information and power transfer,” submitted to IEEE Tarns. on Comm. (Avaiable: http://arxiv.org/abs/1105.4999).
 [15] S. G. Foss and S. A. Zuyev, “On a Voronoi aggregative process related to a bivariate Poisson process,” Advances in Applied Prob., vol. 28, pp. 965–981, Dec. 1996.
 [16] F. Baccelli, M. Klein, M. Lebourges, and S. Zuyev, “Stochastic geometry and architecture of communication networks,” Telecomm. Systems, vol. 7, no. 13, pp. 209–227, 1997.
 [17] J. G. Andrews, F. Baccelli, and R. K. Ganti, “A tractable approach to coverage and rate in cellular networks,” IEEE Trans. on Comm., vol. 59, pp. 3122–3134, Nov. 2011.
 [18] T. D. Novlan, H. S. Dhillon, and J. G. Andrews, “Analytical modeling of uplink cellular networks,” submitted to IEEE Tarns. on Wireless Comm. (Avaiable: http://arxiv.org/abs/1203.1304).
 [19] M. Haenggi, J. G. Andrews, F. Baccelli, O. Dousse, and M. Franceschetti, “Stochastic geometry and random graphs for the analysis and design of wireless networks,” IEEE Journal on Selected Areas in Comm., vol. 27, pp. 1029–1046, Jul. 2009.
 [20] D. Stoyan, W. S. Kendall, and J. Mecke, Stochastic Gemoetry and its Applications. Wiley, 2nd ed., 1995.
 [21] S. Verdu, Multiuser Detection. Cambridge, UK: Cambridge, 1998.
 [22] S. Luo, R. Zhang, and T. J. Lim, “Optimal savethentransmit protocol for energy harvesting wireless transmitters,” submitted to IEEE Tarns. on Wireless Comm. (Avaiable: http://arxiv.org/abs/1204.1240).
 [23] F. Baccelli, B. Blaszczyszyn, and P. Muhlethaler, “An ALOHA protocol for multihop mobile wireless networks,” IEEE Trans. on Information Theory, vol. 52, pp. 421–36, Feb. 2006.
 [24] J. F. C. Kingman, Poisson processes. Oxford University Press, 1993.
 [25] M. Haenggi, “On distances in uniformly random networks,” IEEE Trans. on Information Theory, vol. 51, pp. 3584–86, Oct. 2005.
 [26] K. Huang, “Spatial throughput of mobile ad hoc networks with energy harvesting,” submitted to IEEE Trans. on Information Theory (Available: http://arxiv.org/abs/1111.5799).
 [27] S. B. Lowen and M. C. Teich, “Powerlaw shot noise,” IEEE Trans. on Information Theory, vol. 36, pp. 1302–1318, Nov. 1990.
 [28] T. S. Rappaport, Wireless Communications Principles and Practice. Upper Saddle River, NJ: Prentice Hall, second ed., 2002.
 [29] S. P. Weber, J. G. Andrews, X. Yang, and G. de Veciana, “Transmission capacity of wireless ad hoc networks with successive interference cancelation,” IEEE Trans. on Information Theory, vol. 53, pp. 2799–2814, Aug. 2007.
a Proof of Lemma 1
To facilitate the proof, redenote as . It is almost surely that given and an arbitrary active mobile , there exists a set of coefficients such that . Consequently, a set of random coefficients can be defined under the constraint that is uniformly distributed in the corresponding Voronoi cell. It follows that
(28) 
Given and , define a function such that for , . Intuitively, expands (or shrink) the Euclidean plane by a factor if (or ). Let denote the mean measure of and define a measure with as . It follows that where denotes the measure (area) of . Since is a measurable function that does not maps distinct points to a single point, applying Mapping Theorem [24, p18] gives that is a homogeneous PPP with density . As a result, the random variable is identically distributed as . Combining this fact and (28) yields the desired result.
B Proof of Lemma 3
For isotropic MPT, using (7) and applying Campbell’s Theorem give
(29)  
(30) 
where (29) results from the stationarity of and using the polarcoordinate system. The corresponding expression of in (19) follows from (30).
Next, consider directed MPT and the corresponding expression for is derived as follows. To facilitate analysis, is rewritten from (31) as
(31) 
For ease of notation, define . Using (31) and using the result in (19),
(32) 
Note that measures the shortest distance between a point in the PPP to a fixed point and have the following probabilitydensity function [25]
(33) 
Given the PDF, it is obtained that
(34) 
Substituting (34) into (32) gives the desired result in (20), completing the proof.
C Proof of Lemma 4
Without loss of generality, assume is located at the origin since according to Slyvnyak’s Theorem,
(35) 
where is given in (7) and (31) for isotropic and directed MPT, respectively. The previous assumption that is at the origin is unnecessary for the proof.
Consider isotropic MPT. Define . Note that any PB with a distance from can supply received raw power higher than . Furthermore, as a result of the assumption . Inspired by the approach in [29] that studies outage probability for mobile ad hoc networks, the probability of no power shortage is expanded as
(36) 
where is a disk centered at and with a radius . It follows that the powershortage probability can be lower bounded as
(37) 
The derived result in in the lemma statement follows.
The result for directed MPT can be obtained similarly. This completes the proof.