Isolation probabilities in dynamic soft random geometric graphs

Isolation probabilities in dynamic soft random geometric graphs

Carl P. Dettmann School of Mathematics, University of Bristol, University Walk, Bristol BS81TW, UK    Orestis Georgiou School of Mathematics, University of Bristol, University Walk, Bristol BS81TW, UK Toshiba Telecommuncations Research Laboratory, 32 Queens Square, Bristol BS1 4ND, United Kingdom
July 3, 2019

We consider soft random geometric graphs, constructed by distributing points (nodes) randomly according to a Poisson Point Process, and forming links between pairs of nodes with a probability that depends on their mutual distance, the “connection function.” Each node has a probability of being isolated depending on the locations of the other nodes; we give analytic expressions for the distribution of isolation probabilities. Keeping the node locations fixed, the links break and reform over time, making a dynamic network; this is a good model of a wireless ad-hoc network with communication channels undergoing rapid fading. We use the above isolation probabilities to investigate the distribution of the time to transmit information to all the nodes, finding good agreement with numerics.

Many kinds of complex networks such as transport, power, social and neuronal networks are spatial in character Barthélemy (2011), that is, the nodes and perhaps also links have a physical location. Geometry structures the network in that the probability of a link between two nodes is related to their mutual distance.

Consider a wireless ad-hoc network where nodes (devices) communicate directly with each other rather than a central router and where their locations may be considered random; examples include sensor Goyal and Tripathy (2012) and vehicular Zeadally et al. (2012) networks and the Internet of Things Reina et al. (2013). In wireless networks the probability of a link decreases with the distance between nodes. As time evolves, the links form a dynamic network Zhang et al. (). The communication channel exhibits rapid fading, so that some time later, the state of the system may be chosen independently with the same distance-dependent probabilities. Here we assume that the nodes remain in fixed locations, at least on the rapid fading timescale. See supplemental material ani () for an animation in a square domain of length with nodes, showing connected components in different colours and pausing when the whole network is connected. The link probability between nodes of mutual distanace is Eq. (1) below with and .

If the link probability is either zero or one everywhere, there is randomness only due to the node locations. This is the case for the original random geometric graph (RGG) model Gilbert (1961), in which nodes connect if and only if their mutual distance is less than a threshold . If the link probability somewhere lies strictly between 0 and 1, there are two sources of randomness, in the node locations and the links. Here, we fix the node locations (“quenched disorder”), and study the randomness due to the links, as in the above dynamic wireless network application. This system has also been studied using an approach based on graph entropy Coon (2016).

We distribute nodes and links according to the following spatial inhomogeneous random graph model Penrose (): Place nodes in space according to a Poisson Point Process (PPP) with intensity measure in -dimensional space ; we usually consider . This means that the number of nodes in a bounded set is Poisson distributed with mean and independent of the number of points in any set disjoint with . Thus the average number of nodes in the whole system is , possibly infinite. The simplest case is where is proportional to Lebesgue measure, that is, where is the (constant) density and is the volume of . In this case, we often replace by a cube with opposite faces identified (a flat torus). Then .

Now we form links between each pair of nodes with locations , , independently with probability . Here we consider soft RGGs, for which where denotes the Euclidean (or in general some other) length and is called the connection function. It is possible with information about node locations and links to construct a connection function for any spatial network, and thus model it as a soft RGG. In practice, though, the link independence assumption may not be accurate. In the case of wireless communication networks, there are detailed theories of the physics of the communication channel leading to a variety of connection functions; see Refs. Durgin et al. (2002); Dettmann and Georgiou (2016); Penrose et al. (2016).

Figure 1: A Poisson Point Process with density . Nodes are coloured by isolation probability, using the Rayleigh connection function, Eq. (1) with and periodic boundary conditions.

