Social-Aware Peer Discovery for Energy Harvesting-Based Device-To-Device Communications

In Device-to-Device (D2D) communications, the first step is to find all of the neighboring peers in the network by performing a peer discovery process. Most previous studies use the social behaviors of the users to adjust the sending rates of the peer discovery messages (i.e., beacons) under the constraint of consumed power for increasing the Peer Discovery Ratio (PDR). However, these studies do not consider the potential for energy harvesting, which allows for the User Equipments (UEs) to procure additional power within charging areas. Accordingly, this paper proposes an Energy-Ratio Rate Decision (ERRD) algorithm that comprises three steps, namely Social Ratio Allocation (SRA), Energy Ratio Allocation (ERA), and Beacon Rate Decision (BRD). The SRA step determines the allocated power quantum for each UE from the total budget power based on the social behavior of the UE. The ERA step then adjusts this allocated power quantum in accordance with the power that is harvested by the UE. Finally, the BRD step computes the beacon rate for the UE based on the adjusted power quantum. The simulation results show that ERRD outperforms the previously-reported Social-Based Grouping (SBG) algorithm by 190% on the PDR for a budget power of one watt and 8% for a budget power of 20 watts.


Introduction
The demands that are placed on wireless communications have exponentially increased in recent years due to the proliferation of User Equipments (UEs) and the unceasing development of new mobile services. Consequently, the lack of spectrum resources has emerged as a significant concern for communication operators. Currently, any UE that wishes to transmit data to another UE must transmit this data via the Base Station (BS) in the cellular network. However, if the distance between the two UEs is sufficiently small, then the potential exists for the UEs to communicate directly, thereby saving cellular spectrum resources. Accordingly, a new communication paradigm designated as Device-to-Device (D2D) communication has been proposed as a means of achieving short-distance transmissions in 5G networks with improved resource utilization efficiency.
In implementing D2D communications, the first, and most critical, step is that of peer discovery, in which the UEs attempt to identify all of their neighbors in the network [1]. The peer discovery process is generally performed using a beacon mechanism and it aims to maximize the Peer Discovery Ratio (PDR), which is defined as the number of peers that are found in the discovery process divided is the first paper to consider the social-aware peer discovery problem for D2D communications in an energy-harvesting environment, (2) regarding the solution orientation, we propose ERRD, which determines the beacon rate according to UE's social behavior and the harvested power quantum to increase its PDR, and (3) regarding the evaluation orientation, some of the simulations are conducted to investigate many important parameters along with significant observations. The remainder of this paper is organized, as follows. Section 2 introduces the related work in the field. Section 3 describes the system model that is considered in the present study and formulates the related problem. Section 4 introduces the ERRD algorithm and describes its detailed operation. Section 5 presents and discusses the simulation results. Finally, Section 6 provides some brief concluding remarks and indicates the intended direction of future research.

Related Works
The problem of peer discovery in D2D communications has attracted considerable attention in the literature. Broadly speaking, existing approaches can be classified as either social-aware or social-unaware (see Table 1). Within each classification, the proposed mechanisms can be further divided as autonomous, network-assisted, or network-controlled. In autonomous mechanisms, the UEs find nearby peers by themselves, and hence the major issue lies in determining efficient methods for broadcasting the beacons with a minimum number of collisions. By contrast, in network-assisted solutions, the BS collects surrounding the information, determines a suitable beacon rate for each UE, and allocates a proper amount of resources to each UE to perform its transmissions. Finally, in network-controlled mechanisms, the BS directly helps the UEs to find peers by locating UEs.

Category
Ref.

