On the Hitting Time Index of Graphs

Abstract

The hitting time index of a simple connected undirected graph G is a recently defined topological descriptor, HT( G ) , which is computed using the expected hitting times of the random walk on G . In this article, we find a new upper bound for HT( G ) , given in terms of the Kirchhoff index, and a new degree-weighted Kirchhoff index. Then we look at HT( G ) when G is edge transitive, both in the regular and non-regular cases. We find closed-form formulas for some families of Cayley graphs, the complete bipartite graphs K m,n , and some subdivision graphs.

Share and Cite:

Palacios, J. and Santamaría, A. (2025) On the Hitting Time Index of Graphs. Open Journal of Applied Sciences, 15, 2464-2478. doi: 10.4236/ojapps.2025.158164.

1. Introduction

Throughout this article, a graph G=( V,E ) shall be taken as a finite simple connected undirected graph with vertex set V={ 1,2,,n } , edge set E , and vertex degrees d 1 , d 2 ,, d n . The reader may consult reference [1] for all pertinent graph theoretical details.

Mathematical Chemistry makes use of these graphs in order to model molecules by identifying the vertices as the atoms and the edges as the atomic bonds between the vertices. A plethora of descriptors, or topological indices, i.e., real-valued functions defined on the domain of all graphs, have been used in trying to capture the physico-chemical properties of the molecules and to classify them according to the values of their indices.

Several of these descriptors are defined in terms of electrical or probabilistic concepts, and these two areas are intimately related as we will point out below. A notable such descriptor is the Kirchhoff index defined in [2] as

K( G )= i<j R ij , (1)

where R ij is the effective resistance across vertices i and j when the graph is considered as an electrical network, and all the edges are endowed with a resistance of 1 ohm. This descriptor was a natural evolution of an earlier one, the Wiener index, introduced in [3] as

W( G )= i<j d( i,j ), (2)

where d( i,j ) is the distance in the graph G between vertices i and j . It is clear that for a tree T one has W( T )=K( T ) and for a general G , W( G )K( G ) .

The simple random walk on G is the Markov chain { X n ,n0 } with state space V and uniform transition probabilities, that is, the walk jumps from vertex i to any of its d i neighboring vertices with probability 1 d i . The hitting time T j of the vertex j is the random variable that counts the number of jumps needed by the random walk to reach the vertex j for the first time:

T j =inf{ n0: X n =j },

and its expected value, when the process is started in state i is denoted by E i T j . Notice that the definition implies E i T i =0 , and this should not be confused with the mean return time to i , E i T i + = 2| E | d i , defined for the random variable T i + =inf{ n1: X n =i } . Reference [4] contains the basic facts about hitting times of Markov chains.

In [5], we showed that there is a close relationship between the Kirchhoff index of a graph G and the expected hitting times for the random walk on G , specifically

K( G )= 1 2| E | i<j ( E i T j + E j T i ), (3)

thus providing probabilistic ideas for the analysis of this index. Reference [6] gives a good introduction to the interplay between electric networks and random walks on graphs. One important result in this regard is the formula for the commute time between two vertices i and j , found in [7] and stating

E i T j + E j T i =2| E | R ij , (4)

which can be simplified, in case there is some symmetry that guarantees both hitting times are equal, to

E i T j =| E | R ij . (5)

We will be using these results below.

A recent probabilistic/electrical index was proposed by Camby et al. in [8], the random walk index, defined as follows: across any two vertices i and j , a battery is placed between them so that a 1 ampere current enters i and exits j . This setup implies a voltage v x ij on all vertices xV , and a potential drop on any edge ( x,y ) given by v x ij v y ij . If one inverts the polarity of the battery, then the potential difference on the edge ( x,y ) is v y ij v x ij , and therefore, to avoid the dependance on the polarity of the battery, the authors consider the absolute value | v x ij v y ij | , and the summation of these values over all edges of the graph, i.e.:

d ^ ij = ( x,y )E | v x ij v y ij |.

The function d ^ defined on the pairs of vertices ij by the value d ^ ij is shown to be a metric on the set of vertices, and the random walk index is defined as

RW( G )= i<j d ^ ij . (6)

