Age of Information in Wireless Powered Networks in Low SNR Region for Future 5G

Wireless powered communication technology has a great potential to power low-power wireless sensor networks and Internet of Things (IoT) for real-time applications in future 5G networks, where age of information (AoI) plays a very important performance metric. This paper studies the system average AoI of a wireless powered network, where a wireless-powered user harvests energy from a wireless power source (WPS) and then transmits data packets to its access point (AP) by using the harvested energy. The user generates data packets with some probability and adopts the first-come-first-served (FCFS) service policy. For such a system, by using the queuing theory and the probability models, we derive a closed-form expression of the system average AoI. We also formulate an optimization problem to minimize the AoI by optimizing the data packet generating probability, and find its solution by simple calculation and search. Simulation results demonstrate the correctness of our obtained analytical results. It also shows that, when the total distance of the two hops is fixed, the system average AoI increases linearly with the increment of the distance of the first hop, and a smaller data packet generating probability should be selected to match a bigger first-hop distance for achieving a smaller system average AoI. Moreover, a smaller data packet size also contributes to a smaller system average AoI.


Introduction
Recently, the widespread use of mobile devices and applications has made the real-time information updating applications such as news, weather forecasting and traffic alert more and more popular [1][2][3]. Timely information updating is also becoming more and more critical for real-time monitoring and control systems, including wireless sensor networks (WSNs) and internet of Things (IoT) for temperature and humidity detection in warehouses [4], safety and state monitoring in industrial production lines, embedded equipments in medical care [5], and road condition detection in automatic droving for future 5G systems [6,7]. The common key point of above-mentioned real-time applications is how to guarantee the freshness of the collected data.
Traditionally, delay and throughput are two important and widely adopted indices to evaluate the system performance of communication networks [8,9]. However, they are insufficient to describe the freshness of the data [10]. Therefore, a new metric, referred to as age of information (AoI), has emerged recently, which is defined to characterize the elapsed time since the last received data was generated [11]. The AoI actually describes the freshness of status updates based on time-varying wireless channel transmissions because it can reflect real world constraints condition that the delivery of a status message requires a nonzero and typically random time in the system [12].
In the past several years, AoI has been analyzed and studied in various queuing systems. For instance, the AoI was analyzed for single-source single-server queues in [11], and for M/M/1 first-come-first-served (FCFS) systems with multiple sources in [13]. In [14], AoI was explored for a multi-class M/G/1 queueing system. In [15], a packet deadline was regarded as a control mechanism to study its impact on the average AoI for M/M/1 queueing systems, and in [16], it was found that the packet waiting in queue was replaced if a new packet arrives service discipline is optimal.
Meanwhile, in order to realize the fresh-data transmission and explore the system performance limit, AoI was minimized for single-hop networks, see e.g., [11,14,17,18], and multi-hop networks, see e.g., [19,20], with different system setups. In [19], it was shown that for general system settings (including arbitrary network topology, packet generation times, packet arrival times, and queue buffer sizes), age-optimality can be achieved. In [20], an energy harvesting (EH) assisted two-hop system was studied, where an EH-enabled source collected measurements from a physical phenomenon and then sent updates to a destination with the help of an EH-enabled relay.
As for EH, it is able to power the device by harvesting energy from external environment [21] and has a great potential to be applied to energy-constrained networks including IoT and WSNs in future 5G [22][23][24][25][26][27][28]. EH technologies can be classified into two categories: the natural energy source-based EH and the radio frequency (RF) signal-based EH [29][30][31][32][33][34]. Compared with natural source, RF signal is easier to control and has less environmental limitations to deploy. Therefore, it is becoming more and more popular. Since RF EH-enabled IoT and sensor networks are expected to be widely employed for industrial control, unmanned driving systems and real-time applications, analyzing AoI performance and designing AoI minimized systems for RF-EH wireless networks have attracted increasing interest, see e.g., [35,36]. In [35], the average AoI of the two-way networks, where the slave node is powered by RF EH supply and the uplink average AoI was analyzed. In [36], the authors studied how to design optimal online status update policies to minimize the long-term average AoI, subject to the energy causality constraint at the sensor.
In this paper, we study the average AoI in a wireless powered network, where a wireless-powered user harvests energy from a wireless power source (WPS) and transmits data packets to its access point (AP). The user generates a data packet with some probability p in each time block, and the generated data are stored in an infinite buffer queue. FCFS service policy is employed, using the queuing theory and the probability model. We derive a closed-form expression of the average AoI for such a system. Numerical results are provided to discuss the system AoI performance. It is observed that there is an optimal p such that the AoI reaches minimal. Therefore, we formulate an optimal problem to minimize the AoI to find the optimal p. Since the problem is non-convex, we use the one-dimensional search to find the optimal p. Moreover, with other parameters being fixed, we find that the minimal average AoI linearly increase with the increment of packet size. In addition, we also analyze the effects of other factors on the average AoI. It is found that, when the total distance of the two hops is fixed, the system average AoI increases linearly with the increment of the distance of the first hop, and a smaller data packet generating probability should be selected to match a bigger first-hop distance for achieving a smaller system average AoI. Moreover, a smaller data packet size also contributes to a smaller system average AoI.
The rest of the paper is organized as follows. In Section 2, we present the system model, including the channel model, the data transmission model and the energy transfer model. In Section 3, we derive an explicit expression of the average AoI. Section 4 discusses the system AoI performance via simulations. Section 5 summarizes the paper with some conclusions.