Type Main Goal Approach
Social-unaware [2] Autonomous Min resource UEs periodically and synchronously send beacons using FDM. [3] Autonomous Min resource UEs use trellis tone modulation multiple-access scheme. [4] Autonomous Quick recovery UEs use a common channel and a group of channels to send beacons. [5] Autonomous LTE compatible UEs listen to SRS channel to identify nearby UEs. [6] Autonomous LTE compatible UEs listen to SRS channel to detect active UEs. [7] Autonomous Min energy UEs determine the beacon rate based on their state. [8] Network-Assisted Min collision BS determines the beacon rate for UEs depending on the number of requests sent by UE. [9] Network-Assisted Min collision UEs authorized to perform discovery contend to transmit beacons.
[10] Network-Assisted Min interference Similar to [9], but interference from cellular UEs imposed on D2D pairs is also considered. [11] Network-Assisted Min resource UEs send preamble to nearby UEs and BS allocates uplink RBs for UEs. [12] Network-Assisted Beacon schedule BS roughly estimates the location of UEs by measuring channel components. [13] Network-Controlled Min resource BSs locate UEs by AOA. [14] Network-Controlled Min energy Wifi scans are first used to determine the UE locations and BS then sends D2D broadcast.
Social-aware [17] Autonomous Trust UEs find trusted UEs. [18] Autonomous Two-hop pairing UEs send requests to trusted UEs, which forward request to all one-hop UEs. [19] Autonomous Hybrid attributes UEs use three key social attributes to construct neighbor lists.
[20] Autonomous Max content delivery Social relationship is used as a weight for D2D pair formation and content sharing. [21] Network-Assisted Max PDR BS determines the beacon rate based on contact rate. The literature contains several autonomous mechanisms for peer discovery using a social-unaware approach [2][3][4][5][6][7]. For example, the FlashLinQ protocol that was proposed in [2] uses frequency division multiplexing (FDM) to propagate the beacons through the network [2]. Notably, the beacons are transmitted both periodically and synchronously, and hence FlashLinQ provides an effective means of estimating the amount of consumed resources and timing the resource consumption to minimize resource snatching. The scheme that is presented in [3] replaces the FDM-based peer discovery process in [2], with a non-orthogonal multiple-access scheme that is referred to as Trellis Tone Modulation Multiple-Access (TTMMA). TTMMA uses single-tone transmissions and achieves long-distance discovery due to its low Peak-to-Average Power Ratio (PAPR). Furthermore, it makes a higher discovery capacity through its use of a non-orthogonal resource assignment mechanism possible. Based on the assumption of a synchronized superframe structure among the UEs, the peer discovery mechanism in [4] uses both a common channel and a group of channels to send beacons subject to certain rules and procedures that are designed to minimize the discovery time. To ensure compatibility with the standard LTE protocol, the UEs contained in [5] perform peer discovery by listening to the Sounding Reference Signal (SRS) channel, being originally designed for data uploading purposes in LTE networks and peer UEs can access it. However, the method is only capable of detecting the presence of peers, i.e., not active peers that are interested in D2D discovery and communications. Accordingly, the mechanism in [6] exploits the standardized uplink signal structure in the SRS channel to accomplish both the detection of the active peers and the identification of their beacons. Finally, in the peer discovery method that was proposed in [7], the UEs reside in one of five different states, namely, Keep Alive, Advertise, Discover, Follow, and Passive, depending on their behavior. The UEs in each state then employ a particular beacon rate that is chosen in advance in such a way as to maximize the power saving in the network.
The literature contains various proposals for network-assisted peer discovery methods [8][9][10][11][12]. In [8], the UEs advertise their presence using a random access mechanism and the BS accepts the D2D requests, allocates resource blocks (RB), and chooses an appropriate beacon rate for each UE, depending on the number of requests that they produce. In the peer discovery method that is proposed in [9], the UEs authorized to perform discovery in a given discovery interval contend to transmit their beacons in a time-frequency multiplexed pool of network-allocated resources. In [9], it is assumed that the transmitted beacons are always successful, i.e., other peers can successfully receive them. However, the method in [10] considers a more realistic network model that is based on the Poisson Point Process (PPP), in which the effect of the Channel State Information (CSI) on the performance of the D2D discovery process is taken into account by considering the interference that is imposed on the D2D pairs by nearby cellular users. In [11], each UE sends a preamble to the nearby UEs via a newly-introduced physical channel, and the UEs that receive this preamble send a corresponding report message to the BS by means of a normal random access procedure. The BS then allocates an uplink RB for each reported preamble, such that the UEs that initially sent the preamble can send a further report message to the BS. Finally, the BS, by comparing their reported preambles, identifies pairs of UEs that are in close proximity to one another. The method in [12] exploits the fact that some of the channel components of the UEs are spatially correlated to enable the BS to make a rough estimate of the UE locations by measuring these components. The BS then schedules the transmissions of the beacons that are sent by nearby UEs, such that nearby UEs transmit their beacons at similar times.
Network-controlled methods have several important advantages over Autonomous and Network-Assisted methods, including low power consumption, reduced interference, and a low beacon transmission cost. Consequently, the authors in [13] proposed a peer discovery method based on the Angle of Arrival (AOA) measurements that were obtained by multiple BSs and further analyzed the performance of network-controlled D2D discovery in random spatial networks. The authors in [14] proposed a centralized novel approach, called ROOMMATEs, which utilizes the ubiquitous WiFi network, which combines with BS for indoor peer discovery. ROOMMATEs is an unsupervised approach that can provide different granularity location information. However, none of the studies in [2][3][4][5][6][7][8][9][10][11][12][13][14] consider the potential for improving the peer discovery performance by exploiting the social behaviors of the UEs in the network.
In fact, the literature contains very few proposals for social-aware peer discovery mechanisms [17][18][19][20][21]. Among those methods that have been proposed, three schemes [17][18][19] adopt an autonomous approach. The method in [17] focuses on the security of the data transmissions and it chooses UEs with high social ties (i.e., high trust) to perform D2D communications. However, by adopting such an approach, it is possible that no UEs may be available for pairing. Consequently, the authors in [18] proposed a two-hop pairing process, in which any UE failing to find a trusted UE with its wanted contents for D2D communications sends a request to all the one-hop neighbors of its trusted UEs. The method in [19] uses three key social attributes, namely the trust degree, the similarity degree between UEs, and the center degree of each UE, to construct a unified metric with which to construct neighbor lists for peer discovery. Reference [20] addresses the content delivery problem that is related to optimization of peer discovery and resource allocation by combining both the social and physical layer information in D2D networks. The social relationship, which is modeled as the probability of selecting similar contents, is used as a weight to characterize the impact of social features on D2D pair formation and content sharing. The peer discovery in the physical layer depends on the received signal power by UEs, while the social relationship is mainly used to optimize the context delivery. The method in [21], which is called Social-Based Grouping (SBG), adopts a network-assisted approach to perform peer discovery. UEs are grouped based on their social feature and centrality, and the UEs in the same group are assigned the same beacon rate. The BS, in accordance with the rate at which each UE contacts the other UEs in the network, determines the beacon rate. The UE with more contacts will send beacons at a higher rate.
Similar to SBG, our study also focuses on the network-assisted approach to perform peer discovery. However, between [21] and our study, there are some big differences: (1) SBG determines the beacon rate with only considering the social feature, while ERRD determines it with not only considering the social feature, but also the harvested power. (2) The UEs in SBG are grouped and the UEs in the same group are assigned the same beacon rate due to its high complexity, while each UE in ERRD has its individual beacon rate due to its low complexity. (3) Since ERRD extra considers the harvested power to determine more suitable beacon rate, its PDR can be significantly better than that of SBG. This can be easily observed in Section 5.

Problem Description
This section commences by introducing the system model and associated notations. The problem statement is then formally defined.

System Model
The interval between the time t at which UE i comes within range of UE j and the time at which it was last within range of UE j, t 0 , is referred to as the D2D contact interval of the two UEs, CI i,j , and it is defined as where . denotes the distance between the two UEs; L i and L j are the locations of UEs i and j, respectively; and, R i,j is the coverage range between the two UEs. Based on the contact interval between the two UEs, the D2D contact rate between them, denoted as λ i,j , is defined as where E[.] denotes the expectation. The average contact rate of UE i, denoted as λ i , is then computed as where N is the total number of UEs in the network. Let CI i,j follow a cumulative distribution function (CDF) of F CI i,j (x), with rate λ i,j . Assume further that F CI i,j (x) is a uniform distribution, and can thus be written as Let the social ratio of UE i be defined as the square root of its average contact rate over the sum of the square root of each average contact rate, i.e., In modeling the energy-harvesting environment, it is assumed that the charging devices (CDs) and BS convert their power into RF signals. Furthermore, any UEs within the coverage of these CDs or the BS acquire this RF signal and then convert it into power through special equipment. The power that is obtained by UE i from energy harvesting, denoted as OP i , is thus calculated as where σ is an energy harvesting efficiency factor that reflects the ability of UE i to change the RF signal into power; v is the path loss exponent that is caused by interference, which increases with an increasing distance; and, L i is the location of UE i. In addition, CDP m , CDR m , and CDL m are the transmitted power, coverage, and location, respectively, of CD m, when UE i lies within its coverage [30][31][32].