Partly inspired by the ideas in [8], we defined in [9] a new probabilistic index, the hitting time index HT( G ) of a graph G , as

HT( G )= i<j D( i,j ), (7)

where D( i,j )=max{ E i T j , E j T i } .

We showed that D( i,j ) is actually a distance on the set of all vertices of G , found some inequalities involving HT( G ) , KG) , W( G ) , and RW( G ) , and computed HT( G ) for some families of graphs. These results were extended in [10], where we found tighter inequalities and closed-form expressions of HT( G ) for more families of graphs. The hitting time index is a good discriminator, in the sense that different graphs will get different values of HT( G ) , unlike the random walk index, which gives the same values to cycle graphs and linear graphs, for example, as mentioned in [9]. Moreover, recently there has been some interest in finding graphs that preserve the value of their Kirchhoff index after the removal of a vertex (see [11]). For example, for the 5-cycle, a removal of a vertex yields a 4-path, and the value of the Kirchhoff index of both graphs is 10. This may be interpreted as a discriminating weakness, to be contrasted with the HT index, which gives a value of 50 to the first graph and a value of 38 to the second one. This issue deserves further research.

A simple way to compute HT( G ) starts by finding the transition probability P of the random walk on G , with entries P( i,j )= 1 d i in case i and j are neighbors, and zero otherwise. We also use the matrix W , all of whose rows are identical to the stationary distribution π= 1 2| E | [ d 1 , d 2 ,, d n ] . Then we find the so-called fundamental matrix

Z= ( IP+W ) 1 ,

where I is the n×n identity matrix. Now, the matrix E of expected hitting times is found via the matrix Z . Its entries are

E( i,j )= E i T j = ( Z( j,j )Z( i,j ) )/ π j .

See [4] for a discussion of the matrices Z and E . Finally, the HT index of the graph under consideration is found by adding ( n1 2 ) terms,

i<j max{ E( i,j ),E( j,i ) }.

In this article, we continue with the study of HT( G ) . Specifically, we prove the inequality

HT( G )2| E |K( G ) K min ( G ),

where K min ( G ) is a new Kirchhoffian index defined by

K min ( G )= i<j R ij min{ d i , d j },

which resembles indices such as the multiplicative and the additive degree-Kirchhoff indices, introduced in [12] and [13], respectively.

We also find closed-form formulas for the values of the HT index in some families of graphs which are edge-transitive. For this concept of symmetry in graphs, we need first the definition of a graph automorphism, which is a bijection of the set of vertices of the graph that preserves adjacencies. Then we define a graph to be edge-transitive if for every pair of edges there is an automorphism that maps one edge onto the other. Thus, for example, the n-star graph is edge-transitive since every edge can be mapped into any other through a suitable rotation. The reader should be aware of the more restrictive definition of arc-transitivity where the automorphism maps any edge in any of its two orientations into any other edge. It is clear that arc-transitive graphs are edge-transitive, while the opposite is not true, as the n-star graph shows: in this graph, every edge can be mapped in one direction, but not in both onto any other. Examples of arc-transitive graphs are the n-cycle graph and the three-dimensional cube, where any edge, in any of its two orientations, can be mapped onto any other edge through rotations and symmetries.

2. A New Upper Bound

We showed in [9] that

HT( G )2| E |K( G ). (8)

This inequality was improved in [10] as follows:

HT( G )2| E |K( G )W( G ), (9)

where the equality is attained in case G= P 2 .

We also showed that for any tree T we have

HT( T )2| E |K( T ) W 2 ( T ), (10)

where W 2 ( G )= i<j d ( i,j ) 2 is the generalized Wiener index with parameter 2. The equality is attained for P 2 and P 3 .

We will prove a new inequality, which is better than (8) and non-comparable to (9), based on the following fact whose proof can be found in [14]:

Lemma 1. For a simple random walk on a graph G and any a,bG we have

E i T j R ij ( 2| E | d j ). (11)

Now we prove the new upper bound in the following.

Proposition 1. For all G we have

HT( G )2| E |K( G ) K min ( G ), (12)

where the equality is attained in case G= P 2 .

Proof. For all i,jG we have

