This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Age-Threshold Slotted ALOHA for Optimizing Information Freshness in Mobile Networks

Fangming Zhao, Nikolaos Pappas, Chuan Ma, Xinghua Sun, Tony Q. S. Quek, and Howard H. Yang    Fangming Zhao, Nikolaos Pappas, Senior Member, IEEE, Chuan Ma, Member, IEEE,
Xinghua Sun, Member, IEEE, Tony Q. S. Quek, Fellow, IEEE, and Howard H. Yang, Member, IEEE
Abstract

We optimize the Age of Information (AoI) in mobile networks using the age-threshold slotted ALOHA (TSA) protocol. The network comprises multiple source-destination pairs, where each source sends a sequence of status update packets to its destination over a shared spectrum. The TSA protocol stipulates that a source node must remain silent until its AoI reaches a predefined threshold, after which the node accesses the radio channel with a certain probability. Using stochastic geometry tools, we derive analytical expressions for the transmission success probability, mean peak AoI, and time-average AoI. Subsequently, we obtain closed-form expressions for the optimal update rate and age threshold that minimize the mean peak and time-average AoI, respectively. In addition, we establish a scaling law for the mean peak AoI and time-average AoI in mobile networks, revealing that the optimal mean peak AoI and time-average AoI increase linearly with the deployment density. Notably, the growth rate of time-average AoI under TSA is half of that under SA. When considering the optimal mean peak AoI, the TSA protocol exhibits comparable performance to the traditional slotted ALOHA protocol. These findings conclusively affirm the advantage of TSA in reducing higher-order AoI, particularly in densely deployed networks.

Index Terms:
Age of information, mobile networks, status updating protocol, age threshold, interference.

I Introduction

Timeliness has emerged as a fundamental prerequisite for Internet-of-Things (IoT) services, encompassing various applications such as autonomous driving, smart factory operations, and intelligent healthcare. Acquiring timely updates from massively deployed sensors is crucial for making prompt and accurate decisions in these applications[1].

To assess the timeliness of delivered messages, the notion of Age of Information (AoI) has been proposed in [2][3], where the authors showed that minimizing AoI is fundamentally different from optimizing conventional metrics such as delay or throughput. As a result, plenty of research efforts have been devoted to ensuring timely information delivery over communication networks. This paper employs the age-threshold slotted ALOHA (TSA), a variant of the classic ALOHA protocol, to enhance AoI in random access networks. Specifically, TSA exploits the AoI information contained in the feedback messages to control each transmitter’s channel access, thereby improving AoI. We study the minimization of AoI through the lens of the signal-to-interference-plus-noise ratio (SINR) model, by optimally adjusting the parameters of TSA according to the network configuration. The resultant scheme has low complexity and is particularly relevant to large-scale wireless systems.

I-A Related Works

In the early stages of AoI analysis for large-scale networks, many researchers adopted the traditional analytical framework of slotted ALOHA. For example, [4] characterized and optimized the average AoI by optimal tuning the channel access probability. [5] optimized the peak AoI using First-Come-First-Serve (FCFS) and Last-Come-First-Serve (LCFS) strategies, showing that with optimal packet arrival rate and channel access probability, the peak AoI scales linearly with the number of access nodes. [6] considered a multi-channel ALOHA network, enhancing AoI performance via joint radio access control and resource allocation. [7] explores the average AoI performance in a heterogeneous random access network where stochastic arrival nodes coexist with generate-at-will nodes. The function of irregular repetition slotted ALOHA (IRSA) on various AoI performance metrics has been demonstrated in [8]. The dynamic frame slotted ALOHA (FSA) has been studied in [9]. These efforts have alleviated the performance bottlenecks of the AoI metric in large-scale networks. However, the primary objective behind the conventional ALOHA-like protocol was to mitigate channel contention, thereby enhancing transmission success probability or throughput. It was not explicitly designed to cater to the distinct attributes of the AoI metric. Hence, it becomes essential to explore random access protocols that closely align with the attributes of the AoI metric.

To design an age-aware random access network, it is essential to (a) mitigate network interference, (b) minimize the waiting time for information packets at the transmitters, and (c) ensure equalization of the update packet reception intervals. A straightforward strategy to address these needs is to prioritize nodes with lower timeliness for updates, allowing nodes with higher information freshness to remain silent temporarily. Based on this idea, several works [10, 11, 12, 13, 14, 15] extended the Slotted ALOHA (SA) into an age-aware variant. Specifically, a threshold-based age-dependent random access protocol was proposed in [10, 11], where each node accesses the channel only when its instantaneous AoI exceeds a predetermined threshold. A distributed transmission strategy was proposed in [12] based on the age gain, defined as the reduction of instantaneous AoI when packets are successfully delivered. An Index-Prioritized Random Access scheme was proposed in [13], where nodes access the channel according to their indices that reflect the urgency of the update. Minislotted threshold ALOHA, designed to minimize the average AoI in the high payload scenario, is investigated in [14]. The dynamic FSA with age gain threshold is considered in [15].

TABLE I: A summary of the literature on the age-threshold aloha policy.
Interference Analysis AoI Performance Analysis Optimal Configuration Performance Limits Scaling law
Decoding model Spatial Attribute Stable Region Numerical Closed-form Numerical Closed-form
[10] Collison model
[11] Collison model
[12] Collison model
[13] Collison model
[14] Collison model
[15] Collison model
[16] SINR model
Proposed Work SINR model

The aforementioned research endeavors predominantly concentrate on the single-cell setting, where multiple nodes transmit data to a central access point. However, it is imperative to consider the ad-hoc network scenario for a more comprehensive understanding. Another noteworthy aspect is that the interference in the above references is modeled by the collision model[17], wherein concurrent transmissions always lead to communication failures. Such a model neglects essential physical elements in wireless channels, such as fading, path loss, and spatial network topology, oversimplifying the co-channel interference effects. In contrast, the SINR model is preferable to the update packet verdict on the receiver end, as it can adequately account for these critical features in the ad-hoc network scenario.

Consequently, a line of recent studies has been carried out, analyzing and optimizing AoI in wireless networks on the basis of the SINR model[18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30]. Specifically, the upper and lower bounds of the time-average AoI in Poisson bipolar networks are provided in [18] based on the favorable/dominant system argument, and the upper bound is improved by reconstruction of the dominant system in [19]. The authors in [20] derived the average peak AoI with different packet arrival patterns, i.e., the time-triggered and event-triggered traffic for large-scale wireless networks, while the authors in [21] and [22] analyze the AoI statistics under different buffer sizes and package management disciplines. The authors in [23] present a comprehensive study on the interplay between throughput and AoI in a cellular-based IoT network. Explicit optimal expression of peak AoI of slotted ALOHA, and its optimal parameters are derived in [24], the distinctions and similarities between the optimal parameters for the peak AoI and average AoI are extensively discussed in [25]. The energy-AoI tradeoff is further explored in [26]. The authors in [27] examine the AoI performance of random access networks under frame-slotted ALOHA protocol. Locally adaptive strategies for optimizing AoI have been studied in [28], [29], and [30], where the channel access probability, update rate, and transmit power at each node are adjusted according to the local communication environment to minimize AoI.

The existing results primarily pertain to the conventional slotted ALOHA protocol, whilst the effects of the age threshold remain unexplored. Although a recent work [16] has investigated the effect of TSA on AoI performance in static Poisson bipolar networks, the analysis involves complicated expressions of the SINR meta distribution, prohibiting it from optimizing TSA to reap the full potential of this protocol. In this paper, we consider the AoI optimization under TSA in a high mobility case, with which we can establish tractable analysis and obtain closed-form expression of the optimal network configuration, AoI performance limits, and scaling law that provide insights into network designs.

I-B Contributions

The main contributions of this paper are summarized below.

  • We develop a mathematical framework for analyzing the impact of the TSA protocol on the AoI performance of random access networks. We establish a fixed-point equation to characterize the transmission success probability, encompassing key features such as the update rate, age threshold, channel fading, deployment density, and co-channel interference.

  • We derive analytical expressions of the mean peak AoI and time-average AoI over the typical link. Our analysis culminates in an optimization strategy to minimize these AoI metrics, achieved by a joint optimization of the update rate and age threshold. Specifically, we investigate the performance limits of the mean peak AoI under TSA protocol and identify the corresponding parameter pairs that enable its minimization. We also propose an alternative iterative algorithm to achieve the optimal parameter configuration that minimizes the time-average AoI. Additionally, we derive closed-form approximations for these optimal parameters, which are asymptotically accurate.

  • We establish a scaling law of AoI in random access networks. Notably, it reveals that although the optimal mean peak AoI and time-average AoI increase linearly with the deployment density, the growth rate of time-average AoI under TSA is only half of that under SA.

  • Considering the inherent instability of ALOHA-like networks, especially under conditions of high interference levels and traffic load, we propose parameter tuning strategies for both mean peak AoI and time-average AoI in the presence of a bistable region, aiming to bolster the robustness of the TSA protocol network.

To better highlight the technical contributions of this paper, we summarize existing literature that dealt with age-threshold policy in Table I.

TABLE II: key notations
Notations Definition Notations Definition
λ\lambda Spatial deployment density Δ^0\hat{\Delta}_{0} Mean peak AoI of typical link
rr Transceiver distance Δ¯0\bar{\Delta}_{0} Time-average AoI of typical link
η\eta Update rate η^|A\hat{\eta}^{*}_{|A} Optimal update rate (Fixed AA, mean peak AoI)
AA Age-threshold A^|η\hat{A}^{*}_{|\eta} Optimal age-threshold (Fixed η\eta, mean peak AoI)
α\alpha Path-loss fading coefficient (A^,η^)(\hat{A}^{*},\hat{\eta}^{*}) Optimal age-threshold &\& update rate (mean peak AoI)
ρ\rho SNR η¯|A\bar{\eta}^{*}_{|A} Optimal update rate (Fixed AA, time-average AoI)
θ\theta Decoding SINR threshold A¯|η\bar{A}^{*}_{|\eta} Optimal update rate (Fixed η\eta, time-average AoI)
psp_{\mathrm{s}} Transmission success probability (A¯,η¯)(\bar{A}^{*},\bar{\eta}^{*}) Optimal age-threshold &\& update rate (time-average AoI)
psLp_{\mathrm{s}}^{L} psp_{\mathrm{s}} in low efficiency state Δ^0\hat{\Delta}_{0}^{*} Optimal mean peak AoI
psHp_{\mathrm{s}}^{H} psp_{\mathrm{s}} in high efficiency state Δ¯0\bar{\Delta}_{0}^{*} Optimal time-average AoI

II System Model

II-A Spatial Configuration

We consider a random access network comprised of source-destination pairs. The source nodes are scattered according to a homogeneous Poisson point process (HPPP) Φs={Xi}i=1\Phi_{\mathrm{s}}=\{X_{i}\}_{i=1}^{\infty} of density λ\lambda. Each source is paired with a destination in the distance rr and oriented in a random direction. According to the displacement theorem[31, 32], locations of the destination nodes constitute an independent HPPP Φd\Phi_{\mathrm{d}} of the same density. We add a receiver at the origin to the point process Φd\Phi_{\mathrm{d}}. We also add its tagged transmitter, denoted by X0X_{0}, to the point process Φs\Phi_{\mathrm{s}}. We refer to this link pair as the typical one. 111This setting is a large-scale analog to the classical model of Random Networks [33], where the fixed distance between each transmitter-receiver pair represents the average value. It is noteworthy that the analysis developed in this article can be extended to investigate the AoI performance in networks with centralized infrastructures where multiple sources can associate with a common destination and are located at random distances from the destinations [34, 35]. Based on Slivnyak’s theorem[31, 32], we can concentrate on analyzing the performance of the typical link. A sequence of status information is generated at each source, encapsulated into information packets, and transmitted over a shared spectrum. When a source node sends out information packets, it transmits at a fixed power. We assume the channel between any pair of nodes is affected by the Rayleigh fading with a unit mean, which varies independently across time slots, and path loss that follows the power law. We also assume the received signal is subjected to white Gaussian thermal noise. The main notations used throughout this article are summarized in Table II.

II-B Temporal Pattern

In this network, time is segmented into equal-length intervals, in which each slot is the minimum duration required to transmit an information packet (and receive feedback from the destination). We assume the network is synchronized[36, 37]. Each source adopts the generate-at-will policy [38] for the status updating. Particularly, if a node decides to transmit, it will generate a new sample at the beginning of the time slot and send the information packet to the destination immediately. By the end of the same time slot, the packet is successfully decoded if the received SINR exceeds a decoding threshold; otherwise, the transmission fails, and in that case, the packet is dropped. The delivery of packets, therefore, incurs a delay of one time slot. Additionally, the positions of the links vary according to the high mobility model [24, 39], i.e., spatial locations of the link pairs are (re)generated independently in each time slot according to the same Poisson bipolar network model. 222By leveraging the framework of [40], one can extend the results obtained in this paper to investigate AoI in a mobile network with spatial correlations over time.

II-C Transmission Protocol

We employ the TSA protocol for status updating. Specifically, every source node stays inactive until its AoI reaches a predefined threshold, denoted by AA, upon which the source node turns on the status updating mode: it generates a fresh sample with probability η\eta at the beginning of the subsequent time slot, encapsulates that information into a packet, and immediately sends it to the destination. If the received SINR surpasses a decoding threshold θ\theta, the transmission succeeds; the receiver then feeds back an ACK to the source, and the age is reset to one.333In this paper, we make the commonly adopted assumption that the ACK feedback is error-free [39, 11, 12]. It is, nonetheless, noteworthy that the analysis can be extended to account for the feedback error of ACK messages. Otherwise, the source will generate a new sample in the next time slot with probability η\eta. This probability is usually referred to as the update rate.

II-D Performance Metric

We put the main focus on the notion of AoI in this paper. As depicted in Fig. 1, AoI grows linearly with time in the absence of new updates at the destination, and it drops to the time elapsed since the generation of the latest packet received. Formally, the age evolution process over the typical link can be written as follows [2, 38]

Δ0(t)=tG0(t),\Delta_{0}(t)=t-G_{0}(t), (1)

where G0(t)G_{0}(t) denotes the generation time of the latest packet successfully delivered over this link.

Refer to caption
Figure 1: An illustration of the evolution of the age of information over the link under consideration in the discrete-time model.

In this paper, we use the mean peak AoI and time-average AoI over the typical link as our performance metrics, defined respectively as [21]

Δ^0=limN1Nn=1NΔ0(T0(n))\displaystyle\hat{\Delta}_{0}=\lim_{N\rightarrow\infty}\frac{1}{N}\sum_{n=1}^{N}\Delta_{0}(T_{0}(n)) (2)

where T0(n)T_{0}(n) denotes the time slot at which the nn-th packet is successfully delivered, and

Δ¯0=limT1Tt=1TΔ0(t).\bar{\Delta}_{0}=\lim_{T\to\infty}\frac{1}{T}\sum_{t=1}^{T}\Delta_{0}(t). (3)

Since the point processes formed by the TSA protocol still exhibit ergodicity and stationarity, the statistical characteristics of AoI performance for the typical link are identical to those of the other links.

III AoI Analysis

III-A SINR at the Typical Receiver

Since all the source nodes use the same spectrum for packet delivery, each transmission would be affected by others owing to co-channel interference. As such, at the typical destination, the received SINR at time slot tt is[41]

SINR0(t)=h00(t)rαi0hi0(t)vi(t)Xiα+ρ1,\mathrm{SINR}_{0}(t)=\frac{h_{00}(t)r^{-\alpha}}{\sum_{i\neq 0}h_{i0}(t)v_{i}(t)||X_{i}||^{-\alpha}+\rho^{-1}}, (4)

where hi0(t)exp(1)h_{i0}(t)\sim\exp(1) represents the small-scale fading between transmitter ii and the typical receiver, α\alpha is the path-loss exponent, \|\cdot\| stands for the Euclidean norm, ρ\rho denotes the signal-to-noise ratio (SNR) between the transmit power of each source node and the noise power, and vi(t)v_{i}(t) is a binary function, where vi(t)=1v_{i}(t)=1 indicates that transmitter ii initiates a packet transmission at time slot tt, and vi(t)=0v_{i}(t)=0 otherwise.

III-B Transmission Success Probability

According to the transmission policy, an information packet is successfully delivered if the received SINR exceeds a decoding threshold θ\theta. Therefore, the probability of successfully transmitting a packet from the typical source node can be written as [32]

ps(t)=(SINR0(t)>θ).p_{\mathrm{s}}(t)=\mathbb{P}\left(\mathrm{SINR}_{0}(t)>\theta\right). (5)

Under the high mobility random walk model, the received SINRi(t)\text{SINR}_{i}(t) of each transmitter ii, ii\in\mathbb{N}, is independent and identically distributed (i.i.d.) across time tt. By symmetry, the transmission success probability is also identical across the transmitters. As such, we drop the time index from ps(t)p_{\mathrm{s}}(t) in the sequel. The following theorem provides a formal characterization of this quantity.

Theorem 1.

The transmission success probability under TSA is given by the following fixed-point equation

ps=exp(λcηr21+Aηpsθrαρ),p_{\mathrm{s}}=\exp\left(-\frac{\lambda c\eta r^{2}}{1+A\eta p_{\mathrm{s}}}-\frac{\theta r^{\alpha}}{\rho}\right), (6)

in which

c=πθ2αΓ(12α)Γ(1+2α),c=\pi\theta^{\frac{2}{\alpha}}\Gamma\left(1-\frac{2}{\alpha}\right)\Gamma\left(1+\frac{2}{\alpha}\right), (7)

where Γ()\Gamma(\cdot) denotes the Gamma function [42].

Proof:

Please see Appendix A. ∎

The expression in (6) accounts for the impact of several key factors, including the update rate, age threshold, path loss, and deployment density, on the transmission success probability. As opposed to other variants of SA, e.g., the FSA [27], the transmission success probability under TSA is given in the form of a fixed-point equation, which steams from the spatial-temporal interactions amongst the source nodes across the network.

The parameter cc in (6) reflects the level of spatial contention [43], measuring the network’s capability of spatial reuse by quantifying how fast the transmission success probability deteriorates when the interferers’ density increases.

Remark 1.

In the static deployment case [16], the transmission success probability of each source varies from node to node and cannot be characterized by Theorem 1. Consequently, it requires deriving the SINR meta-distribution to characterize the link performance. Although this approach provides a fine-grained perspective for network analysis, it often leads to complex analytical expressions that hinder the exploration of AoI limits. In contrast, the mobile network model considered by this paper mitigates the influence of spatial correlations among the nodes, enabling tractable results that provide useful network design guidelines, which will be detailed in the subsequent sections.

The solution(s) to fixed-point equation (6) can be obtained via fixed-point iterations. Specifically, we denote by ps,1(0,1]p_{\mathrm{s},1}\in(0,1] and ps,np_{\mathrm{s},n} the initial point and estimation of transmission success probability after the nnth iteration, respectively. The iteration process can be expressed as follows

ps,n+1=exp(λcηr21+Aηps,nθrαρ).p_{\mathrm{s},n+1}=\exp\left(-\frac{\lambda c\eta r^{2}}{1+A\eta p_{\mathrm{s},n}}-\frac{\theta r^{\alpha}}{\rho}\right). (8)

While the convergence of the fixed-point iteration is immune to the initial value, depending on the network configuration, different starting points may lead to different solutions. Therefore, we investigate the conditions of different distributions of roots of (6) in the following.

Corollary 1.

The fixed-point equation in (6) has three distinct roots, 0<psL<psM<psH<10<p_{\mathrm{s}}^{L}<p_{\mathrm{s}}^{M}<p_{\mathrm{s}}^{H}<1, if

λcr2>4η\lambda cr^{2}>\frac{4}{\eta} (9)

and

Al<A<Ah,A_{l}<A<A_{h}, (10)

where AlA_{l} and AhA_{h} are given respectively as follows

Al=λcr2(12+141λcr2η)1ηexp(θrαρ112+141λcr2η),\displaystyle A_{l}=\frac{\lambda cr^{2}\left(\frac{1}{2}+\sqrt{\frac{1}{4}-\frac{1}{\lambda cr^{2}\eta}}\right)-\frac{1}{\eta}}{\exp\left(-\frac{\theta r^{\alpha}}{\rho}-\frac{1}{\frac{1}{2}+\sqrt{\frac{1}{4}-\frac{1}{\lambda cr^{2}\eta}}}\right)}, (11)
Ah=λcr2(12141λcr2η)1ηexp(θrαρ112141λcr2η).\displaystyle A_{h}=\frac{\lambda cr^{2}\left(\frac{1}{2}-\sqrt{\frac{1}{4}-\frac{1}{\lambda cr^{2}\eta}}\right)-\frac{1}{\eta}}{\exp\left(-\frac{\theta r^{\alpha}}{\rho}-\frac{1}{\frac{1}{2}-\sqrt{\frac{1}{4}-\frac{1}{\lambda cr^{2}\eta}}}\right)}. (12)

When A=AlA=A_{l} or A=AhA=A_{h}, (6) has two different solutions; otherwise, (6) has a unique solution.

Proof:

Please see Appendix B. ∎

Notably, in the scenario that (6) has three distinct roots, if ps,0(0,psM)p_{\mathrm{s},0}\in(0,p_{\mathrm{s}}^{M}) or ps,0(psM,1]p_{\mathrm{s},0}\in(p_{\mathrm{s}}^{M},1], the fixed-point iteration (8) results in psLp_{\mathrm{s}}^{L} and psHp_{\mathrm{s}}^{H}, respectively. However, regardless of the initial value of ps,0p_{\mathrm{s},0}, there is no guarantee from (8) that the ultimate solution would be psMp_{\mathrm{s}}^{M}.

Refer to caption
Figure 2: Approximate trajectory of ps,np_{\mathrm{s},n}. Spatial deployment density λ=0.15\lambda=0.15, update rate η=1\eta=1, decoding threshold θ=0dB\theta=0~{}\text{dB}, SNR ρ=20dB\rho=20~{}\text{dB}, path-loss exponent α=3.8\alpha=3.8, and transceiver distance r=3r=3.

