Outage performance and energy efficiency optimization of wireless‑powered millimeter‑wave sensor networks

variety

the sensor can receive a continuous power supply without having to replace the battery [2,3].Although solar energy, kinetic energy, and floor vibration can be used to power the sensors, the dynamic nature of these resources, as well as a lack of prior knowledge of energy distribution will bring great difficulties to the protocol design of the systems.In contrast, harvesting energy from ambient radio frequency (RF) signals presents a potentially more efficient approach owing to its wireless nature, ready availability, low cost, and the feasibility of small form factor implementation [4].This method proves particularly suitable for wireless energy transfer (WPT) applications in scenarios such as smart homes, forest fire monitoring, and indoor IoT systems [5][6][7].With the growing demand of RF energy harvesting, commercialized products, such as Powercaster and Cota systems have been entered the market [8].
Although WPT techniques promise tremendous benefits, some key challenges must be resolved before they are deployed in the WSNs.Firstly, because the great attenuation of the electromagnetic waves with the increase of distance, only a portion of the energy radiated by an energy transmitter can be harvested by an energy receiver [9].Secondly, in WPT, the sensitivity of the current energy receivers is as high as − 10 dBm, while the sensitivity of the information receiver is only − 60 dBm [10].At present, there are three main approaches to solve the above problems: i) forming the RF energy into narrow beams; ii) keeping the distance between transmitter and receiver as smaller as possible in the design of network architecture; iii) adopting massive MIMO technique to enhance spectral efficiency and energy efficiency (EE) [11].
n the fifth-generation (5G) communications, the implementation of cost-effective, low-power, and ultra-dense cell systems has the potential to enhance network performance in terms of energy efficiency and load balancing.Due to the physical properties of mmWave communication, i.e.,,,, a frequency band ranging from 30 to 300 GHz and corresponding wavelengths from 10 to 1 mm [12], mmWave technology can leverage larger quantities of available spectrum.Additionally, it can utilize massive antenna arrays to facilitate highly directional transmission and reception.This capability is achievable at higher frequencies owing to the smaller size of antenna components [13].Furthermore, mmWave communications can be realized in the range of 150-200 m in the urban environments [14].As a result, mmWave is considered as the prime candidate for 5G systems, while mmWave can offer all three solutions to enhance WPT performance [15].Therefore, the above factors ensure the availability and feasibility of WPT by mmWave in 5G networks.