D( i,j )=max{ E i T j , E j T i }max{ R ij ( 2| E | d j ), R ji ( 2| E | d i ) } = R ij ( 2| E |min{ d i , d j } ).

And thus

HT( G )= i<j D( i,j ) i<j R ij ( 2| E |min{ d i , d j } ) =2| E |K( G ) i<j R ij min{ d i , d j }=2| E |K( G ) K min ( G ).

Obviously, our new (12) is always better than (21). It is also better than (9) in case

W( G )< K min ( G ).

This happens, for example, in the case of the complete graph K n , where W( K n )=( n 2 )< K min ( K n )= ( n1 ) 2 , for n>2 . Also, for the n-cycle, for any two vertices i,j such that d( i,j )=k , the corresponding term R ij min{ d i , d j } in K min ( G ) has a value 2 k( nk ) n k because k n 2 , and the inequality is strict for k< n 2 . Also, for any tree T

W( T )=K( T ) K min ( T ),

and the inequality is strict if there is at least a pair i,jT such that neither i nor j are leaves, i.e., the inequality is strict for all trees other than the n-star graph, n2 .

On the other hand, it is easy to produce examples of graphs G where K min ( G )<W( G ) : take G n to be the n-star graph with an additional edge connecting any two leaves (in other words, G n is a triangle with n3 pendant vertices attached to the same vertex of the triangle). Then all summands in the computations of W( G n ) and K min ( G n ) are the same except for the vertices in the triangle, which contribute a total of 3 for W( G n ) and a total of 4 for K min ( G n ) , and the 2( n3 ) pairs of vertices i,j , where i is either of the vertices with degree 2 of the triangle and j is any pendant vertex. Each of these pairs contributes a distance 2 for W( G n ) , but only an effective resistance 2 3 +1 for K min ( G n ) . Therefore, 4+2( 2 3 +1 )( n3 )<3+4( n3 ) and thus K min ( G n )<W( G n ) .

This discussion shows then that the bounds (9) and (12) are noncomparable.

Another observation: (12) is better than (9) for all trees other than the n-star graph, n2 , but it seems to be worse than (10) for any tree other than the 2-star, though this conjecture needs a formal proof.

3. Edge-Transitive Graphs

As discussed in [15], an edge-transitive graph is either regular, or it is bipartite, where the partition consists of the set V 1 of all those vertices with the larger degree M , and the set V 2 of those vertices with the smaller degree m . The regular case is perhaps of lesser interest because for such case, the HT index is equal to the Kirchhoff index multiplied by the number of edges of the graph, and potentially we could find HT( G ) through the value K( G ) , which has been studied extensively. Still, we discuss this case first.

As was shown in [15], for a regular edge-transitive graph we have

Lemma 2.

E a T b =n1, (13)

whenever ( a,b )E . Also,

E a T b = E b T a , (14)

whenever d( a,b ) is even, where d( a,b ) is the distance in the graph between vertices a and b .

Immediate implications of (13) and (14) are collected in the following.

Corollary 1. For regular edge-transitive graphs we have

D( a,b )=n1, (15)

whenever ( a,b )E .

D( a,b )=| E | R ab , (16)

whenever d( a,b ) is even.

Corollary 1 gives a substantial amount of information on how to compute the HT index, but not all that is needed. For d( a,b )>2 , it is not clear how to get E a T b if d( a,b ) is odd, and even when d( a,b ) is even, it is not evident what the value of R ab may be. Let diam( G )=max{ d( a,b ):a,bV } . A first simple case where we can identify hitting times E a T b other than the case when d( a,b )=1 is given in the following

Proposition 2. If diam( G )=2 then E a T b =n , whenever d( a,b )=2 and therefore

HT( G )= n 2 ( n 2 nd ).

Proof. By conditioning on the first jump of the walk, it is clear that when d( a,b )=2 we have E a T b =1+ E c T b , where c and b are neighbors. But we know by the previous lemma that E c T b =n1 , so E a T b =n . Now, counting the pairs of vertices at distances 1 and 2, we have

HT( G )= nd 2 ( n1 )+[ ( n 2 ) nd 2 ]n= n 2 ( n 2 nd ) (17)