Problem Statement
Before formally defining the problem statement, the used notations are listed in Table 2. As shown, the notations fall into six categories that relate to the system, range, power, location, contact, and beacon, respectively. The system parameters define the number of UEs in the network and the total peer discovery time, respectively, while the range parameters describe the coverage of the devices. The power parameters define the transmission powers of the devices and the various power variables that are used in the ERRD model. The location parameters define the positions of the devices. The contact parameters describe the contact behaviors of the UEs. Finally, the beacon parameters define the beacon rates of the UEs. The energy harvesting technology that is considered in this study is wireless power transfer (WPT) via radio frequency (RF), as shown in Figure 1 [23]. The basic idea is that the RF transmitter transmits radio signal towards the receiving antenna at the desired frequency and power level. The RF receiver then applies electromagnetic radiation to charge the battery, that is, the receiving antenna receives the traveling signal and the rectifier converts the alternating current (AC) to direct current (DC) to charge the battery. The energy carrier can be located at 2.4 GHz or 5 GHz frequency band under the considerations that these bands are internationally reserved for Industrial, Scientific, and Medical (ISM) purposes. Energy harvesting can explore sufficient power sources in significant radio coverage by increasing the transmitting power under the regulation of the government.  This study attempts to determine the suitable beacon rates of all UEs to explore most peers in a D2D network under the constraints of budget power, the limitation of overall power consumption. This study further introduces a new scenario, i.e., energy harvesting, which previous studies did not consider. The energy harvesting enables UEs within charging areas to procure additional power from surrounding energy sources. Thus, under considering the UE may be out-of-energy and the study determines the suitable beacon rates of all UEs to maximize PDR according to the budget power and harvested power.
The problem that we investigated is formally described, as follows. First, for a beacon rate of UE i equal to μi, the peer discovery ratio (PDR) of the UE can be calculated as [33] The second term denotes the missing probability and the probability UE i cannot detect the contact. In other words, it is the probability that the UE i cannot find other peers.
According to the arrival contact rate of each UE, the aim of the peer discovery process is to maximize the total PDR, as This study attempts to determine the suitable beacon rates of all UEs to explore most peers in a D2D network under the constraints of budget power, the limitation of overall power consumption. This study further introduces a new scenario, i.e., energy harvesting, which previous studies did not consider. The energy harvesting enables UEs within charging areas to procure additional power from surrounding energy sources. Thus, under considering the UE may be out-of-energy and the study determines the suitable beacon rates of all UEs to maximize PDR according to the budget power and harvested power.
The problem that we investigated is formally described, as follows. First, for a beacon rate of UE i equal to µ i , the peer discovery ratio (PDR) of the UE can be calculated as [33] The second term denotes the missing probability and the probability UE i cannot detect the contact. In other words, it is the probability that the UE i cannot find other peers.
According to the arrival contact rate of each UE, the aim of the peer discovery process is to maximize the total PDR, as where TP is the transmission power used to send a beacon message and BP is the total budget power quantum of the network. However, if a UE is out-of-energy, it cannot find any peers in the network and any other peers cannot find it. Thus, in performing the peer discovery process, the objective given in Equation (8) should be modified, as follows: where P i is the power reserved for peer discovery of UE i. Differing from Equation (8), the UE cannot find other peers and it cannot be found by any other peers when P i ≤ 0 in Equation (9), because it is a more reasonable condition. Also note that the constraint denotes that the total power consumption must be less than the budget power quantum. Thus, the problem statement is formally given, as follows.
Input: the parameters marked input in Table 2. Output: the beacon rate of UE i, µ i , ∀i.

Solutions
Before introducing the ERRD algorithm in detail, the overall operations of UEs and BS are first described. The peer discovery is conducted per time period T. At the beginning of time period, each UE will send its ID and its harvested power quantum to BS. Thus, BS will know which UEs in its coverage and immediately extracts their social features from the database. According to UEs' social features and the harvested power quantum, and the pre-determined budget power quantum, BS can run ERRD to determine the beacon rate for each UE, and then sends the determined beacon rate to the corresponding UE. After receiving this rate, the UE will send the beacons accordingly. Since each UE only calculates the harvested power quantum, according to Equation (8), its operation is very simple. Therefore, below, we focus on describing the algorithm in BS, i.e., ERRD.
In D2D communications, the beacon rate is proportional to the amount of consumed power. Hence, the ERRD algorithm that is proposed in this study first virtually allocates a proper power quantum to each UE and then determines the corresponding beacon rate according to this allocated power quantum. As described in the following sub-sections, ERRD comprises three steps, namely Social Ratio Allocation (SRA), Energy Ratio Allocation (ERA), and Beacon Rate Decision (BRD). SRA first allocates the budget power among the UEs, depending on their social ratios. For each UE, ERA then adjusts this allocated power quantum according to the amount of energy that is harvested by the UE. Finally, BRD computes the beacon rate of each UE, depending on the adjusted allocated power quantum.
Below, the concept of ERRD is first explained. As ERRD is composed of three steps, we describe SRA, ERA, and BRD in sequence in Section 4.1. After describing the concept of ERRD, we formally exhibit the pseudo code of ERRD in Section 4.2. Finally, an example to illustrate the overall ERRD operation is given in Section 4.3.

Social Ratio Allocation
Let AP be the total allocated power quantum and GP i be the amount of power that is allocated to UE i. Since the peer discovery process is subject to the constraint that the total consumed power must be less than or equal to budget power BP, AP is initialized as BP. According to [21], when F CI i,j (x) has a uniform distribution (as shown in Equation (4)), the maximum PDR is achieved when the beacon rate of each UE is set proportional to its social ratio, i.e., µ i ∝ r i . Thus, the initially allocated power quantum for every UE is given by GP i = AP × r i .

Energy Ratio Allocation
The ERA step is the most critical step in ERRD. In order to properly explain the step, it is appropriate to introduce the following intuitive thought regarding the approach for allocating the power quantum, depending on the amount of power that is obtained from energy harvesting. The smart approach employed by ERRD is then introduced. •

