Isolation probabilities in dynamic soft random geometric graphs

Isolation probabilities in dynamic soft random geometric graphs


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 [1], 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 [2] and vehicular [3] networks and the Internet of Things [4]. In wireless networks the probability of a link decreases with the distance between nodes. As time evolves, the links form a dynamic network [5]. 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 [6] 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. (Equation 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 [7], 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 [8].

We distribute nodes and links according to the following spatial inhomogeneous random graph model [9]: 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. [10].

A Poisson Point Process with density \rho=2. Nodes are coloured by isolation probability, using the Rayleigh connection function, Eq. () with \eta=2 and periodic boundary conditions.
A Poisson Point Process with density . Nodes are coloured by isolation probability, using the Rayleigh connection function, Eq. () 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. ?. 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. [13] 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 [11]. 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. (Equation 1),


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

which is called the Roman harmonic number [14].

For not necessarily an integer, we can use Eq. (Equation 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.

The pdf of P_{iso}, by direct simulation in d=2 (points) and using Eqs. (,  ) with \alpha=100 (black solid lines). Top: \eta=2. Bottom: \rho=0.5. For \eta=2 and x>0.6, the first term of Eq. () is shown as a coloured line.
The pdf of , by direct simulation in (points) and using Eqs. (, ) with (black solid lines). Top: . Bottom: . For and , the first term of Eq. () is shown as a coloured line.
The pdf of P_{iso}, by direct simulation in d=2 (points) and using Eqs. (,  ) with \alpha=100 (black solid lines). Top: \eta=2. Bottom: \rho=0.5. For \eta=2 and x>0.6, the first term of Eq. () is shown as a coloured line.
The pdf of , by direct simulation in (points) and using Eqs. (, ) with (black solid lines). Top: . Bottom: . For and , the first term of Eq. () 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 [15], 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. ?. 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. (Equation 7) yields

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

Alternatively, we can take the limit . If increases proportional to , we see from Eq. (Equation 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. (Equation 6):

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

Probability that no nodes will remain isolated for time T. Here, d=\eta=2, L=20, 10^4 configurations: points. Theoretical curves, Eq. ().
Probability that no nodes will remain isolated for time . Here, , , configurations: points. Theoretical curves, Eq. ().

Eq. (Equation 10) has been confirmed by numerical simulation; see Fig. ?. 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. (Equation 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. [16]. 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.


  1. bibitemNoStop
  2. in (, ) pp. bibitemNoStop
  3. bibitemNoStop
  4. in (, ) pp. bibitemNoStop
  5. bibitemNoStop
  6. bibitemNoStop
  7. bibitemNoStop
  8. in ()bibitemNoStop
  9. bibitemNoStop
  10. bibitemNoStop
  11. bibitemNoStop
  12. bibitemNoStop
  13. (, )bibitemNoStop
  14. bibitemNoStop
  15. bibitemNoStop
  16. bibitemNoStop
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 minumum 40 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