Related work
In communication systems, WPT can be realized by utilizing the following two architectures, i.e.,,,, wireless-powered communication networks (WPCNs) and simultaneous wireless information and power transfer (SWIPT) [16,17].For WPCNs, the harvestthen-transmit (HTT) protocol was proposed with a general multi-user wireless-powered interference channel in [18], where the communication time is divided into two phases, i.e.,,,, wireless energy harvesting (WEH) and wireless information transmission (WIT).
The key of SWIPT is that both information and energy can be extracted simultaneously from the same received RF signals via time-switching (TS) scheme or power-splitting (PS) scheme [19].Two SWIPT cooperative spectrum sharing strategies were proposed aiming to improve the spectral and energy efficiencies for cognitive IoT networks in 6G communication systems, where IoT devices can harvest energy from RF signals and can be acted as orthogonal frequency division multiple (OFDM) relay to access the primary spectrum [20].Multicast communication of a satellite and aerial-integrated network with rate-splitting multiple access was investigated to satisfy the demands of massive IoT devices accessing [21].
At present, a large number of WPT related studies have been carried out in microwave networks (below 6 GHz).In [22], the authors developed a WPT-based WIT scheme in relay-based CRNs, where a secondary network can share the spectrum and harvest energy to assist the data transmission of the primary network.The authors in [23] considered a multiple relay-assisted cooperative cognitive radio system, where the relays first harvested RF energy from the primary transmitter and then adopted decode-andforward (DF) method to transmit data.In addition, three best relay selection schemes were also investigated.The authors of [24] considered a two-way amplify-and-forward (AF) relay network, where two nodes exchange information with assistance of an energy harvesting MIMO relay.The authors further optimized the power allocation ratios, precoding matrix, and RF phase shifting factor to maximize the sum-rate of the considered network [25].investigated the wireless-powered relay network with harvest-receiveforward time-switching scheme, where the relay can random harvest RF energy from dedicated energy signal and other ambient resources.The maximal offline throughput in different scenarios was considered by optimizing the time fraction and reception rate.A novel architecture design for the green WSNs in smart cities was considered in [26], where an efficient energy and information transfer protocol was proposed to enhance the performance of transmission rate.In addition, beamforming design is other key factor for WEH and WIT.In [11], the authors considered a hybrid satellite-terrestrial relay network to optimally joint design beamforming vectors.The authors of [27] proposed a joint optimization design for a NOMA-based satellite-terrestrial integrated networks, which can share the mmWave spectrum with a cellular network.
Although the research on mmWave communications and related topics have been extensive, the investigation on the combination of mmWave and WPT is still limited.The authors of [28] proposed an analytical framework to investigate capability of mmWave WPT with large-scale antennas, i.e.,,,, energy coverage probability, average harvested energy, and achievable rate, by utilizing stochastic geometry method.In [29], the authors analyzed the performance of WPT in a large-scale mmWave tactical network, in which a person equipped with WEH nodes were powered by another person equipped with power beacon (PB) nodes.Taking human as potential blockages, the energy coverage probability and overall success probability were derived.A directional power transfer scheme based on beamforming was proposed in [30] for a large-scale network, where each PB can steer its beam toward the direction with maximum receiving power during the training phase.The aforementioned architectures for microwave networks have been investigated in the area of mmWave communications.A SWPITenabled mmWave network is considered in [31].A corresponding optimal PS scheme was proposed to improve the system performance in terms of reducing communication interruption and information loss.The authors in [32] considered a mmWave MIMO-NOMA system with SWIPT and studied maximizing the achievable sum-rate and total energy efficiency.Furthermore, the secrecy performance of SWIPT-enable mmWave ultra-dense networks was investigated with a stochastic geometry framework [33].In order to both improve the system security and reduce energy consumption, the authors in [34] proposed an alternative optimization algorithm to address the problem.In [35], the authors represented a secrecy-energy efficient hybrid beamforming design for mmWave-based satellite-terrestrial integrated network, where the sub-optimal solutions with two robust beamforming schemes were obtained to solve the considered optimization problem.Moreover, the multibeam satellite system with secure energy efficient beamforming was investigated in [36], where an alternating optimization scheme was applied to address the secrecy energy efficiency problem.

Motivation and contribution
With the further development of 5G technologies, IoT systems based on mmWave sensor networks will be implemented everywhere.This paper investigates a typical mmWave sensor network where each sensor can harvest RF energy in downlink and then use the harvested energy for data transmission in uplink.This model can provide an effective solution for sensors that are inconvenient for battery replacement or continuous power supply in IoT systems, such as intelligent transportation systems and human data monitoring.Unlike [31,32], which either only focuses on deriving close-form expressions of system performance or only optimizes the energy efficiency of the system.In order to further study the wireless-powered mmWave sensor networks, we first theoretically analyze the system performance, which aims to derive the closed-form solution of the beam outage probability within a stochastic geometry framework random in uplink after the sensors harvest RF energy in downlink.Then, according to the proposed transmission protocol, the transmission slots in downlink and uplink are jointly optimized to guarantee the maximum energy efficiency in uplink.The main contributions can be summarized as follows: • We investigate a wireless-powered mmWave sensor network within a cellular, in which the hybrid access point (HAP) is located at the central point of the cellular network and the sensors are spatially distributed around it.The HTT protocol is adopted, where the whole transmission slot is divided into two sub-slots.In the first sub-slot, the HAP transmits RF energy signal to multiple sensors through the downlink.In the second sub-slot, these sensors use the harvested energy to transmit their individual signals to the HAP under the TDMA mode.According to the proposed transmission protocol, an analytical expression of the system beam outage probability is derived by using stochastic geometry theory.In the meanwhile, for two special cases, the approximate closed expressions of the system beam outage probability are further derived.• In order to maximize the energy efficiency of the system, it can be achieved by optimizing the transmission slot allocations of WPT in downlink and WIT in uplink.Therefore, we propose a corresponding optimization problem, which can be can be converted into its parametric form.Then, the new formed optimization problem can be solved by adopting the Lagrange dual multipliers method with corresponding KKT conditions.Finally, the Dinkelbach algorithm is adopted to iteratively obtain the global optimal timeslot allocation.
The numerical results provide useful insights that the beam outage probability of the considered system is significantly impacted by various parameters, such as central angle of service sector, the number of antennas, fading severity parameters, density of homogeneous Poisson point process (HPPP), etc.Moreover, the simulation results also reveal that the proposed algorithm for optimal timeslot allocation can effectively improve the performance of the system compared to the corresponding ergodic analyses in terms of energy efficiency and sum achievable rate.
The rest of this paper is organized as follows.In Section II, network topology, channel model, and transmission protocol are presented.Section III derives the beam outage probability of the system and optimal time allocation algorithm is also proposed in this section.Section IV contains simulation results.Section V summarizes this paper.