Intuitive thought
The SRA step in the ERRD algorithm allocates the budget power quantum among the UEs based on their social ratios. However, some UEs can obtain power from energy harvesting, and hence their obtained power quantum OP i may exceed the allocated power quantum GP i . In this case, the UEs can utilize the obtained power quantum, rather than the originally allocated power quantum to perform their beacon transmissions. The unused portion of the allocated power quantum can then be returned to AP for the re-allocation to other UEs. Conversely, if OP i is less than GP i , i.e., the UE only acquires little (or no) energy via harvesting, the UE uses the original allocated power quantum GP i and returns its obtained power quantum OP i to AP. Combining these two cases, the unused power quantum of UE i is given by the minimum of OP i and GP i , i.e., min(OP i , GP i ).
Since AP and GP i are dependent, it is necessary to obtain them alternately while using an iterative approach. Let x k denote the value of x in the k-th iteration and k = 0 denote the initial value. Therefore, the initial value of AP, AP 0 , is set as the budget power quantum, BP, and it is allocated to each UE in accordance with r i , i.e., P 0 i = AP 0 × r i . Any unused power quantum must be returned to AP. Consequently, AP is equal to the sum of BP and the returned power quantum of all the UEs, i.e., Once AP k+1 is calculated, the new allocated power quantum of each UE can be obtained as GP k+1 i = AP k+1 × r i , and used to calculate AP k+2 accordingly. As shown in Lemma 1, AP is non-decreasing as the number of iterations increases. Proof.
Intuitively, AP k can be iteratively calculated until no further change in its value is obtained, i.e., AP k+1 = AP k . However, while such an approach is technically feasible, it requires many iterations to converge, causing high complexity, and it is hence impractical for real-world peer discovery applications. Consequently, the following smart approach is proposed instead.

•
Smart approach According to the relationship between OP i and GP i , ERRD classifies the UEs into two groups, namely G 1 when OP i < GP i and G 2 when OP i ≥ GP i . Therefore, Equation (10) can be re-formulated as Lemma 2: In each iteration, the UEs in G 2 may be shifted to G 1 , but the UEs in G 1 cannot be shifted to G 2 . Proof.

1.
From Lemma 1, AP is non-decreasing and GP i is also non-decreasing, i.e., GP k+1 Since OP i is fixed, the condition OP i < GP k i implies that OP i < GP k+1 i . Therefore, any member in G 1 will not be shifted to G 2 . 3.
Consequently, in each iteration, only the UEs in G 2 may be shifted to G 1 , but no UEs in G 1 will be shifted to G 2 .
From Lemma 2, the ERA process seeks to shift any UEs belonging to G 2 to G 1 if possible in each iteration in order to reduce the convergence time. However, to achieve this, two issues must first be addressed, namely (1) which UE should be considered first to be shifted and (2) whether this UE can actually be shifted.
ERA determines the answer to the first issue by inspecting the ratio of the obtained power quantum over the corresponding social ratio, i.e., OPR i = OP i /r i . A smaller value of OPR i implies that OP i is more likely to be less than GP i . In other words, UE i is more likely to belong to G 1 . Thus, in the k-th iteration, according to the current grouping, UE m k , whose OPR m k is the smallest among those of all the UEs in G k 2 , is chosen as the pivot UE, and it is most likely to be shifted from G k 2 to G k 1 , as described in the following lemma.

Lemma 3:
The pivot UE m k among G k 2 is most likely be shifted to G k 1 .

2.
In the k-th iteration, as AP k is fixed, UE m k has the smallest OPR in G 2 , i.e., and hence OP m k is most likely to be less than GP k m k . Consequently, UE m k is most likely to be shifted from G k 2 to G k 1 .
ERA considers whether or not this UE can actually be shifted after determining the pivot UE. To achieve this, a virtual critical point, defined as OP m k = GP m k , is considered. At this critical point, the virtual total allocated power quantum, VAP k , obtained using m k as the baseline, is given as The virtual allocated power quantum for m k is then obtained as VGP k m k =VAP k × r m k . The virtual critical point represents the threshold between G k 2 and G k 1 . Thus, checking whether VGP k m k is more than OP m k provides an efficient means of determining whether UE m k belongs to G k 2 or G k 1 . If VGP k m k is larger than OP m k , the pivot UE m k should be shifted from G k 2 to G k 1 and the iteration process should continue to the next round. Otherwise, the pivot UE m k should remain in G k 2 . As the pivot UE m k among G k 2 is the most likely be shifted to G k 1 , the other UEs belonging to G k 2 also remain in G k 2 . In other words, the grouping process is complete and no further changes in the memberships of G 1 and G 2 are required.
Once the grouping process is finished (in iteration *), GP * m * can be computed as After solving Equation (13), GP * m * can be obtained as AP * can then be computed as AP * = GP * m * /r m * . In addition, each GP * i can be computed with AP * . Finally, the allocated power quantum for each UE is set as the maximum between GP * i and OP i .

Beacon Rate Decision
In the BRD step, since OP i is considered, the sum of all the allocated power quanta is equal to the budget power quantum plus all the obtained power quanta, as shown in Lemma 4. Proof.

1.
For UE i in G 2 , the allocated power quantum is GP * i . As the UE has an obtained power quantum of OP i , the allocated power quantum occupying the budget power quantum is equal to GP * i − OP i . Thus, it follows that Due to this value exceeding the BP quantum, the final allocated power quantum should be normalized through multiplication by R, i.e., the ratio of BP over the sum of the allocated power quanta. That is, Finally, the µ i of each UE i is computed as the allocated power quantum multiplied by R and divided by TP.

Pseudo Code
Algorithm 1 shows the pseudo code of the ERRD algorithm, where lines 1-7 correspond to the SRA step, lines 8-24 describe the ERA step, and lines 25-29 are the BRD step. In the SRA step, lines 1-7 initialize the variables and compute r i , GP i , and OPR i for all the UEs. In the ERA step, lines 8-12 classify the UEs into two groups. In particular, UEs with an OP i greater than or equal to GP i are shifted into group G 2 . Lines 15-17 then select UE m whose OPR m is the smallest among those of all the UEs in G 2 as the pivot UE and compute VAP and VGP m accordingly. Line 18 checks whether the pivot UE meets the condition that OP m is larger than VGP m . If the condition holds, the final AP is computed in lines 19 and 20. Otherwise, lines 14-23 are repeated iteratively until this condition is reached. Finally, in the BRD step, R is calculated in line 25. Line 27 then computes GP i based on the final AP obtained in the ERA process and line 28 normalizes the allocated power quantum and computes the normalized µ i .

Algorithm 1 Energy Ratio Rate Allocation Algorithm
if OP i ≥ GP i then 10 :