Examples of graphs included in the previous proposition are the complete bipartite graphs K n,n and the Cayley unitary graphs Cay( Z p k , U p k ) , for p prime (see [15] for details on the latter). A more involved example, not covered in proposition 1, is the unitary Cayley graph Cay( Z 2p , U 2p ) , for p prime, that has diameter 3, and their hitting times depend only on the distances 1, 2 or 3 between their vertices, with values

E T 1 =2p1,E T 2 = 2 ( p1 ) 2 p2 ,E T 3 = p( 2p3 ) p2 . (18)

Now we simply count pairs of vertices at distances 1, 2, 3. For distance 1 we have | V |d 2 such pairs, where d is the common degree, resulting in p 2 p pairs. For distance 3, we use theorem 3.3 in [15], and notice that the vertex 0 only has the vertex p at distance 3, and since this happens for all vertices, the total number of pairs at distance 3 is p . A simple subtraction yields that the number of pairs of vertices at distance 2 is p 2 p . Putting these facts and (18) together yields the following.

Proposition 3. For any prime p we have

HT( Cay( Z 2p , U 2p ) )=( p 2 p )( 2p1 )+( p 2 p ) 2 ( p1 ) 2 p2 + p 2 ( 2p3 ) ( p2 ) .

As a corollary, since for Cay( Z 2p , U 2p ) we have | E |=p( p1 ) , we can also give a closed-form formula for the Kirchhoff index of these graphs:

Corollary 2. For any prime p we have

K( Cay( Z 2p , U 2p ) )=( 2p1 )+ 2 ( p1 ) 2 p2 + p( 2p3 ) ( p2 )( p1 ) .

The Kirchhoff index of Cayley graphs has been studied in some depth, usually through the calculation of Laplacian eigenvalues (see [16] [17]) so the alternative approach presented here may be of some interest. Our calculations, based on probabilistic concepts and a symmetry condition, are generally simpler than those involving eigenvalues, and lead to closed form formulas which are also relatively simple. Now we will be interested in the non-regular bipartite case. It was shown in [15] that

Lemma 3. Whenever ( a,b )E , a V 1 and b V 2 we have

E a T b =2| V 2 |1, (19)

and

E b T a =2| V 1 |1. (20)

Equation (14) also applies in this case of non-regular edge-transitive graphs.

We start applying the previous result to the complete bipartite graph K m,n . Without loss of generality, let us assume that mn . Then, the partition of V is V 1 ={ vV:d( v )=n } and V 2 ={ vV:d( v )=m } . Clearly | V 1 |=m and | V 2 |=n . By Lemma 2, we have, for a V 1 ,b V 2 :

E a T b =2n1and E b T a =2m1.

Thus, for a V 1 ,b V 2 :

D( a,b )=2m1. (21)

The graph K m,n has diameter 2, so in order to find all the remaining different values of the expected hitting times, we only have to find E a T b when both a,b V 1 and when both a,b V 2 . Let us start with a,b V 1 . Because the graph is bipartite and has diameter 2, and conditioning on the first jump, it is clear that

E a T b =1+ E c T b ,

for c V 2 , i.e.,

E a T b =1+2m1=2m.

For the reverse expected hitting time we get the same result, E b T a =2m , so that for a,b V 1 we get

D( a,b )=2m. (22)

A similar argument yields, for a,b V 2 that

E a T b =1+2n1=2n,

and that in this case

D( a,b )=2n. (23)

Putting together (21), (22) and (23), we obtain the following.

Proposition 4. For m,n1 we have

HT( K m,n )= m 2 ( m1 )+ n 2 ( n1 )+nm( 2max{ m,n }1 ). (24)

Proof. We first assume mn . There are m( m1 ) 2 entries in the matrix E that we must include in the computation of HT( K m,n ) with value equal to 2m . Likewise, we must include n( n1 ) 2 entries with the value 2n , and finally we must include nm entries with the value 2m1 . Adding, we obtain the desired result, with max{ m,n }=m . If we assume nm we get the same result except that now max{ m,n }=n