Methods/experimental
In this paper, we mainly study the outage probability when the signal is transmitted by energy-constrained sensors in a wireless-powered millimeter-wave sensor network (WP-MWSN), and optimize the downlink and uplink transmission slot allocation to improve the EE of sensors in uplink transmission.The beam outage probability is calculated by obtaining the cumulative distribution function (CDF) of the achievable rate at the sensor.Then, by constructing the Lagrange function corresponding to the optimization problem and the KKT condition, Dinkelhach algorithm is utilized to iteratively obtain the optimal allocations of downlink and uplink transmission slots.Finally, Matlab simulation is used to demonstrate the performance of the derived outage probability and verify whether it matches with the corresponding Monte Carlo simulation.In addition, Matlab simulation is further used to verify the advantages of optimized downlink and uplink slot allocation in improving the EE of the system.

Network topology
As shown in Fig. 1, we consider a wireless-powered millimeter-wave sensor network (WP-MWSN) consisting of one HAP and multiple randomly deployed sensors in a single cell, which is modeled as a disc with radius R. The HAP is plugged into a constant power supply and equipped with a uniform linear array of M antennas.All sensors harvest energy from the HAP's transmitted energy signals and equipped with a small mmWave antenna array, where several power dividers are used to combine sub-arrays into one channel for signal reception and transmission [37].As a result, the small antenna array on the sensor can be regarded as an antenna unit.We assume that the HAP is located at the origin of the disc.The spatial distribution of sensors is followed as a homogeneous Poisson point process (HPPP) ={x k } , k ≥ 1 , in Euclidean plane with density , where x k is defined as the coordinates of the k-th sensor.
The protocol on harvest-then-transmit (HTT) is followed in this paper.According to the TDMA scheme, each transmission interval T is divided into two sub-intervals for downlink and uplink, respectively.Sensors first preform WEH in the downlink and then send data in the uplink towards WIT.The received noise at all receivers is additive white Gaussian noise (AWGN) with zero mean and average power σ 2 0 .

Channel model
Measurement findings showed that mmWave links undergo different channel conditions, which are characterized by a line of sight (LOS) link and several non-line of sight (NLOS) links [38].Therefore, the mmWave link between the HAP and k-th sensor can be modeled as where g k,L and δ k,L denote the complex gain and normalized direction of the link with k-th sensor for LOS path, respectively; g k,NL and δ k,NL represent the complex gain and normalized direction of the link with k-th sensor for NLOS path, respectively; β L and β NL represent the path-loss exponents for LOS path and NLOS path, respectively; d k represents the Euclidean distance from the HAP to the k-th sensor.In addition, a(δ k ) is the array steering vector, which can be expressed as where (•) T indicates the transpose of matrix, δ k =2d sin (φ k ) ω denotes the normalized direction of the LOS link with the k-th sensor, where φ k ∈ −π 2, π 2 is the physical angle of the k-th sensor, d and ω denote the distance between two adjacent antennas and the carrier wavelength, respectively.Furthermore, each mmWave channel has different propagation property since it characterized by a strong LOS link as if in free space and several weak NLOS links.In the mmWave channel, the path-loss of NLOS component is much larger than the path-loss of LOS component, which because that the power of g k,L associated with NLOS link is lower 20 dB than the complex gain associated with LOS link [39].Therefore, the effect of NLOS links can be neglected, i.e.,,,, g l k,NL = 0 , ∀l , and only LOS component is consid- ered for mmWave links.The mmWave channel response of the link between the HAP and k-th sensor can be simplified as (1) where g k and δ k represent the complex gain and normalized direction of the link with k- th sensor for LOS path, respectively, with g k assuming to be a complex normal random variable, i.e.,,,, g k ∼ CN (0, 1) .It should be noted that the analytical method in this paper can be used to analyze the outage performance of NLOS user, where the probability density function (PDF) of the channel gain between the HAP and NLOS user, and then substituting it into the expression.