Here, we use one of the simplest models: We assume Rayleigh fading, corresponding to diffuse scattering of the signal, which leads to an exponentially distributed channel gain . The signal power decays as where is called the path loss exponent. Free propagation gives the inverse square law , whilst more cluttered environments have a faster decay of the signal, leading to larger measured values of . A link may be made if the signal to noise ratio, proportional to , reaches a given threshold, leading to the connection probability


for some constant that determines the length scale; we measure length in these units and so take hereafter. Observe that the limit gives the original RGG model.

In order to understand transmission of information through a dynamic network, we must first analyse the instantaneous isolation probability of a node, that is, the probability that it has no links. This will determined by the locations of nearby nodes; see Fig. 1. Considering all nodes together, there is a distribution of isolation probabilities.

The isolation probability of a node at in some configuration of the PPP is

We note that in a PPP, the distribution of points found by conditioning on a node at a fixed position is unaffected (ie Palm distribution of a PPP is the same PPP); see Ref. Haenggi (2012) for the theory of PPPs. To find the distribution of , we use the probability generating functional (PGF)

for arbitrary function where the first equality is the definition and second follows for a PPP. The function needs to satisfy some mild conditions, for example (a) , or (b) and as the case here. In particular,

which is the connectivity mass, important for understanding the overall (multihop) connection probability of an ad-hoc network when  Dettmann and Georgiou (2016). However, we can also find further information about the distribution of , namely for the th moment is


If the PPP has constant density we find

independent of the location of the node. Here, is the total (solid) angle in dimensions, namely for respectively. Now, for Rayleigh fading, we have from Eq. (1),




For integer we can expand the parentheses to yield a finite sum


which is called the Roman harmonic number Roman (1992).

For not necessarily an integer, we can use Eq. (4): Numerical integration provides an efficient and stable means of calculation, whilst asymptotically expanding the integral for large gives

where is the Euler constant. Thus we have


where is the volume of the unit ball in dimensions.

When , that is, , we have a further simplification


where is the digamma function, and we have used its standard expansion for large argument. For integer

is the usual harmonic number.

Figure 2: The pdf of , by direct simulation in (points) and using Eqs. (35 8) with (black solid lines). Top: . Bottom: . For and , the first term of Eq. (9) is shown as a coloured line.

We can now attempt to extract the probability density function (pdf) of , which we will denote as with , from these moments. For a general distribution on a finite interval this is called the Hausdorff moment problem, and the solution is unique if it exists. For the most general numerical approach we follow Mnatsakanov Mnatsakanov (2008), who gives for a general pdf determined from integer moments , an approximation depending on a positive integer parameter :


The function depends on only through , and is hence piecewise constant for any fixed . It converges to the correct function as . It is possible to use this to get a numerical approximation to , using high precision arithmetic to overcome problems from cancellations; see Fig. 2. We see that is singular at or or both; when it is almost symmetrical at . It is never quite symmetrical: For we must have and then the third central moment is .

In the case of , we can write

giving a representation as an inverse Mellin transform

This is however intractable either analytically or numerically.

Still for , the large asymptotics does however give information on the behaviour of near , the distribution of highly isolated nodes. Making an ansatz

multiplying by and integrating gives

which by comparison with Eq. (7) yields


For moderate it is more accurate to keep just the first term, as shown in Fig. 2.

Alternatively, we can take the limit . If increases proportional to , we see from Eq. (3) that the only effect is to change the effective density . If is constant, we find , so that with

and so corresponds to a distribution that is sharply peaked at .

Now, we return to the problem of a dynamic network, assuming , constant density and neglecting boundary effects. The network chooses links anew each time . In order to ensure information can reach every node in the network, we need to ensure that no node is isolated for the considered time interval where is the number of time steps. At each time step, the isolation probabilities of the nodes are a PPP on with intensity . The probability that a node with is isolated for consecutive timesteps is simply . Denoting the event that none of the nodes are isolated during these timesteps by we can again use the PGF (noting that the number of nodes is almost surely finite):