Formula (24) is symmetric in n and m , as it should, because HT( K m,n )=HT( K n,m ) . For the particular case m=n we obtain HT( K n,n )= n 2 ( 4n3 ) . (this could also be found using (17) above.) On the other hand, the Kirchhoff index is K( K n.n )=4n3 , as can be seen in [18]. Therefore, in this case we have HT( K n,n )=| E |K( K n,n ) , which holds, as was pointed out in [9], because K n,n is walk-regular. Another interesting particular case is the star graph, S m = K m1,1 , for which we get

HT( K m1,1 )= (m1) 2 +( m1 )( 2m3 )=( m1 )( m 2 m1 ),

which coincides with the value for HT( S m ) obtained in [9]. So far, the graphs studied had diameter either 2 or 3. Now we will look at some graphs with diameter 4. One simple way to produce non-regular edge-transitive graphs is to introduce subdivision graphs: given a graph G=( V,E ) , its subdivision graph G =( V , E ) is built by splitting every edge of G with an extra vertex having degree 2. It is clear then that G is bipartite, with one set of the partition being the 2-degree vertices, and the other set the original vertices in V . Moreover, | V |=| V |+| E | and | E |=2| E | . We will denote by T i ' the hitting time of i on the subdivision graph G . Then we have:

Lemma 4. If G is arc-transitive then G is edge-transitive.

Proof. If ϕ is the automorphism of G that maps (either way) edges into edges, then the obvious extension of ϕ on V will also map any edge of G into any edge of G except maybe for the need of the application of a symmetry

Lemma 2 allows us to conclude that

Proposition 5. For iV,j V V we have

E i T j =2| E |1, (25)

and

E j T i =2| V |1, (26)

and therefore

D( i,j )=2| E |1. (27)

For the next proposition, we will need the following result, which can be found in [19]:

Lemma 5. Multiplying by the same factor r the resistance of every resistor in an electric network implies that the effective resistance between any two nodes of the network is multiplied by the same factor r .

Proposition 6. For i,jV , ( i,j )E , we have

E i T j = E j T i =D( i,j )=4( | V |1 ), (28)

More generally, for i,jV , ( i,j )E , we have

E i T j = E j T i =D( i,j )=4 E i T j . (29)

Proof. In G , all pairs of vertices in V are at an even distance, therefore (14) applies and

E i T j =| E | R ij =4| E | R ij =4 E i T j ,

proving (29). In particular, if d( i,j )=2 in V then d( i,j )=1 in V , and thus E i T j =| V |1 , proving (28).

Here we have used the fact that the introduction of the vertices splitting all edges amounts to doubling the value of all the resistances between neighboring vertices in V , implying, by the lemma, the doubling of any effective resistance between vertices in V , and thus R ij =2 R ij for i,jV

We will apply these formulas to a couple of examples. Let us first consider the complete graph K n =( V,E ) and its subdivision K n =( V , E ) . Then we have

Proposition 7. The expected hitting time E i T j in K n equals

(i) n 2 n1 if iV,j V V and d( i,j )=1 ;

(ii) 2n1 if the roles of i and j in (i) are reversed;

(iii) 4( n1 ) if i,jV ;

(iv) 4n3 if i V V,jV , and d( i,j )=3 ;

(v) n 2 +n3 if the roles of i and j in (iv) are reversed;

(vi) n 2 1 if i,j V V and d( i,j )=2 ;

(vii) n 2 +n2 if i,j V V and d( i,j )=4 .

Proof. (i), (ii) and (iii) are direct consequences of propositions 4 and 5. For (iv) we condition on the first jump of the walk and obtain

E i T j =1+ 1 2 E a T j + 1 2 E b T j ,

where a and b are the two neighbors of i . Now both a,bV and d( a,j )=d( b,j )=2 , so we can apply (iii) and get

E i T j =1+ E a T j =1+4( n1 ).

Now for (v), (vi) and (vii) we obtain a 3×3 linear system, again conditioning on the first jump of the walk. If i,j V V with d( i,j )=2 we have

E i T j =1+ 1 2 E a T j + 1 2 E b T j ,

with a,bV neighbors of i , d( a,j )=1 , d( b,j )=3. By (i), E a T j = n 2 n1 so we can write