Energy harvesting and information transmission
Based on the HTT protocol, the whole transmission interval T can be divided into WEH and WIT sub-intervals, which is shown in Fig. 2. Assuming K sensors existed in the service sector, all sensors first harvest energy in the downlink during the WEH sub-interval, then the WIT sub-interval of uplink is divided into K slots and all sensors transmit their individual signal to HAP with TDMA mode by utilizing the whole harvested energy of each sensor.
During WEH sub-interval α 0 T , HAP sends energy beamforming vector to serve multiple sensors with a normalized direction δ .The energy beam is given by where a δ can be obtained from (2) with δ being distributed uniformly in normalized direction in [−1, 1] .As shown in Fig. 1, the sensors are located in the fan-shaped sector with normalized direction δ with central angle 2δ , which can be obtain through a low- complexity training process or the Global Positioning System (GPS).Then, the effective channel gain of the k-th sensor is given by where (3) 2 Frame structure of TDMA-based HTT protocol for the considered WP-MWSN denotes the Fejér kernel [40].To simplify the analysis, we set δ=0 and the received RF power at the k-th sensor can be expressed as where P HAP is constant transmit power of the HAP.Then, the amount of harvested energy at the k-th sensor based on the linear WEH model can be calculated as where η ∈ (0, 1] represents the energy conversion efficiency.The harvested energy from the AWGN is negligible.Note that many studies focus on non-linear WEH since it can better reflect the saturation effect of energy harvesting.However, we consider the lowpower mmWave communication systems, in which the transmit power of mmWave base station is relatively low, so the rectification process can be approximately expressed as linear [41].Therefore, linear WEH model is adopted in this paper. During the WIT sub-interval (1 − α 0 )T , with the specified time-slot, K sensors will transfer information to HAP using all harvested energy with TDMA mode.As a result, the achievable throughput of the k-th sensor at the HAP is given by where α k T denotes the allocated time duration of WIT for the k-th sensor, h ′ k and g ′ k are the channel and complex gain between the k-th sensor and HAP, respectively.

Beam outage probability and optimal time allocation
In this section, we commence with an analysis of the beam outage probability within the considered system.Subsequently, we determine the optimal time allocation for WEH and WIT involving K sensors.The focus of this analysis is on a representative sensor within the service sector; thus, for the sake of simplicity in description, we omit the subscripts used in the previous section.Throughout the following discussions, we set T = 1.To facilitate a clearer expression of the beam outage probability, we establish a scenario where the transmission duration in the uplink for each user is identical.This choice is motivated by the difficulty in obtaining the probability density function for any arbitrary timeslot assignment.