A pictorial example is provided in Fig. 2, where we set the parameters that satisfy λcr2>4η\lambda cr^{2}>\frac{4}{\eta} and showcase six distinct curves, arranged based on the value of the age threshold. Specifically, they are: (i)(i) A=0A=0, (ii)(ii) A=10(0,Al)A=10\in(0,A_{l}), (iii)(iii) A=Al31A=A_{l}\approx 31, (iv)(iv) A=50(Al,Ah)A=50\in(A_{l},A_{h}), (v)(v) A=Ah135A=A_{h}\approx 135, and (vi)(vi) A=150(Ah,+)A=150\in(A_{h},+\infty). Let us take the case of A=50(Al,Ah)A=50\in(A_{l},A_{h}) as an example. In this case, the region of ps,np_{\mathrm{s},n} has four possible scenarios: (1) 0<ps,n<psL0<p_{\mathrm{s},n}<p_{\mathrm{s}}^{L}; (2) psL<ps,n<psMp_{\mathrm{s}}^{L}<p_{\mathrm{s},n}<p_{\mathrm{s}}^{M}; (3) psM<ps,n<psHp_{\mathrm{s}}^{M}<p_{\mathrm{s},n}<p_{\mathrm{s}}^{H}; and (4) psH<ps,n<1p_{\mathrm{s}}^{H}<p_{\mathrm{s},n}<1. In the first and second scenarios, ps,np_{\mathrm{s},n} will gradually approach psLp_{\mathrm{s}}^{L}. In contrast, ps,np_{\mathrm{s},n} tends toward psHp_{\mathrm{s}}^{H} in the third and fourth cases. As such, in this example, there are two candidate points, namely, psLp_{\mathrm{s}}^{L} and psHp_{\mathrm{s}}^{H}, that ps,np_{\mathrm{s},n} can converge to. In a similar vein, we can show that ps,np_{\mathrm{s},n} only converges to psHp_{\mathrm{s}}^{H} (or psLp_{\mathrm{s}}^{L}) in other cases.

In the terminology of random network analysis [44], psLp_{\mathrm{s}}^{L} and psHp_{\mathrm{s}}^{H} are steady-state points of the fixed-point equation (6) since they are attainable by controlling the input to the fixed-point iteration in (8). In contrast, psMp_{\mathrm{s}}^{M} is not a steady-state point in this context.

To this end, we define the following region of network parameter configurations under which the fixed-point equation of transmission success probability has steady-state solutions, and we term these ranges the stable regions.

  • Bistable region:444In the context of ALOHA-like random access networks, the bi-stable characteristic refers to a phenomenon where the system’s performance (measured by, e.g., transmission success probability) can exist in one of two statistically stable equilibrium points under certain conditions[45, 44, 46, 47]. The bi-stable region is the unstable region in practice, since the system will easily converge to the state where every transmitter has a small transmission success probability. In this case, data transmissions often fail, significantly deteriorating the performance of real-time applications.

    ={(λ,θ,r,α,ρ,A,η)|λcr2>4η,Al<A<Ah},\mathcal{B}=\left\{\left(\lambda,\theta,r,\alpha,\rho,A,\eta\right)\Big{|}\lambda cr^{2}>\frac{4}{\eta},A_{l}<A<A_{h}\right\},

    in which the network has two steady-state points psLp_{\mathrm{s}}^{L} and psHp_{\mathrm{s}}^{H}.

  • Monostable region: =psHpsL=¯\mathcal{M}=\mathcal{M}_{p_{\mathrm{s}}^{H}}\cup\mathcal{M}_{p_{\mathrm{s}}^{L}}=\bar{\mathcal{B}} (¯\bar{\mathcal{B}} is the complement of set \mathcal{B}); in this case, the network has one steady-state point, which can be further partitioned into the following.

    (i) High-efficiency region:

    psH={(λ,θ,r,α,ρ,A,η)|λcr24ηandλcr2>4η,AAh},\mathcal{M}_{p_{\mathrm{s}}^{H}}=\left\{\left(\lambda,\theta,r,\alpha,\rho,A,\eta\right)\big{|}\lambda cr^{2}\leq\frac{4}{\eta}~{}\text{and}~{}\lambda cr^{2}>\frac{4}{\eta},A\geq A_{h}\right\},

    in which the iteration converges to psHp_{\mathrm{s}}^{H}.

    (ii) Low-efficiency region:

    psL={(λ,θ,r,α,ρ,A,η)|λcr2>4η,AAl},\mathcal{M}_{p_{\mathrm{s}}^{L}}=\left\{\left(\lambda,\theta,r,\alpha,\rho,A,\eta\right)\big{|}\lambda cr^{2}>\frac{4}{\eta},A\leq A_{l}\right\},

    in which the iteration converges to psLp_{\mathrm{s}}^{L}.

The following corollary summarizes the properties of the steady-state points in terms of the configuration parameters.

Corollary 2.

The steady-state transmission success probabilities psLp_{\mathrm{s}}^{L} and psHp_{\mathrm{s}}^{H}are increasing functions of the age threshold AA, while they decrease in terms of the spatial deployment density λ\lambda and update rate η\eta.

Proof:

Please see Appendix C. ∎

Corollary 3.

The fixed-point iteration in (8) always converges with a linear convergence rate.

Proof:

Please see Appendix D. ∎

It is intuitive that in densely deployed networks, the transmission success probability will benefit (more) from the TSA protocol due to its capability in mitigating channel contention. This can also be observed from Fig. 2 by comparing the value of transmission success probability at A=0A=0 with those at A>0A>0.

III-C AoI Statistics

Using the transmission success probability, we can derive closed-form expressions for the mean peak AoI and time-average AoI as follows.

Theorem 2.

Under the TSA protocol, the mean peak AoI over the typical link is

Δ^0=A+1ηps,\hat{\Delta}_{0}=A+\frac{1}{\eta p_{\mathrm{s}}}, (13)

and the time-average AoI is

Δ¯0=A+12+1ηpsA+12(1+Aηps),\bar{\Delta}_{0}=\frac{A+1}{2}+\frac{1}{\eta p_{\mathrm{s}}}-\frac{A+1}{2\left(1+A\eta p_{\mathrm{s}}\right)}, (14)

where psp_{\mathrm{s}} is given in (6).

Proof:

Please see Appendix E. ∎

Leveraging the above theorem, we can investigate different regimes of the age threshold to demystify its influence on the AoI performance. Firstly, if we set the age threshold AA to be zero, the results in Theorem 2 reduce to the mean peak AoI and time-average AoI in large-scale networks under the SA protocol[10], i.e.,

Δ^0SA=Δ¯0SA=exp(λcr2η+θrαρ1)η.\hat{\Delta}^{\mathrm{SA}}_{0}=\bar{\Delta}^{\mathrm{SA}}_{0}=\frac{\exp\left(\lambda cr^{2}\eta+\theta r^{\alpha}\rho^{-1}\right)}{\eta}. (15)

Consequently, we can compare the AoI achieved under the SA and TSA protocols to identify the operating regimes where TSA outperforms SA.

We first compare the mean peak AoI attained under the SA and TSA protocols for A>0A>0, which yields

Δ^0SAΔ^0TSA=1ηpsSA(A+1ηpsTSA)\displaystyle\hat{\Delta}_{0}^{\mathrm{SA}}-\hat{\Delta}_{0}^{\mathrm{TSA}}=\frac{1}{\eta p_{\mathrm{s}}^{\mathrm{SA}}}-\left(A+\frac{1}{\eta p_{\mathrm{s}}^{\mathrm{TSA}}}\right)
=(exp(λcr2η)exp(λcr2η1+AηpsTSA))ηexp(ρ/θrα)A\displaystyle=\frac{\left(\exp\left(\lambda cr^{2}\eta\right)-\exp\left(\frac{\lambda cr^{2}\eta}{1+A\eta p_{\mathrm{s}}^{\mathrm{TSA}}}\right)\right)}{\eta\exp\left(\rho/\theta r^{\alpha}\right)}-A
=(a)(λcr2(AηpsTSA1+AηpsTSA)exp(λcr2ηξp+θrαρ))A\displaystyle\overset{(a)}{=}\left(\lambda cr^{2}\left(\frac{A\eta p_{\mathrm{s}}^{\mathrm{TSA}}}{1+A\eta p_{\mathrm{s}}^{\mathrm{TSA}}}\right)\exp\left(\lambda cr^{2}\eta\xi_{\mathrm{p}}+\frac{\theta r^{\alpha}}{\rho}\right)\right)-A
λcr2(Aη1+AηpsTSA)A,\displaystyle\geq{\lambda cr^{2}\left(\frac{A\eta}{1+A\eta p_{\mathrm{s}}^{\mathrm{TSA}}}\right)}-A, (16)

where (a)(a) follows from the mean value theorem, in which ξp[11+Aηps,1]\xi_{\mathrm{p}}\in[\frac{1}{1+A\eta p_{\mathrm{s}}},1]. As such, when λcr21+AηpsTSAη\lambda cr^{2}\geq\frac{1+A\eta p_{\mathrm{s}}^{\mathrm{TSA}}}{\eta}, we have

Δ^0SAΔ^0TSA=λcr2(Aη1+AηpsTSA)A0.\hat{\Delta}_{0}^{\mathrm{SA}}-\hat{\Delta}_{0}^{\mathrm{TSA}}=\lambda cr^{2}\left(\frac{A\eta}{1+A\eta p_{\mathrm{s}}^{\mathrm{TSA}}}\right)-A\geq 0. (17)

Following similar steps, we can compare the time-average AoI attained under the SA and TSA protocols as follows:

Δ¯0SAΔ¯0TSA>1ηpsSA(A+12+1ηpsTSA)\displaystyle\bar{\Delta}_{0}^{\mathrm{SA}}-\bar{\Delta}_{0}^{\mathrm{TSA}}>\frac{1}{\eta p_{\mathrm{s}}^{\mathrm{SA}}}-\left(\frac{A+1}{2}+\frac{1}{\eta p_{\mathrm{s}}^{\mathrm{TSA}}}\right)
=(exp(λcr2η)exp(λcr2η1+AηpsTSA))ηexp(ρ/θrα)A+12\displaystyle=\frac{\left(\exp\left(\lambda cr^{2}\eta\right)-\exp\left(\frac{\lambda cr^{2}\eta}{1+A\eta p_{\mathrm{s}}^{\mathrm{TSA}}}\right)\right)}{\eta\exp\left(\rho/\theta r^{\alpha}\right)}-\frac{A+1}{2}
=(a)(λcr2(AηpsTSA1+AηpsTSA)exp(λcr2ηξa+θrαρ))A+12\displaystyle\overset{(a)}{=}\left(\lambda cr^{2}\left(\frac{A\eta p_{\mathrm{s}}^{\mathrm{TSA}}}{1+A\eta p_{\mathrm{s}}^{\mathrm{TSA}}}\right)\exp\left(\lambda cr^{2}\eta\xi_{\mathrm{a}}+\frac{\theta r^{\alpha}}{\rho}\right)\right)-\frac{A+1}{2}
λcr2(Aη1+AηpsTSA)A+12,\displaystyle\geq{\lambda cr^{2}\left(\frac{A\eta}{1+A\eta p_{\mathrm{s}}^{\mathrm{TSA}}}\right)}-{\vphantom{\lambda cr^{2}\left(\frac{A\eta}{1+A\eta p_{\mathrm{s}}^{\mathrm{TSA}}}\right)}\frac{A+1}{2}}, (18)

where (a)(a) follows from the mean value theorem, in which ξa[11+Aηps,1]\xi_{\mathrm{a}}\in[\frac{1}{1+A\eta p_{\mathrm{s}}},1]. When λcr2(1+A)(1+AηpsTSA)2Aη\lambda cr^{2}\geq\frac{(1+A)(1+A\eta p_{\mathrm{s}}^{\mathrm{TSA}})}{2A\eta}, we have

Δ¯0SAΔ¯0TSA>λcr2(Aη1+AηpsTSA)A+120.\bar{\Delta}_{0}^{\mathrm{SA}}-\bar{\Delta}_{0}^{\mathrm{TSA}}>\lambda cr^{2}\left(\frac{A\eta}{1+A\eta p_{\mathrm{s}}^{\mathrm{TSA}}}\right)-\frac{A+1}{2}\geq 0. (19)

The inequalities (17) and (19) demonstrate that the TSA has the potential of outperforming SA in reducing AoI.

However, if we raise the age threshold to an excessively large value, i.e., AA\to\infty, then we have

Δ^0A+,\displaystyle\hat{\Delta}_{0}\sim A\rightarrow+\infty, (20)
Δ¯0A2+,\displaystyle\bar{\Delta}_{0}\sim\frac{A}{2}\rightarrow+\infty, (21)

implying that a large age threshold will be detrimental to the AoI.

To this end, it is apparent that the age threshold, as well as the update rate, shall be adequately tuned to reveal the full potential of TSA. We will detail the approaches to obtain the optimal parameters in the following sections. For ease of exposition, we start our derivations without considering the influence of the bi-stable region. In Section VI, we will demonstrate how to adjust the parameters to ensure network stability while concurrently optimizing the AoI performance in the presence of a bi-stable region.

IV Mean Peak AoI Optimization

This section optimizes the mean peak AoI by adjusting the age threshold and update rate. The optimization problem can be formulated as

minη,AΔ^0s.t.η(0,1],A.\begin{split}&\underset{\eta,A}{\min}\quad\hat{\Delta}_{0}\\ &~{}\mathrm{s.t.}\quad\eta\in(0,1],\\ &\quad~{}~{}~{}~{}~{}A\in\mathbb{N}.\end{split} (22)

The age threshold is usually set as an integer, leading to a mixed integer programming problem in (22), which is non-trivial to solve. As such, we first convert the problem into a continuous one by relaxing the age threshold to a continuous variable and then round the result to an integer. More precisely, the relaxed optimization problem can be written as the following

minη,AΔ^0s.t.η(0,1],A+.\begin{split}&\underset{\eta,A}{\min}\quad\hat{\Delta}_{0}\\ &~{}\mathrm{s.t.}\quad\eta\in(0,1],\\ &\quad~{}~{}~{}~{}~{}A\in\mathbb{R^{+}}.\end{split} (23)

This optimization problem can be decomposed into two subproblems: 1) optimizing the update rate η\eta given an age threshold AA, and 2) optimizing the age threshold AA given an update rate η\eta.

In the following, we will first investigate the solutions to these subproblems respectively to obtain a comprehensive understanding of the differences and similarities between the optimal parameters. Then, we solve (23) by jointly tuning the update rate and age threshold.

IV-A Optimizing the Update Rate

The following theorem presents the optimal update rate for a given age threshold.

Lemma 1.

Given an age threshold AA, the optimal mean peak AoI Δ^0η=η^A\hat{\Delta}_{0}^{\eta=\hat{\eta}^{*}_{A}} is given by