E i T j =1+ 1 2 ( n 2 n1 )+ 1 2 E b T j . (30)

Also,

E b T j =1+ 2 n1 E i T j + n3 n1 E c T j , (31)

where c V V and d( c,j )=4 . Finally

E c T j =1+ E b T j . (32)

Solving for E i T j , E b T j and E c T j in the system determined by equations (30), (31) and (32) finishes the proof.

With the results in the previous proposition, we can prove the following:

Proposition 8. For n4 we have

HT( K n )=n( n1 )( n 2 n1 )+4( n1 )( n 2 )+n( n1 2 )( n 2 1 ) +[ ( ( n 2 ) 2 )n( n1 2 ) ]( n 2 +n2 ) +[ ( n+( n 2 ) 2 )( ( n 2 ) 2 )3( n 2 ) ]( n 2 +n3 ). (33)

= n( n1 )( n 4 +4 n 3 5 n 2 4 ) 8 . (34)

Proof. Every edge in the original graph K n generates two pairs of vertices at distance 1 in K n , so there are n( n1 ) such pars of vertices, where the larger hitting time between them is n 2 n1 , according to (i) and (ii) of the previous theorem. That takes care of the first summand of formula (33).

As far as pairs of vertices at distance 2 are concerned, there are ( n 2 ) such pairs generated by all possible pairs of vertices in the original K n . The hitting time for these is 4( n1 ) , given in (iii). In addition, we need to consider pairs of subdivision points of incident edges of the original graph K n which are at distance 2 in K n . The number of incident edges in K n is easily found by looking at a fixed vertex, which generates ( n1 2 ) such incident pairs of edges, and then multiplying by the n possible choices of the pivot vertex, for a total of n( n1 2 ) pairs of subdivision points at distance 2 in K n , whose hitting time between them is n 2 1 by (vi). This justifies the second and third summands in (33).

Next we look at vertices at distance 4 in K n . These are determined by all pairs of non-incident pairs of edges in K n , and their number is equal to the total number of pairs of edges minus the number of those which are incident, i.e., ( ( n 2 ) 2 )n( n1 2 ) . For these, the value of the hitting time between them, either way, is the same, given by (vii), as n 2 +n2 . This is the argument for the fourth summand of (33).

Finally, pairs of vertices in K n at distance 3 are found by subtracting from the grand total ( n+( n 2 ) 2 ) , all other pairs of vertices found in this proof. For these, the larger hitting times are n 2 +n3 as (iv) and (v) show. This justifies the last summand in (33). Adding then all the products of possible hitting times and pairs of vertices where these values occur, we get (33), and then some algebra simplifies this formula to the final expression (34)

The next example is the complete bipartite graph K n,n =( V,E ) and its subdivision graph K n,n =( V , E ) . Conditioning, it is easy to see that in K n,n we have E i T j =2n1 if i and j belong to different sets of the partition, and E i T j =2n if i and j belong to the same partition set. Then we have:

Proposition 9. The expected hitting time E i T j in K n,n equals

(i) 2 n 2 1 if iV,j V V and d( i,j )=1 ;

(ii) 4n1 if the roles of i and j in (i) are reversed;

(iii) 4( 2n1 ) if i,jV , and d( i,j )=2 ;

(iv) 8n if i,jV , and d( i,j )=4 ;

(v) 8n1 if i V V,jV and d( i,j )=3 ;

(vi) 2 n 2 +4n1 if the roles of i and j in (v) are reversed;

(vii) 2 n 2 +2n if i,j V V and d( i,j )=2 ;

(viii) 2 n 2 +4n if i,j V V and d( i,j )=4 .

Proof. (i) and (ii) are immediately derived from proposition 4; (iii) and (iv) follow from proposition 5. Now for (v) we condition on the first jump of the walk to obtain

E i T j =1+ 1 2 E a T j + 1 2 E b T j , (35)

with a,bV the neighbors of i and d( a,j )=4 , d( b,j )=2 . We then use (iii) and (iv) in (35) to conclude

E i T j =1+ 1 2 4( 2n1 )+ 1 2 8n=8n1.

For the three remaining cases, we condition on the first jump in order to obtain a 3×3 linear system. Indeed, if i,j V v and d( i,j )=2 we have