Beam outage probability
The outage probability is a vital performance measure to characterize the performance of wireless networks.Mathematically, the beam outage probability is defined as the probability that the achievable rate of a sensor is below a given threshold r.In order to insight the closed-form expression of the beam outage probability for the considered system, we simplify to set that the remaining WIT duration is divided equally for the K sensors.It should be noted that if the sub-timeslot allocation of each (7) sensor is different, then the probability density function (PDF) of the sub-timeslot allocation needs to be obtained in advance for calculating the exactly expression of beam outage probability.Since the optimization goal of the system needs to be determined according to the actual application, the PDF of sub-timeslot allocation is difficult to obtain.Thus, we equally divided the uplink transmission timeslot for each sensor.Considering K sensors located at the service sector to transmit their individual signal to the HAP, the expression of beam outage probability is given by where r denotes the target transmission rate.In order to derive the closed-form expression of (10), we first derive the cumulative distribution function (CDF) of a typical sensor's achievable rate, which is shown as the following Lemma.
Lemma 1: Assume there are K sensors within the service sector, the CDF of achievable rate for a typical sensor at the HAP can be expressed as where Proof: Please refer to Appendix.
It should be noted that it is difficult to obtain the closed-form expression of (11) since an integral part existed.In order to obtain further insights from Lemma 1, we simplify some components of the derived expression to obtain approximate results of (11) as follows: 1) Special Case 1: This case corresponds to a high SNR scenario with a small number of sensors, i.e.,,,, the SNR of transmit power for HAP is high.Then, we have A → 0 .On the basis of the fact that exp −y ≈ 1 − y for y → 0 , (11) can be rewritten as 2) Special Case 2: In this case, in addition to considering the case of high SNR and small number of sensors, we also consider the case where the central angle 2δ of ser- vice sector is small, i.e.,,,,δ → 0 .Since sin y ≈ y for y → 0 and sinc y =sin y y , we can get F M (φ) ≈ Msinc πφM 2 .Furthermore, based on the fact that sinc y ≈ 1 − y 2 6 and 1 − y 2 ≈ 1 − 2y for y → 0 , the Fejér kernel can be approxi- mately expressed as (10) dφ.
Substituting ( 13) into (12) and utilizing the fact that 1 − y −1 ≈ 1 + y for y → 0 , the result of ( 12) can be further approximated as Then, the beam outage probability is shown as following proposition.
Proposition 1: Suppose there are K sensors within the service sector, the beam outage probability for the considered system is given by.
where µ= δπR 2 is the area of the service sector with central angle 2δ , F AR (•) is given in (11).
Then, the probability that the area of the service sector with density contains K sensors is given as µ K exp (−µ) K ![42].By taking the expectation with respect to the amount of sensors K, the proposition can be proved.

Optimal time allocation
When the system beam outage probability is analyzed above, the transmission timeslot in uplink is evenly divided to each sensor for convenience of calculation.However, in practical applications, uplink sub-timeslots can be allocated more reasonably according to different application requirements.In order to further study how to improve the system performance, we aim to maximize the energy efficiency (EE) of the system, and realize the efficient green communications by optimizing the allocations of WET subtimeslot in downlink and WIT sub-timeslots based on TDMA operation.Therefore, we aim to maximize the EE of the considered system by optimizing the allocation of WIT sub-intervals for sensors in the service sector.The EE is defined as where R sum and E T denote the system throughput with K sensors and total energy con- sumption, respectively, which are given by (13) Note that in the TDMA mode, the HAP can rely on channel reciprocity to acquire full CSI of uplink when the uplink and downlink share the same band [43].According to the channel reciprocity, the result of effective complex gain for uplink g ′ k 2 in ( 9) is equal to g k
In (20), α 0 and {α k } represent the downlink energy transmission duration and uplink data transmission duration for each sensor, respectively.Therefore, the constraint C1 in (20) indicates that the sum of α 0 and {α k } is less than or equal to a whole transmission interval.In the constraint C2, we set α 0 > 0 and α k ≥ 0 , where α k = 0 occurs when the k-th sensor does not harvest RF energy.To obtain the optimal pair of (α 0 , {α k }) , we first prove that the objective function in ( 20) is a concave function.Since the constraint C1 α 0 + K k=1 α k ≤ T illustrating the parameters between the α 0 and {α k } are coupled, we thus need to the con- cavity of the subject function when α 0 is fixed.Correspondingly, the second derivative of the subjective function can be equivalently expressed as which means that optimization problem (P1) is actually a concave optimization problem [44].Therefore, the optimal solution can be obtained by transferring the objective function in (P1) into parametric form and then adopting the Lagrange dual multipliers method.The parametric form of objective function in optimization problem (P1) can be rewritten as where θ denotes a non-negative parameter.For a given θ , the Lagrangian function of optimization problem (P1) is given by ( 19) where ξ is Lagrangian coefficient.Based on the complementary slackness, the optimal lagrange multipliers are zero for constrain C2.We thus omit the corresponding lagrangian multipliers in (23).According to the KKT conditions of lagrangian function (23), we have respectively.From (25), result of ξ can be obtained and then substituting it into (24), we have where � k = α 0 α k .Therefore, all k s ∀k ∈ {1, 2, . . ., K } can be obtained from (27).For (26), we have T =α 0 + K k=1 α k .The optimal solution for pair of (α 0 , {α k }) can be derived as follows Finally, optimal θ and globally optimal solution for the pair of α * 0 , α * k can be decided by adopting Dinkelbach algorithm with iterating, which is shown in Algorithm 1. ( 23)