: end for
Note that the pseudo code of ERRD is run in BS at the beginning of each time period T. All allocating, returning, and re-allocating power quanta from the BS to the UEs are virtually calculated in BS, rather than the real transfer between BS and UEs. Therefore, the communication between BS and UEs happens when each UE sends its ID and its harvested power quantum to BS at the beginning of each time period, and BS sends the determined beacon rate to the corresponding UE after it has executed ERRD. The

Illustrative Example
The following discussions present an illustrative example to demonstrate the detailed operational steps of ERRD. It is assumed that the network contains five UEs with average social contacts, λ i , of 1, 4, 9, 16, and 25, respectively. It is further assumed that the OP i values of the five UEs are 5, 4, 5, 3, and 1, respectively. Finally, the total BP is assumed to be 15 and TP is set as 1.
The SRA step first computes r i in accordance with the average social contacts, λ i , i.e., r 1 = Thus, r 1 , r 2 , r 3 , r 4 , and r 5 are obtained as 1 15 , 2 15 . 3 15 , 4 15 , and 5 15 , respectively (see Table 3). For each UE, GP i is then computed as the product of AP, which is initialized as BP, and r i , i.e., GP i = AP × r i . In other words, GP 1 , GP 2 , GP 3 , GP 4 and GP 5 are obtained as 1, 2, 3, 4, and 5, respectively. In addition, OPR i is calculated as OP i /r i . For example, OPR 1 = 5 1 15 = 75. The UEs are then classified into two groups, namely G 1 or G 2 , by comparing OP i with GP i . In the present example, OP 1 , OP 2 , and OP 3 are greater than GP 1 , GP 2 , and GP 3 , respectively, while OP 4 and OP 5 are less than GP 4 and GP 5 , respectively. Consequently, group G 1 is determined to be {4, 5}, while group G 2 is obtained as {1, 2, 3}. Based on the results that are presented in Table 3, the ERA step selects UE 3 as the pivot UE, since the value of OPR 3 is the smallest among all of the OPR values in G 2, To check whether the grouping process is finished, ERA uses the ratio of r i over r 3 to compute VAP as BP The VAP result is then used to compute VGP 3 as 29 × 3 15 = 29 5 . OP 3 is less than VGP 3 , that is, the virtually allocated power quantum exceeds the obtained power quantum. Hence, UE 3 is shifted from group G 2 to group G 1 . The procedure that is described above is iteratively repeated until the grouping process is complete. For the present example, G * 1 is obtained as {3, 4, 5} and G * 2 , as {1, 2}. Once the grouping process is finished, GP * 2 is GP * m and in accordance with Equation (14)  =30 (line 20). Finally, the values of GP * i are computed using the determined value of AP * . Choosing the larger value between OP i and GP * i for each UE, the total consumed power quantum is obtained as 5 + 4 + 6 + 8 + 10 = 33, which is equal to the sum of BP and all OP i (15+5+4+5+3+1). As the total consumed power quantum is larger than BP, normalization by R = 15 33 = 5 11 is performed and used to compute µ i . For example, µ 1 = 5 × 5 11 × 1 TP = 25 11 . Table 4 shows the final results for all the UEs.

Performance Evaluation
The PDR performance of the proposed ERRD algorithm was compared with that of the previously reported Social-Based Grouping (SBG) algorithm [21]. In SBG, the UEs are grouped based on their social feature: centrality and the UEs in the same group are assigned the same beacon rate. In the present simulations, SBG classifies the UEs into three groups.