E i T j =1+ 1 2 E a T j + 1 2 E b T j ,

where a,bV are the neighbors of i and d( a,j )=1 , d( b,j )=3 . Using (i) we get then

E i T j =1+ 1 2 ( 2 n 2 1 )+ 1 2 E b T j . (36)

Also,

E b T j =1+ 1 n E i T j + n1 n E c T j , (37)

where c V V and d( c,j )=4 . Finally,

E c T j =1+ E b T j . (38)

Solving the system determined by (36), (37) and (38) finishes the proof.

With the results of the previous theorem we can prove the following

Proposition 10. For n2

HT( K n,n )=2 n 2 ( 2 n 2 1 )+4( 2n1 ) n 2 +2n( n 2 )( 2 n 2 +2n ) +[ ( n 2 2 )2n( n 2 ) ]( 2 n 2 +4n )+8 n 2 ( n1 ) +[ ( 2n+ n 2 2 )( n 2 2 )3 n 2 n( n1 ) ]( 2 n 2 +4n1 ) (39)

= n 6 +6 n 5 +5 n 4 +6 n 3 2 n 2 12 n 2 . (40)

Proof. Every edge in the original graph K n,n generates two pairs of vertices at distance 1 in K n,n , so there are 2 n 2 such pars of vertices, where the larger hitting time between them is 2 n 2 1 , according to (i) and (ii) of the previous theorem. That takes care of the first summand in (39).

As for pairs of vertices at distance 2, we see that there are n 2 such pairs generated by the endpoints of all edges in the original K n,n . The hitting time for these is 4( 2n1 ) , given in (iii). In addition, we need to consider pairs of subdivision points of incident edges of the original graph K n,n which are at distance 2 in K n . The number of incident edges in K n,n is easily found by looking at a fixed vertex, which generates ( n 2 ) such incident pairs of edges, and then multiplying by the 2n possible choices of the pivot vertex, for a total of 2n( n 2 ) pairs of subdivision points at distance 2 in K n,n , whose hitting time between them is 2 n 2 +2n by (vii). This justifies the second and third summands in (39).

Next we look at vertices at distance 4 in K n,n . Some are determined by all pairs of non-incident pairs of edges in K n,n , and their number is equal to the total number of pairs of edges minus the number of those which are incident, i.e., ( n 2 2 )2n( n 2 ) . For these, the value of the hitting time between them, either way, is the same, given by (viii), as 2 n 2 +4n . Also, pairs of points in the same part of the partition (of which there are 2( n 2 )=n( n1 ) ) are at distance 2 in K n,n and at distance 4 in K n,n , and the hitting time between them is 8n , according to (iv). This is the argument for the fourth and fifth summands in (39).

Finally, pairs of vertices in K n,n at distance 3 are found by subtracting from the grand total ( 2n+ n 2 2 ) , all other pairs of vertices found in this proof. For these,

the larger hitting times are n 2 +4n1 as (v) and (vi) show. This justifies the last of (39). Then some algebra yields the final expression (40).

Remark 1. For n=2 , our formula produces HT( K 2,2 )=336 , and we have that K 2,2 = C 8 , the 8-cycle. It is known that HT( C n )=| E |K( C 8 )= 1 12 n 2 ( n 2 1 ) , and for n=8 this yields HT( C 8 )=336 , in agreement with the value found for HT( K 2,2 ) . In [9] the value of HT( C n ) is wrongly reported as twice its correct value.

Remark 2. Formula (29) is shown in [20], for general graphs, using a spectral approach. That approach could have been used to find the expected hitting times on our subdivision graphs in terms of the hitting times on the original graphs, but we prefer the simple methods used here.

4. Conclusions

In this article, we have given a new upper bound for the hitting time index in terms of the Kirchhoff index and a new degree-weighted Kirchhoff index. This upper bound is better than or non-comparable to other bounds found in the literature. We have also worked with some families of graphs with diameters 2, 3 and 4, endowed with edge-transitivity, a form of symmetry weaker than edge-transitivity, that allows us to obtain, with probabilistic arguments, closed-form formulas for the hitting time index of these graphs. Our results deepen the understanding of the fairly recent HT index, showing a new relationship with other Kirchhoffian indices and enlarging the universe of graphs whose HT indices can be found with closed-form formulas.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding this article.