Numerical results
In this section, we present some simulation results to verify the accuracy of the theoretical analyses and demonstrate the availability for the proposed optimal time allocation algorithm.Table 1 lists the necessary parameters for numerical results.

Beam outage probability
Figure 3 depicts the beam outage probability versus the transmit SNR of HAP for various target rates.With the increase of the transmit SNR of the HAP, the beam outage probability of the system is improved because the sensors in the service sector can harvest more energy for data transmission.In addition, a higher target rate will result in a higher beam outage probability since the links are more difficult to satisfy the transmission requirement.All analytical results are in good agreement with the Monte Carlo simulations, which verify the accuracy of the derived expression.
Figure 4 illustrates the beam outage probability with respect to the transmit SNR of HAP for the different values of central angles δ .Similar to the previous analysis, with the increase of the HAP's transmit SNR, the beam outage performance can be improved.It can be seen from the figure that the beam outage probability is decreased with a small value of central angle δ , which is due to the fact that the sensors can obtain larger chan- nel gains in a narrower sector as they are closer to the normalized direction δ of HAP.Furthermore, the exact analytical results are well aligned with the simulation results.In Fig. 5, we investigate the beam outage probability versus the number of antennas at the HAP for different value of δ .We can further insight that the central angle δ significantly affects the outcome of beam outage probability by increasing the num- ber of HAP's antennas.For small δ , the beam outage probability is first decreased and then increased to a certain value with an increase of M. This is due to fact that when the number of antennas is increased, the beam formed by the antenna can adequately  cover the sensors in the service sector.However, if the number of antennas continues to increase, even if the generated beamforming can provide more energy to the sensors, the beamwidth will become narrow and can only serve sensors close to the normalized direction.Therefore, the beam outage probability shows an upward trend.Moreover, a larger center angle implies that the HAP needs to serve a larger area.In this scenario, increasing the number of antennas will reduce the coverage of the HAP, which results the system's beam outage probability further increasing.For a service sector with large central angle in mmWave communication systems, the more the antennas, the narrower the beamwidth, which will reduce the number of served sensors and beam outage performance.In contrast, a small number of antennas can achieve a relatively better beam outage performance.The analytical results also coincide exactly with simulation results.
Figure 6 demonstrates the beam outage probability with respect to WEH sub-interval α 0 .From the general trend, with the increase of α 0 , the system's beam outage probability first decreases and then increases.This is due to the sensors being able to harvest more RF energy with the increase of α 0 , which can improve the data transmission rate and beam outage probability to some extent.If α 0 is too long, the sub-interval for data trans- mission will be reduced and the beam outage performance would suffer significantly.As expected, a higher HAP's transmit power can deliver more RF energy for the sensors at a given WEH sub-interval, thus reducing the beam outage probability of the system.Meanwhile, the beam outage performance of the system decreases with the increase of the density , which because that the TDMA scheme is utilized so that each sensor will be allocated less WIT sub-interval as the number of sensors increases.Similarly, the numerical results agree well with the simulation results.
Figure 7 plots the system's beam outage probability versus target rates for different central angles δ .As the target rate increases, the probability that the HAP can suc- cessfully decode data will be decreased, so that the beam outage probability increases.As expected, within a smaller coverage range, the sensors will be located closer to the center angle, so that the sensors can receive energy and transmit data more effectively.
As can be seen, the simulation points are perfectly aligned with the corresponding analytical curves.