Dataset and Environment
The simulations were performed using the Infocom06 user mobility trace [34], which consisted of the D2D communication contacts of 98 individuals that were recorded over the IEEE Infocom Conference in 2006. The first half of the dataset was used to calculate the beacon rates from the contact rates of the UEs, while the second half was used to evaluate and compare the performance of the two schemes (ERRD and SBG). The dataset contains no information regarding the actual physical locations of the users and the BS. Consequently, in performing the simulations, an artificial environment for energy harvesting was created, with dimensions of 600 × 600 m 2 . The simulation field was partitioned into a 3 × 3 grid containing a BS with a 10-watt power in the center and eight CDs with a five-watt power distributed around the outside (see Figure 2). The UEs were uniformly deployed in the 3 × 3 grid initially and their locations were then randomly moved as the simulations proceeded to simulate the mobility of UEs. Thus, the UEs obtained harvested power from different CDs (or the BS) at different points in the simulation process.
rates of the UEs, while the second half was used to evaluate and compare the performance of the two schemes (ERRD and SBG). The dataset contains no information regarding the actual physical locations of the users and the BS. Consequently, in performing the simulations, an artificial environment for energy harvesting was created, with dimensions of 600 × 600 m 2 . The simulation field was partitioned into a 3 × 3 grid containing a BS with a 10-watt power in the center and eight CDs with a five-watt power distributed around the outside (see Figure 2). The UEs were uniformly deployed in the 3 × 3 grid initially and their locations were then randomly moved as the simulations proceeded to simulate the mobility of UEs. Thus, the UEs obtained harvested power from different CDs (or the BS) at different points in the simulation process. The adopted dataset is the same as the dataset used in [21], because it is the most popular realistic dataset used in D2D communications, although it does not include the actual locations of UEs. To the best of our knowledge, there is no typical dataset that includes UE contacts and UE locations. Thus, in this study, we adopted the dataset in [21] and simulated the UE locations and their mobility by ourselves.
Each simulation was run for 1000 seconds. The plotted simulation results were then computed as the average value obtained over 100 simulations that were performed under identical conditions. The default environmental parameters were set as shown in Table 5. In accordance with the Infocom06 dataset, the number of UEs was set as N = 98. Moreover, the path loss exponent, v, of the BS and CD RF transmissions was set as 2, the energy harvesting efficiency factor, σ, was set as 0.7, and the transmission power of the beacon messages was set as 20 mW. The BP of the peer discovery process (with a duration of 1000 s) was set as 10 W. Finally, the mobility (moving speed) of the UEs was set as 1 under the assumption that the mobility in the Infocom06 dataset is 1. In general, as the moving speed of the UEs increases, the contact interval between them reduces, and vice versa. Thus, for a mobility value that is equal to 2, the contact interval between the UEs is equal to half that in the Infocom06 dataset.
Some assumptions are made for simplifying the simulations: (1) The harvested power quantum is according to Equation (6), although this quantum will be affected by interferences or other factors in a real environment. (2) The locations of UEs are limited in this area 600 × 600 m 2 . That is, if a UE moves outside this area, then its location will be randomly located within this area. (3) The beacon is perfectly transmitted, i.e., it will not encounter any collision or be interfered by other noises.
The simulations compared the performance of the ERRD algorithm with that of SBG under different settings of the BP, CDP, and mobility parameters. We compared ERRD and SBG, rather than other solutions because of two points. (1) As described in Section 2, ERRD and SBG belong to the same type: social-aware and network-assisted, but other solutions belong to different types. (2) As SBG used social features to adjust the beacon rate, it always outperforms other solutions without social-awareness. The evidence was exhibited in [21]. The UEs were classified into two types to facilitate the comparison between the two schemes, namely those with an increased beacon rate (IBR) and those with a decreased beacon rate (DBR), respectively. In the former case, the UEs using The adopted dataset is the same as the dataset used in [21], because it is the most popular realistic dataset used in D2D communications, although it does not include the actual locations of UEs. To the best of our knowledge, there is no typical dataset that includes UE contacts and UE locations. Thus, in this study, we adopted the dataset in [21] and simulated the UE locations and their mobility by ourselves.
Each simulation was run for 1000 seconds. The plotted simulation results were then computed as the average value obtained over 100 simulations that were performed under identical conditions. The default environmental parameters were set as shown in Table 5. In accordance with the Infocom06 dataset, the number of UEs was set as N = 98. Moreover, the path loss exponent, v, of the BS and CD RF transmissions was set as 2, the energy harvesting efficiency factor, σ, was set as 0.7, and the transmission power of the beacon messages was set as 20 mW. The BP of the peer discovery process (with a duration of 1000 s) was set as 10 W. Finally, the mobility (moving speed) of the UEs was set as 1 under the assumption that the mobility in the Infocom06 dataset is 1. In general, as the moving speed of the UEs increases, the contact interval between them reduces, and vice versa. Thus, for a mobility value that is equal to 2, the contact interval between the UEs is equal to half that in the Infocom06 dataset.
Some assumptions are made for simplifying the simulations: (1) The harvested power quantum is according to Equation (6), although this quantum will be affected by interferences or other factors in a real environment. (2) The locations of UEs are limited in this area 600 × 600 m 2 . That is, if a UE moves outside this area, then its location will be randomly located within this area. (3) The beacon is perfectly transmitted, i.e., it will not encounter any collision or be interfered by other noises.
The simulations compared the performance of the ERRD algorithm with that of SBG under different settings of the BP, CDP, and mobility parameters. We compared ERRD and SBG, rather than other solutions because of two points. (1) As described in Section 2, ERRD and SBG belong to the same type: social-aware and network-assisted, but other solutions belong to different types. (2) As SBG used social features to adjust the beacon rate, it always outperforms other solutions without social-awareness. The evidence was exhibited in [21]. The UEs were classified into two types to facilitate the comparison between the two schemes, namely those with an increased beacon rate (IBR) and those with a decreased beacon rate (DBR), respectively. In the former case, the UEs using ERRD had a higher beacon rate than those using SBG, while, in the latter case, the beacon rate of the UEs using ERRD was lower than that of those using SBG.  Figure 3 shows the PDR that was obtained under the two schemes for various values of the BP in the range of 0-20 W. For both schemes, the PDR exhibits a logarithmic-type increase as BP increases. This result is reasonable, since, intuitively, as the BP initially increases, the beacon rates of the UEs also increase, and hence a greater number of peers can be found. However, as the BP continues to increase, the number of originally unfound peers reduces, and hence the improvement in the PDR also reduces. The contact intervals in the Infocom06 dataset are not uniformly distributed and some of the intervals are extremely short. Thus, if the beacon rate is assigned a very high value in an attempt to increase the number of discovered peers, the consumed BP significantly increases. Therefore, in practical implementations, the BP should be set in such a way as to achieve a satisfactory tradeoff between the PDR and the consumed BP. ERRD had a higher beacon rate than those using SBG, while, in the latter case, the beacon rate of the UEs using ERRD was lower than that of those using SBG.  Figure 3 shows the PDR that was obtained under the two schemes for various values of the BP in the range of 0-20 W. For both schemes, the PDR exhibits a logarithmic-type increase as BP increases. This result is reasonable, since, intuitively, as the BP initially increases, the beacon rates of the UEs also increase, and hence a greater number of peers can be found. However, as the BP continues to increase, the number of originally unfound peers reduces, and hence the improvement in the PDR also reduces. The contact intervals in the Infocom06 dataset are not uniformly distributed and some of the intervals are extremely short. Thus, if the beacon rate is assigned a very high value in an attempt to increase the number of discovered peers, the consumed BP significantly increases. Therefore, in practical implementations, the BP should be set in such a way as to achieve a satisfactory tradeoff between the PDR and the consumed BP. For a given value of the BP, ERRD consistently achieves a higher PDR than SBG. For example, the PDR achieved under ERRD for a BP of 1 W is around 190% higher than that obtained under SBG, while for a BP of 20 W, the PDR is approximately 8% higher. This performance improvement can be attributed to two main factors. First, the IBR UEs using ERRD can send more beacons than those using SBG, and hence achieve a higher PDR. Second, although the DBR UEs using ERRD send fewer beacons than those using SBG, they have a longer lifetime. Consequently, the UEs can still achieve a For a given value of the BP, ERRD consistently achieves a higher PDR than SBG. For example, the PDR achieved under ERRD for a BP of 1 W is around 190% higher than that obtained under SBG, while for a BP of 20 W, the PDR is approximately 8% higher. This performance improvement can be attributed to two main factors. First, the IBR UEs using ERRD can send more beacons than those using SBG, and hence achieve a higher PDR. Second, although the DBR UEs using ERRD send fewer beacons than those using SBG, they have a longer lifetime. Consequently, the UEs can still achieve a higher PDR than those using SBG. However, the performance improvement that is offered by ERRD reduces with an increasing BP. This finding is reasonable, since, under a large BP, the amount of harvested energy is relatively smaller than that directly allocated by the BS. In other words, for most UEs, OP i < GP i , and hence the value of GP i under ERRD is similar to that under SBG. Consequently, the PDR performance gap between the two schemes reduces.