System Model
We consider a wireless powered network, as depicted in Figure 1, where a user (e.g., a sensor node) desires to transmit its data to its AP (e.g., a sink node). Since the user is lack of energy, it has to harvest energy from a WPS that deployed in the system and used to charge the wireless devices via wireless power transfer. Our considered network model can be regarded as a basic component of complex networks. For example, when multiple nodes are deployed, by using time division multiple access(TDMA) or frequency division multiple access (FDMA), the complex network can be decomposed into multiple point-to-point networks [37][38][39], each of which is like our considered model. It is assumed that the system works in a discrete time manner. That is, the time is divided into blocks with equal interval T b . The time period from the epoch n to the epoch n + 1 is referred to as the time block n. Block Rayleigh fading channel model from WPS to the user and from the user to its AP is assumed, so that the channel coefficient can be considered as a constant in each time block, and it may vary from one time block to the next for information transmission. Let h 1 [n] and h 2 [n] be the channel coefficients of the links from the WPS to the user and from the user to the AP, respectively, associated with time block n. The corresponding power gains |h 1 | 2 and |h 2 | 2 follow the exponential distribution, which can be expressed by and where λ 1 and λ 2 are the exponential distribution parameters. In each block, the user generates a data packet with size of δ bits randomly with a certain probability p, and the generated data packets are first stored in a data buffer, and then transmitted with the FCFS policy. Denote the distance between the WPS and the user and between the user and the AP to be d 1 and d 2 , respectively. Let the transmit power of the WPS be P w . If in time block n, the user performs energy harvesting, the energy received at the user from the WPS in time block n is given by where η is the energy transfer efficiency and µ = η . α is the pass loss factor. Since the energy transfer efficiency η is less than one and the received power is also relatively small due to pass loss, the user may take several blocks to harvest and accumulate energy to complete a block of transmission. The energy accumulated at the user within j blocks is According to Equation (1), |h 1 | 2 follows exponential distribution with parameter λ 1 , so e j follows Erlang(j, λ 1 µP w T b ) distribution, i.e., where . It is assumed that the energy transmission and data transmission are over orthogonal frequency bands. Let the transmit power of the user be P u . Assume that the received signals are suffered from additive white Gaussian noise. If in time block n, the user performs data transmission, the data size can be delivered in block n is (see Endnote [40]-which refers to Reference [35,41]) where B is the system bandwidth and N 0 = n 0 d α 2 with n 0 denoting the noise spectral density. For such a system, our goal is to analyze its average AoI performance in fading channels.

AOI Analysis
In IoT and sensor networks, the transmit power of the devices is usually very low, so that the received energy from the transmitted signals of WPS at the user is relatively very small. Therefore, we analyze the system AoI performance in the low SNR regime, where c[n] can be approximated by

Preliminary Analysis
Denote the number of time blocks required to complete the transmission of one data packet with size of δ to be N I . Define the probability of successfully completing a packet transmission with j time blocks to be p n i j . When the user accumulates sufficient energy, we obtain the following theoretical results. Proposition 1. The probability of successfully completing a packet transmission with j time blocks is where j = 1, 2, 3, · · · , n and θ = λ 2 N 0 δ Proof. According to the definition of p n i j and [42], it can be expressed that Since |h 2 | 2 follows exponential distribution, from Equation (6), c[n] obeys exponential distribution. Thus, the cumulative distribution function (CDF) of c[n] is . Therefore, we arrive at Proposition 1.