Energy efficiency analysis
In this subsection, we demonstrate the energy efficiency (EE) of the system and the sum achievable rates with proposed algorithm, which are also compared to the ergodic EE and achievable sum-rate in [45].Figure 8 depicts the changes of system's EE and sum achievable rate versus HAP's transmit SNR for different numbers of sensors.As can be seen from the figure, with the increase of HAP's transmit SNR, the sum achievable rate of the system gradually increases, while the EE is gradually decreased.This is because although the sum achievable rate of the sensors in uplink transmission is improved, the correspondingly consumed energy consumed is also increased, which ultimately results in a decrease of the system's EE.In the meanwhile, with the increase of the number of sensors, the considered HTT transmission protocol can provide WEH and WIT services for more sensors, but the timeslot allocations of energy harvesting and the data transmission will be decreased, thus reducing both the EE and sum achievable rate reachability of the system.In addition, compared with the ergodic performance analysis of the multi-user uplink transmission mmWave WPCN system considered in [45], both the EE and sum achievable rate of the system can be concurrently improved by proposed jointly optimizing the uplink and downlink transmission timeslot allocations.For example, it can be calculated from the figure that when K = 7, compared with the comparison scheme [45], the proposed optimal timeslot allocation algorithm can increase the sum achievable rate and EE of the system by up to 25% and 26.4%, respectively.
In Fig. 9, the EE and sum achievable rates are evaluated versus HAP's transmit SNR for different center angles.Similar to the analysis of Fig. 8, with the increase of HAP's transmit SNR, the sum achievable rate of the system is increased, while the EE is decreased.The figure also shows that the value of center angle δ has a great impact on the system's performance.When δ is small, the sensors are more concentrated in the near of the Fig. 8 Energy efficiency and sum achievable rate versus HAP's transmit SNR for different numbers of sensors HAP's normalized direction of energy beam, so that the sensors in the service sector can harvest more energy during WEH phase for future WIT, which can promote the system sum achievable rate.Likewise, compared with ergodic performance analysis in [45], both the sum achievable rate and EE based on the proposed timeslot allocation algorithm can be maximized by 12% and 13.8%, respectively, when taking δ=0.1 as an example.

Conclusion and discussion
In this paper, we investigate a typical wireless-powered mmWave sensor network, where a transmission interval is divided in to two sub-intervals.For the first sub-interval, a hybrid access point utilizes energy beamforming to simultaneously power the sensors within the service sector.During the remaining sub-interval, sensors transmit their data to HAP successively by using the harvested energy.According to the stochastic geometry theory, the exact expression of the system's beam outage probability is derived.The approximately closed-form expressions of beam outage probability are obtained based on some reasonable assumptions.Then, in order to maximize the system's energy efficiency, an optimal time allocation algorithm is proposed.The optimization problem can be solved by a proposed Dinkelbach algorithm, which includes parametric form transforming and KKT conditions.Extensive simulation results reveal that system's parameters significantly affect the beam outage probability, while our proposed algorithm can improve the system's energy efficiency and sum achievable rate compared to the corresponding ergodic analyses.
On the basis of this paper, in the future, we continue to carry out further in-depth research, where a mmWave-based wireless-powered communication network will be considered with combination of TDMA and FDMA operations.In order to improve the energy efficiency of the system, a joint optimization problem of multiple resource allocation is investigated and a two-stage design method is proposed to obtain optimal

Appendix
Proof of Lemma 1 The CDF of achievable rate for a typical sensor at the HAP can be written as Pr {ν < r} .In addition, we assume that there is reciprocity between channels h k and h ′ k [46].Let X = g 2 , which is exponential-distributed random variable with PDF g X (x) shown as follows where denotes channel average power.Substituting ( 8) into ( 9), we thus have where (a) follows the fact that g 2 ≥ 0 .Then, (30) can be further written as where A = (1 − α 0 )σ 2 0 2 rK / (1−α 0 ) − 1 (K α 0 ηP HAP ) , (b) follows from the fact that g 2 is exponential distribution, (c) follows by considering the spatial distribution over all possible distances and angles.Moreover, the K sensors are deployed in service sector following HPPP and they are independently and identically distributed (i.i.d) points.Therefore, the PDF of distance for sensors can be obtained as By substituting (31) into (30), the result of (30) can be rewritten as (29)

Fig. 3
Fig. 3 Beam outage probability versus HAP's transmit SNR for different target rates

Fig. 4
Fig. 4 Beam outage probability versus HAP's transmits SNR for different central angles

Fig. 6 Fig. 7
Fig. 6 Beam outage probability versus WEH sub-interval for different HAP's transmit SNRs and various density of HPPP

Fig. 9
Fig.9 Energy efficiency and sum achievable rate versus HAP's transmit SNR for different central angles

Table 1
Lists of necessary parameters