Effect of BP
Some intermediate results can be observed to better understand the differences between ERRD and SBG. First, note that the overall power consumptions for ERRD and SBG are the same as the overall power consumption is always limited by the budget power quantum; so, showing this value is unnecessary. However, the power consumption of each UE is quite different for these two approaches. Thus, we observe two intermediate results: Number of out-of-energy UEs (NoE) and the coefficient of variation (CV) of power in UEs. The NoE represents the number of UEs that power quantum reserved for peer discovery has been exhausted, that is, the power consumption exceeds the value of the power quantum reserved for peer discovery plus the harvested power quantum. The out-of-energy UEs cannot find other peers and cannot be found by other peers. The CV represents the power distribution among the UEs. The lower CV means a better balance of UE battery power. Figure 4 shows the NoE and CV that were obtained under the two schemes for various values of the BP. Observing this figure, NoE increases as BP exceeds a threshold. When BP increases, the UE can send more beacons. Therefore, the probability that a UE exhausts its battery power becomes larger, resulting in the increase of NoE. The out-of-energy UE appears when BP is 6 for SBG, while the out-of-energy UE appears when BP is 16 for ERRD. This is because SBG only considers the social ratio to send the beacons, while ERRD not only considers the social ratio, but also harvested power. When a UE has low harvested power and a high social ratio, ERRD reduces its beacon rate to reduce the probability of out-of-energy. However, in this case, SBG still lets this UE send beacons at a high rate, so it is very likely to be out-of-energy. On the other hand, CVs for ERRD and SBG almost linearly increase as BP increases, but SBG has a sharper slope. SBG generates many out-of-energy UEs and many UEs having much power, as they obtain much harvested power and send few beacons, resulting in more unbalanced battery power distribution among UEs, i.e., a higher CV. higher PDR than those using SBG. However, the performance improvement that is offered by ERRD reduces with an increasing BP. This finding is reasonable, since, under a large BP, the amount of harvested energy is relatively smaller than that directly allocated by the BS. In other words, for most UEs, OPi < GPi, and hence the value of GPi under ERRD is similar to that under SBG. Consequently, the PDR performance gap between the two schemes reduces. Some intermediate results can be observed to better understand the differences between ERRD and SBG. First, note that the overall power consumptions for ERRD and SBG are the same as the overall power consumption is always limited by the budget power quantum; so, showing this value is unnecessary. However, the power consumption of each UE is quite different for these two approaches. Thus, we observe two intermediate results: Number of out-of-energy UEs (NoE) and the coefficient of variation (CV) of power in UEs. The NoE represents the number of UEs that power quantum reserved for peer discovery has been exhausted, that is, the power consumption exceeds the value of the power quantum reserved for peer discovery plus the harvested power quantum. The out-of-energy UEs cannot find other peers and cannot be found by other peers. The CV represents the power distribution among the UEs. The lower CV means a better balance of UE battery power. Figure 4 shows the NoE and CV that were obtained under the two schemes for various values of the BP. Observing this figure, NoE increases as BP exceeds a threshold. When BP increases, the UE can send more beacons. Therefore, the probability that a UE exhausts its battery power becomes larger, resulting in the increase of NoE. The out-of-energy UE appears when BP is 6 for SBG, while the out-of-energy UE appears when BP is 16 for ERRD. This is because SBG only considers the social ratio to send the beacons, while ERRD not only considers the social ratio, but also harvested power. When a UE has low harvested power and a high social ratio, ERRD reduces its beacon rate to reduce the probability of out-of-energy. However, in this case, SBG still lets this UE send beacons at a high rate, so it is very likely to be out-of-energy. On the other hand, CVs for ERRD and SBG almost linearly increase as BP increases, but SBG has a sharper slope. SBG generates many out-of-energy UEs and many UEs having much power, as they obtain much harvested power and send few beacons, resulting in more unbalanced battery power distribution among UEs, i.e., a higher CV.  Figure 5 shows the variation of the PDR with the CDP. In general, a larger CDP indicates that more UEs can harvest energy, or individual UEs can acquire a greater amount of energy. However, the total consumed power in the network is limited to BP, irrespective of the value assigned to CDP. In other words, for a larger CDP, ERRD cannot consume more power than BP, but can only allocate the power quantum to each UE more precisely, such that all of the UEs can more efficiently send their beacons. As shown in Figure 5, the PDR obtained under ERRD rapidly increases as the CDP   Figure 5 shows the variation of the PDR with the CDP. In general, a larger CDP indicates that more UEs can harvest energy, or individual UEs can acquire a greater amount of energy. However, the total consumed power in the network is limited to BP, irrespective of the value assigned to CDP. In other words, for a larger CDP, ERRD cannot consume more power than BP, but can only allocate the power quantum to each UE more precisely, such that all of the UEs can more efficiently send their beacons. As shown in Figure 5, the PDR obtained under ERRD rapidly increases as the CDP first rises since the beacon rate decision made by ERRD reduces the probability that the UEs run out of energy. For a CDP value greater than 2, almost none of the UEs are out of energy. In this case, the PDR slightly increases, since the IBR UEs can find a greater number of peers. However, the PDR is not guaranteed to continuously increase when the number of UEs with an energy-harvesting capability exceeds a certain threshold. Therefore, in implementing the ERRD algorithm, a threshold should be set, whereby when the number of UEs with an energy-harvesting capability is greater than this threshold, ERRD should let some of the UEs store the harvested energy in their batteries, rather than expending it on beacon transmissions in order to provide power for other applications. first rises since the beacon rate decision made by ERRD reduces the probability that the UEs run out of energy. For a CDP value greater than 2, almost none of the UEs are out of energy. In this case, the PDR slightly increases, since the IBR UEs can find a greater number of peers. However, the PDR is not guaranteed to continuously increase when the number of UEs with an energy-harvesting capability exceeds a certain threshold. Therefore, in implementing the ERRD algorithm, a threshold should be set, whereby when the number of UEs with an energy-harvesting capability is greater than this threshold, ERRD should let some of the UEs store the harvested energy in their batteries, rather than expending it on beacon transmissions in order to provide power for other applications. For the SBG scheme, the PDR also increases with an increasing CDP. However, it increases at a slower rate than under ERRD, since, even though SBG does not explicitly consider the power obtained from energy harvesting, the number of out-of-energy UEs still decreases as the number of UEs having an energy-harvesting capability increases. Figure 6 shows the NoE and CV that were obtained under the two schemes for various values of the CDP. Observing this figure, NoE decreases as CDP increases, because UEs can harvest more power. Therefore, the probability that a UE exhausts its battery power becomes lesser, resulting in the decrease of NoE. However, the NoE of SBG is significantly larger than that of ERRD, because of the reasons that are described in Figure 4. On the other hand, CVs for ERRD and SBG almost linearly decrease as CDP increases. When CDP is small, the UE in the BS coverage still can harvest much power, but the UEs in the CD coverage only harvest less power, resulting in a larger CV. Similar to Figure 4, we can see SBG has a sharper slope of CV than ERRD. For the SBG scheme, the PDR also increases with an increasing CDP. However, it increases at a slower rate than under ERRD, since, even though SBG does not explicitly consider the power obtained from energy harvesting, the number of out-of-energy UEs still decreases as the number of UEs having an energy-harvesting capability increases. Figure 6 shows the NoE and CV that were obtained under the two schemes for various values of the CDP. Observing this figure, NoE decreases as CDP increases, because UEs can harvest more power. Therefore, the probability that a UE exhausts its battery power becomes lesser, resulting in the decrease of NoE. However, the NoE of SBG is significantly larger than that of ERRD, because of the reasons that are described in Figure 4. On the other hand, CVs for ERRD and SBG almost linearly decrease as CDP increases. When CDP is small, the UE in the BS coverage still can harvest much power, but the UEs in the CD coverage only harvest less power, resulting in a larger CV. Similar to Figure 4, we can see SBG has a sharper slope of CV than ERRD. Figure 7 shows the effect of the UE mobility on the PDR under the two schemes. As expected, the PDR reduces with an increasing mobility for both schemes, since, as the UEs move more rapidly, the contact rate between them increases, and hence the peers are less easily found for a given beacon rate. Nonetheless, the ERRD algorithm consistently outperforms SBG by around 10% for all the values of the UE mobility. Although the duration available for energy harvesting from one particular CD reduces as the UE mobility increases, the chance of harvesting energy from the other CDs increases. As a result, the UE mobility has no significant effect on the amount of harvested energy under the ERRD scheme. A similar tendency also occurs under the SBG scheme. Consequently, the performance advantage of ERRD over SBG is maintained, irrespective of the value of the UE mobility.   Figure 7 shows the effect of the UE mobility on the PDR under the two schemes. As expected, the PDR reduces with an increasing mobility for both schemes, since, as the UEs move more rapidly, the contact rate between them increases, and hence the peers are less easily found for a given beacon rate. Nonetheless, the ERRD algorithm consistently outperforms SBG by around 10% for all the values of the UE mobility. Although the duration available for energy harvesting from one particular CD reduces as the UE mobility increases, the chance of harvesting energy from the other CDs increases. As a result, the UE mobility has no significant effect on the amount of harvested energy under the ERRD scheme. A similar tendency also occurs under the SBG scheme. Consequently, the performance advantage of ERRD over SBG is maintained, irrespective of the value of the UE mobility.