Proposition 2.
The probability generating function (PGF) and the expectation of N I are and respectively.
Proof. In terms of the definition of PGF, we have By the property of PGF, we have In our model, the harvested energy is allowed to be used for information transmission in the same block. As long as the accumulated energy at the user is greater than or equal to P u T b , it performed a block of transmission. Suppose the number of blocks that need to harvest energy for the user to complete a block of transmission is N E . Since the energy transfer efficiency η is less than one and the received power is also relatively small due to pass loss, the user may take several blocks to harvest energy to complete a block of transmission, i.e., N E ≥ 1. As a result, the transmitter start information delivering at time block j, if Therefore, the actual time required to complete a block of transmission is N E = j. Without loss of generality, we suppose P w = βP u , with β > 0. In terms of Equations (4) and (5), the probability Pr{N E = j} of N E = j can be given by Let the service time of the kth data packet is T (S) k , which represents the time for harvesting energy and the time for transmitting data. Lemma 1. The probability that j time blocks are needed to actually transfer a data packet is Proof. When the user accumulates sufficient energy, in terms of Equation (9), the average time to transmit a data packet is 1 + θ. Since the user needs N E time blocks to harvest energy to perform a block of transmission, that is, T (S) k contains 1 + θ independent and identically distributed variables N E (i), according to Equation (10), one can derive Equation (11).

Proposition 3. The probability generating function (PGF) of T
and the expectation of T Proof. In terms of the definition of PGF, we have

Description of System Average AoI
The AoI is defined as the time elapsed since the last received packet was generated [11], which is used to measure the freshness of information at the destination. We assume that the most recently received data packet at AP in block n was generated at the time U[n]. Therefore, the AoI in block n, can be expressed by Figure 2 illustrates a sample evolution of AoI versus time blocks with initial age ∆ 0 , i.e., ∆[0] = ∆ 0 . Since time is discrete, the AoI is constant within each block and varies from a block to next block. Let n k denote the generation time of the kth data packet, and n k denote the time when the kth data packet is completely transmitted. The number of blocks from the generation of a data packet to the completion of the transmission is T which is called the system time. The interval time between the generation of data packet k − 1 and data packet k is set as I k . That is, Let T It is observed that the AoI increases linearly in time and is reset to a smaller value when a data packet is received. That is, at n k , the AoI is reset to ∆[n] = n k − n k . Over a period of N blocks where K data packets are delivered, the average AoI is defined as As illustrated in Figure 2, the average AoI of system can be calculated as the average area of the blue graphic Q k , i.e.,∆ = lim where K + 1 . In the following, we shall discuss how to derive a closed-form exprtession of the average AoI∆.

Closed-Form Expression of∆
From Figure 2, one can see that the area of Q 0 and R K are limited, with large enough N, Q 0 N and R k N are close to zero. Therefore,∆ can be approximatively give bȳ Moreover, one can see that the area of Q k is the difference of the area of the large red triangle minus the area of the small black triangle, as illustrated in Figure 3. That is, Now, we begin to drive the explicit expressions of E (I k ), E I 2 k and E I k T (T) K as follows.

Expressions of E (I k ) and E I 2 k
Since p is the generation rate of the data packet at the user, i.e., data generating with probability p in each block. Therefore, the inter arrival time I k follows geometric distribution. I k = j means that in the j-th block data packet is successfully generated, but was not generated in the previous j − 1 consecutive blocks. Thus, As a result, and

Expressions of E I k T (T) K
Let Y k be the number of data packets generated during the system time of the kth data packet. It is also a random variable. To describe it, we define {Y k } as an independent and identically distributed random sequence. The probability of generating j data packets during the service time of a data packet can be given by p where j = 0, 1, 2, · · · and ∑ ∞ j=0 p y j = 1. Let L k be the number of data packets in the queue before servicing the kth data packet. Thus, we have By regarding {L k } as a state of the system, it can be expressed by a Markov chain. Its state transition diagram is shown in Figure 4, where the numbers in the figure represents the different status (i.e., value of L k ), and the probability on the curve represents the transition probability from one state to the next. The transition probability is obtained as follows:

Proposition 4.
The probability distribution and the PGF of Y are respectively given by, and Proof. See Appendix A.
Proposition 6. The PGF of the system time of each data packet is given by Proof. See Appendix C.

Lemma 2. The average system time is
Proof. By the property of PGF, we can obtain Recall Equation ( * ), i.e., the system time of data packet k. Note that, when the kth data packet is generated at the user, if the (k − 1)th data packet has completed the service, T ). (27) Proposition 7. For the considered system, Proof. See Appendix D.

Proof. Combining Equations
Furthermore, from Theorem 1, we can get the following two corollaries.
Proof. Let p = , where > 0 is a small enough number. Equation (29) can be approximately written as Since we arrive at Corollary 1.

The Minimal Average AoI of the System
Our objective is to minimize the average AoI. The maximum data rate that can make the data packets queue stable is . Hence, we can formulate the problem as follows: Since it is difficult to theoretically prove that the problem 1 is a convex or non-convex, we cannot directly obtain the optimal solution. Thus, we use the one-dimensional search to find the optimal p, and give simulation results to illustrate that there exists an optimal solution for Problem 1 in the next section (see Endnote [44]).