Δ^0η=η^|A={λcr2exp(1+θrαρ1),ifλcr2>1+Aps,A+(ps,)1,otherwise\hat{\Delta}_{0}^{\eta=\hat{\eta}^{*}_{|A}}=\begin{cases}\lambda cr^{2}\exp\left(1+\theta r^{\alpha}\rho^{-1}\right),&\text{if}~{}\lambda cr^{2}>1+Ap_{s,*}\\ A+\left(p_{\mathrm{s},*}\right)^{-1},&\text{otherwise}\end{cases} (24)

which is achieved when the update rate η\eta is set to be

η^|A={1λcr2Aexp(1θrαρ1),ifλcr2>1+Aps,1,otherwise\hat{\eta}^{*}_{|A}=\begin{cases}\dfrac{1}{\lambda cr^{2}-A\exp\left(-1-\theta r^{\alpha}\rho^{-1}\right)},&\text{if}~{}\lambda cr^{2}>1+Ap_{s,*}\\ 1,&\text{otherwise}\end{cases} (25)

where ps,p_{\mathrm{s},*} is the root of the following equation

ps,=exp(λcr21+Aps,θrαρ).p_{\mathrm{s},*}=\exp\left(-\frac{\lambda cr^{2}}{1+Ap_{\mathrm{s},*}}-\frac{\theta r^{\alpha}}{\rho}\right). (26)
Proof:

Please see Appendix F. ∎

Lemma 1 shows that the optimal update rate shall be tuned to the maximum, i.e., η^|A=1\hat{\eta}^{*}_{|A}=1, when the λcr21+Aps,\lambda cr^{2}\leq 1+Ap_{\mathrm{s},*}. This is because when the nodes are sparsely deployed, each receiver has a high chance of successfully receiving updates. As such, the AoI performance can be effectively improved by increasing update rates. In contrast, when the nodes are deployed densely, each source shall decrease its update rate to circumvent channel contention and promote timely information delivery across the network.

IV-B Optimizing the Age Threshold

The following theorem characterizes the optimal age threshold for a given update rate.

Lemma 2.

Given an update rate η\eta, the optimal mean peak AoI Δ^0A=A^η\hat{\Delta}_{0}^{A=\hat{A}^{*}_{\eta}} is given by

Δ^0A=A^|η={λcr2exp(1+θrαρ),ifλcr2>1η1ηexp(λcr2η+θrαρ),otherwise\hat{\Delta}_{0}^{A=\hat{A}^{*}_{|\eta}}=\begin{cases}\lambda cr^{2}\exp\left(1+\frac{\theta r^{\alpha}}{\rho}\right),&\text{if}~{}~{}\lambda cr^{2}>\dfrac{1}{\eta}\\ \frac{1}{\eta}\exp\left(\lambda cr^{2}\eta+\frac{\theta r^{\alpha}}{\rho}\right),&\text{otherwise}\end{cases} (27)

which is achieved when the age threshold AA is set to be

A^|η={(λcr21η)exp(1+θrαρ),ifλcr2>1η0,otherwise.\hat{A}^{*}_{|\eta}=\begin{cases}\left(\lambda cr^{2}-\dfrac{1}{\eta}\right)\exp\left(1+\dfrac{\theta r^{\alpha}}{\rho}\right),&\text{if}~{}~{}\lambda cr^{2}>\dfrac{1}{\eta}\\ 0,&\text{otherwise}.\end{cases} (28)
Proof:

Please see Appendix G. ∎

Lemma 2 indicates that when the interference level is high, i.e., λcr2η>1\lambda cr^{2}\eta>1, imposing an age threshold to the ALOHA-like protocol is beneficial for reducing mean peak AoI. Conversely, when the interference level is low, the conventional SA protocol is preferable in decreasing the mean peak AoI.

IV-C Joint Optimization

By integrating results from the above, the theorem below presents a closed-form result for the optimal pair of the update rate and age threshold that minimizes the mean peak AoI.

Theorem 3.

The optimal mean peak AoI is given by

Δ^0={λcr2exp(1+θrαρ1),ifλcr2>1exp(λcr2+θrαρ1),otherwise\hat{\Delta}_{0}^{*}=\begin{cases}\lambda cr^{2}\exp\left(1+\theta r^{\alpha}\rho^{-1}\right),&\text{if}~{}~{}\lambda cr^{2}>1\\ \exp\left(\lambda cr^{2}+\theta r^{\alpha}\rho^{-1}\right),&\text{otherwise}\end{cases} (29)

which is achieved when the age threshold and update rate together satisfy

(A^,η^)={(A1,η1),ifλcr2>1(0,1),otherwise(\hat{A}^{*},\hat{\eta}^{*})=\begin{cases}\left(A_{1},\eta_{1}\right),&\text{if}~{}~{}\lambda cr^{2}>1\\ \left(0,1\right),&\text{otherwise}\end{cases} (30)

where (A1,η1)\left(A_{1},\eta_{1}\right) satisfy the following

A1=(λcr21η1)exp(1+θrαρ)A_{1}=\left(\lambda cr^{2}-\frac{1}{\eta_{1}}\right)\exp\left(1+\frac{\theta r^{\alpha}}{\rho}\right) (31)

and η1\eta_{1} shall be confined in the range η1[1λcr2,1]\eta_{1}\in\left[\frac{1}{\lambda cr^{2}},1\right].

Proof:

Please see Appendix H. ∎

From (31), we observe that the optimal age threshold increases with the optimal update rate to mitigate the induced interference. Moreover, when λcr21\lambda cr^{2}\leq 1, the optimal TSA configuration reduces to the conventional SA. In addition to these, a (perhaps) more important message conveyed by Theorem 3 is that by appropriately adjusting the update rate, SA can attain the same optimal mean peak AoI as TSA. In other words, TSA and SA are equivalent in minimizing the mean peak AoI.

Remark 2.

The optimal tuning of the update rate and the age threshold depends on the statistical traffic information (such as the spatial deployment density), instead of the real-time number of active nodes requesting transmission or prior experience. Such parameters are collected when the network is initially deployed.

V Time-Average AoI Optimization

This section presents the optimization of time-average AoI by solving the following optimization problem

minη,AΔ¯0s.t.η(0,1],A.\begin{split}&\underset{\eta,A}{\min}\quad\bar{\Delta}_{0}\\ &~{}\mathrm{s.t.}\quad\eta\in(0,1],\\ &\quad~{}~{}~{}~{}~{}A\in\mathbb{N}.\end{split} (32)

Following similar approaches in the previous section, we relax the mixed integer programming problem (32) to a continuous one. Then, we decompose it into two subproblems and solve each one respectively. Finally, we integrate the whole procedure to establish a joint optimization. The concrete steps are detailed below.

V-A Optimizing the Update Rate

First, we optimize the update rate under a fixed age threshold, where the structural result is summarized as follows.

Lemma 3.

Given an age threshold AA, the optimal time-average AoI Δ¯0η=η¯|A\bar{\Delta}_{0}^{\eta=\bar{\eta}^{*}_{|A}} can be achieved by setting the update rate as

η¯|A={1λcr2Aexp(1θrαρ1),ifλcr2>1+Aps,;1,otherwise,\bar{\eta}^{*}_{|A}=\begin{cases}\dfrac{1}{\lambda cr^{2}-A\exp\left(-1-\theta r^{\alpha}\rho^{-1}\right)},&\text{if}~{}\lambda cr^{2}>1+Ap_{\mathrm{s},*};\\ 1,&\text{otherwise},\end{cases} (33)

where ps,p_{s,*} is given by the solution to the following

ps,=exp(λcr21+Aps,θrαρ).p_{\mathrm{s},*}=\exp\left(-\frac{\lambda cr^{2}}{1+Ap_{\mathrm{s},*}}-\dfrac{\theta r^{\alpha}}{\rho}\right). (34)

In this case, the time-average AoI is

Δ¯0η=η¯|A={2(λcr2exp(1+θrαρ1))2+A(A+1)2λcr2exp(1+θrαρ1)A,ifλcr2>1+Aps,;A+12+A12(1+Aps,)+1ps,(1+Aps,),otherwise.\bar{\Delta}_{0}^{\eta=\bar{\eta}^{*}_{|A}}=\begin{cases}\dfrac{2\left(\lambda cr^{2}\exp\left(1+\theta r^{\alpha}\rho^{-1}\right)\right)^{2}+A(A+1)}{2\lambda cr^{2}\exp\left(1+\theta r^{\alpha}\rho^{-1}\right)}-A,\\ ~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}\text{if}~{}\lambda cr^{2}>1+Ap_{\mathrm{s},*};\\ \\ \dfrac{A+1}{2}+\dfrac{A-1}{2(1+Ap_{\mathrm{s},*})}+\dfrac{1}{p_{\mathrm{s},*}(1+Ap_{\mathrm{s},*})},\\ ~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}\text{otherwise}.\end{cases} (35)
Proof:

Please see the Appendix I. ∎

The condition that λcr21+Aps,\lambda cr^{2}\leq 1+Ap_{\mathrm{s},*} can be expressed as

ps,exp(1θrαρ1).p_{\mathrm{s},*}\geq\exp\left(-1-\theta r^{\alpha}\rho^{-1}\right). (36)

This inequality indicates that if the transmission success probability is lower than exp(1θrαρ1)\exp\left(-1-\theta r^{\alpha}\rho^{-1}\right) when η=1\eta=1, the update rate should be reduced to ensure the transmission success probability stabilizes at this value to optimize the time-average AoI. Because the age threshold always contributes to the transmission success probability, the optimal update rate goes up with the increase of the age threshold.

A comparison between (25) and (33) reveals that for a given age threshold, the optimal update rates that minimize the mean peak AoI and time-average AoI admit the same structure. This finding seems surprising, as the two performance metrics have very different analytical expressions (cf. Theorem 1). And the reason boils down to that, in essence, the core of optimizing the update rate for mean peak AoI and time-average AoI lies at maximizing the spatial throughput ηps\eta p_{\mathrm{s}}.

V-B Optimizing the Age Threshold

In a similar vein, we explore the optimal structure of the age threshold under a fixed update rate. The following result provides a systematic approach to achieving this.

Lemma 4.

Given an update rate η\eta, the optimal time-average AoI Δ¯0A=A¯|η\bar{\Delta}_{0}^{A=\bar{A}^{*}_{|\eta}} can be achieved by setting the age threshold as

A¯|η={A1,ifλcr2>1η𝕎0(η2exp(θrαρ));0,otherwise,\bar{A}^{*}_{|\eta}=\begin{cases}A_{1},&\text{if}~{}~{}\lambda cr^{2}>\frac{1}{\eta}\mathbb{W}_{0}\left(\frac{\eta}{2}\exp\left(-\frac{\theta r^{\alpha}}{\rho}\right)\right);\\ 0,&\text{otherwise},\end{cases} (37)

where 𝕎0()\mathbb{W}_{0}(\cdot) is the principal branch of Lambert W function [48], and A1A_{1} is the single non-zero root of the following equation

(1+A1ηps)(ηexp(λcr2η1+A1ηpsθrαρ)1)λcr2η(1+A1ηps)2+(1+A1ηps)3=λcr2η.\begin{split}&\left(1+A_{1}\eta p_{\mathrm{s}}\right)\left(\eta\exp\left(\frac{-\lambda cr^{2}\eta}{1+A_{1}\eta p_{\mathrm{s}}}-\frac{\theta r^{\alpha}}{\rho}\right)-1\right)\\ &\quad\quad-\lambda cr^{2}\eta\left(1+A_{1}\eta p_{\mathrm{s}}\right)^{2}+\left(1+A_{1}\eta p_{\mathrm{s}}\right)^{3}=\lambda cr^{2}\eta.\end{split} (38)
Proof:

Please see the Appendix J. ∎

Remark 3.

The computational complexity (in terms of the number of iterations) in solving (38) can be reduced by first computing (1+A1ηps)(1+A_{1}\eta p_{\mathrm{s}}) and then calculate psp_{\mathrm{s}}, instead of directly solving (6) for the transmission success probability. The optimal age threshold A1A_{1} can be obtained by combining the numerical evaluations of (1+A1ηps)(1+A_{1}\eta p_{\mathrm{s}}) and psp_{\mathrm{s}}.

According to Lemma 4, we note that the optimal age threshold is A¯|η=0\bar{A}^{*}_{|\eta}=0 when

λcr21η𝕎0(η2exp(θrαρ)),\displaystyle\lambda cr^{2}\leq\frac{1}{\eta}\mathbb{W}_{0}\left(\frac{\eta}{2}\exp\left(-\frac{\theta r^{\alpha}}{\rho}\right)\right), (39)

indicating that when network interference is mild, one shall not impose any waiting duration at the source nodes but employ an SA-like protocol for the status updating. In contrast, TSA protocol is preferable when interference is severe, as it reduces concurrent transmissions, thus mitigating potential contentions. As a result, the nodes with urgent updating demands will benefit from a better transmission environment, which facilitates age performance.

The condition in (39) can be equivalently written as

ηln(2λcr2exp(θrαρ1))λcr2.\eta\leq\frac{-\ln\left(2\lambda cr^{2}\exp\left(\theta r^{\alpha}\rho^{-1}\right)\right)}{\lambda cr^{2}}. (40)

This implies that TSA outperforms SA in terms of AoI when the update frequency is relatively high. Furthermore, when λcr2>12exp(θrαρ1)\lambda cr^{2}>\tfrac{1}{2}\exp(-\theta r^{\alpha}\rho^{-1}), TSA protocol consistently achieves a better time-average AoI than SA. Notably, these observations are in line with those drawn in [12], which are concluded based on a collision channel model.

While it is possible to solve equation (38) numerically, obtaining a closed-form solution (which could provide design insights) is challenging. In that respect, we leverage (tight) upper and lower bounds of the time-average AoI to facilitate our optimization. Specifically, we can derive a lower bound to the average AoI as

Δ¯0LB=Aηps+12ηps+12ηps(1+Aηps),\displaystyle\bar{\Delta}^{\text{LB}}_{0}=\frac{A\eta p_{\mathrm{s}}+1}{2\eta p_{\mathrm{s}}}+\frac{1}{2\eta p_{\mathrm{s}}(1+A\eta p_{\mathrm{s}})}, (41)

and an upper bound as

Δ¯0UB=Aηps+12ηps+12ηps(1+Aηps)+12.\displaystyle\bar{\Delta}_{0}^{\text{UB}}=\frac{A\eta p_{\mathrm{s}}+1}{2\eta p_{\mathrm{s}}}+\frac{1}{2\eta p_{\mathrm{s}}(1+A\eta p_{\mathrm{s}})}+\frac{1}{2}. (42)

Then, we take a derivative of Δ¯0LB\bar{\Delta}^{\text{LB}}_{0} (and/or Δ¯0UB\bar{\Delta}^{\text{UB}}_{0}) with respect to AA and assign it to zero, which gives

(1+Aηps)3λcr2η(1+Aηps)2(1+Aηps)λcr2η=0.(1\!+\!A\eta p_{\mathrm{s}})^{3}-\lambda cr^{2}\eta(1\!+\!A\eta p_{\mathrm{s}})^{2}-(1\!+\!A\eta p_{\mathrm{s}})-\lambda cr^{2}\eta=0. (43)

By solving this equation, we obtain a closed-form expression for the suboptimal age threshold

A~|η=Ω23+(λcr2η3)Ω13+(λcr2η)2+33ηΩ13exp(3λcr2ηΩ13Ω23+λcr2ηΩ13+(λcr2η)2+3θrαρ),\tilde{A}^{*}_{|\eta}=\frac{\Omega^{\frac{2}{3}}+(\lambda cr^{2}\eta-3)\Omega^{\frac{1}{3}}+(\lambda cr^{2}\eta)^{2}+3}{3\eta\Omega^{\frac{1}{3}}\exp\!\left(\dfrac{-3\lambda cr^{2}\eta\Omega^{\frac{1}{3}}}{\Omega^{\frac{2}{3}}\!+\!\lambda cr^{2}\eta\Omega^{\frac{1}{3}}\!+\!(\lambda cr^{2}\eta)^{2}\!+\!3}-\dfrac{\theta r^{\alpha}}{\rho}\right)}, (44)

where

Ω=(λcr2η)3+18λcr2η+33(λcr2η)4+11(λcr2η)21.\Omega=(\lambda cr^{2}\eta)^{3}+18\lambda cr^{2}\eta+3\sqrt{3}\sqrt{(\lambda cr^{2}\eta)^{4}+11(\lambda cr^{2}\eta)^{2}-1}. (45)

Correspondingly, we have a closed-form expression, given in (46), for the time-average AoI under the suboptimal threshold. It is noteworthy that the difference between the bounds and the time-average AoI is at most 1212(1+Aηps)[0,12)\frac{1}{2}-\frac{1}{2(1+A\eta p_{\mathrm{s}})}\in[0,\frac{1}{2}), hence the suboptimal solution differs from the optimal one by (at most) a constant. Indeed, by comparing (14) and (42), we can see that the time-average AoI approaches the upper bound when the age threshold AA becomes large. Therefore, the solution given in (44) is asymptotically optimal.

Remark 4.

Based on (28) and (44), one can show that the optimal age threshold that minimizes the time-average AoI is larger than that for minimizing the mean peak AoI, since optimizing the former metric necessitates reducing the variance in the status update intervals. This observation indicates that the TSA parameters need to be carefully configured when optimizing different AoI metrics.

Δ¯0A=A~|η=(Ω23+λcr2ηΩ13+(λcr2η)2+3)2+9Ω23(1ηexp(3λcr2ηΩ13Ω23+λcr2ηΩ13+(λcr2η)2+3θrαρ))6ηΩ13(Ω23+λcr2ηΩ13+(λcr2η)2+3)exp(3λcr2ηΩ13Ω23+λcr2ηΩ13+(λcr2η)2+3θrαρ)+12.\displaystyle\bar{\Delta}_{0}^{A=\tilde{A}^{*}_{|\eta}}=\frac{\left(\Omega^{\frac{2}{3}}+\lambda cr^{2}\eta\Omega^{\frac{1}{3}}+(\lambda cr^{2}\eta)^{2}+3\right)^{2}\!\!+9\Omega^{\frac{2}{3}}\left(1-\eta\exp\left(\dfrac{-3\lambda cr^{2}\eta\Omega^{\frac{1}{3}}}{\Omega^{\frac{2}{3}}\!+\!\lambda cr^{2}\eta\Omega^{\frac{1}{3}}\!+\!(\lambda cr^{2}\eta)^{2}+3}-\dfrac{\theta r^{\alpha}}{\rho}\right)\right)}{6\eta\Omega^{\frac{1}{3}}\left(\Omega^{\frac{2}{3}}\!+\!\lambda cr^{2}\eta\Omega^{\frac{1}{3}}\!+\!(\lambda cr^{2}\eta)^{2}\!+\!3\right)\exp\left(\dfrac{-3\lambda cr^{2}\eta\Omega^{\frac{1}{3}}}{\Omega^{\frac{2}{3}}+\lambda cr^{2}\eta\Omega^{\frac{1}{3}}+(\lambda cr^{2}\eta)^{2}+3}-\dfrac{\theta r^{\alpha}}{\rho}\right)}+\frac{1}{2}. (46)

 

V-C Joint Optimization

The previous subsections have explored the optimization of the update rate (resp. age threshold) given the age threshold (resp. update rate). This part will investigate how to jointly adjust the update rate and age threshold to further reduce the time-average AoI.

Algorithm 1 Iterative algorithm for AoI optimization
1:Input λ\lambda, rr, θ\theta, ρ\rho, α\alpha, nn, ϵ\epsilon
2:Initialize i=1i=1 and choose η¯1\bar{\eta}^{*}_{1} from (0,1] at random
3:repeat
4:***Optimize age threshold under a fixed update rate***
5:     if λcr2η¯i>𝕎0(η¯i2exp(θrαρ1))\lambda cr^{2}\bar{\eta}^{*}_{i}>\mathbb{W}_{0}\left(\frac{\bar{\eta}^{*}_{i}}{2}\exp\left(-\theta r^{\alpha}\rho^{-1}\right)\right) then
6:         Use Bisection method to search for ψ=1+A¯iη¯ips\psi=1+\bar{A}^{*}_{i}\bar{\eta}^{*}_{i}p_{\mathrm{s}} from (38), with error accuracy ϵ\epsilon and in the interval [1,2eλcr2][1,2e\lambda cr^{2}]
7:         Compute ps=exp(λcr2ψ1θrαρ1)p_{\mathrm{s}}=\exp\left(-\lambda cr^{2}\psi^{-1}-\theta r^{\alpha}\rho^{-1}\right)
8:         Compute A¯i=(ψ1)/η¯ips\bar{A}^{*}_{i}=(\psi-1)/{\bar{\eta}^{*}_{i}p_{\mathrm{s}}}
9:         Compute Δ¯0,i\bar{\Delta}_{0,i} based on A¯i\bar{A}^{*}_{i}, η¯i\bar{\eta}^{*}_{i}, psp_{\mathrm{s}}, and (14)
10:     else
11:         Let A¯i=0\bar{A}^{*}_{i}=0
12:         Compute Δ¯0,i=exp(λcr2η¯i+θrαρ1)/η¯i\bar{\Delta}_{0,i}=\exp\left(\lambda cr^{2}\bar{\eta}^{*}_{i}+\theta r^{\alpha}\rho^{-1}\right)/\bar{\eta}^{*}_{i}      
13:     ii+1i\leftarrow i+1
14:***Optimize update rate under a fixed age threshold***
15:     Compute ps,p_{\mathrm{s},*} base on η=1\eta=1, A¯i1\bar{A}^{*}_{i-1}, and (6)
16:     if λcr2>1+A¯i1ps,\lambda cr^{2}>1+\bar{A}^{*}_{i-1}p_{\mathrm{s},*} then
17:         Compute η¯i\bar{\eta}^{*}_{i} based on A¯i1\bar{A}^{*}_{i-1}, and (33)
18:         Compute Δ¯0,i\bar{\Delta}_{0,i} based on A¯i1\bar{A}^{*}_{i-1}, and (35)
19:     else
20:         Let η¯i=1\bar{\eta}^{*}_{i}=1
21:         Compute Δ¯0,i\bar{\Delta}_{0,i} based on A¯i1\bar{A}^{*}_{i-1}, and (35)      
22:until |Δ¯0,iΔ¯0,i1|<ϵ|\bar{\Delta}_{0,i}-\bar{\Delta}_{0,i-1}|<\epsilon
23:Assign A¯A¯i1\bar{A}^{*}\leftarrow\bar{A}^{*}_{i-1}, η¯η¯i\bar{\eta}^{*}\leftarrow\bar{\eta}^{*}_{i} and Δ¯0Δ¯0,i\bar{\Delta}_{0}^{*}\leftarrow\bar{\Delta}_{0,i}
24:return A¯\bar{A}^{*}, η¯\bar{\eta}^{*}, Δ¯0\bar{\Delta}_{0}^{*}

Since the age threshold and update rate have a composite influence on the time-average AoI (cf. (6) and (14)), it is challenging to characterize the optimal values of these factors explicitly. As such, based on Lemmas 3 and 4, we propose an iterative approach, presented in Algorithm 1, to jointly optimize the age threshold and update rate. Specifically, let ii denote the iteration index, while η¯i\bar{\eta}^{*}_{i} and A¯i\bar{A}^{*}_{i} represent the update rate and age threshold in the iith iteration, respectively. As detailed in Algorithm 1, we initialize η¯1\bar{\eta}^{*}_{1} by a randomly chosen value from (0,1](0,1], and then obtain A¯1\bar{A}^{*}_{1} using Lemma 4. Based on η¯1\bar{\eta}^{*}_{1} and Lemma 3, we further calculate A¯2\bar{A}^{*}_{2} and η¯2\bar{\eta}^{*}_{2}. Subsequently, we can compute each update of η¯i\bar{\eta}_{i}^{*} and A¯i\bar{A}^{*}_{i} for i2i\geq 2. The iterations repeat until reaching the termination condition |Δ¯0,iΔ¯0,i1|<ϵ|\bar{\Delta}_{0,i}-\bar{\Delta}_{0,i-1}|<\epsilon, where ϵ\epsilon is an error accuracy.

V-C1 Convergence and Optimality Analysis of the Proposed Algorithms

Upon dividing the problem (32), in Lemmas 3 and 4, we have demonstrated that each of the sub-problems, i.e., (V.A) and (V.B), possess a unique optimal value. Moreover, the objective function is unimodular in AA and η\eta (detailed in Appendix H and I). Consequently, obtaining the optimal solution for each sub-problem in the iterations enables the alternative iteration algorithm to converge to the global optimum.

V-C2 Complexity Analysis of the Proposed Algorithms

The sub-problem of fixing ηi\eta_{i} and computing A¯i\bar{A}^{*}_{i} based on explicit expression (33) has a complexity of O(1)O(1). The sub-problem of fixing A¯i1\bar{A}^{*}_{i-1} and computing η¯i\bar{\eta}^{*}_{i} has demonstrated to have a unique optimal solution in Lemma 3. By utilizing bisection search and the method outlined in Remark 3, the complexity of this sub-problem is O(log22eλcr2ϵ)O(\log_{2}\frac{2e\lambda cr^{2}}{\epsilon}). The maximum iteration count is nn. Consequently, the proposed algorithms have a computational complexity of O(nlog22eλcr2ϵ)O(n\log_{2}\frac{2e\lambda cr^{2}}{\epsilon}).

Furthermore, we can derive a structural form for the age threshold and update rate pair by leveraging the analytical expression in (46), which is asymptotically optimal.

Theorem 4.

The optimal time-average AoI Δ¯0\bar{\Delta}_{0}^{*} is achieved when the age threshold and the update rate are jointly set as

(A¯,η¯)={(A,1),ifλcr2>𝕎0(12exp(θrαρ));(0,1),otherwise,\left(\bar{A}^{*},\bar{\eta}^{*}\right)=\begin{cases}(A_{*},1),&\text{if}~{}\lambda cr^{2}>\mathbb{W}_{0}\left(\frac{1}{2}\exp\left(-\frac{\theta r^{\alpha}}{\rho}\right)\right);\\ \left(0,1\right),&\text{otherwise},\end{cases} (47)

where a closed-form approximation of AA_{*} is given by

AΩ23+(λcr23)Ω13+(λcr2)2+33Ω13exp(3λcr2Ω13Ω23+λcr2Ω13+(λcr2)2+3θrαρ),A_{*}\approx\frac{\Omega_{*}^{\frac{2}{3}}+(\lambda cr^{2}-3)\Omega_{*}^{\frac{1}{3}}+(\lambda cr^{2})^{2}+3}{3\Omega_{*}^{\frac{1}{3}}\exp\left(\dfrac{-3\lambda cr^{2}\Omega_{*}^{\frac{1}{3}}}{\Omega_{*}^{\frac{2}{3}}+\lambda cr^{2}\Omega_{*}^{\frac{1}{3}}+(\lambda cr^{2})^{2}+3}-\dfrac{\theta r^{\alpha}}{\rho}\right)}, (48)

in which

Ω=(λcr2)3+18λcr2+33(λcr2)4+11(λcr2)21.\Omega_{*}=(\lambda cr^{2})^{3}+18\lambda cr^{2}+3\sqrt{3}\sqrt{(\lambda cr^{2})^{4}+11(\lambda cr^{2})^{2}-1}. (49)

The corresponding time-average AoI is given by (50).

Δ~0={(Ω23+λcr2Ω13+(λcr2)2+3)2+9Ω23(1exp(3λcr2Ω13Ω23+λcr2Ω13+(λcr2)2+3θrαρ))6Ω13(Ω23+λcr2Ω13+(λcr2)2+3)exp(3λcr2Ω13Ω23+λcr2Ω13+(λcr2)2+3θrαρ)+12,ifλcr2>𝕎0(exp(θrαρ)2),exp(λcr2+θrαρ1),otherwise.\tilde{\Delta}_{0}^{*}=\begin{cases}\dfrac{\left(\Omega_{*}^{\frac{2}{3}}+\lambda cr^{2}\Omega_{*}^{\frac{1}{3}}+(\lambda cr^{2})^{2}+3\right)^{2}\!\!+9\Omega_{*}^{\frac{2}{3}}\left(1-\exp\left(\dfrac{-3\lambda cr^{2}\Omega_{*}^{\frac{1}{3}}}{\Omega_{*}^{\frac{2}{3}}\!+\!\lambda cr^{2}\Omega_{*}^{\frac{1}{3}}\!+\!(\lambda cr^{2})^{2}+3}-\dfrac{\theta r^{\alpha}}{\rho}\right)\right)}{6\Omega_{*}^{\frac{1}{3}}\left(\Omega_{*}^{\frac{2}{3}}\!+\!\lambda cr^{2}\Omega_{*}^{\frac{1}{3}}\!+\!(\lambda cr^{2})^{2}\!+\!3\right)\exp\left(\dfrac{-3\lambda cr^{2}\Omega_{*}^{\frac{1}{3}}}{\Omega_{*}^{\frac{2}{3}}+\lambda cr^{2}\Omega_{*}^{\frac{1}{3}}+(\lambda cr^{2})^{2}+3}-\dfrac{\theta r^{\alpha}}{\rho}\right)}+\dfrac{1}{2},&\text{if}~{}\lambda cr^{2}>\mathbb{W}_{0}\left(\frac{\exp\left(-\frac{\theta r^{\alpha}}{\rho}\right)}{2}\right),\\ \exp(\lambda cr^{2}+\theta r^{\alpha}\rho^{-1}),&\text{otherwise}.\end{cases} (50)

 

Proof:

Please see the Appendix K. ∎

We can observe from the above analysis that at the optimal operating region of time-average AoI, the update rate of TSA always equals one while the age threshold varies in accordance with the interference level.

Theorem 4 also allows us to investigate the scaling property of AoI under TSA. More precisely, when we densify the infrastructure (i.e., let λ\lambda\rightarrow\infty), the network becomes interference-limited (namely, the thermal noise is negligible), and the AoI scales as follows.

Corollary 4.

The optimal time-average AoI under TSA protocol scales with the deployment density as

limλcr2Δ¯0,TSAλcr2=e2,\lim_{\lambda cr^{2}\to\infty}\frac{\bar{\Delta}_{0}^{*,\mathrm{TSA}}}{\lambda cr^{2}}=\frac{e}{2}, (51)

in which the optimal age threshold A{A}^{*} satisfies

limλcr2A¯λcr2=e.\lim_{\lambda cr^{2}\to\infty}\frac{\bar{A}^{*}}{\lambda cr^{2}}=e. (52)

In contrast, the optimal time-average AoI under SA protocol scales as555 The optimal mean peak AoI obeys the same scaling law under TSA and SA, and is identical to that of time-average AoI under the SA protocol.

limλcr2Δ¯0,SAλcr2=e.\lim_{\lambda cr^{2}\to\infty}\frac{\bar{\Delta}_{0}^{*,\mathrm{SA}}}{\lambda cr^{2}}=e. (53)

This corollary unveils that the optimal mean peak AoI and time-average AoI under TSA and SA both scale linearly with the interference level, but TSA can reduce the increasing rate of time-average AoI to half of that under conventional SA. Note that the results in Corollary 4 are obtained from the SINR model, and they coincide with the conclusions in [12] developed under the collision model.

Remark 5.

The TSA protocol retains the low computational complexity as well as low signaling overhead characteristics of ALOHA-like protocols. By introducing local AoI comparison (with the age threshold) during the channel contention period and utilizing the AoI information of the receiver carried in the ACK feedback, it forms a closed-loop control system that significantly enhances the network AoI performance.

VI AoI optimization with bi-stable Region

In the previous derivations, we concentrate on optimizing the TSA protocol while putting aside the root distribution of the fixed-point equation in (6). According to the discussions in Section III-B, if λcr2η>4\lambda cr^{2}\eta>4, the solutions to (6) fall in the bi-stable region, i.e., the network can stabilize in two distinct states, typically a high-efficiency state (i.e., psHp_{\mathrm{s}}^{H}) and a low-efficiency state (i.e., psLp_{\mathrm{s}}^{L}), depending on the initial conditions and the dynamics of the system. This section aims to address such an issue.

First of all, we can assess the difference between the two candidates steady state psHp_{\mathrm{s}}^{H} and psLp_{\mathrm{s}}^{L} according to the region described in Section III-B. By examining the cases A=AlA=A_{l} and A=AhA=A_{h}, the feasible region for psHp_{\mathrm{s}}^{H} and psLp_{\mathrm{s}}^{L} can be determined as follows:

psH(exp(θrαρ112+141λcr2η),1),\displaystyle p_{\mathrm{s}}^{H}\in\left(\exp\left(-\tfrac{\theta r^{\alpha}}{\rho}-\tfrac{1}{\frac{1}{2}+\sqrt{\frac{1}{4}-\frac{1}{\lambda cr^{2}\eta}}}\right),~{}1\right), (54)
psL(0,exp(θrαρ112141λcr2η)).\displaystyle p_{\mathrm{s}}^{L}\in\left(0,~{}\exp\left(-\tfrac{\theta r^{\alpha}}{\rho}-\tfrac{1}{\frac{1}{2}-\sqrt{\frac{1}{4}-\frac{1}{\lambda cr^{2}\eta}}}\right)\right). (55)

By taking a ratio between them, as

psHpsL>exp(θrαρ112+141λcr2η)exp(θrαρ112141λcr2η)=exp((λcr2η)24λcr2η).\begin{split}\frac{p_{\mathrm{s}}^{H}}{p_{\mathrm{s}}^{L}}&>\tfrac{\exp\left(-\tfrac{\theta r^{\alpha}}{\rho}-\tfrac{1}{\frac{1}{2}+\sqrt{\frac{1}{4}-\frac{1}{\lambda cr^{2}\eta}}}\right)}{\exp\left(-\tfrac{\theta r^{\alpha}}{\rho}-\tfrac{1}{\frac{1}{2}-\sqrt{\frac{1}{4}-\frac{1}{\lambda cr^{2}\eta}}}\right)}\\ &=\exp\left(\sqrt{(\lambda cr^{2}\eta)^{2}-4\lambda cr^{2}\eta}\right).\end{split} (56)

This inequality demonstrates that in the bi-stable region, if the steady-state point shifts from psHp_{\mathrm{s}}^{H} to psLp_{\mathrm{s}}^{L}, the performance loss increases exponentially with product of the interference level and update rate, i.e., λcr2η\lambda cr^{2}\eta, if the parameters are not properly tuned. For example, at λcr2η=6\lambda cr^{2}\eta=6, the performance loss can be more than 3030 times.

In light of the above challenge, we propose parameter tuning strategies for both mean peak AoI and time-average AoI to enhance the robustness against the bi-stage region situation.

Fig. 3 illustrates an example of the mono-stable and bi-stable regions to facilitate explanations of the optimization process. It also highlights instances where the optimization results from previous discussions fall within these two regions.

Refer to caption
Figure 3: Bistable region \mathcal{B} and monostable region \mathcal{M}.

VI-A Mean Peak AoI

For the mean peak AoI optimization, the optimal parameters pairs (31) intersect with the boundary line of the bi-stable region AhA_{h}. Then, the optimal parameter pairs are bisected by the curve AhA_{h}, with one portion lying in the mono-stable region and the other in the bi-stable region. Therefore, to achieve the optimal mean peak AoI, we only need to select optimal parameters pairs that fall within the mono-stable region. In short, the performance limits of mean peak AoI can be achieved even when the bi-stable region exists; only the tuning region for optimal parameter settings is limited.

The region of optimal update rate η1\eta_{1} when optimal mean peak AoI is achieved in (31) can be expressed as

η1[1λcr2,min{1,4.3509λcr2}].\eta_{1}\in\left[\frac{1}{\lambda cr^{2}},\min\left\{1,\frac{4.3509}{\lambda cr^{2}}\right\}\right]. (57)

The number 4.35094.3509 is the approximate value of λcr2η\lambda cr^{2}\eta when curve (31) intersects with AhA_{h}.

VI-B Time-Average AoI

For the time-average AoI optimization, the optimal performance is influenced by the bi-stable region, as the optimal operating point resides within it. Then, the suboptimal time-average AoI is achieved when A=AhA=A_{h}, which is given by

minη[4λcr2,1]Ah+12+1ηpsA=AhAh+12(1+AhηpsA=Ah).\min_{\eta\in\left[\frac{4}{\lambda cr^{2}},1\right]}~{}~{}\frac{A_{h}+1}{2}+\frac{1}{\eta p^{A=A_{h}}_{\mathrm{s}}}-\frac{A_{h}+1}{2\left(1+A_{h}\eta p^{A=A_{h}}_{\mathrm{s}}\right)}. (58)

This outcome can be intuitively interpreted as follows:

First, due to the unstable nature of the bi-stable region, the network parameters must be fine-tuned to ensure that the system predominantly operates within the monostable region, avoiding the bi-stable region. Therefore, when the bi-stable region exists, the suboptimal parameters should be chosen in the boundary between the mono-stable and bi-stable regions, i.e., AlA_{l} or AhA_{h} (the boundary line falls within the monostable region).

Then, based on Corollary 2, we can determine that, for a given update rate, the transmission success probability on the AlA_{l} curve is lower than the undesired steady-state point within the bi-stable region, and also lower than that on the AhA_{h}, which characterizes an extremely adverse data transmission environment. Therefore, the suboptimal parameters shall be chosen on the boundary of AhA_{h} rather than AlA_{l}.

Finally, upon investigating the time-average AoI performance along the boundary line AhA_{h}, we found that the time-average AoI decreases first and then increases as the update rate η\eta increases, which means the time-average AoI has a unique optimal value in the curve of AhA_{h}. Then, the following gradient information will be incorporated, i.e.,

Δ¯A=Ahη=12Ahη(11ηpsA=Ah(1+AhηpsA=Ah)2)\displaystyle\tfrac{\partial{\bar{\Delta}^{A=A_{h}}}}{\partial{\eta}}=\tfrac{1}{2}\tfrac{\partial{A_{h}}}{\partial{\eta}}\left(1-\tfrac{1-\eta p_{\mathrm{s}}^{A=A_{h}}}{\left(1+A_{h}\eta p_{\mathrm{s}}^{A=A_{h}}\right)^{2}}\right) (59)
(psA=Ah+ηpsA=Ahη)(1(ηpsA=Ah)2Ah(Ah+1)2(1+AhηpsA=Ah)2),\displaystyle-\left(p_{\mathrm{s}}^{A=A_{h}}+\eta\tfrac{\partial{p_{\mathrm{s}}^{A=A_{h}}}}{\partial{\eta}}\right)\left(\tfrac{1}{\left(\eta p_{\mathrm{s}}^{A=A_{h}}\right)^{2}}-\tfrac{A_{h}(A_{h}+1)}{2\left(1+A_{h}\eta p_{\mathrm{s}}^{A=A_{h}}\right)^{2}}\right),

where the gradient Ahη\frac{\partial{A_{h}}}{\partial{\eta}} and psA=Ahη\frac{\partial{p_{\mathrm{s}}^{A=A_{h}}}}{\partial{\eta}} are given as follows

Ahη=1η2exp(112141λcr2η)exp(θrαρ1),\displaystyle\tfrac{\partial{A_{h}}}{\partial{\eta}}=\tfrac{1}{\eta^{2}}\exp{\left(\tfrac{1}{\tfrac{1}{2}-\sqrt{\tfrac{1}{4}-\tfrac{1}{\lambda cr^{2}\eta}}}\right)}\exp{\left(\theta r^{\alpha}\rho^{-1}\right)}, (60)
psA=Ahη=λcr2psA=Ahλcr2Ahη2psA=Ah(1+AhηpsA=Ah)2\displaystyle\tfrac{\partial{p_{\mathrm{s}}^{A=A_{h}}}}{\partial{\eta}}=\tfrac{\lambda cr^{2}p_{\mathrm{s}}^{A=A_{h}}}{\lambda cr^{2}A_{h}\eta^{2}p_{\mathrm{s}}^{A=A_{h}}-\left(1+A_{h}\eta p_{\mathrm{s}}^{A=A_{h}}\right)^{2}} (61)
×(1exp(λcr2η1+AhηpsA=Ah)exp(112141λcr2η)),\displaystyle~{}\times\left(1-\exp\left(\tfrac{-\lambda cr^{2}\eta}{1+A_{h}\eta p_{\mathrm{s}}^{A=A_{h}}}\right)\exp{\left(\tfrac{1}{\tfrac{1}{2}-\sqrt{\tfrac{1}{4}-\tfrac{1}{\lambda cr^{2}\eta}}}\right)}\right),

the corresponding update rate in optimization problem (58) can be solved effectively by numerical method.

VII Numerical Results

In this section, we conduct simulations to verify the accuracy of the developed analysis and then use the numerical evaluations to explore the impact of different network parameters on the AoI performance. Specifically, at the beginning of each simulation run, the locations of source-destination pairs are scattered over a 100×100100\times 100 square area of unit size according to independent PPPs. We add the typical link to the network where the receiver is located at the center of the area. In each time slot, the locations of the links are rearranged independently according to the same PPPs except for the typical one. Each simulation lasts for 10610^{6} time slots. We record the AoI statistics for each time slot of the typical link and then average them to obtain the mean peak and time-average AoI performance.

Refer to caption
Figure 4: Transmission success probability versus age threshold. Network parameters are set as: η=1\eta=1, λ{0.001,0.003,0.005,0.01}\lambda\in\{0.001,0.003,0.005,0.01\}, θ=0dB\theta=0~{}\text{dB}, ρ=20dB\rho=20~{}\text{dB}, α=3.8\alpha=3.8, and r=3r=3.
Refer to caption
(a) A=0A=0 and A=10A=10.
Refer to caption
(b) A=20A=20 and A=50A=50.
Figure 5: AoI performance versus update rate under different age thresholds. Network parameters are set as: λ=0.02\lambda=0.02, θ=0\theta=0 dB, ρ=20dB\rho=20~{}\text{dB}, α=3.8\alpha=3.8, and r=4r=4.

VII-A Transmission Success Probability

Fig. 4 plots the transmission success probability as a function of the age threshold under a variety of spatial deployment densities. We notice that the analytical results and simulations are almost indistinguishable, confirming our derivations in Theorem 1. Moreover, under the same configuration of other network parameters, the larger the spatial deployment density, the worse the transmission success probability due to the surging co-channel interference. On the other hand, the transmission success probability increases with the age threshold, demonstrating the effectiveness of the TSA protocol in interference mitigation, where the gain becomes particularly significant when the spatial deployment density is large.

VII-B Information Freshness versus Update Rate

Fig. 5 depicts the mean peak AoI and time-average AoI as functions of the update rate under different age thresholds. We can see from this figure that an optimal update rate exists that minimizes the mean peak AoI and time-average AoI, whereas the optimal value of the update rate is dependent on the age threshold. When the age threshold is relatively large, the optimal update rate also tends to be higher. This observation corroborates the interchangeable roles of both in mitigating channel contention.

Furthermore, we observe that the performance limits achieved by the TSA protocol in optimizing the mean peak AoI are equivalent to those achieved by the SA protocol, further substantiating this perspective. Yet, we also observe that integrating an age threshold into the SA protocol is beneficial for reducing time-average AoI, while such a threshold must be adjusted appropriately. This is attributed to the fact that, compared with adjusting the update rate, modifying the age threshold not only mitigates channel conflicts but also serves to equalize the intervals of update packet receptions. Consequently, it offers more pronounced gains for higher-order AoI metrics. Another noteworthy observation is that the optimal update rate (with a fixed age threshold) that realizes both the optimal mean peak AoI and the optimal time-average AoI is identical, whatever η|A<1\eta^{*}_{|A}<1 in Fig. 5(a) or η|A=1\eta^{*}_{|A}=1 in Fig. 5(b), corroborating the precision of our theoretical findings.

VII-C Information Freshness versus Age Threshold

Next, Fig. 6 plots the mean peak AoI and time-average AoI as functions of the age threshold under different spatial deployment densities. In this figure, the update rate is set to one (which is often optimal, according to Theorem 4). We can see that the solution given in (44) is almost indistinguishable from the optimality, confirming the accuracy of our analysis. Additionally, we note that the optimal age threshold increases with the spatial deployment density, owing to the more stringent requirement for interference management. Besides, the age threshold to achieve the optimal time-average AoI is higher than that for the optimal mean peak AoI because of the requirement for minimal fluctuations in the data update intervals.

Refer to caption
Figure 6: AoI performance versus age threshold under different deployment densities. Network parameters are set as: η=1\eta=1, θ=0.5\theta=0.5, ρ=20dB\rho=20~{}\text{dB}, α=3.8\alpha=3.8, and r=5r=5.
Refer to caption
Figure 7: TSA performance gain. Network parameters are set as: θ=0\theta=0 dB, ρ=20dB\rho=20~{}\text{dB}, α=3.8\alpha=3.8, and r=3r=3.
Refer to caption
(a) The first moment of the update interval
Refer to caption
(b) The second moment of the update interval
Figure 8: TSA mechanism analysis. Network parameters are set as: θ=0\theta=0 dB, ρ=20dB\rho=20~{}\text{dB}, α=3.8\alpha=3.8, and r=3r=3.

VII-D Scaling Law

Refer to caption
Figure 9: The variance of AoI versus age threshold under different deployment densities. Network parameters are set as: θ=0\theta=0 dB, ρ=20dB\rho=20~{}\text{dB}, α=3.8\alpha=3.8, r=3r=3, and η=1\eta=1.

Fig. 7 illustrates the performance gain of TSA over SA by drawing the optimal time-average AoI as a function of spatial deployment density. Specifically, at each value of λ\lambda, we calculate the optimal age threshold and update rate pair using Algorithm 1 and then compute the optimal time-average AoI. We also plot the results obtained from Theorem 4 in this figure. We can see that the optimal mean peak AoI and time-average AoI under SA and TSA increases almost linearly with the spatial deployment density, whilst the time-average AoI attained by TSA is approximately half of that under SA. This gain is mainly attributed to the advantage of TSA in (aa) equalizing the status updating interval and (bb) mitigating mutual interference, hence benefiting the AoI. Such observation also aligns with the theoretical analysis in Corollary 4. We also notice that the gap between the optimal AoI and that in Theorem 4 quickly diminishes, confirming the accuracy of our theoretical results.

VII-E Higher-Order AoI Statistics

To demystify the essential factors of TSA in reducing AoI, Fig. 8 compares the first and second moments of the update interval (the time interval between two consecutive successful transmissions of the updates) under TSA and SA, where the age threshold and update rate are configured in the AoI-optimal regimes. Fig. 8(a) shows that when the protocol configurations are adjusted on mean peak AoI optimal, TSA and SA result in the same average length of the update interval, leading to equal mean peak AoI. In contrast, if the TSA configuration is set toward optimal time-average AoI, it increases the average length of the update interval (and hence the mean peak AoI), but such a difference decreases as spatial deployment density increases. Fig. 8(b) unveils that the gain from TSA in improving time-average AoI mainly stems from controlling the fluctuation of the update interval and thus maintaining a smaller second moment of the update interval.

In addition to the average of AoI, we plot the variance of this metric [27] in Fig. 9, quantifying the deviation of the AoI performance at each link pair compared to the averaged result of the network. This figure shows that compared to SA (i.e., when A=0A=0), TSA can substantially reduce the AoI variance with an adequately set age threshold (e.g., A=10A=10 when λ=0.05\lambda=0.05, which results in a reduction at the level of an order of magnitude). The reason can be attributed to the fact that the age threshold implicitly stipulates that the sources in good communication environments, namely, those located in sparse areas, shall remain silent for most of the time, since they have a relatively high transmission success probability; in contrast, the source nodes situated in a congested spatial area will be activated more frequently for packet transmissions, because their transmissions often fail. As a result, it reduces the dispersion of AoI across the network, which, in turn, enhances the fairness amongst the source nodes.

VIII Conclusion

In this paper, we have developed an analytical framework to minimize AoI in random access networks by optimizing the TSA status updating protocol. We have derived analytical expressions for the transmission success probability, the mean peak AoI, and the time-average AoI of a typical link, accounting for effects from the update rate, age threshold, channel fading, and interference. Furthermore, we have obtained closed-form solutions for the optimal update rate (resp. age threshold) given a fixed age threshold (resp. update rate); we have also established a closed-form structure for the (almost) optimal update rate and age threshold pair. Our study showed that at the optimal operating point for time-average AoI, the update rate should be set to one while the age threshold varies according to the interference level. The optimization results have also revealed that under the same network configuration, the mean peak AoI and time-average AoI have an identical optimal update rate, yet the optimal threshold for time-average AoI is higher than that for mean peak AoI.

Furthemore, our analysis has identified the regime in which TSA outperforms SA in reducing AoI. Specifically, for the mean peak AoI, the TSA protocol exhibits comparable performance to the traditional slotted ALOHA protocol due to the interchangeable roles of tuning update rate and age threshold in mitigating channel contention. For the time-average AoI, TSA protocol outperforms SA due to its ability to equalize the intervals of update packet reception. The analysis has also disclosed the AoI scaling property, showing that the optimal mean peak AoI and time-average AoI under TSA and SA grow linearly with the deployment density; however, TSA can reduce the growth rate of time-average AoI to half of that under SA, thus, being efficient from an operation perspective.

It should be emphasized that while this study predominantly centers on the generate-at-will model, the proposed framework can also accommodate stochastic arrival patterns. Incorporating stochastic arrivals, exploring the intricate interplay among the age threshold, packet arrival, and channel access in TSA networks presents a compelling topic of discussion. Additionally, research has shown that the TSA protocol offers predictable gains for higher-order AoI metrics. However, the optimal operating point may not necessarily align perfectly with the metrics discussed in this paper. Therefore, utilizing the framework presented in this article to analyze and optimize other high-order AoI metrics with linear growth, as well as AoI metrics associated with nonlinear penalty functions, emerges as a significant direction for subsequent investigations. Moreover, enhancing the efficiency of IoT access in terms of operational costs, such as energy consumption and network lifetime, represents another promising avenue for investigation.

Appendix A Proof of Theorem 1

The transmission success probability can be derived as

(SINR0>θ)=(h00rαi0hi0viXiα+ρ1>θ)\displaystyle\mathbb{P}(\text{SINR}_{0}>\theta)=\mathbb{P}\left(\frac{h_{00}r^{-\alpha}}{\sum_{i\neq 0}h_{i0}v_{i}||X_{i}||^{-\alpha}+\rho^{-1}}>\theta\right)
=(h00>θrα(i0hi0viXiα+ρ1))\displaystyle=\mathbb{P}\left(h_{00}>\theta r^{\alpha}\left(\sum_{i\neq 0}h_{i0}v_{i}||X_{i}||^{-\alpha}+\rho^{-1}\right)\right)
=exp(θrαρ1)i0exp(θrαhi0viXiα)\displaystyle=\exp\left(-\theta r^{\alpha}\rho^{-1}\right)\prod_{i\neq 0}\exp\left(-\theta r^{\alpha}h_{i0}v_{i}||X_{i}||^{-\alpha}\right)
=(a)exp(θrαρ1)i0(11+θrαviXiα)\displaystyle\overset{(a)}{=}\exp\left(-\theta r^{\alpha}\rho^{-1}\right)\prod_{i\neq 0}\left(\frac{1}{1+\theta r^{\alpha}v_{i}||X_{i}||^{-\alpha}}\right)
(b)eθrαρ1i0(1η(Δi(t)>A)+η(Δi(t)>A)1+θrαXiα)\displaystyle\overset{(b)}{\approx}e^{-\theta r^{\alpha}\rho^{-1}}\prod_{i\neq 0}\left(1-\eta\mathbb{P}(\Delta_{i}(t)>A)+\frac{\eta\mathbb{P}(\Delta_{i}(t)>A)}{1+\theta r^{\alpha}||X_{i}||^{-\alpha}}\right)
=eθrαρ1i0(1η(Δi(t)>A))1+Xiα/(θrα)),\displaystyle=e^{-\theta r^{\alpha}\rho^{-1}}\prod_{i\neq 0}\left(1-\frac{\eta\mathbb{P}(\Delta_{i}(t)>A))}{1+||X_{i}||^{\alpha}/\left(\theta r^{\alpha}\right)}\right), (62)

where (a) follows since hijh_{ij} are i.i.d. random variables satisfy hijexp(1)h_{ij}\sim\exp\left(1\right), (b) follows as viv_{i} are independent of each other, and (vi=1)=η(Δi(t)>A)\mathbb{P}(v_{i}=1)=\eta\mathbb{P}(\Delta_{i}(t)>A). The probability that the AoI Δi(t)\Delta_{i}(t) exceed threshold AA can be expressed as

(Δi(t)>A)=1ηps,iA+1ηps,i=11+Aηps,i,\mathbb{P}(\Delta_{i}(t)>A)=\frac{\frac{1}{\eta p_{\mathrm{s},i}}}{A+\frac{1}{\eta p_{\mathrm{s},i}}}=\frac{1}{1+A\eta p_{\mathrm{s},i}}, (63)

With the mobility model, the correlation of the received SINR of each transmitter becomes i.i.d. Thus, we drop the ii index. The spatiotemproal interaction among the queues can be captured by a identification psp_{\mathrm{s}}, which can be expressed as

ps\displaystyle p_{\mathrm{s}} =(SINR0>θ)\displaystyle=\mathbb{P}(\text{SINR}_{0}>\theta)
=eθrαρ1i0(1η1+Aηps1+Xiα/(θrα))\displaystyle=e^{-\theta r^{\alpha}\rho^{-1}}\prod_{i\neq 0}\left(1-\frac{\frac{\eta}{1+A\eta p_{\mathrm{s}}}}{1+||X_{i}||^{\alpha}/(\theta r^{\alpha})}\right)
=eθrαρ1×\displaystyle=e^{-\theta r^{\alpha}\rho^{-1}}\times
exp(λx2[1(1η1+Aηps1+Xiα/(θrα))]𝑑x)\displaystyle~{}~{}~{}\exp\left(-\lambda\int_{x\in\mathbb{R}^{2}}\left[1-\left(1-\frac{\frac{\eta}{1+A\eta p_{\mathrm{s}}}}{1+||X_{i}||^{\alpha}/(\theta r^{\alpha})}\right)\right]dx\right)
=eθrαρ1exp(2πλα0η1+Aηps1+u/θrαu2α1𝑑u)\displaystyle=e^{-\theta r^{\alpha}\rho^{-1}}\exp\left(-\frac{2\pi\lambda}{\alpha}\int_{0}^{\infty}\frac{\frac{\eta}{1+A\eta p_{\mathrm{s}}}}{1+u/\theta r^{\alpha}}u^{\frac{2}{\alpha}-1}du\right)
=exp(λπηr2θ2αΓ(12α)Γ(1+2α)1+Aηpsθrαρ).\displaystyle=\exp\left(-\frac{\lambda\pi\eta r^{2}\theta^{\frac{2}{\alpha}}\Gamma\left(1-\frac{2}{\alpha}\right)\Gamma\left(1+\frac{2}{\alpha}\right)}{1+A\eta p_{\mathrm{s}}}-\frac{\theta r^{\alpha}}{\rho}\right). (64)

The transmission success probability psp_{\mathrm{s}} can then be expressed as (6) in Theorem 1.

Appendix B Root analysis of (6): Corollary 1

We begin by constructing the following auxiliary function

f(ps)=lnpsMN+psK,f(p_{\mathrm{s}})=-\ln p_{\mathrm{s}}-\frac{M}{N+p_{\mathrm{s}}}-K, (65)

where M=λcr2AM=\frac{\lambda cr^{2}}{A}, N=1AηN=\frac{1}{A\eta}, and K=θrbρ1K=\theta r^{b}\rho^{-1}. It is evident that f(ps)=0f(p_{\mathrm{s}})=0 is the equivalent transformation of the fixed-point equation (6), and can be employed to analyze the number of non-zero root for simplicity. Subsequently, we take a derivative of f(ps)f(p_{\mathrm{s}}), that is

f(ps)=ϕ(ps)ps(N+ps)2.f^{{}^{\prime}}(p_{\mathrm{s}})=\frac{\phi(p_{\mathrm{s}})}{p_{\mathrm{s}}(N+p_{\mathrm{s}})^{2}}. (66)

The numerator of f(ps)f^{{}^{\prime}}(p_{\mathrm{s}}), i.e., ϕ(ps)\phi(p_{\mathrm{s}}) is a quadratic function, which can be expressed as

ϕ(ps)=(ps+NM2)2+M24MN.\phi(p_{\mathrm{s}})=-\left(p_{\mathrm{s}}+N-\frac{M}{2}\right)^{2}+\frac{M^{2}}{4}-MN. (67)

The following lemma demonstrates that the number of non-zero roots of (66) is determined by the zero-points of roots of (67).

Lemma 5.

f(ps)f(p_{\mathrm{s}}) exhibits three non-zero roots of 0<psL<psM<psH<10<p_{\mathrm{s}}^{L}<p_{\mathrm{s}}^{M}<p_{\mathrm{s}}^{H}<1, if and only if ϕ(ps)=0\phi(p_{\mathrm{s}})=0 has two non-zero roots 0<ps,1<ps,2<10<p^{{}^{\prime}}_{\mathrm{s},1}<p^{{}^{\prime}}_{\mathrm{s},2}<1 with f(ps,1)<0f(p^{{}^{\prime}}_{\mathrm{s},1})<0 and f(ps,2)>0f(p^{{}^{\prime}}_{\mathrm{s},2})>0; in the cases where f(ps,1)=0f(p^{{}^{\prime}}_{\mathrm{s},1})=0 or f(ps,2)=0f(p^{{}^{\prime}}_{\mathrm{s},2})=0; f(ps)=0f(p_{\mathrm{s}})=0 has two non-zero roots, otherwise, f(ps)=0f(p_{\mathrm{s}})=0 will solely have one non-zero root 0<ps<10<p_{\mathrm{s}}<1.

Proof:

Due to the fact limps0f(ps)>0\underset{p_{\mathrm{s}}\to 0}{\lim}f(p_{\mathrm{s}})>0, f(1)=MN+1K<0f(1)=-\frac{M}{N+1}-K<0, f(ps)f(p_{\mathrm{s}}) is a continuous function, according to the zero-point theorem, f(ps)f(p_{\mathrm{s}}) possesses a non-zero root at least. Furthermore, the denominator of the function f(ps)f^{{}^{\prime}}(p_{\mathrm{s}}) is greater than zero for ps(0,1]p_{\mathrm{s}}\in(0,1]. Consequently, the number of non-zero roots of ϕ(ps)=0\phi(p_{\mathrm{s}})=0 for ps(0,1]p_{\mathrm{s}}\in(0,1] determines the number of non-zero roots of f(ps)=0f^{{}^{\prime}}(p_{\mathrm{s}})=0 for ps(0,1]p_{\mathrm{s}}\in(0,1], and we will only consider the root of ϕ(ps)\phi(p_{\mathrm{s}}) in the following scenarios.

(1) Assume that ϕ(ps)=0\phi(p_{\mathrm{s}})=0 has no non-zero root when ps(0,1]p_{\mathrm{s}}\in(0,1], it follows that in this case ϕ(ps)<0\phi(p_{\mathrm{s}})<0 for ps(0,1]p_{\mathrm{s}}\in(0,1]. Consequently, f(ps)<0f^{{}^{\prime}}(p_{\mathrm{s}})<0 and f(ps)f(p_{\mathrm{s}}) decreases monotonically when ps(0,1]p_{\mathrm{s}}\in(0,1], Hence, f(ps)f(p_{\mathrm{s}}) has only one non-zero root in this scenario.

(2) Assume that ϕ(ps)=0\phi(p_{\mathrm{s}})=0 has one non-zero root 0<ps10<p_{\mathrm{s}}^{{}^{\prime}}\leq 1, we have the following scenarios:

(2.a) If ϕ(ps)=0\phi(p_{\mathrm{s}})=0 has two roots when psp_{\mathrm{s}}\in\mathbb{R}, and one of them in range ps(0,1)p_{\mathrm{s}}\in(0,1), we denote the root in (0,1)(0,1) as psp_{\mathrm{s}}^{{}^{\prime}}. Consequently, when ps(0,ps)p_{\mathrm{s}}\in(0,p_{\mathrm{s}}^{{}^{\prime}}), f(ps)f(p_{\mathrm{s}}) monotonically decrease and ps(ps,1)p_{\mathrm{s}}\in(p_{\mathrm{s}}^{{}^{\prime}},1), f(ps)f(p_{\mathrm{s}}) monotonically increase. Since we have f(1)<0f(1)<0, thus f(ps)f(p_{\mathrm{s}}) only has one non-zero root.

(2.b) If ϕ(ps)=0\phi(p_{\mathrm{s}})=0 has one root when psp_{\mathrm{s}}\in\mathbb{R}, and this root in range ps(0,1)p_{\mathrm{s}}\in(0,1), then f(ps)f(p_{\mathrm{s}}) monotonically decrease when ps(0,1)p_{\mathrm{s}}\in(0,1), so f(ps)f(p_{\mathrm{s}}) only has one non-zero root. Combining two cases, when ϕ(ps)=0\phi(p_{\mathrm{s}})=0 has one non-zero root, f(ps)=0f(p_{\mathrm{s}})=0 only has one non-zero root.

Combining these two cases, when ϕ(ps)=0\phi(p_{\mathrm{s}})=0 has one non-zero root, f(ps)=0f(p_{\mathrm{s}})=0 only has one non-zero root.

(3) Assume that ϕ(ps)=0\phi(p_{\mathrm{s}})=0 has two non-zero roots 0<ps,1<ps,2<10<p^{{}^{\prime}}_{\mathrm{s},1}<p^{{}^{\prime}}_{\mathrm{s},2}<1, then we have

(3.a) If ps,2=1p_{\mathrm{s},2}^{{}^{\prime}}=1, then ϕ(ps)<0\phi(p_{\mathrm{s}})<0 for ps(0,ps,1)p_{\mathrm{s}}\in(0,p_{\mathrm{s},1}^{{}^{\prime}}) and ϕ(ps)>0\phi(p_{\mathrm{s}})>0 for ps(ps,1,1)p_{\mathrm{s}}\in(p_{\mathrm{s},1}^{{}^{\prime}},1). Consequently, f(ps)<0f^{{}^{\prime}}(p_{\mathrm{s}})<0 for ps(0,ps,1)p_{\mathrm{s}}\in(0,p_{\mathrm{s},1}^{{}^{\prime}}) and f(ps)>0f^{{}^{\prime}}(p_{\mathrm{s}})>0 for ps(ps,1,1)p_{\mathrm{s}}\in(p_{\mathrm{s},1}^{{}^{\prime}},1), indicating that f(ps)f(p_{\mathrm{s}}) monotonically decreases for ps(0,ps,1)p_{\mathrm{s}}\in(0,p_{\mathrm{s},1}^{{}^{\prime}}), and increases for ps(ps,1,1]p_{\mathrm{s}}\in(p_{\mathrm{s},1}^{{}^{\prime}},1]. Since f(1)<0f(1)<0, we can conclude that in this case, f(ps)=0f(p_{\mathrm{s}})=0 only one non-zero root 0<psH<10<p_{\mathrm{s}}^{H}<1.

(3.b) If ps,2<1p^{{}^{\prime}}_{\mathrm{s},2}<1, then ϕ(ps)<0\phi(p_{\mathrm{s}})<0 for ps(0,ps,1)(ps,2,1)p_{\mathrm{s}}\in(0,p^{{}^{\prime}}_{\mathrm{s},1})\cup(p^{{}^{\prime}}_{\mathrm{s},2},1), and ϕ(ps)>0\phi(p_{\mathrm{s}})>0 for ps(ps,1,ps,2)p_{\mathrm{s}}\in(p^{{}^{\prime}}_{\mathrm{s},1},p^{{}^{\prime}}_{\mathrm{s},2}). As a result, f(ps)<0f^{{}^{\prime}}(p_{\mathrm{s}})<0 for ps(0,ps,1)(ps,2,1)p_{\mathrm{s}}\in(0,p^{{}^{\prime}}_{\mathrm{s},1})\cup(p^{{}^{\prime}}_{\mathrm{s},2},1), and f(ps)>0f^{{}^{\prime}}(p_{\mathrm{s}})>0 for ps(ps,1,ps,2)p_{\mathrm{s}}\in(p^{{}^{\prime}}_{\mathrm{s},1},p^{{}^{\prime}}_{\mathrm{s},2}), indicating that f(ps)f(p_{\mathrm{s}}) monotonically decreases for ps(0,ps,1)(ps,2,1)p_{\mathrm{s}}\in(0,p^{{}^{\prime}}_{\mathrm{s},1})\cup(p^{{}^{\prime}}_{\mathrm{s},2},1), and increases for p(p1,p2)p\in(p^{{}^{\prime}}_{1},p^{{}^{\prime}}_{2}).Then, we have if f(ps,1)>0f(p_{\mathrm{s},1}^{{}^{\prime}})>0 or f(ps,2)<0f(p_{\mathrm{s},2}^{{}^{\prime}})<0, f(ps)=0f(p_{\mathrm{s}})=0 has one root 0<psH10<p_{\mathrm{s}}^{H}\leq 1; when f(ps,1)=0f(p^{{}^{\prime}}_{\mathrm{s},1})=0 or f(ps,2)=0f(p^{{}^{\prime}}_{\mathrm{s},2})=0; f(ps)=0f(p_{\mathrm{s}})=0 has two non-zero roots; otherwise, f(ps)=0f(p_{\mathrm{s}})=0 has three non-zero roots 0<psL<psM<psH10<p_{\mathrm{s}}^{L}<p_{\mathrm{s}}^{M}<p_{\mathrm{s}}^{H}\leq 1 in which f(ps,1)<0f(p_{\mathrm{s},1}^{{}^{\prime}})<0 and f(ps,2)>0f(p_{\mathrm{s},2}^{{}^{\prime}})>0. ∎

Lemma 4 further concludes the specific condition for ϕ(ps)=0\phi(p_{\mathrm{s}})=0 has two non-zero roots 0<ps,1<ps,2<10<p^{{}^{\prime}}_{\mathrm{s},1}<p^{{}^{\prime}}_{\mathrm{s},2}<1 with f(ps,1)<0f(p^{{}^{\prime}}_{\mathrm{s},1})<0 and f(ps,2)>0f(p^{{}^{\prime}}_{\mathrm{s},2})>0.

Lemma 6.

ϕ(ps)=0\phi(p_{\mathrm{s}})=0 has two non-zero roots 0<ps,1<ps,2<10<p^{{}^{\prime}}_{\mathrm{s},1}<p^{{}^{\prime}}_{\mathrm{s},2}<1 with f(ps,1)<0f(p^{{}^{\prime}}_{\mathrm{s},1})<0 and f(ps,2)>0f(p^{{}^{\prime}}_{\mathrm{s},2})>0 if and only if λcr2>4η\lambda cr^{2}>\frac{4}{\eta}, and Al<A<AhA_{l}<A<A_{h}, where AlA_{l} and AhA_{h} are respectively given by

Al=λcr2(12+141λcr2η)1ηexp(θrαρ112+141λcr2η);\displaystyle A_{l}=\frac{\lambda cr^{2}\left(\frac{1}{2}+\sqrt{\frac{1}{4}-\frac{1}{\lambda cr^{2}\eta}}\right)-\frac{1}{\eta}}{\exp\left(-\frac{\theta r^{\alpha}}{\rho}-\frac{1}{\frac{1}{2}+\sqrt{\frac{1}{4}-\frac{1}{\lambda cr^{2}\eta}}}\right)}; (68)
Ah=λcr2(12141λcr2η)1ηexp(θrαρ112141λcr2η).\displaystyle A_{h}=\frac{\lambda cr^{2}\left(\frac{1}{2}-\sqrt{\frac{1}{4}-\frac{1}{\lambda cr^{2}\eta}}\right)-\frac{1}{\eta}}{\exp\left(-\frac{\theta r^{\alpha}}{\rho}-\frac{1}{\frac{1}{2}-\sqrt{\frac{1}{4}-\frac{1}{\lambda cr^{2}\eta}}}\right)}. (69)
Proof:

Considering the properties of quadratic function, ϕ(ps)=0\phi(p_{\mathrm{s}})=0 has two non-zero roots when limps0ϕ(ps)<0\underset{p_{\mathrm{s}}\to 0}{\lim}~{}\phi(p_{\mathrm{s}})<0, limps1ϕ(ps)<0\underset{p_{\mathrm{s}}\to 1}{\lim}~{}\phi(p_{\mathrm{s}})<0 and the peak value of ϕ(ps)\phi(p_{\mathrm{s}}) is larger than zero, that is:

limps0ϕ(ps)=N2<0\displaystyle~{}\underset{p_{\mathrm{s}}\to 0}{\lim}~{}\phi(p_{\mathrm{s}})=-N^{2}<0 (70)
limps1ϕ(ps)=(1+NM2)2+M24MN<0\displaystyle~{}\underset{p_{\mathrm{s}}\to 1}{\lim}~{}\phi(p_{\mathrm{s}})=-\left(1+N-\frac{M}{2}\right)^{2}+\frac{M^{2}}{4}-MN<0 (71)
maxps(0,1)ϕ(ps)=M24MN>0\displaystyle\underset{p_{\mathrm{s}}\in(0,1)}{\text{max}}~{}\phi(p_{\mathrm{s}})=\frac{M^{2}}{4}-MN>0 (72)

Combining (70), (71) and (72), we have

4η<λcr2<min{A(1Aη+1)2,2A+2η}\frac{4}{\eta}<\lambda cr^{2}<\min\left\{A\left(\frac{1}{A\eta}+1\right)^{2},2A+\frac{2}{\eta}\right\} (73)

For 4η<2η+2A\frac{4}{\eta}<\frac{2}{\eta}+2A, we have 1Aη2A<0\frac{1}{A\eta^{2}}-A<0 . Then,

A(1Aη+1)2(2A+2η)=1Aη2A<0A\left(\frac{1}{A\eta}+1\right)^{2}-\left(2A+\frac{2}{\eta}\right)=\frac{1}{A\eta^{2}}-A<0 (74)

which means (73) can be simplified as

4η<λcr2<A(1Aη+1)2\frac{4}{\eta}<\lambda cr^{2}<A\left(\frac{1}{A\eta}+1\right)^{2} (75)

Then, we consider the condition that f(ps,1)<0f(p^{{}^{\prime}}_{\mathrm{s},1})<0 and f(ps,2)>0f(p^{{}^{\prime}}_{\mathrm{s},2})>0, i.e.,

f(ps,1)\displaystyle f(p^{{}^{\prime}}_{\mathrm{s},1}) =MM2M24MNln(M2NM24MN)\displaystyle=-\frac{M}{\frac{{M}}{2}-\sqrt{\frac{{M}^{2}}{4}-MN}}-\ln\left(\frac{M}{2}-N-\sqrt{\frac{{M}^{2}}{4}-MN}\right)
θrαρ1<0,\displaystyle-\theta r^{\alpha}\rho^{-1}<0, (76)

and

f(ps,2)\displaystyle f(p^{{}^{\prime}}_{\mathrm{s},2}) =MM2+M24MNln(M2N+M24MN)\displaystyle=-\frac{M}{\frac{M}{2}+\sqrt{\frac{{M}^{2}}{4}-MN}}-\ln\left(\frac{M}{2}-N+\sqrt{\frac{{M}^{2}}{4}-MN}\right)
θrαρ1>0,\displaystyle-\theta r^{\alpha}\rho^{-1}>0, (77)

By solving the inequalities (B) and (B), we then have

A>Al=λcr2(12+141λcr2η)1ηexp(θrαρ112+141λcr2η)A>A_{l}=\frac{\lambda cr^{2}\left(\frac{1}{2}+\sqrt{\frac{1}{4}-\frac{1}{\lambda cr^{2}\eta}}\right)-\frac{1}{\eta}}{\exp\left(-\frac{\theta r^{\alpha}}{\rho}-\frac{1}{\frac{1}{2}+\sqrt{\frac{1}{4}-\frac{1}{\lambda cr^{2}\eta}}}\right)} (78)

and

A<Ah=λcr2(12141λcr2η)1ηexp(θrαρ112141λcr2η).A<A_{h}=\frac{\lambda cr^{2}\left(\frac{1}{2}-\sqrt{\frac{1}{4}-\frac{1}{\lambda cr^{2}\eta}}\right)-\frac{1}{\eta}}{\exp\left(-\frac{\theta r^{\alpha}}{\rho}-\frac{1}{\frac{1}{2}-\sqrt{\frac{1}{4}-\frac{1}{\lambda cr^{2}\eta}}}\right)}. (79)

Then, by solving the inequality λcR2<A(1Aη+1)2\lambda cR^{2}<A(\frac{1}{A\eta}+1)^{2}, we have

A>λcr2(12+141λcr2η)1η,A>\lambda cr^{2}\left(\frac{1}{2}+\sqrt{\frac{1}{4}-\frac{1}{\lambda cr^{2}\eta}}\right)-\frac{1}{\eta}, (80)

or

A<λcr2(12141λcr2η)1ηA<\lambda cr^{2}\left(\frac{1}{2}-\sqrt{\frac{1}{4}-\frac{1}{\lambda cr^{2}\eta}}\right)-\frac{1}{\eta} (81)

When the condition Al<A<AhA_{l}<A<A_{h} is satisfied, the inequality (80) must be satisfied, and the inequality (81) will not be satisfied. Then, the condition (75) can be simplified as λcr2>4η\lambda cr^{2}>\frac{4}{\eta}. ∎

Consequently, Corollary 1 can be proved by combining Lemma 5 and lemma 6.

Appendix C Proof of Monotony of psLp^{L}_{\mathrm{s}} and psHp^{H}_{\mathrm{s}}: Corollary 2

By taking a partial derivative of psp_{\mathrm{s}} with respect to AA, η\eta, and λ\lambda, we have

psA=λcr2η2ps2(1+Aηps)2λcr2Aη2ps=λcr2ps2A21ϕ(ps);\displaystyle\frac{\partial{p_{\mathrm{s}}}}{\partial{A}}=\frac{\lambda cr^{2}\eta^{2}p_{\mathrm{s}}^{2}}{(1+A\eta p_{\mathrm{s}})^{2}-\lambda cr^{2}A\eta^{2}p_{\mathrm{s}}}=\frac{-\lambda cr^{2}p_{\mathrm{s}}^{2}}{A^{2}}\frac{1}{\phi(p_{\mathrm{s}})}; (82)
psη=λcr2ps(1+Aηps)2λcr2Aη2ps=λcr2psA2η21ϕ(ps);\displaystyle\frac{\partial{p_{\mathrm{s}}}}{\partial{\eta}}=\frac{-\lambda cr^{2}p_{\mathrm{s}}}{(1+A\eta p_{\mathrm{s}})^{2}-\lambda cr^{2}A\eta^{2}p_{\mathrm{s}}}=\frac{\lambda cr^{2}p_{\mathrm{s}}}{A^{2}\eta^{2}}\frac{1}{\phi(p_{\mathrm{s}})}; (83)
psλ=cr2ηps(1+Aηps)(1+Aηps)2λcr2Aη2ps=cr2ηps(1+Aηps)A2η21ϕ(ps),\displaystyle\frac{\partial{p_{\mathrm{s}}}}{\partial{\lambda}}=\frac{-cr^{2}\eta p_{\mathrm{s}}(1+A\eta p_{\mathrm{s}})}{(1+A\eta p_{\mathrm{s}})^{2}-\lambda cr^{2}A\eta^{2}p_{\mathrm{s}}}=\frac{cr^{2}\eta p_{\mathrm{s}}(1+A\eta p_{\mathrm{s}})}{A^{2}\eta^{2}}\frac{1}{\phi(p_{\mathrm{s}})}, (84)

where ϕ(ps)\phi(p_{\mathrm{s}}) is consistence with equation (67). Then we prove ϕ(psL)\phi(p_{\mathrm{s}}^{L}) and ϕ(psH)\phi(p_{\mathrm{s}}^{H}) both smaller than zero.

  • (1)

    If ϕ(ps)=0\phi(p_{\mathrm{s}})=0 has no root for ps(0,1]p_{\mathrm{s}}\in(0,1], then ϕ(ps)<0\phi(p_{\mathrm{s}})<0 for ps(0,1]p_{\mathrm{s}}\in(0,1]. In this case, (6) has one-zero root psHp_{\mathrm{s}}^{H}, and we have ϕ(psH)<0\phi(p_{\mathrm{s}}^{H})<0.

  • (2)

    If ϕ(ps)=0\phi(p_{\mathrm{s}})=0 has one root for 0<ps,110<p^{{}^{\prime}}_{\mathrm{s},1}\leq 1, we have:

    • (2.a)

      If ϕ(ps)=0\phi(p_{\mathrm{s}})=0 has two roots when psp_{\mathrm{s}}\in\mathbb{R}, and one of them in range ps(0,1]p_{\mathrm{s}}\in(0,1], then ϕ(ps)<0\phi(p_{\mathrm{s}})<0 for ps(0,ps,1)p_{\mathrm{s}}\in(0,p^{{}^{\prime}}_{\mathrm{s},1}) and ϕ(ps)>0\phi(p_{\mathrm{s}})>0 for ps(ps,1,1)p_{\mathrm{s}}\in(p^{{}^{\prime}}_{\mathrm{s},1},1). The fixed point equation (6) has one-zero root psHp_{\mathrm{s}}^{H}, and psH<ps,1p_{\mathrm{s}}^{H}<p^{{}^{\prime}}_{\mathrm{s},1}, we then have ϕ(psH)<0\phi(p_{\mathrm{s}}^{H})<0.

    • (2.b)

      If ϕ(ps)=0\phi(p_{\mathrm{s}})=0 has one root when psp_{\mathrm{s}}\in\mathbb{R}, and this root in range ps(0,1)p_{\mathrm{s}}\in(0,1), then ϕ(ps)<0\phi(p_{\mathrm{s}})<0 for ps(0,1]p_{\mathrm{s}}\in(0,1]. (6) has one-zero root ϕ(psH)\phi(p_{\mathrm{s}}^{H}). We have ϕ(psH)<0\phi(p_{\mathrm{s}}^{H})<0.

  • (3)

    If ϕ(ps)=0\phi(p_{\mathrm{s}})=0 has two roots that 0<ps,1<ps,2<10<p^{{}^{\prime}}_{\mathrm{s},1}<p^{{}^{\prime}}_{\mathrm{s},2}<1, two scenarios as follows:

    • (3.a)

      If ps,2=1p^{{}^{\prime}}_{\mathrm{s},2}=1, then ϕ(ps)<0\phi(p_{\mathrm{s}})<0 for ps(0,ps,1)p_{\mathrm{s}}\in(0,p^{{}^{\prime}}_{\mathrm{s},1}) and ϕ(ps)>0\phi(p_{\mathrm{s}})>0 for ps(ps,1,1)p_{\mathrm{s}}\in(p^{{}^{\prime}}_{\mathrm{s},1},1), (6) has one root psH<ps,1p_{\mathrm{s}}^{H}<p^{{}^{\prime}}_{\mathrm{s},1}. We then have ϕ(psH)<0\phi(p^{H}_{\mathrm{s}})<0.

    • (3.b)

      If ps,2<1p^{{}^{\prime}}_{\mathrm{s},2}<1, then ϕ(ps)<0\phi(p_{\mathrm{s}})<0 for ps(0,ps,1)(ps,2,1)p_{\mathrm{s}}\in(0,p^{{}^{\prime}}_{\mathrm{s},1})\cup(p^{{}^{\prime}}_{\mathrm{s},2},1), and ϕ(ps)>0\phi(p_{\mathrm{s}})>0 for ps(ps,1,ps,2)p_{\mathrm{s}}\in(p^{{}^{\prime}}_{\mathrm{s},1},p^{{}^{\prime}}_{\mathrm{s},2}). In this case, (6) may has one root psH(0,ps,1)(ps,2,1)p_{\mathrm{s}}^{H}\in(0,p^{{}^{\prime}}_{\mathrm{s},1})\cup(p^{{}^{\prime}}_{\mathrm{s},2},1), or three roots ϕ(psL)<ps,1<ϕ(psM)<ps,2<ϕ(psH)\phi(p_{\mathrm{s}}^{L})<p^{{}^{\prime}}_{\mathrm{s},1}<\phi(p_{\mathrm{s}}^{M})<p^{{}^{\prime}}_{\mathrm{s},2}<\phi(p_{\mathrm{s}}^{H}). We then have ϕ(psL)<0\phi(p_{\mathrm{s}}^{L})<0 and ϕ(psH)<0\phi(p_{\mathrm{s}}^{H})<0.

In summary, ϕ(psL)<0\phi(p_{\mathrm{s}}^{L})<0 and ϕ(psH)<0\phi(p_{\mathrm{s}}^{H})<0 can be proved by combining above cases. Thus, it can be obtained from (82) – (84) that psA>0\frac{\partial{p_{\mathrm{s}}}}{\partial{A}}>0, psη<0\frac{\partial{p_{\mathrm{s}}}}{\partial{\eta}}<0 and psλ<0\frac{\partial{p_{\mathrm{s}}}}{\partial{\lambda}}<0 at both the steady states psLp_{\mathrm{s}}^{L} and psHp_{\mathrm{s}}^{H}. Therefore, the denominator of derivation satisfied the following inequality:

(1+Aηps)2λcr2Aη2ps>0.(1+A\eta p_{\mathrm{s}})^{2}-\lambda cr^{2}A\eta^{2}p_{\mathrm{s}}>0. (85)

Appendix D Convergence of (8): Corollary 3

By examining the cases A=AlA=A_{l} and A=AhA=A_{h}, the feasible region for psHp_{\mathrm{s}}^{H} and psLp_{\mathrm{s}}^{L} can be determined as follows:

psH(exp(θrαρ112+141λcr2η),1),\displaystyle p_{\mathrm{s}}^{H}\in\left(\exp\left(-\tfrac{\theta r^{\alpha}}{\rho}-\tfrac{1}{\frac{1}{2}+\sqrt{\frac{1}{4}-\frac{1}{\lambda cr^{2}\eta}}}\right),~{}1\right), (86)
psL(0,exp(θrαρ112141λcr2η)).\displaystyle p_{\mathrm{s}}^{L}\in\left(0,~{}\exp\left(-\tfrac{\theta r^{\alpha}}{\rho}-\tfrac{1}{\frac{1}{2}-\sqrt{\frac{1}{4}-\frac{1}{\lambda cr^{2}\eta}}}\right)\right). (87)

Without loss of generality, we consider psHp_{\mathrm{s}}^{H}, which exists when (1) λcr2η4\lambda cr^{2}\eta\leq 4 and (2) λcr2η>4\lambda cr^{2}\eta>4 and A>AlA>A_{l}. For A>0A>0, we assume the initial iterative value ps,1H=1p_{\mathrm{s},1}^{H}=1. Then, we prove the sequence of {ps,1H,ps,2H,,ps,nH,ps,n+1H,}\{p_{\mathrm{s},1}^{H},p_{\mathrm{s},2}^{H},...,p_{\mathrm{s},n}^{H},p_{\mathrm{s},n+1}^{H},...\} is a decreasing sequence. For ease of exposition, let us denote by

f(psH)=exp(λcr2η1+AηpsHθrαρ).f(p_{\mathrm{s}}^{H})=\exp\left(-\frac{\lambda cr^{2}\eta}{1+A\eta p_{\mathrm{s}}^{H}}-\frac{\theta r^{\alpha}}{\rho}\right). (88)

Note that the function f(psH)f(p^{H}_{\mathrm{s}}) increases as the increasing psHp^{H}_{\mathrm{s}}. We then use the mathematical indication method, which is

ps,2H\displaystyle p_{\mathrm{s},2}^{H} =exp(λcr2η1+Aηps,1Hθrαρ)\displaystyle=\exp\left(-\frac{\lambda cr^{2}\eta}{1+A\eta p_{\mathrm{s},1}^{H}}-\frac{\theta r^{\alpha}}{\rho}\right)
=exp(λcr2η1+Aηθrαρ)<1=ps,1H,\displaystyle=\exp\left(-\frac{\lambda cr^{2}\eta}{1+A\eta}-\frac{\theta r^{\alpha}}{\rho}\right)<1=p_{\mathrm{s},1}^{H}, (89)

and

ps,3H\displaystyle p_{\mathrm{s},3}^{H} =exp(λcr2η1+Aηps,2Hθrαρ)\displaystyle=\exp\left(-\frac{\lambda cr^{2}\eta}{1+A\eta p_{\mathrm{s},2}^{H}}-\frac{\theta r^{\alpha}}{\rho}\right)
<exp(λcr2η1+Aηps,1Hθrαρ)=ps,2H.\displaystyle<\exp\left(-\frac{\lambda cr^{2}\eta}{1+A\eta p_{\mathrm{s},1}^{H}}-\frac{\theta r^{\alpha}}{\rho}\right)=p_{\mathrm{s},2}^{H}. (90)

Then, suppose ps,nH<ps,n1Hp_{\mathrm{s},n}^{H}<p_{\mathrm{s},n-1}^{H} holds for nn, we have

ps,n+1H\displaystyle p_{\mathrm{s},n+1}^{H} =exp(λcr2η1+Aηps,nHθrαρ)\displaystyle=\exp\left(-\frac{\lambda cr^{2}\eta}{1+A\eta p_{\mathrm{s},n}^{H}}-\frac{\theta r^{\alpha}}{\rho}\right) (91)
<exp(λcr2η1+Aηps,n1Hθrαρ)=ps,nH.\displaystyle<\exp\left(-\frac{\lambda cr^{2}\eta}{1+A\eta p_{\mathrm{s},n-1}^{H}}-\frac{\theta r^{\alpha}}{\rho}\right)=p_{\mathrm{s},n}^{H}. (92)

As such, the sequence of {ps,1H,ps,2H,,ps,nH,ps,n+1H,}\{p_{\mathrm{s},1}^{H},p_{\mathrm{s},2}^{H},...,p_{\mathrm{s},n}^{H},p_{\mathrm{s},n+1}^{H},...\} is a decrease sequence when ps,1H=1p_{\mathrm{s},1}^{H}=1.

With a similar approach, we can prove that the root psHp_{\mathrm{s}}^{H} is a lower bound of the sequence {ps,1H,ps,2H,,ps,nH,ps,n+1H,}\{p_{\mathrm{s},1}^{H},p_{\mathrm{s},2}^{H},...,p_{\mathrm{s},n}^{H},p_{\mathrm{s},n+1}^{H},...\}. First, we have psH<ps,1H=1p_{\mathrm{s}}^{H}<p_{\mathrm{s},1}^{H}=1. Then, we have

ps,2H\displaystyle p_{\mathrm{s},2}^{H} =exp(λcr2η1+Aηps,1Hθrαρ)\displaystyle=\exp\left(-\frac{\lambda cr^{2}\eta}{1+A\eta p_{\mathrm{s},1}^{H}}-\frac{\theta r^{\alpha}}{\rho}\right)
>exp(λcr2η1+AηpsHθrαρ)=psH\displaystyle>\exp\left(-\frac{\lambda cr^{2}\eta}{1+A\eta p_{\mathrm{s}}^{H}}-\frac{\theta r^{\alpha}}{\rho}\right)=p_{\mathrm{s}}^{H} (93)

Then, suppose ps,nH>psHp_{\mathrm{s},n}^{H}>p_{\mathrm{s}}^{H} holds for nn, we have

ps,n+1H\displaystyle p_{\mathrm{s},n+1}^{H} =exp(λcr2η1+Aηps,nHθrαρ)\displaystyle=\exp\left(-\frac{\lambda cr^{2}\eta}{1+A\eta p_{\mathrm{s},n}^{H}}-\frac{\theta r^{\alpha}}{\rho}\right)
>exp(λcr2η1+AηpsHθrαρ)=psH.\displaystyle>\exp\left(-\frac{\lambda cr^{2}\eta}{1+A\eta p_{\mathrm{s}}^{H}}-\frac{\theta r^{\alpha}}{\rho}\right)=p_{\mathrm{s}}^{H}. (94)

This completes the proof that psHp_{\mathrm{s}}^{H} is a lower bound of the sequence {ps,1H,ps,2H,,ps,nH,ps,n+1H,}\{p_{\mathrm{s},1}^{H},p_{\mathrm{s},2}^{H},...,p_{\mathrm{s},n}^{H},p_{\mathrm{s},n+1}^{H},...\}.

The proof of convergence of psHp_{\mathrm{s}}^{H} and its convergence rate is shown as follows:
(1) When λcr2η4\lambda cr^{2}\eta\leq 4, the steady-state point psHp_{\mathrm{s}}^{H} exists, and we have

|ps,n+1Hps,nH||ps,nHps,n1H|\displaystyle\frac{|p_{\mathrm{s},n+1}^{H}-p_{\mathrm{s},n}^{H}|}{\big{|}p_{\mathrm{s},n}^{H}-p_{\mathrm{s},n-1}^{H}\big{|}}
=exp(θrαρ)|exp(λcr2η1+Aηps,nH)exp(λcr2η1+Aηps,n1H)||ps,nHps,n1H|\displaystyle=\exp{\left(-\frac{\theta r^{\alpha}}{\rho}\right)}\frac{|\exp{\left(-\frac{\lambda cr^{2}\eta}{1+A\eta p_{\mathrm{s},n}^{H}}\right)}-\exp{\left(-\frac{\lambda cr^{2}\eta}{1+A\eta p_{\mathrm{s},n-1}^{H}}\right)}|}{|p_{\mathrm{s},n}^{H}-p_{\mathrm{s},n-1}^{H}|}
=(a)|ps,nHps,n1H|exp(λcr2η1+Aηξθrαρ)λcr2Aη2(1+Aηξ)2|ps,nHps,n1H|\displaystyle\overset{(a)}{=}\frac{|p_{\mathrm{s},n}^{H}-p_{\mathrm{s},n-1}^{H}|\exp\left(-\frac{\lambda cr^{2}\eta}{1+A\eta\xi}-\frac{\theta r^{\alpha}}{\rho}\right)\frac{\lambda cr^{2}A\eta^{2}}{\left(1+A\eta\xi\right)^{2}}}{|p_{\mathrm{s},n}^{H}-p_{\mathrm{s},n-1}^{H}|}
<(b)λcr2Aη2ξ(1+Aηξ)2=λcr2Aη2A2η2ξ+1ξ+2Aη\displaystyle\overset{(b)}{<}\frac{\lambda cr^{2}A\eta^{2}\xi}{\left(1+A\eta\xi\right)^{2}}=\frac{\lambda cr^{2}A\eta^{2}}{A^{2}\eta^{2}\xi+\frac{1}{\xi}+2A\eta}
(c)λcr2η41.\displaystyle\overset{(c)}{\leq}\frac{\lambda cr^{2}\eta}{4}\leq 1. (95)

(2) When λcr2η>4\lambda cr^{2}\eta>4 and A>AlA>A_{l}, the steady-state point psHp_{\mathrm{s}}^{H} exists, we have

|ps,n+1Hps,nH||ps,nHps,n1H|=(a)exp(λcr2η1+Aηξθrαρ)λcr2Aη2(1+Aηξ)2<(b)λcr2Aη2ξ(1+Aηξ)2(d)1,\begin{split}&\frac{|p_{\mathrm{s},n+1}^{H}-p_{\mathrm{s},n}^{H}|}{\big{|}p_{\mathrm{s},n}^{H}-p_{\mathrm{s},n-1}^{H}\big{|}}\\ &\overset{(a)}{=}\exp\left(-\frac{\lambda cr^{2}\eta}{1+A\eta\xi}-\frac{\theta r^{\alpha}}{\rho}\right)\frac{\lambda cr^{2}A\eta^{2}}{\left(1+A\eta\xi\right)^{2}}\\ &\overset{(b)}{<}\frac{\lambda cr^{2}A\eta^{2}\xi}{\left(1+A\eta\xi\right)^{2}}\overset{(d)}{\leq}1,\end{split} (96)

where (a) follows the mean theorem, (b) follows the ξ>exp(λcr2η1+Aηξθrαρ)\xi>\exp\left(-\frac{\lambda cr^{2}\eta}{1+A\eta\xi}-\frac{\theta r^{\alpha}}{\rho}\right) when ξ(ps,n+1H,ps,nH)\xi\in\left(p_{\mathrm{s},n+1}^{H},p_{\mathrm{s},n}^{H}\right), the last inequality (c) follows the fundamental inequality, and (d) follows A>AlA>A_{l} and ξ>ps,n+1H>psH>exp(112+141λcr2ηθrαρ)\xi>p_{\mathrm{s},n+1}^{H}>p_{\mathrm{s}}^{H}>\exp\left(-\frac{1}{\frac{1}{2}+\sqrt{\frac{1}{4}-\frac{1}{\lambda cr^{2}\eta}}}-\tfrac{\theta r^{\alpha}}{\rho}\right).

The convergence of psLp^{L}_{\mathrm{s}} can be shown via similar steps with initial iterative point ps,1L=exp(112141λcr2ηθrαρ)p_{\mathrm{s},1}^{L}=\exp\left(-\frac{1}{\frac{1}{2}-\sqrt{\frac{1}{4}-\frac{1}{\lambda cr^{2}\eta}}}-\tfrac{\theta r^{\alpha}}{\rho}\right), which is upper bound of the psLp^{L}_{\mathrm{s}}.

Then, the fixed-point iteration in (8) always converges with a linear convergence rate.

Appendix E AoI performance derivation: Theorem 2

Let us denote InI_{n} as the time interval between two consecutive transmission attempts over the typical link and JkJ_{k} the waiting time at the destination between successful receptions of the kk-th and (k+1k+1)-th updates, respectively. We have the following relationship between InI_{n} and JkJ_{k} according to the age threshold-based transmission protocol:

Jk=A+n=1NInJ_{k}=A+\sum_{n=1}^{N}I_{n} (97)

where NN is a random variable that represents the number of attempts between two successful transmissions. We further introduce a variable LkL_{k}, which represents the area under the AoI evolution curve across the kk-th successful update, as follows:

Lk=j=1Jkj=12Jk(Jk+1).L_{k}=\sum_{j=1}^{J_{k}}j=\frac{1}{2}J_{k}(J_{k}+1). (98)

Then, we can calculate the mean peak AoI and time average AoI over the typical link as:

Δ^0=𝔼[Jk],\hat{\Delta}_{0}=\mathbb{E}[J_{k}], (99)

and

Δ¯0=𝔼[Lk]𝔼[Jk]=12+𝔼[Jk2]2𝔼[Jk].\bar{\Delta}_{0}=\frac{\mathbb{E}[L_{k}]}{\mathbb{E}[J_{k}]}=\frac{1}{2}+\frac{\mathbb{E}[J_{k}^{2}]}{2\mathbb{E}[J_{k}]}. (100)

Since the transmission success probability psp_{\mathrm{s}} is i.i.d over the time, we can compute 𝔼[Jk]\mathbb{E}[J_{k}] and 𝔼[Jk2]\mathbb{E}[J_{k}^{2}] respectively as the following:

𝔼[Jk]=A+n=1𝔼[nps(1ps)n1]=A+𝔼[In]ps,\mathbb{E}[J_{k}]=A+\sum_{n=1}^{\infty}\mathbb{E}[np_{\mathrm{s}}(1-p_{\mathrm{s}})^{n-1}]=A+\frac{\mathbb{E}[I_{n}]}{p_{\mathrm{s}}}, (101)

and

𝔼[Jk2]=A2+2A𝔼[n=1NIn]+𝔼[(n=1NIn)2]=A2+2A𝔼[In]+𝔼[In2]ps+2(1ps)(𝔼[In])2ps2.\begin{split}\mathbb{E}[J_{k}^{2}]&=A^{2}+2A\mathbb{E}\left[\sum_{n=1}^{N}I_{n}\right]+\mathbb{E}\left[\left(\sum_{n=1}^{N}I_{n}\right)^{2}\right]\\ &=A^{2}+\frac{2A\mathbb{E}[I_{n}]+\mathbb{E}[I_{n}^{2}]}{p_{\mathrm{s}}}+\frac{2(1-p_{\mathrm{s}})\left(\mathbb{E}[I_{n}]\right)^{2}}{p_{\mathrm{s}}^{2}}.\end{split} (102)

when a source node is allowed to transmit, it generates new updates with probability η\eta independently over time, we have

𝔼[In]=1η,\displaystyle\mathbb{E}[I_{n}]=\frac{1}{\eta}, (103)
𝔼[In2]=2ηη2.\displaystyle\mathbb{E}[I_{n}^{2}]=\frac{2-\eta}{\eta^{2}}. (104)

The age performance is derived by combining (99), (100), (101), (102), (103) and (104).

Appendix F Proof of lemma 1

By substituting (6) into (13), and take a partial derivative of Δ^0\hat{\Delta}_{0} with respect to η\eta, we have

Δ^0η=1η2ps+λcr2ηps((1+Aηps)2λcr2Aη2ps).\frac{\partial{\hat{\Delta}_{0}}}{\partial{\eta}}=-\frac{1}{\eta^{2}p_{\mathrm{s}}}+\frac{\lambda cr^{2}}{\eta p_{\mathrm{s}}((1+A\eta p_{\mathrm{s}})^{2}-\lambda cr^{2}A\eta^{2}p_{\mathrm{s}})}. (105)

At the two extreme operating points of η\eta (i.e., η=0\eta=0 and η=1\eta=1), we have

limη0Δ^0η<0,\lim_{\eta\to 0}\frac{\partial{\hat{\Delta}_{0}}}{\partial{\eta}}<0, (106)

and

limη1Δ^0η=1ps+λcr2(1+Aps)2λcr2Aps.\lim_{\eta\to 1}\frac{\partial{\hat{\Delta}_{0}}}{\partial{\eta}}=-\frac{1}{p_{\mathrm{s}}}+\frac{\lambda cr^{2}}{(1+Ap_{\mathrm{s}})^{2}-\lambda cr^{2}Ap_{\mathrm{s}}}. (107)

when λcr2>1+Aps,\lambda cr^{2}>1+Ap_{\mathrm{s},{*}}, we have limη1Δ^0η>0\underset{{\eta\to 1}}{\lim}\frac{\partial{\hat{\Delta}_{0}}}{\partial{\eta}}>0, the mean peak AoI can then be optimized when η(0,1)\eta\in(0,1). By combining Δ^0η=0\frac{\partial{\hat{\Delta}_{0}}}{\partial{\eta}}=0 and (6), the optimal age threshold in mean peak AoI optimization can be obtained as

η^|A=1λcr2Aexp(1θrαρ1),\hat{\eta}^{*}_{|A}=\dfrac{1}{\lambda cr^{2}-A\exp\left(-1-\theta r^{\alpha}\rho^{-1}\right)}, (108)

The optimal mean peak AoI can be obtained by substituting (118) into (13).

When λcr21+Aps,\lambda cr^{2}\leq 1+Ap_{\mathrm{s},{*}}, on the other hand, the optimal age threshold is given by η=1\eta=1, and the corresponding optimal mean peak AoI can be obtained by combining η=1\eta=1 and (13).

Appendix G Proof of lemma 2

By substituting (6) into (13), and take a partial derivative of Δ^0\hat{\Delta}_{0} with respect to AA, we have

Δ^0A=1λcr2η(1+Aηps)2λcr2Aη2ps.\frac{\partial{\hat{\Delta}_{0}}}{\partial{A}}=1-\frac{\lambda cr^{2}\eta}{(1+A\eta p_{\mathrm{s}})^{2}-\lambda cr^{2}A\eta^{2}p_{\mathrm{s}}}. (109)

At the two extreme operating points of AA (i.e., A=0A=0 and AA\to\infty), we have

limA0Δ^0A=1λcr2η,\lim_{A\to 0}\frac{\partial{\hat{\Delta}_{0}}}{\partial{A}}=1-\lambda cr^{2}\eta, (110)

and

limAΔ^0A=1,\lim_{A\to\infty}\frac{\partial{\hat{\Delta}_{0}}}{\partial{A}}=1, (111)

when λcr2>1η\lambda cr^{2}>\frac{1}{\eta}, we have limA0Δ^0A<0\underset{{A\to 0}}{\lim}\frac{\partial{\hat{\Delta}_{0}}}{\partial{A}}<0, the mean peak AoI can then be optimized when A(0,)A\in(0,\infty). By combining Δ^0A=0\frac{\partial{\hat{\Delta}_{0}}}{\partial{A}}=0 and (6), the optimal age threshold in mean peak AoI optimization can be obtained as

A^|η=(λcr21η)exp(1+θrαρ1).\hat{A}^{*}_{|\eta}=\left(\lambda cr^{2}-\dfrac{1}{\eta}\right)\exp\left(1+\theta r^{\alpha}\rho^{-1}\right). (112)

The optimal mean peak AoI can be obtained by substituting (112) into (13).

On the other hand, when λcr21η\lambda cr^{2}\leq\frac{1}{\eta}, the optimal age threshold is given by A=0A=0, and the corresponding optimal mean peak AoI is obtained by combining A=0A=0 and (13).

Appendix H Proof of Theorem 3

Combining equation (6), Δ^0η=0\frac{\partial{\hat{\Delta}_{0}}}{\partial{\eta}}=0 and Δ^0A=0\frac{\partial{\hat{\Delta}_{0}}}{\partial{A}}=0, we have

ps=exp(1θrαρ1).p_{\mathrm{s}}=\exp\left(-1-\theta r^{\alpha}\rho^{-1}\right). (113)

Then, the optimal mean peak AoI can be achieved when (A^,η^)(\hat{A}^{*},\hat{\eta}^{*}) satisfy

A^=(λcr21η^)exp(1+θrαρ1).\hat{A}^{*}=\left(\lambda cr^{2}-\frac{1}{\hat{\eta}^{*}}\right)\exp\left(1+\theta r^{\alpha}\rho^{-1}\right). (114)

Due to the fact A0A\geq 0, the optimal parameter pair (A>0,η)(A>0,\eta) exists when the condition λcr2>min{1η}=1\lambda cr^{2}>\min\{\frac{1}{\eta}\}=1 is satisfied; Otherwise, the optimal parameter setting is (A=0,η=1)(A=0,\eta=1). The corresponding optimal mean peak AoI can be obtained by substitute (A,η)(A,\eta) of above two situation in (13).

Appendix I Proof of Lemma 3

By substituting (6) into (14) and taking a partial derivative of Δ¯0\bar{\Delta}_{0} with respect to η\eta, we have

Δ¯0η=(psλcr2ηps(1+Aηps)2λcr2Aη2ps)(A2(1+Aηps)21η2ps2).\frac{\partial{\bar{\Delta}_{0}}}{\partial{\eta}}=\left(p_{\mathrm{s}}-\frac{\lambda cr^{2}\eta p_{\mathrm{s}}}{\left(1+A\eta p_{\mathrm{s}}\right)^{2}-\lambda cr^{2}A\eta^{2}p_{\mathrm{s}}}\right)\left(\frac{A}{2(1+A\eta p_{\mathrm{s}})^{2}}-\frac{1}{\eta^{2}p_{\mathrm{s}}^{2}}\right). (115)

At the two extreme operating points of η\eta (i.e., η=0\eta=0 and η=1\eta=1), we have limη0Δ¯η\underset{{\eta\to 0}}{\lim}\frac{\partial{\bar{\Delta}}}{\partial{\eta}}\to-\infty, and

limη1Δ¯0η\displaystyle\lim_{\eta\to 1}\frac{\partial{\bar{\Delta}_{0}}}{\partial{\eta}} (116)
=(ps,λcr2ps,(1+Aps,)2λcr2Aps,)(A2(1+Aps,)21ps,2).\displaystyle=\left(p_{\mathrm{s},*}-\frac{\lambda cr^{2}p_{\mathrm{s},*}}{\left(1+Ap_{\mathrm{s},*}\right)^{2}-\lambda cr^{2}Ap_{\mathrm{s},*}}\right)\left(\frac{A}{2(1+Ap_{\mathrm{s},*})^{2}}-\frac{1}{p_{\mathrm{s},*}^{2}}\right).

We notice the last term of (116) satisfies

A2(1+Aps,)21(ps,)2=Aps,224Aps,2A2ps,22(1+Aps,)2ps,2<0.\frac{A}{2(1+Ap_{\mathrm{s},*})^{2}}-\frac{1}{(p_{\mathrm{s},*})^{2}}=\frac{Ap_{\mathrm{s},*}^{2}-2-4Ap_{\mathrm{s},*}-2A^{2}p_{\mathrm{s},*}^{2}}{2(1+Ap_{\mathrm{s},*})^{2}p_{\mathrm{s},*}^{2}}<0. (117)

Therefore, if the first term of (116) is negative (this is achieved when λcr2>1+Aps,\lambda cr^{2}>1+Ap_{s,*}), the time average AoI can be optimized by a specific η(0,1)\eta\in(0,1). As such, by solving the equation Δ¯0η=0\frac{\partial{\bar{\Delta}_{0}}}{\partial{\eta}}=0, the optimal update rate is obtained as:

η¯|A=1λcr2Aexp(1θrαρ1).\bar{\eta}^{*}_{|A}=\frac{1}{\lambda cr^{2}-A\exp\left(-1-\theta r^{\alpha}\rho^{-1}\right)}. (118)

Then, the optimal time-average AoI follows by substituting (118) into (14). On the other hand, if λcr21+Aps,\lambda cr^{2}\leq 1+Ap_{s,*}, the optimal update rate is achieved when η=1\eta=1, and the optimal time-average AoI is obtained by combining η=1\eta=1 and (14).

Appendix J Proof of Lemma 4

By substituting (6) into (14), and taking a partial derivative of Δ¯0\bar{\Delta}_{0} with respect to AA, we have

Δ¯0A=A2η2ps2+2Aηps+ηps2(1+Aηps)2λcr2η(2+4Aηps+A2η2ps2Aη2ps2)2(1+Aηps)2((1+Aηps)2λcr2Aη2ps).\begin{split}\frac{\partial{\bar{\Delta}_{0}}}{\partial{A}}&=\frac{A^{2}\eta^{2}p_{\mathrm{s}}^{2}+2A\eta p_{\mathrm{s}}+\eta p_{\mathrm{s}}}{2\left(1+A\eta p_{\mathrm{s}}\right)^{2}}\\ &-\frac{\lambda cr^{2}\eta\left(2+4A\eta p_{\mathrm{s}}+A^{2}\eta^{2}p_{\mathrm{s}}^{2}-A\eta^{2}p_{\mathrm{s}}^{2}\right)}{2\left(1+A\eta p_{\mathrm{s}}\right)^{2}\left(\left(1+A\eta p_{\mathrm{s}}\right)^{2}-\lambda cr^{2}A\eta^{2}p_{\mathrm{s}}\right)}.\end{split} (119)

The two extreme operating points of the age threshold lead to limA0Δ¯0A=12ηpsA=0λcr2η\underset{{A\to 0}}{\lim}\frac{\partial{\bar{\Delta}_{0}}}{\partial{A}}=\frac{1}{2}\eta p_{\mathrm{s}}^{A=0}-\lambda cr^{2}\eta and limAΔ¯0A=12\underset{A\to\infty}{\lim}\frac{\partial{\bar{\Delta}_{0}}}{\partial{A}}=\frac{1}{2}, respectively.

When λcr2η>𝕎0(η2exp(θrαρ1))\lambda cr^{2}\eta>\mathbb{W}_{0}\left(\frac{\eta}{2}\exp\left(-\theta r^{\alpha}\rho^{-1}\right)\right), we have limA0Δ¯0A<0\underset{{A\to 0}}{\lim}\frac{\partial{\bar{\Delta}_{0}}}{\partial{A}}<0. Since Δ¯0A\frac{\partial{\bar{\Delta}_{0}}}{\partial{A}} is continuous in AA, Intermediate Value Theorem assures that Δ¯0A=0\frac{\partial{\bar{\Delta}_{0}}}{\partial{A}}=0 has solutions in this case. As such, by combining Δ¯0A=0\frac{\partial{\bar{\Delta}_{0}}}{\partial{A}}=0 and (6), the optimal age threshold can be obtained from solving the following equation

(1+Aηps)3λcr2η(1+Aηps)2+(ηps1)(1+Aηps)λcr2η=0.(1+A\eta p_{\mathrm{s}})^{3}-\lambda cr^{2}\eta(1+A\eta p_{\mathrm{s}})^{2}+(\eta p_{\mathrm{s}}-1)(1+A\eta p_{\mathrm{s}})-\lambda cr^{2}\eta=0. (120)

(38) follows by combining (6) and (120).

The above analysis provides the conditions for the existence of roots. Next, we show that (38) (resp. (120)) has a unique non-zero root. We denote by u=1+Aηpsu=1+A\eta p_{\mathrm{s}} and introduce an auxiliary function as follows

f(u)=u3λcr2ηu2+(ηexp(θrαρ)exp(λcr2ηu)1)uλcr2η.\begin{split}&f(u)=u^{3}-\lambda cr^{2}\eta u^{2}\\ &\quad+\left(\eta\exp\left(\tfrac{\theta r^{\alpha}}{\rho}\right)\exp\left(-\tfrac{\lambda cr^{2}\eta}{u}\right)-1\right)u-\lambda cr^{2}\eta.\end{split} (121)

Since f(u)=0f(u)=0 is equivalent to (38), it can be used to analyze the number of roots.

(1.1) When 𝕎0(η2exp(θrαρ1))<λcr2η1\mathbb{W}_{0}\left(\frac{\eta}{2}\exp\left(-\theta r^{\alpha}\rho^{-1}\right)\right)<\lambda cr^{2}\eta\leq 1, the first derivative of f(u)f(u) can be expressed as follows:

f(u)u=3u22λcr2ηu1+ηexp(θrαρ)exp(λcr2ηu)(1+λcr2ηu).\begin{split}&\frac{\partial{f(u)}}{\partial{u}}=3u^{2}-2\lambda cr^{2}\eta u-1\\ &\quad~{}~{}~{}~{}+\eta\exp\left(-\tfrac{\theta r^{\alpha}}{\rho}\right)\exp\left(-\tfrac{\lambda cr^{2}\eta}{u}\right)\left(1+\tfrac{\lambda cr^{2}\eta}{u}\right).\end{split} (122)

Because 3u22λcr2ηu103u^{2}-2\lambda cr^{2}\eta u-1\geq 0 for u1u\geq 1 and the last term of f(u)u\frac{\partial{f(u)}}{\partial{u}} is positive for u[1,)u\in[1,\infty). We have f(u)u>0\frac{\partial{f(u)}}{\partial{u}}>0 for u[1,)u\in[1,\infty). Additionally, since f(u)f(u) is monotonically increasing for u[1,)u\in[1,\infty), f(u)=0f(u)=0 has at most one root in this region.

(1.2) When 1<λcr2η31<\lambda cr^{2}\eta\leq 3, the second derivative of f(u)f(u) satisfies

2f(u)u2=6u2λcr2η+ηexp(θrαρ1)exp(λcr2ηu)(λcr2η)2u3>0.\begin{split}&\frac{\partial^{2}{f(u)}}{\partial{u^{2}}}=6u-2\lambda cr^{2}\eta\\ &\quad~{}~{}+\eta\exp\left(-\theta r^{\alpha}\rho^{-1}\right)\exp\left(-\tfrac{\lambda cr^{2}\eta}{u}\right)\tfrac{(\lambda cr^{2}\eta)^{2}}{u^{3}}>0.\end{split} (123)

Therefore, the term f(u)u\frac{\partial f(u)}{\partial u} increases for u[1,)u\in[1,\infty). We denote by v=λcr2ηv=\lambda cr^{2}\eta and let g(v)=minu[1,)f(u)u=f(u)u|u=1g(v)=\underset{{u\in[1,\infty)}}{\min}\frac{\partial{f(u)}}{\partial{u}}=\frac{\partial{f(u)}}{\partial{u}}|_{u=1}, which can be expressed as

g(v)=22v+ηexp(θrαρ1)exp(v)(1+v).\begin{split}&g(v)=2-2v+\eta\exp\left(-\theta r^{\alpha}\rho^{-1}\right)\exp\left(-v\right)(1+v).\end{split} (124)

The first derivative of g(v)g(v) satisfies g(v)v<0\frac{\partial g(v)}{\partial v}<0, indicating that g(v)g(v) decreases monotonically for v=λcr2η(1,3)v=\lambda cr^{2}\eta\in(1,3). Additionally, notice that g(1)>0g(1)>0 and g(3)<0g(3)<0, g(v)=0g(v)=0 must have a root in (1,3)(1,3), which we denote by ψ\psi. Consequently, if if λcr2(1,ψ]\lambda cr^{2}\in(1,\psi], f(u)f(u) increases for u(1,)u\in(1,\infty); if λcr2(ψ,3]\lambda cr^{2}\in(\psi,3], f(u)f(u) decreases first, and then increases with uu. Moreover, since f(1)<0f(1)<0 and limuf(u)>0\underset{u\to\infty}{\lim}f(u)>0, f(u)f(u) has a unique root in this case.

(1.3) When λcr2η>3\lambda cr^{2}\eta>3, the third derivative of f(u)f(u) satisfy

3f(u)u3=6+ηexp(θrαρ1)×exp(λcr2ηu)(λcr2η)2(λcr2η3u)u5>6+(λcr2η)3u3exp(λcr2ηu)(1uu2)>6(λcr2η)3u3exp(λcr2ηu)>0.\begin{split}&\frac{\partial^{3}{f(u)}}{\partial{u^{3}}}=6+\eta\exp\left(-\theta r^{\alpha}\rho^{-1}\right)\times\\ &~{}~{}~{}~{}\exp\left(-\frac{\lambda cr^{2}\eta}{u}\right)\frac{(\lambda cr^{2}\eta)^{2}(\lambda cr^{2}\eta-3u)}{u^{5}}\\ &>6+\frac{(\lambda cr^{2}\eta)^{3}}{u^{3}}\exp\left(-\frac{\lambda cr^{2}\eta}{u}\right)\left(\frac{1-u}{u^{2}}\right)\\ &>6-\frac{(\lambda cr^{2}\eta)^{3}}{u^{3}}\exp\left(-\frac{\lambda cr^{2}\eta}{u}\right)>0.\end{split} (125)

Then, the second derivative 2f(u)u2\frac{\partial^{2}{f(u)}}{\partial{u^{2}}} increases for u[1,)u\in[1,\infty).

(a) If 2f(u)u2|u=1<0\frac{\partial^{2}{f(u)}}{\partial{u^{2}}}|_{u=1}<0, f(u)u\frac{\partial{f(u)}}{\partial{u}} decreases first, then increases for u[1,)u\in[1,\infty), and f(u)u|u=1<0\frac{\partial{f(u)}}{\partial{u}}|_{u=1}<0, limuf(u)u>0\underset{u\to\infty}{\lim}\frac{\partial{f(u)}}{\partial{u}}>0. Then, f(u)f(u) decreases first and then increases versus uu. Since f(1)<0f(1)<0 and limuf(u)>0\underset{u\to\infty}{\lim}f(u)>0, f(u)f(u) has a root at most.

(b) If 2f(u)u2|u=10\frac{\partial^{2}{f(u)}}{\partial{u^{2}}}|_{u=1}\geq 0, f(u)u\frac{\partial{f(u)}}{\partial{u}} monotonically increases for u[1,)u\in[1,\infty), and f(u)u|u=1<0\frac{\partial{f(u)}}{\partial{u}}|_{u=1}<0, limuf(u)u>0\underset{u\to\infty}{\lim}\frac{\partial{f(u)}}{\partial{u}}>0. Then, f(u)f(u) decreases first and then increases versus uu. Since f(1)<0f(1)<0 and limuf(u)>0\underset{u\to\infty}{\lim}f(u)>0, f(u)=0f(u)=0 has a single root.

The discussions in (1.1), (1.2), and (1.3) confirm that (120) has a unique root, and the optimal age threshold A1A_{1} is thus also unique. The results further indicate that when λcr2η>𝕎0(η2exp(θrαρ1))\lambda cr^{2}\eta>\mathbb{W}_{0}\left(\frac{\eta}{2}\exp\left(-\theta r^{\alpha}\rho^{-1}\right)\right), Δ¯0\bar{\Delta}_{0} decreases first, then increases with AA, the solution of (38) is the global minimum.

When λcr2η𝕎0(η2exp(θrαρ1))\lambda cr^{2}\eta\leq\mathbb{W}_{0}\left(\frac{\eta}{2}\exp\left(-\theta r^{\alpha}\rho^{-1}\right)\right), Δ¯0\bar{\Delta}_{0} monotonically increases for A[0,+)A\in[0,+\infty). The optimal age threshold is thus given by A=0A=0, and the corresponding optimal time-average AoI can be obtained by combining A=0A=0 and (14).

Appendix K Proof of Theorem 4

We denote by Δ¯0=min{η}Δ¯0A=A~|η\bar{\Delta}_{0}^{*}=\underset{\{\eta\}}{\min}~{}\bar{\Delta}_{0}^{A=\tilde{A}^{*}_{|\eta}}. On the one hand, for λcr21η𝕎0(η2exp(θrαρ))\lambda cr^{2}\leq\frac{1}{\eta}\mathbb{W}_{0}\left(\frac{\eta}{2}\exp\left(-\frac{\theta r^{\alpha}}{\rho}\right)\right), we have

Δ¯0A=A|η=1ηexp(λcr2η+θrαρ1).\bar{\Delta}_{0}^{A=A^{*}_{|\eta}}=\frac{1}{\eta}\exp\left(\lambda cr^{2}\eta+\theta r^{\alpha}\rho^{-1}\right). (126)

Then, we take a derivative of Δ¯0A=A|η\bar{\Delta}_{0}^{A=A^{*}_{|\eta}} with respect to η\eta, which gives

Δ¯0A=A|ηη=λcr2η1η2exp(λcr2η+θrαρ1).\tfrac{\partial{\bar{\Delta}_{0}^{A=A^{*}_{|\eta}}}}{\partial{\eta}}=\frac{\lambda cr^{2}\eta-1}{\eta^{2}}\exp\left(\lambda cr^{2}\eta+\theta r^{\alpha}\rho^{-1}\right). (127)

Following this result, we note that if

λcr2<1η𝕎0(η2exp(θrαρ)),\lambda cr^{2}<\frac{1}{\eta}\mathbb{W}_{0}\left(\frac{\eta}{2}\exp\left(-\frac{\theta r^{\alpha}}{\rho}\right)\right), (128)

Due to the fact,

1η𝕎0(η2exp(θrαρ))<1η𝕎0(12)<1η,\frac{1}{\eta}\mathbb{W}_{0}\left(\frac{\eta}{2}\exp\left(-\frac{\theta r^{\alpha}}{\rho}\right)\right)<\frac{1}{\eta}\mathbb{W}_{0}\left(\frac{1}{2}\right)<\frac{1}{\eta}, (129)

there is Δ¯0A=A|ηη<0\frac{\partial{\bar{\Delta}_{0}^{A=A^{*}_{|\eta}}}}{\partial{\eta}}<0, namely Δ¯0A=A|η\bar{\Delta}_{0}^{A=A^{*}_{|\eta}} decreases with the update rate η\eta in this region.

On the other hand, when λcr2>1η𝕎0(η2exp(θrαρ))\lambda cr^{2}>\frac{1}{\eta}\mathbb{W}_{0}\left(\frac{\eta}{2}\exp\left(-\frac{\theta r^{\alpha}}{\rho}\right)\right), Δ¯0A=A¯|η\bar{\Delta}_{0}^{A=\bar{A}^{*}_{|\eta}} can be explicitly approximated by (46), and we have

Δ¯0A=A~|ηη=3((Ω13)G(η)Ω13G(η))2(G(η))23Ω13exp(3λcr2ηΩ13G(η)+θrαρ)η2G(η),\tfrac{\partial{\bar{\Delta}_{0}^{A=\tilde{A}^{*}_{|\eta}}}}{\partial{\eta}}=\tfrac{-3\left((\Omega^{\frac{1}{3}})^{{}^{\prime}}G(\eta)-\Omega^{\frac{1}{3}}G^{{}^{\prime}}(\eta)\right)}{2(G(\eta))^{2}}\\ -\tfrac{3\Omega^{\frac{1}{3}}\exp\left(\frac{3\lambda cr^{2}\eta\Omega^{\frac{1}{3}}}{G(\eta)}+\frac{\theta r^{\alpha}}{\rho}\right)}{\eta^{2}G(\eta)}, (130)

where G(η)=Ω23+λcr2ηΩ13+(λcr2η)2+3G(\eta)=\Omega^{\frac{2}{3}}+\lambda cr^{2}\eta\Omega^{\frac{1}{3}}+(\lambda cr^{2}\eta)^{2}+3 and Ω\Omega is given in (45). Since the exponential term is greater than one, then

Δ¯0A=A~|ηη3η((Ω13)ηG(η)Ω13ηG(η))2η2(G(η))26ηΩ13G(η)2η2(G(η))2\displaystyle\dfrac{\partial{\bar{\Delta}_{0}^{A=\tilde{A}^{*}_{|\eta}}}}{\partial{\eta}}\leq-\frac{3\eta((\Omega^{\frac{1}{3}})^{{}^{\prime}}\eta G(\eta)-\Omega^{\frac{1}{3}}\eta G^{{}^{\prime}}(\eta))}{2\eta^{2}(G(\eta))^{2}}-\frac{6\eta\Omega^{\frac{1}{3}}G(\eta)}{2\eta^{2}(G(\eta))^{2}}
=(a)3(G(η)ηΩ13G(η)(Ω13+η(Ω13)))2η(G(η))23Ω132ηG(η)<(b)0,\displaystyle\quad~{}\overset{(a)}{=}\frac{3\left(G^{{}^{\prime}}(\eta)\eta\Omega^{\frac{1}{3}}-G(\eta)(\Omega^{\frac{1}{3}}+\eta(\Omega^{\frac{1}{3}})^{{}^{\prime}})\right)}{2\eta(G(\eta))^{2}}-\frac{3\Omega^{\frac{1}{3}}}{2\eta G(\eta)}\overset{(b)}{<}0, (131)

where (aa) follows algebraic manipulations. We further clarify the second inequality (bb) in the following.

Since G(η)G(\eta), Ω\Omega, and η\eta are positive, we only need to show that G(η)ηΩ13G(η)(Ω13+η(Ω13))<0G^{{}^{\prime}}(\eta)\eta\Omega^{\frac{1}{3}}-G(\eta)(\Omega^{\frac{1}{3}}+\eta(\Omega^{\frac{1}{3}})^{{}^{\prime}})<0. We expand this term as follows:

G(η)ηΩ13G(η)(Ω13+η(Ω13))=η2(λcr21)2Ω13+Ωη3Ω23(Ω23(λcr2η)23)Ω(3η2)Ω13.\begin{split}&G^{{}^{\prime}}(\eta)\eta\Omega^{\frac{1}{3}}-G(\eta)(\Omega^{\frac{1}{3}}+\eta(\Omega^{\frac{1}{3}})^{{}^{\prime}})=-\eta^{2}(\lambda cr^{2}-1)^{2}\Omega^{\frac{1}{3}}\\ &\quad\quad~{}~{}+\frac{\Omega^{{}^{\prime}}\eta}{3\Omega^{\frac{2}{3}}}\left(\Omega^{\frac{2}{3}}-(\lambda cr^{2}\eta)^{2}-3\right)-\Omega-(3-\eta^{2})\Omega^{\frac{1}{3}}.\end{split} (132)

In (132), the sign of terms in the second line requires further clarification. The Ωη\Omega^{{}^{\prime}}\eta in this term can be expressed as

Ωη=Ω+2(λcr2η)3+33(λcr2η)4+1(λcr2η)4+11(λcr2η)21,\Omega^{{}^{\prime}}\eta=\Omega+2(\lambda cr^{2}\eta)^{3}+3\sqrt{3}\frac{(\lambda cr^{2}\eta)^{4}+1}{\sqrt{(\lambda cr^{2}\eta)^{4}+11(\lambda cr^{2}\eta)^{2}-1}}, (133)

which satisfy Ωη<3Ω\Omega^{{}^{\prime}}\eta<3\Omega. Then, when Ω23(λcr2η)23>0\Omega^{\frac{2}{3}}-(\lambda cr^{2}\eta)^{2}-3>0, the terms in second line satisfy the following inequality

Ωη3Ω23(Ω23(λcr2)2η23)Ω(3η2)(Ω13)<Ω13(Ω23(λcr2)2η23)Ω(3η2)(Ω13)<0.\begin{split}&\quad~{}~{}\frac{\Omega^{{}^{\prime}}\eta}{3\Omega^{\frac{2}{3}}}\left(\Omega^{\frac{2}{3}}-(\lambda cr^{2})^{2}\eta^{2}-3\right)-\Omega-(3-\eta^{2})(\Omega^{\frac{1}{3}})\\ &~{}<\Omega^{\frac{1}{3}}\left(\Omega^{\frac{2}{3}}-(\lambda cr^{2})^{2}\eta^{2}-3\right)-\Omega-(3-\eta^{2})(\Omega^{\frac{1}{3}})<0.\end{split} (134)

Then, G(η)ηΩ13G(η)(Ω13+η(Ω13))<0G^{{}^{\prime}}(\eta)\eta\Omega^{\frac{1}{3}}-G(\eta)(\Omega^{\frac{1}{3}}+\eta(\Omega^{\frac{1}{3}})^{{}^{\prime}})<0. On the other hand, when Ω23(λcr2η)23<0\Omega^{\frac{2}{3}}-(\lambda cr^{2}\eta)^{2}-3<0, due to Ωη>0\Omega^{{}^{\prime}}\eta>0, the condition G(η)ηΩ13G(η)(Ω13+η(Ω13))<0G^{{}^{\prime}}(\eta)\eta\Omega^{\frac{1}{3}}-G(\eta)(\Omega^{\frac{1}{3}}+\eta(\Omega^{\frac{1}{3}})^{{}^{\prime}})<0 is still satisfied. Therefore, Δ¯0A~|ηη<0\frac{\partial{\bar{\Delta}_{0}^{\tilde{A}^{*}_{|\eta}}}}{\partial{\eta}}<0, which means that Δ¯0A=A|η\bar{\Delta}_{0}^{A=A^{*}_{|\eta}} decreases monotonically with update rate η\eta in this region.

Therefore, Δ¯0=min{η}Δ¯0A=A~|η\bar{\Delta}_{0}^{*}=\underset{\{\eta\}}{\min}~{}\bar{\Delta}_{0}^{A=\tilde{A}^{*}_{|\eta}} when η=1\eta=1, and Theorem 4 can be obtained by combining Lemma 4, η=1\eta=1 and (44).

References

  • [1] N. Pappas, M. A. Abd-Elmagid, B. Zhou, W. Saad, and H. S. Dhillon, Age of Information: Foundations and Applications. Cambridge University Press, 2023.
  • [2] S. Kaul, R. Yates, and M. Gruteser, “Real-time status: How often should one update?,” in Proc. IEEE INFOCOM, Orlando, FL, USA, pp. 2731–2735, Mar. 2012.
  • [3] R. D. Yates and S. K. Kaul, “The age of information: Real-time status updating by multiple sources,” IEEE Trans. Inf. Theory, vol. 65, no. 3, pp. 1807–1827, Sep. 2019.
  • [4] R. D. Yates and S. K. Kaul, “Status updates over unreliable multiaccess channels,” in Proc. IEEE Int. Symp. Inf. Theory (ISIT), Aachen, Germany, pp. 331–335, Jun. 2017.
  • [5] W. Zhan, D. Wu, X. Sun, Z. Guo, P. Liu, and J. Liu, “Peak age of information optimization of slotted Aloha: FCFS versus LCFS,” IEEE Trans. Netw. Sci. Eng., pp. 1–13, May 2023.
  • [6] B. Yu, Y. Cai, and D. Wu, “Joint access control and resource allocation for short-packet-based mmtc in status update systems,” IEEE J. Sel. Areas Commun., vol. 39, no. 3, pp. 851–865, Aug. 2021.
  • [7] M. Salimnejad and N. Pappas, “On the age of information in a two-user multiple access setup,” Entropy, vol. 24, no. 4, pp. 542–563, Apr. 2022.
  • [8] A. Munari, “Modern random access: An age of information perspective on irregular repetition slotted ALOHA,” IEEE Trans. Commun., vol. 69, no. 6, pp. 3572–3585, Feb. 2021.
  • [9] F. Vázquez-Gallego, J. Alonso-Zarate, and L. Alonso, “Reservation dynamic frame slotted-ALOHA for wireless M2M networks with energy harvesting,” in Proc. IEEE ICC, London, UK, pp. 5985–5991, Jun. 2015.
  • [10] H. Chen, Y. Gu, and S. C. Liew, “Age-of-information dependent random access for massive IoT networks,” in Proc. IEEE INFOCOM Workshop, Toronto, Canada, pp. 930–935, Jul. 2020.
  • [11] O. T. Yavascan and E. Uysal, “Analysis of slotted ALOHA with an age threshold,” IEEE J. Sel. Areas Commun., vol. 39, no. 5, pp. 1456–1470, May 2021.
  • [12] X. Chen, K. Gatsis, H. Hassani, and S. S. Bidokhti, “Age of information in random access channels,” IEEE Trans. Inf. Theory, vol. 68, no. 10, pp. 6548–6568, Jun. 2022.
  • [13] J. Sun, Z. Jiang, B. Krishnamachari, S. Zhou, and Z. Niu, “Closed-form Whittle’s index-enabled random access for timely status update,” IEEE Trans. Commun., vol. 68, no. 3, pp. 1538–1551, Mar. 2020.
  • [14] M. Ahmetoglu, O. T. Yavascan, and E. Uysal, “Mista: An age-optimized slotted ALOHA protocol,” IEEE Internet of Things J., vol. 9, no. 17, pp. 15484–15496, May 2022.
  • [15] M. Moradian, A. Dadlani, A. Khonsari, and H. Tabassum, “Age-aware dynamic frame slotted ALOHA for machine-type communications,” IEEE Trans. Commun., vol. 72, no. 5, pp. 2639–2654, May 2024.
  • [16] H. H. Yang, N. Pappas, T. Q. S. Quek, and M. Haenggi, “Analysis of the age of information in age-threshold slotted ALOHA,” in Int. Symposium Modeling and Optimization in Mobile, Ad Hoc, and Wireless Networks (WiOpt), Singapore, pp. 366–373, Aug. 2023.
  • [17] J. Massey and P. Mathys, “The collision channel without feedback,” IEEE Trans. Inf. Theory, vol. 31, no. 2, pp. 192–204, Mar. 1985.
  • [18] Y. Hu, Y. Zhong, and W. Zhang, “Age of information in Poisson networks,” in Proc. Int. Conf. Wireless Commun. Signal Process. (WCSP) Hangzhou, China, pp. 1–6, Oct. 2018.
  • [19] P. D. Mankar, M. A. Abd-Elmagid, and H. S. Dhillon, “Spatial distribution of the mean peak age of information in wireless networks,” IEEE Trans. Wireless Commun., vol. 20, no. 7, pp. 4465–4479, Feb. 2021.
  • [20] M. Emara, H. Elsawy, and G. Bauch, “A spatiotemporal model for peak AoI in uplink IoT networks: Time versus event-triggered traffic,” IEEE Internet of Things J., vol. 7, no. 8, pp. 6762–6777, Aug. 2020.
  • [21] H. H. Yang, C. Xu, X. Wang, D. Feng, and T. Q. S. Quek, “Understanding age of information in large-scale wireless networks,” IEEE Trans. Wireless Commun., vol. 20, no. 5, pp. 3196–3210, May 2021.
  • [22] H. H. Yang, A. Arafa, T. Q. S. Quek, and H. V. Poor, “Spatiotemporal analysis for age of information in random access networks under last-come first-serve with replacement protocol,” IEEE Trans. Wireless Commun., vol. 21, no. 4, pp. 2813–2829, Sep. 2022.
  • [23] P. D. Mankar, Z. Chen, M. A. Abd-Elmagid, N. Pappas, and H. S. Dhillon, “Throughput and age of information in a cellular-based IoT network,” IEEE Trans. Wireless Commun., vol. 20, no. 12, pp. 8248–8263, Jun. 2021.
  • [24] X. Sun, F. Zhao, H. H. Yang, W. Zhan, X. Wang, and T. Q. S. Quek, “Optimizing age of information in random-access Poisson networks,” IEEE Internet of Things J., vol. 9, no. 9, pp. 6816–6829, Sep. 2022.
  • [25] F. Zhao, X. Sun, W. Zhan, X. Wang, and X. Chen, “Information freshness in random-access Poisson network: Average AoI versus peak AoI,” in Proc. IEEE VTC2022-Fall, London, United Kingdom, pp. 1–6, Sep. 2022.
  • [26] F. Zhao, X. Sun, W. Zhan, X. Wang, J. Gong, and X. Chen, “Age-energy tradeoff in random-access Poisson networks,” IEEE Trans. Green Commun. Netw., vol. 6, no. 4, pp. 2055–2072, Dec. 2022.
  • [27] Z. Yue, H. H. Yang, M. Zhang, and N. Pappas, “Age of information under frame slotted ALOHA-based status updating protocol,” IEEE J. Sel. Areas Commun., vol. 41, no. 7, pp. 2071–2089, May 2023.
  • [28] H. H. Yang, A. Arafa, T. Q. S. Quek, and H. V. Poor, “Optimizing information freshness in wireless networks: A stochastic geometry approach,” IEEE Trans. Mobile Comput., vol. 20, no. 6, pp. 2269–2280, Jun. 2021.
  • [29] H. H. Yang, M. Song, C. Xu, X. Wang, and T. Q. S. Quek, “Locally adaptive status updating for optimizing age of information in Poisson networks,” IEEE Trans. Mobile Comput., vol. 22, no. 12, pp. 7343–7354, Dec. 2023.
  • [30] M. Song, H. H. Yang, H. Shan, J. Lee, and T. Q. S. Quek, “Age of information in wireless networks: Spatiotemporal analysis and locally adaptive power control,” IEEE Trans. Mobile Comput., vol. 22, no. 6, pp. 3123–3136, Dec. 2023.
  • [31] F. Baccelli and B. Blaszczyszyn, Stochastic Geometry and Wireless Networks: Volume I Theory. Boston, MA, USA: Now Foundations and Trends, 2009.
  • [32] M. Haenggi, Stochastic Geometry for Wireless Networks. Cambridge University Press, 2012.
  • [33] P. Gupta and P. Kumar, “The capacity of wireless networks,” IEEE Trans. Inf. Theory, vol. 46, no. 2, pp. 388–404, Mar. 2000.
  • [34] M. Haenggi, “On distances in uniformly random networks,” IEEE Trans. Inf. Theory, vol. 51, no. 10, pp. 3584–3586, Oct. 2005.
  • [35] S. S. Kalamkar and M. Haenggi, “The spatial outage capacity of wireless networks,” IEEE Trans. Wireless Commun., vol. 17, no. 6, pp. 3709–3722, Jun. 2018.
  • [36] M. Hyder and K. Mahata, “Zadoff–chu sequence design for random access initial uplink synchronization in lte-like systems,” IEEE Trans. Wireless Commun., vol. 16, no. 1, pp. 503–511, Jan. 2017.
  • [37] H. Chougrani, S. Kisseleff, W. A. Martins, and S. Chatzinotas, “Nb-iot random access for nonterrestrial networks: Preamble detection and uplink synchronization,” IEEE Internet of Things J., vol. 9, no. 16, pp. 14913–14927, Aug. 2022.
  • [38] Y. Sun, E. Uysal-Biyikoglu, R. D. Yates, C. E. Koksal, and N. B. Shroff, “Update or wait: How to keep your data fresh,” IEEE Trans. Inf. Theory, vol. 63, no. 11, pp. 7492–7508, Nov. 2017.
  • [39] K. Stamatiou and M. Haenggi, “Random-access Poisson networks: Stability and delay,” IEEE Commun. Lett., vol. 14, no. 11, pp. 1035–1037, Nov. 2010.
  • [40] Z. Gong and M. Haenggi, “The local delay in mobile Poisson networks,” IEEE Trans. Wireless Commun., vol. 12, no. 9, pp. 4766–4777, Sep. 2013.
  • [41] H. H. Yang, T. Q. S. Quek, and H. Vincent Poor, “A unified framework for sinr analysis in Poisson networks with traffic dynamics,” IEEE Trans. Commun., vol. 69, no. 1, pp. 326–339, Jan. 2021.
  • [42] M. Haenggi, “The meta distribution of the SIR in Poisson bipolar and cellular networks,” IEEE Trans. Wireless Commun., vol. 15, no. 4, pp. 2577–2589, Dec. 2016.
  • [43] M. Haenggi, “The local delay in Poisson networks,” IEEE Trans. Inf. Theory, vol. 59, no. 3, pp. 1788–1802, Nov. 2013.
  • [44] L. Dai, “Toward a coherent theory of CSMA and ALOHA,” IEEE Trans. Wireless Commun., vol. 12, no. 7, pp. 3428–3444, Jun. 2013.
  • [45] A. Carleial and M. Hellman, “Bistable behavior of ALOHA-type systems,” IEEE Trans. Commun., vol. 23, no. 4, pp. 401–410, Apr. 1975.
  • [46] W. Zhan and L. Dai, “Massive random access of machine-to-machine communications in lte networks: Modeling and throughput optimization,” IEEE Trans. Wireless Commun., vol. 17, no. 4, pp. 2771–2785, Feb. 2018.
  • [47] L. Dai, “A theoretical framework for random access: Stability regions and transmission control,” IEEE/ACM Trans. Networking, vol. 30, no. 5, pp. 2173–2200, Oct. 2022.
  • [48] R. M. Corless, G. H. Gonnet, D. E. G. Hare, D. J. Jeffrey and D. H. Knuth, “On the Lambert W function,” Adv. Comput. Math., vol. 5, pp. 329–359, Dec. 1996.