Conclusions and Future Works
This study has proposed a social-aware peer discovery mechanism, designated as ERRD, for energy harvesting-based D2D communications. ERRD first allocates the budget power of the network among the UEs, depending on their social ratios and then adjusts the allocated power quantum of each UE, depending on its harvested energy. Finally, ERRD sets the beacon rates of the   Figure 7 shows the effect of the UE mobility on the PDR under the two schemes. As expected, the PDR reduces with an increasing mobility for both schemes, since, as the UEs move more rapidly, the contact rate between them increases, and hence the peers are less easily found for a given beacon rate. Nonetheless, the ERRD algorithm consistently outperforms SBG by around 10% for all the values of the UE mobility. Although the duration available for energy harvesting from one particular CD reduces as the UE mobility increases, the chance of harvesting energy from the other CDs increases. As a result, the UE mobility has no significant effect on the amount of harvested energy under the ERRD scheme. A similar tendency also occurs under the SBG scheme. Consequently, the performance advantage of ERRD over SBG is maintained, irrespective of the value of the UE mobility.

Conclusions and Future Works
This study has proposed a social-aware peer discovery mechanism, designated as ERRD, for energy harvesting-based D2D communications. ERRD first allocates the budget power of the network among the UEs, depending on their social ratios and then adjusts the allocated power quantum of each UE, depending on its harvested energy. Finally, ERRD sets the beacon rates of the

Conclusions and Future Works
This study has proposed a social-aware peer discovery mechanism, designated as ERRD, for energy harvesting-based D2D communications. ERRD first allocates the budget power of the network among the UEs, depending on their social ratios and then adjusts the allocated power quantum of each UE, depending on its harvested energy. Finally, ERRD sets the beacon rates of the UEs based on their adjusted power quanta subject to the constraint that the total power allocated to the UEs may not exceed the BP quantum of the network. ERRD improves the PDR performance by increasing the beacon rates of the high-energy-harvesting UEs, thereby increasing the number of peers that they can discover; and, reducing the beacon rates of the low-energy-harvesting UEs, thereby extending their lifetimes and prolonging the period for which they can participate in the peer discovery process. The simulation results have shown that ERRD outperforms the SBG scheme that is reported in the literature by around 8-190%, depending on the BP quantum of the network. The PDR performance of ERRD does not significantly increase as the BP increases beyond 10 W or the CDP exceeds 2. However, in mobile environments, ERRD retains a 10% performance advantage over SBG, irrespective of the moving speed of the UEs.
Currently, 5G D2D communications and energy-harvesting devices are not so popular to give a realistic example of use case. Therefore, the paper is more research-oriented, rather than system-oriented. However, we believe that D2D communications and energy-harvesting devices will become more popular in the future. The scenario that is considered in the study will actually happen and our proposed solution, ERRD, can be applied.
A complete D2D communication not only includes peer discovery, but also mode selection and resource allocation. Future studies will aim to establish a more comprehensive energy harvesting-based D2D communication model that includes mode selection and resource allocation. In its current form, ERRD considers the past behavior (sociality) and present condition (energy harvesting) of the UEs, but it does not consider the impact of the surrounding information, such as the number of peers in the network or the interference. Therefore, in future studies, ERRD will be extended to take such information into account in order to obtain a more robust estimation of the most appropriate beacon rate for D2D communications. Finally, the performance of ERRD in a realistic environment should be further investigated, as our study proves the outperformance of ERRD by simulations. For example, a UE who has a social-network account and, at the same time, can harvest energy can enhance how much of PDR. To observe this, ERRD should be realistically implemented in BS. We will pay the efforts on this implementation in the future.

Conflicts of Interest:
The authors declare no conflict of interest.