Averaging over configurations of the PPP,

But the integral is just which we calculated in Eq. (6):


from which we find that the time to ensure all nodes are connected at least once is

Figure 3: Probability that no nodes will remain isolated for time . Here, , , configurations: points. Theoretical curves, Eq. (10).

Eq. (10) has been confirmed by numerical simulation; see Fig. 3. Thus, for low density the required time grows as a stretched exponential, controlled by the path loss exponent . When , it reduces to simply , with the probability distribution determined by the highly isolated nodes as in Eq. (9).

Strictly speaking all our results for isolation probabilities apply to one dimensional networks. However, in this case, transmission of information is limited by large gaps, rather than nodes that are likely to be isolated. In the original RGG, transmission can occur if and only if there are no gaps larger than the link range ; see Ref. Devroye (1981). For the soft RGG, it is quite likely that link may be made between nodes that are not directly adjacent to the gaps, and estimating the probability of connectivity, even at a single point in time, remains an open problem.

In conclusion, we have investigated the distribution of isolation probabilities in quenched soft random geometric graphs. This has allowed an analysis of the performance of a dynamic soft RGG model of wireless ad-hoc networks with fixed nodes and rapid channel fading. We obtained explicit formulas for the probability that no node will be isolated for time steps, with good numerical agreement. In contrast to networks with mobile nodes, the transmission of information is greatly hindered by extremes of the quenched disorder, namely highly isolated nodes. In the future it would be interesting to consider boundary effects, and other nonuniform node distributions, which are characteristic of many spatial complex networks.

The authors would like to thank the directors of Toshiba Telecommunications Research Laboratory and the EPSRC (grant EP/N002458/1) for their support. They are grateful to Justin Coon for helpful discussions.


  • Barthélemy (2011) M. Barthélemy, Physics Reports 499, 1 (2011).
  • Goyal and Tripathy (2012) D. Goyal and M. R. Tripathy, in 2012 Second International Conference on Advanced Computing & Communication Technologies (IEEE, 2012) pp. 474–480.
  • Zeadally et al. (2012) S. Zeadally, R. Hunt, Y.-S. Chen, A. Irwin,  and A. Hassan, Telecommunication Systems 50, 217 (2012).
  • Reina et al. (2013) D. G. Reina, S. L. Toral, F. Barrero, N. Bessis,  and E. Asimakopoulou, in Internet of Things and Inter-Cooperative Computational Technologies for Collective Intelligence (Springer, 2013) pp. 89–113.
  • (5) X. Zhang, C. Moore,  and M. Newman, “Random graph models for dynamic networks,” arXiv:1607.07570.
  • (6) Supplemental material, file animated_graph.gif.
  • Gilbert (1961) E. N. Gilbert, J. Soc. Indust. Appl. Math. 9, 533 (1961).
  • Coon (2016) J. P. Coon, in IEEE Globecom 2016 (accepted) (2016).
  • (9) M. D. Penrose, “Inhomogeneous random graphs, isolated vertices, and Poisson approximation,” arXiv:1507.07132.
  • Durgin et al. (2002) G. D. Durgin, T. S. Rappaport,  and D. A. De Wolf, IEEE Transactions on Communications 50, 1005 (2002).
  • Dettmann and Georgiou (2016) C. P. Dettmann and O. Georgiou, Physical Review E 93, 032313 (2016).
  • Penrose et al. (2016) M. D. Penrose et al., The Annals of Applied Probability 26, 986 (2016).
  • Haenggi (2012) M. Haenggi, Stochastic geometry for wireless networks (Cambridge University Press, 2012).
  • Roman (1992) S. Roman, The American mathematical monthly 99, 641 (1992).
  • Mnatsakanov (2008) R. M. Mnatsakanov, Statistics & Probability Letters 78, 1869 (2008).
  • Devroye (1981) L. Devroye, The Annals of Probability 9, 860 (1981).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description