References

[1] Buckley, F. and Harary, F. (1990) Distance in Graphs. Addison-Wesley.
[2] Bonchev, D., Balaban, A.T., Liu, X. and Klein, D.J. (1994) Molecular Cyclicity and Centricity of Polycyclic Graphs. I. Cyclicity Based on Resistance Distances or Reciprocal Distances. International Journal of Quantum Chemistry, 50, 1-20.[CrossRef]
[3] Wiener, H. (1947) Structural Determination of Paraffin Boiling Points. Journal of the American Chemical Society, 69, 17-20.[CrossRef] [PubMed]
[4] Grinstead, M. and Snell, J.L. (1997) Introduction to Probability. 2nd Edition, American Mathematical Society.
[5] Palacios, J.L. (2000) Resistance Distance in Graphs and Random Walks. International Journal of Quantum Chemistry, 81, 29-33.[CrossRef]
[6] Doyle, P. and Snell, J. (1984) Random Walks and Electric Networks. American Mathematical Society.[CrossRef]
[7] Chandra, A.K., Raghavan, P., Ruzzo, W.L. and Smolensky, R. (1989) The Electrical Resistance of a Graph Captures Its Commute and Cover Times. Proceedings of the Twenty-First Annual ACM Symposium on Theory of ComputingSTOC’89, Seattle, 14-17 May 1989, 574-586.[CrossRef]
[8] Camby, E., Caporossi, G., Paiva, M.H.M. and Segatto, M.E.V. (2017) Expected Distance Based on Random Walks. Journal of Mathematical Chemistry, 56, 618-629.[CrossRef]
[9] Palacios, J.L. (2024) A New Probabilistic Molecular Index. Match Communications in Mathematical and in Computer Chemistry, 93, 499-507.[CrossRef]
[10] Palacios, J.L. and Santamarı́a, A. (2025) Some New Results for the Hitting Time Index. Match Communications in Mathematical and in Computer Chemistry, 94, 427-445.[CrossRef]
[11] Ergotić, S.M. and Gottwald, K.K. (2025) Graphs That Preserve Kirchhoff Index after the Removal of a Vertex. Boletín de la Sociedad Matemática Mexicana, 31, Article No. 88.[CrossRef]
[12] Chen, H. and Zhang, F. (2007) Resistance Distance and the Normalized Laplacian Spectrum. Discrete Applied Mathematics, 155, 654-661.[CrossRef]
[13] Gutman, I., Feng, L. and Yu, L. (2012) Degree Resistance Distance of Unicyclic Graphs. Transactions on Combinatorics, 1, 27-40.
[14] Palacios, J.L., Gómez, E. and Del Río, M. (2014) Hitting Times of Walks on Graphs through Voltages. Journal of Probability, 2014, Article ID: 852481.
[15] Palacios, J. and Renom, J. (1998) Random Walks on Edge Transitive Graphs. Statistics & Probability Letters, 37, 29-34.[CrossRef]
[16] Gao, X., Luo, Y. and Liu, W. (2011) Resistance Distances and the Kirchhoff Index in Cayley Graphs. Discrete Applied Mathematics, 159, 2050-2057.[CrossRef]
[17] Wang, Y., Zhang, Q. and Yuan, K. (2025) Kirchhoff Indices of Cayley Graphs on . Applied Mathematics and Computation, 489, Article ID: 129132.[CrossRef]
[18] Bapat, R.B., Karimi, M. and Liu, J. (2017) Kirchhoff Index and Degree Kirchhoff Index of Complete Multipartite Graphs. Discrete Applied Mathematics, 232, 41-49.[CrossRef]
[19] Palacios, J.L. (2003) A Probabilistic Proof of an Electric Network Property. Ciencia, 11, 252-253.
[20] Chen, H. (2017) Hitting Times for Random Walks on Subdivision and Triangulation Graphs. Linear and Multilinear Algebra, 66, 117-130.[CrossRef]

Copyright © 2026 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.