Fact 1. The average AoI converges to a constant with the increment of P u .
Although it is difficult to prove Fact 1 and the converged constant is also hard to derive mathematically, we may analyze it from the physical perspective. It is known that with the increment of P u , the information rate of the second hop is increased, which may decrease the system AoI. However, due to the Shannon capacity theorem, the information rate of the second hop cannot be increased infinitely, so the system AoI cannot be decreased infinitely, which must converge to some value with the increment of P u .

Numerical Results
In this section, numerical simulations are conducted to discuss the system AoI performance. The simulation parameters are set according to [35]. Specifically, the distances between the WPS and the user and between the user and the AP, are set to unit distance, i.e., d 1 = d 2 = 1 meter. The pass loss factor is α = 2. The transmit power of the user is set as P u = 0.01 W, and β = 1, i.e., P w = P u . The size of data packet is δ = 8 bits. The energy transmission efficiency is η = 0.6, the system bandwidth is B = 1 MHz, the block length is set as T b = 10 −3 s, and the noise spectral density is n 0 = 4 × 10 −7 . The parameters of the channel power gain |h 1 | 2 and |h 2 | 2 are set as λ 1 = 3 and λ 2 = 3, respectively.
As shown in Figure 5, the analytical results match the simulation ones well, which validates the correctness of our analytical results. Moreover, when p is close to 0, the average AoI is very large because the interval of data generation I k is very large in this case. When p approaches p max , the average AoI also increases to be infinite because, in this case, the waiting time in the queue becomes very long. This observation is consistent with Corollary 1 and 2. It is also observed that there exists a unique optimal p such that the system average AoI achieves a minimum.q To further discuss the system performance, we fix D = d 1 + d 2 to be two meters and then increase d 1 gradually, as illustrated in Figure 6. The simulation results are plotted in a 3D figure as shown in Figure 7. One can observe that with, the increment of d 1 , the system average AoI increases linearly, and the smaller p should be selected for the bigger d 1 in order to achieve the system minimal average AoI. Figure 8 plots the system average AoI versus data generation probability p and data size δ in a 3D figure. One can see that the minimal average AoI increases linearly with the increment of packet size δ, and the smaller p should be selected for the larger δ. In order to clearly show this, for each given data size δ, we mark the corresponding minimum AoI and it is plotted versus δ in Figure 9, which can be considered as the lower bound of the system average AoI. This observation also implies that, to keep the data fresh, the smaller packet size is preferred.  Figure 10 shows the system minimal average AoI versus P u with different configurations of T b . It is shown that, for a given block length T b , the minimum AoI decreases with increment of P u and then tends to be flat. That is, the average AoI converges to a constant with the increment of P u , which is consistent with the result in Fact 1. Moreover, the larger the the block length T b is, the smaller the minimal AoI of the system is.

Conclusions
In this paper, we studied the system average AoI of a wireless powered communication network, where a wireless powered user harvests energy from a WPS and then transmits data packets to its AP by using the harvested energy. By using the queuing theory and some typical probability models on the channel fading, we derived a closed-form expression of the system average AoI. Some interesting results are obtained, and there exists an optimal generating probability of sensor data (user) p such that the AoI reaches a minimum. Simulation results also show that, when the total distance of the two hops is fixed, the system AoI increases linearly with the increment of the distance of the first hop, and a smaller data packet generating probability should be selected to match a bigger first-hop distance for achieving a smaller system average AoI. Moreover, a smaller data packet size of users will result in contributing to a smaller system average AoI.

Abbreviations
The following abbreviations are used in this manuscript: For the case Y k = j ≥ 1, we have In terms of the definition of PGF, the PGF of Y k is
The two sides of the equation are multiplied by z j , we can obtain where G Y (z) is given by Equation (A3). In order to calculate G L k (z), Equation (A4) can be expressed as According to the properties of the PGF, we have lim z→1 G L k (z) = π 0 1 − p − (1+θ)λ 1 βµ p = 1.

Appendix C
Proof. In the first-come first-served queue, the queue length at the packet departure time being transmitted is equal to the number of arriving data packets during the system time of the leaving data packet. Thus, where j = 0, 1, 2, · · · , andj ≤ i.
According to the definition of PGF, we have Combining formulas G L k (z) in Proposition 5 and G L k (s) = G T (1 − p − ps) in Equation (A6), and replacing (1 − p − ps) with z, we can get G T (z).

Appendix D
Proof. In order to calculate E(I k T (W) k ), we define an auxiliary function as where the last two { . =} obtained by Lemma 3 [42]. G T (z) and E(T (T) k−1 ) are given by Equations (25) and (26). Thus, we have where G T ((1 − p)z) and G T ((1 − p)z) can be calculated by Equation (25).