Hitting time for random walks on the Sierpinski network and the half Sierpinski network

We consider the unbiased random walk on the Sierpinski network (S_n◦N) and the half Sierpinski network (HS_n◦N), where n is the generation. Different from the existing works on the Sierpinski gasket, S_n◦N is generated by the nested method and HS_n◦N is half of S_n◦N based on the vertical cutting of the symmetry axis. We study the hitting time on S_n◦N and HS_n◦N. According to the complete symmetry and structural properties of S_n◦N, we derive the exact expressions of the hitting time on the nth generation of S_n◦N and HS_n◦N. The curves of the hitting time for the two networks are almost consistent when n is large enough. The result indicates that the diffusion efficiency of HS_n◦N has not changed greatly compared with S_n◦N at a large scale.

1 Introduction

In recent decades, complex networks have attracted great interest in the scientific community [1,2] and are considered valuable tools for describing real-world systems in the nature and society [3]. They also have a significant impact on mobility patterns, network sampling, community detection, signal transmission, virus spreading, epidemic control, and link prediction [4,5]. In recent years, complex networks have not only been studied in the field of mathematics due to their intersectionality and complexity but also more scholars in other disciplines have begun to pay attention to them, which involves system science, statistical physics, computer and information science, etc. [6–9] Common analysis methods and tools include the graph theory, combinatorial mathematics, matrix theory, probability theory, and stochastic process.

A fractal is a rough or fragmented geometric shape that can be divided into multiple parts, and each part is approximately a reduced-version copy of the whole. According to this definition, a fractal feature is a property known as self-similarity, which means that a fractal has self-similarity [10–13]. In recent years, fractals have attracted a surge of attention in various scientific fields [14,15]. The fractal theory [16,17] has always been a very popular and active theory. This is due to the self-similar structure that exists and the crucial impact of the idea of fractals on a large variety of scientific disciplines, such as molecular biology, pharmaceutical chemistry, optics, economics, and ecology [18]. In addition, many complex networks and real networks are generated by the self-similar fractal network [19–22]. There are many classical fractal models, such as the Cantor set, the Koch curve, and the Sierpinski gasket [23,24]. These structures have obviously become the focus of research, and many potential characterizations have been discovered. In 25–36, the network evolved from the Sierpinski carpet and some vital properties of the Sierpinski gasket were considered.

It has a great theoretical and practical significance to study the hitting time of random walks, which is the mean of the first-passage time of a random walker starting from any site on the network to a trap (a perfect absorber). It can characterize various other dynamical processes taking place on the network, such as mobility patterns and virus spreading. For example, in the communication and information industry, the hitting time of a random walk model can be used to study and simulate information transmission and latency, data collection, quantification and the prediction of communication and search costs, etc. In the field of biology, random walk models can be used to study and describe the spread of infectious diseases and metabolic fluxes among organisms. In the computer industry, the study of characteristics results in community detection, computer vision, collaborative recommendation, and image segmentation. It also describes the diffusion efficiency of different networks. Kozak and Balakrishnan [37,38] studied random walks on a two-dimensional Sierpinski gasket, three-dimensional Sierpinski tower, and d-dimensional Sierpinski model and gave the analytical expression of the hitting time. Wu [39] studied the random walk on the half Sierpinski gasket and gave the formula for the hitting time. Qi [40] got the expression of the hitting time for several absorbing random walks on Sierpinski graphs and hierarchical graphs.

In this paper, we study random walks and discuss the hitting time on the Sierpinski network model (S_n◦N) and the half Sierpinski network model (HS_n◦N) obtained by vertically cutting S_n◦N based on the symmetry axis. Compared with S_n◦N, the global self-similarity structure of HS_n◦N is lost, and only the local self-similarity is preserved. We expect to obtain the exact expression of the hitting time on HS_n◦N from the complete symmetry of S_n◦N and show whether there is any effect on the diffusion efficiency of the cut network. The mathematical combination method and fractal theory are applied. The remainder of this paper is organized as follows: in Section 2, we introduce our models. In Section 3, we have a detailed calculation of the hitting time of S_n◦N. Also, the detailed calculation of the hitting time of HS_n◦N is given in Section 4. In the last section, we draw the conclusion.

2 Preliminaries

In this section, we introduce the structure of S_n◦N and HS_n◦N and some concepts about several types of random walks, which will be used in the following.

2.1 The structure of S_n◦N and HS_n◦N

Actually, S_n◦N can be constructed in a nested manner with a self-similar structure. We separately define a triangle and a 3-regular graph as S and N. As shown in Figure 1, we can make three copies of the 3-regular graph N and then embed N into S by nodes to obtain the initial generation, which is recorded as S₁◦N. Making three copies of S₁◦N and then embedding the initial generation S₁◦N into S to produce the second generation is known as S(S₁◦N). For the convenience of the following description, the second generation is abbreviated as S₂◦N. We get the nth generation S(S_n−1◦N) by repeating the aforementioned process, making three copies of S_n−1◦N and then embedding the (n−1)th generation S_n−1◦N into S, which is abbreviated as S_n◦N. Specifically, S_n◦N is divided into three parts marked as $S_{n-1}^{(i)}◦N$ with i = 1, 2, and 3, according to the structure and iteration method. The upper half of Figure 2 shows the first two generations of S_n◦N. It should be noted that all nodes will be labeled sequentially from the top to the bottom by the site index i. The corner node 1 is the trap node, also called the target node. Otherwise, the corner node 1, the left-hand corner node, and the right-hand corner node of the bottom row on S_n◦N are represented by the set A. Therefore, it is convenient to refer to these three corner nodes as 1, L, and R, respectively. Also, the hitting time is represented by $T_1^{(n)}$, $T_L^{(n)}$, and $T_R^{(n)}$, where n is the generation.

FIGURE 1

FIGURE 1. Generation process of the initial generation S₁◦N.

FIGURE 2

FIGURE 2. First two generations of S_n◦N and HS_n◦N for n = 1 and 2.

As shown in the lower half of Figure 2, HS_n◦N is obtained by cutting the corresponding S_n◦N along the vertical symmetry axis. The cutting method of the network cannot equally divide all the nodes because vertices on the partition line will be retained during segmentation.

From the aforementioned construction, we can easily derive the total number of nodes on S_n◦N and HS_n◦N to be

$N_n=\frac{5\cdot 3^n+3}{2},(1)$

$\begin{array}{cc}\hfill H_n& =\frac{N_n-\left(n+2\right)}{2}+\left(n+2\right)\hfill \\ \hfill & =\frac{5\cdot 3^n+2n+7}{4}.\hfill \end{array}(2)$

The three connecting vertices C₁, C₂, and C₃, which are named after connecting the three regions $S_{n-1}^{(i)}◦N(i=\mathrm{1,2,3})$ on S_n◦N, are marked separately by the numbers as follows: $C_1=\frac{5\cdot 3^{n-1}+3-2^n}{2}$, $C_2=C_1+2^{n-1}=\frac{5\cdot 3^{n-1}+3}{2}$, and $C_3=C_2+5\cdot 3^{n-1}=\frac{5\cdot 3^n+3-2^n}{2}$. According to the position of the three connecting nodes relative to the trap node 1, we divide them into two categories, denoted by sets ${\mathrm{\Omega }}_n^1$ and ${\mathrm{\Omega }}_n^1$, i.e.,

${\mathrm{\Omega }}_n^1=\left\{C_1,C_2\right\},(3)$

${\mathrm{\Omega }}_n^2=\left\{C_3\right\}.(4)$

In addition, we use ${\mathrm{\Omega }}_n^f$ to denote the set of n non-trap nodes on the vertical secant line, except site 2.

2.2 Calculation of intermediate quantities

The hitting time is the mean of the first-passage time of a random walker starting from any site on the network to a trap node. In order to determine the hitting time for S_n◦N, we introduce the following intermediate quantities, all of which are the probability of a Markov chain on S_n◦N ending at corner node 1. In the process, the chain starts at a certain vertex and stops whenever it visits any of the three corner nodes of S_n◦N, where n is the generation.

p_n: The starting state is a corner node other than site 1. At least one transition is performed.

p₁(n): The starting state is a special vertex belonging to ${\mathrm{\Omega }}_n^1$.

p₂(n): The starting state is a special vertex belonging to ${\mathrm{\Omega }}_n^2$.

$\left\{\begin{array}{c}{p}_{n+1}=p_n{p}_2\left(n+1\right)+p_n{p}_1\left(n+1\right),\hfill \\ p_1\left(n+1\right)=\frac{1}{2}\left[p_n+\left(1-p_n\right)p_1\left(n+1\right)\right]+\frac{1}{2}\left[p_n{p}_2\left(n+1\right)+\left(1-2p_n\right)p_1\left(n+1\right)\right],\hfill \\ p_2\left(n+1\right)=\frac{1}{2}\left[p_n{p}_1\left(n+1\right)+\left(1-2p_n\right)p_2\left(n+1\right)\right]+\frac{1}{2}\left[p_n{p}_1\left(n+1\right)+\left(1-2p_n\right)p_2\left(n+1\right)\right].\hfill \end{array}\right.(5)$

The first equality of the equation group is explained. A random walk in S_n+1◦N, starting from the corner node L and ending whenever the walker reaches any three corner nodes of S_n+1◦N, is limited to jumping at least one step. By definition, the walker stops at the corner node R or 1 with probability p_n+1. In addition, the walker must first reach a connecting vertex C₁ or C₃ in order to reach one of the corner nodes R or 1. By definition and symmetry, the probability of reaching each of the connecting vertices C₁ or C₃ is p_n, where the connecting vertex C₁ belongs to ${\mathrm{\Omega }}_{n+1}^1$ and the other connecting vertex C₃ belongs to ${\mathrm{\Omega }}_{n+1}^2$, and the walker has the remaining probability 1−2p_n returning to the corner node L, where it stops walking. Therefore, the first equality is obtained.

We next certify the second equality in the equation system (Eq. 5). Considering a random walk in S_n+1◦N, which starts from the connecting vertex C₁ and ends whenever the walker reaches any three corner nodes of S_n+1◦N, it jumps at least one step. By definition, the probability that the walker stops at corner node 1 is p₁(n + 1). The probability of getting in $S_n^{(1)}◦N$ and $S_n^{(2)}◦N$ is $\frac{1}{2}$. If it enters $S_n^{(1)}◦N$, it has a probability of p_n to arrive at corner node 1, where the walker stops jumping. Also, the walker reaches the other two corner nodes of $S_n^{(1)}◦N$ with probability 1−p_n, i.e., connecting vertices C₁ and C₂ belonging to ${\mathrm{\Omega }}_{n+1}^1$. As a result, we can write the first item on the right side of the second equality. Similarly, if it enters $S_n^{(2)}◦N$, it reaches the corner node L or the connecting vertex C₃ with the same probability p_n. When it reaches the corner node L, it stops walking. When it reaches the connecting vertex C₃, it will continue to walk to the target node. Then, the walker has a probability 1−2p_n of returning to the connecting vertex C₁. Therefore, we can write the second item on the right side of the second equality. From these analyses, the aforementioned equality is obtained.

Ultimately, we testify to the third equality of the system of Eq. 5. A random walk in S_n+1◦N, starting from the connecting vertex C₃ and ending whenever the walker reaches any three corner nodes of S_n+1◦N, jumps at least one step. By definition, the probability that the walker stops at corner node 1 is p₂(n + 1). The probability of getting into $S_n^{(2)}◦N$ and $S_n^{(3)}◦N$ is $\frac{1}{2}$. If it enters $S_n^{(2)}◦N$, the probability of the walker reaching the corner node L is p_n, and the walker will stop walking at this point. In addition, the probability of the walker reaching the connecting vertex C₁ is p_n, and the probability of returning to the starting site C₃ is 1−2p_n. As a result, we can draw up the primary item on the right side of the last equality in the system of Eq. 5. Due to the complete symmetry of S_n+1◦N, the random walk on $S_n^{(3)}◦N$ is the same as that of $S_n^{(2)}◦N$. Therefore, the last equation is true.

It is easy to know that the initial value $p_1=\frac{4}{15}$. Through simplification, the final solution is obtained:

$\left\{\begin{array}{c}{p}_n=\frac{4}{15}{\left(\frac{3}{5}\right)}^{n-1},\hfill \\ p_1\left(n\right)=\frac{2}{5},\hfill \\ p_2\left(n\right)=\frac{1}{5}.\hfill \end{array}\right.(6)$

We next define the corresponding hitting time. Corner nodes 1, L, and R are represented by set A, where site 1 is set as the trap node; connecting vertices C₁, C₂, and C₃ are indicated by set I.

T_L→1(n): The hitting time from the corner node L to the corner node 1 in the nth generation.

T_L→A(n): The hitting time from the corner node L to any node in set A in the nth generation.

T_I→A(n): The hitting time from any vertex in set I to any node in set A in the nth generation.

$\left\{\begin{array}{c}{T}_{L\to 1}\left(n\right)=T_{L\to A}\left(n\right)+\left(1-p_n\right)T_{L\to 1}\left(n\right),\hfill \\ T_{L\to A}\left(n+1\right)=T_{L\to A}\left(n\right)+2p_n{T}_{I\to A}\left(n+1\right),\hfill \\ T_{I\to A}\left(n+1\right)=\frac{1}{2}\left[T_{L\to A}\left(n\right)+\left(1-p_n\right)T_{I\to A}\left(n+1\right)\right]+\frac{1}{2}\left[T_{L\to A}\left(n\right)+\left(1-p_n\right)T_{I\to A}\left(n+1\right)\right].\hfill \end{array}\right.(7)$

We now attest to the first equality in the equation system (Eq. 7). T_L→1(n) is the hitting time for a random walk in S_n◦N starting from the corner node L and ending at the trap node 1. According to the structure of S_n◦N, the walker must first reach any of the three corner nodes of S_n◦N in order to reach trap node 1, taking expected timesteps T_L→A(n). In such a process, the probability of reaching destination node 1 is p_n, where the walker stops jumping. Also, the probability of reaching corner nodes R and L is 1−p_n, from which the walker must continue to bounce T_L→1(n) steps to reach trap node 1. Thus, we get the first expression.

Subsequently, we prove the second equality in the equation system (Eq. 7). T_L→A(n + 1) is the hitting time for a random walk starting from the corner node L to any of the three corner nodes of S_n+1◦N for the first time, under the limitation that the walker jumps at least one step. In order to reach any of corner nodes in S_n+1◦N, the walker starting from the corner node L must first visit one of the corner nodes in S_n◦N. It is worth noting that these points belong to $S_n^{(2)}◦N$. This process is expected of T_L→A(n) timesteps. In this process, there is the same probability of reaching the connecting vertices C₁ and C₃, which is p_n. The probability of returning to the corner node L is 1−2p_n, stopping at this point. If the walker reaches any of the connecting vertices C₁ and C₃, it will continue to jump T_I→A(n + 1) steps to visit any of the corner nodes of S_n+1◦N.

Ultimately, we testify to the third equality of the system in Eq. 7. Consider a random walk, which starts from any one of the connecting nodes C₁, C₂, and C₃ and ends whenever it reaches any corner nodes on S_n+1◦N. By definition, the hitting time is T_I→A(n + 1). We now study the random walk from the connecting vertex C₃, which may perform the following two processes. Both processes happen with the probability of $\frac{1}{2}$. In the first process, the walker goes inside $S_n^{(2)}◦N$, taking T_L→A(n) timesteps to reach the three corner nodes of $S_n^{(2)}◦N$. The probability of reaching the corner node L is p_n, where the walking process is over. Also, the probability of reaching the other two corner nodes belonging to $S_n^{(2)}◦N$ is 1−p_n, i.e., the connecting vertices C₁ and C₃. The walker needs to take further T_I→A(n + 1) timesteps before being absorbed. According to the symmetry of S_n◦N, the second process is the same as the first one.

By substituting the value of $p_n=\frac{4}{15}{(\frac{3}{5})}^{n-1}$ into the aforementioned equation and simplifying it. In addition, it is easy to know T_L→A(1) = 4, so we can get the following:

$\left\{\begin{array}{cc}\hfill T_{L\to 1}\left(n\right)& =3\cdot 5^n,\hfill \\ \hfill T_{L\to A}\left(n\right)& =4\cdot 3^{n-1},\hfill \\ \hfill T_{I\to A}\left(n\right)& =3\cdot 5^{n-1}.\hfill \end{array}\right.(8)$

3 Formula of hitting time on S_n◦N

The Sierpinski network in any given generation n, using arrays (a, b, c) to represent a piece of interior points, for instance, (2, 5, 6) or (10, 16, 17), as shown in Figure 3; let (I₁, J₁, K₁) label the three vertices of the smallest triangle containing the three interior points (a, b, c), simultaneously, for instance, (1, 7, 9) or (7, 20, 22). According to the numerical results, the hitting time can be rewritten as follows:

$T_a\left(n\right)+T_b\left(n\right)+T_c\left(n\right)=T_{I_1}\left(n\right)+T_{J_1}\left(n\right)+T_{K_1}\left(n\right)+9.(9)$

FIGURE 3

FIGURE 3. Generation of the n = 3 Sierpinski network.

Let (i₁, j₁, k₁) label the three sites, which are one of any minimum size 1 lacunary triangle on S_n◦N. For example, (3, 4, 8) or (12, 13, 21); (I₁, J₁, K₁) also label the three vertices of the triangle containing (i₁, j₁, k₁) as its central lacunary region, such as (1, 7, 9) or (7, 20, 22), referring to Figure 3. According to the numerical results, it is easy to see the following:

$T_{i_1}\left(n\right)+T_{j_1}\left(n\right)+T_{k_1}\left(n\right)=T_{I_1}\left(n\right)+T_{J_1}\left(n\right)+T_{K_1}\left(n\right)+9\cdot 5^0.(10)$

In the same way, we now let (i₂, j₂, k₂) denote the three vertices of a lacunary region of size 2 in the network, such as (7, 9, 22) or (35, 37, 63); (I₂, J₂, K₂) label the three vertices of the triangle containing (i₂, j₂, k₂) as its central lacunary region, such as (1, 20, 24) or (20, 61, 65). It then follows the scaling derived previously that is as follows:

$T_{i_2}\left(n\right)+T_{j_2}\left(n\right)+T_{k_2}\left(n\right)=T_{I_2}\left(n\right)+T_{J_2}\left(n\right)+T_{K_2}\left(n\right)+9\cdot 5^1.(11)$

Therefore, moving up the hierarchy, if (i_r, j_r, k_r) are the sites demarcating a lacunary triangle of size r in the ascending order of size, starting from the smallest in size 1, and if (I_r, J_r, K_r) label the vertices of the triangle with (i_r, j_r, k_r) as the central lacunary region, then

$T_{i_r}\left(n\right)+T_{j_r}\left(n\right)+T_{k_r}\left(n\right)=T_{I_r}\left(n\right)+T_{J_r}\left(n\right)+T_{K_r}\left(n\right)+9\cdot 5^{r-1}.(12)$

The foregoing suggests how the hitting time T_total(n) may be computed for the arbitrary n. This is carried out by suitably regrouping the terms in the sum ${\sum }_{i=2}^{N_n}{T}_i(n)$ and systematically and repeatedly using Eq. 12 as one moves upward through triangles of increasing sizes. It needs some combinatorics and involves the enumeration of the number of lacunary triangles of each size in the network. The final result can be expressed entirely in terms of the known numerical factors and the combination $\left(T_1(n)+T_L(n)+T_R(n)\right)$. Since T₁(n) ≡ 0, while T_L(n) = T_R(n) = 3 ⋅ 5ⁿ, this leads directly to the desired expression for T_total(n). As an illustration of the procedure, consider the network of generation n = 3, with N_n = 69. Dropping for a moment the generation superscript for brevity and with L = 61 and R = 69, as shown in Figure 3, we obtain the following:

$T_{\mathit{total}}\left(3\right)=\sum _{i=2}^{69}{T}_i\left(3\right)=3\left(T_1+T_L+T_R\right)+3^2\cdot 9+3^2\cdot 9\cdot 5^0+5\left\{3^0\left(T_1+T_L+T_R\right)+3^1\cdot 9\cdot 5^1+3^1\left[\left(T_1+T_L+T_R\right)+3^0\cdot 9\cdot 5^2\right]\right\}.(13)$

Thus, T_total(3) has been recast in terms of the sum (T₁ + T_L + T_R) of the hitting time from the three primary sites. The meaning of 3² of the second term on the right side of the equation is that there are nine pieces of interior points in S₃◦N. It should be noted that the modulus 3² of the third term represents the number of regions of size 1, the coefficient 3¹ of the factor multiplying 9 ⋅ 5¹ is the number of lacunary triangles of size 2, and 3⁰ of the factor multiplying 9 ⋅ 5² is the number of lacunary triangles of size 3 on the n = 3 network.

We may now carry out a similar procedure for the case of general n. The analog of Eq. 13 yields the following:

$\begin{array}{ll}\hfill T_{\mathit{total}}\left(n\right)& =\sum _{i=2}^{N_n}{T}_i\left(n\right)\hfill \\ \hfill & =\left(3+5\cdot \sum _{m=0}^{n-2}{3}^m\right)\left(T_1\left(n\right)+T_L\left(n\right)+T_R\left(n\right)\right)+2\cdot 3^{n-1}\cdot 9+\left(3^{n-2}\cdot 9\cdot 5\right)\sum _{m=1}^{n-1}{5}^m.\hfill \end{array}(14)$

The numerical value is substituted to get the following result:

$T_{\mathit{total}}\left(n\right)=\frac{25}{4}\cdot 5^n\cdot 3^n+3\cdot 5^n-\frac{1}{4}\cdot 3^n.(15)$

Therefore, the hitting time on S_n◦N is as follows:

$\begin{array}{ll}\hfill \bar{T}\left(n\right)& =\frac{1}{N_n-1}{T}_{\mathit{total}}\left(n\right)\hfill \\ \hfill & =\frac{1}{N_n-1}\sum _{i=2}^{N_n}{T}_i\left(n\right)\hfill \\ \hfill & =\frac{25\cdot 5^n\cdot 3^n+12\cdot 5^n-3^n}{2\left(3^n\cdot 5+1\right)}.\hfill \end{array}(16)$

4 Formula of hitting time on HS_n◦N

We divide the random walk into two processes on S_n◦N and HS_n◦N, except node 2. The first process is that a walker starts from the starting point to node 3 or 4, which is the sum of the hitting time called T_g(n). The second procedure is a random walk from node 3 or 4 to the trap node, which is called T₃(n). In addition, T₂(n) denotes the hitting time from site 2 to trap node 1.

Similarly, the first process of HS_n◦N is a random walk from any site to node 3, and the hitting time of this process is recorded as H_g(n). The second process is recorded as H₃(n), and H₂(n) represents the hitting time from site 2 to trap node 1 on HS_n◦N. Thence, we have the following equations:

$\left\{\begin{array}{cc}\hfill T_{\mathit{total}}\left(n\right)& =T_g\left(n\right)+\left(N_n-2\right)T_3\left(n\right)+T_2\left(n\right),\hfill \\ \hfill H_{\mathit{total}}\left(n\right)& =H_g\left(n\right)+\left(H_n-2\right)H_3\left(n\right)+H_2\left(n\right).\hfill \end{array}\right.(17)$

Let T_r(n) be the mean of the first return time for a random walker starting from node 3 in the first process of HS_n◦N. Then, T_r(n) can also be the mean of the first return time to node 3 or 4 in the first process of S_n◦N. Therefore, according to the structure of HS_n◦N and S_n◦N, we have the following relationship:

$\left\{\begin{array}{cc}\hfill T_3\left(n\right)& =\frac{1}{6}+\frac{1}{6}\left[1+T_2\left(n\right)\right]+\frac{1}{6}\left[1+T_3\left(n\right)\right]+\frac{3}{6}\left[T_r\left(n\right)+T_3\left(n\right)\right],\hfill \\ \hfill H_3\left(n\right)& =\frac{1}{5}+\frac{1}{5}\left[1+H_2\left(n\right)\right]+\frac{3}{5}\left[T_r\left(n\right)+H_3\left(n\right)\right].\hfill \end{array}\right.(18)$

Moreover, from the numerical results, we have T₂(n) = 3 ⋅ 3ⁿ and H₂(n) = 2 ⋅ 3ⁿ. So, we get the following relation:

$T_3\left(n\right)-H_3\left(n\right)=\frac{1}{2}\left(3^n+1\right).$

Because the nodes on the cut line are retained, the hitting time T_i(n) from any node i on the secant line to the goal node 1 is also divided into two parts. The first process is a random walk from node i to node 3 or 4, which is denoted as T_i→3,4(n). The hitting time of the second process is denoted as T₃(n), which represents a random walk from site 3 to destination node 1. Therefore, T_i(n) can be rewritten as follows:

$T_i\left(n\right)=T_{i\to \mathrm{3,4}}\left(n\right)+T_3\left(n\right).(19)$

Consequently, the formula for the hitting time of nodes that belong to ${\mathrm{\Omega }}_n^f$ is as follows:

$\begin{array}{ll}\hfill T^f\left(n\right)& =\sum _{i\in {\mathrm{\Omega }}_n^f}\left(T_{i\to \mathrm{3,4}}\left(n\right)+T_3\left(n\right)\right)\hfill \\ \hfill & ={\bar{T}}^f\left(n\right)+n\cdot T_3\left(n\right).\hfill \end{array}(20)$

In S_n◦N, in addition to the nodes on the secant line, other nodes are evenly divided on both the sides in the segmentation process. HS_n◦N contains the nodes on the left half of the secant line and the nodes in set ${\mathrm{\Omega }}_n^f$. Therefore, it can be obtained by the following:

$H_g\left(n\right)=\frac{1}{2}\left[T_g\left(n\right)-{\bar{T}}^f\left(n\right)\right]+{\bar{T}}^f\left(n\right).(21)$

Eqs 17 and 20 are inserted into the aforementioned formula to get the following solution:

$H_g\left(n\right)=\frac{1}{2}\left[T_{\mathit{total}}\left(n\right)-\left(N_n+n-2\right)T_3\left(n\right)-T_2\left(n\right)+T^f\left(n\right)\right].(22)$

Then, substituting the aforementioned expression (Eq. 22) into the second formula in the equation group (Eq. 17), we get the following:

$\begin{array}{ll}\hfill H_{\mathit{total}}\left(n\right)& =H_g\left(n\right)+\left(H_n-2\right)H_3\left(n\right)+H_2\left(n\right)=\frac{1}{2}\left[T_{\mathit{total}}\left(n\right)-\left(N_n+n-2\right)T_3\left(n\right)-T_2\left(n\right)+T^f\left(n\right)\right]\hfill \\ \hfill & +\left(H_n-2\right)H_3\left(n\right)+H_2\left(n\right)=\frac{1}{2}\left[T_{\mathit{total}}\left(n\right)-\frac{3^n+1}{2}\left(N_n+n-2\right)+3^n+T^f\left(n\right)\right].\hfill \end{array}(23)$

We focus on solving the expression of T^f(n). In order to facilitate the calculation, we rewrite the three connecting vertices C₁, C₂, and C₃ of S_n◦N as d_n, e_n, and f_n, respectively. According to the structural characteristics of S_n◦N, the sites d_n and e_n in the nth generation can be labeled as the corner nodes L_n−1 and R_n−1 in the (n−1)th generation. The labeling method of nodes on S_n◦N can be seen in Figure 4.

FIGURE 4

FIGURE 4. Labeling method of nodes on S_n◦N.

If corner nodes 1, L_n, and R_n are set as goal nodes, then the walker starting from site f_y(y = 1, …, n) is captured by the trap set A after 3 ⋅ 5^y−1 random walks on average in S_n◦N, which can be seen from the analysis of the equation group (Eq. 7) and results. In addition, the arriving node L_y or R_y has a probability of $\frac{2}{5}$. Thus, $T_{f_y}(n)(y=1,\dots ,n)$ can be denoted as follows:

$\begin{array}{ll}\hfill T_{f_y}\left(n\right)& =3\cdot 5^{y-1}+\frac{2}{5}{T}_{L_y}\left(n\right)+\frac{2}{5}{T}_{R_y}\left(n\right)\hfill \\ \hfill & =3\cdot 5^{y-1}+\frac{4}{5}{T}_{L_y}\left(n\right)\hfill \\ \hfill & =3\cdot 5^{y-1}+\frac{4}{5}{T}_{d_{y+1}}\left(n\right).\hfill \end{array}$

Similarly, starting from the point d_y+1(y = 1,…, n) and getting to the site L_y or R_y has a probability of $\frac{2}{5}$ after 3 ⋅ 5^y random walks on average. Then, the walker is captured by the trap set A. In addition, there are $\frac{2}{5}$ and $\frac{1}{5}$ probabilities reaching nodes L_y+1 and R_y+1, respectively. Therefore, the aforementioned equation can be expressed as follows:

$\begin{array}{ll}\hfill T_{f_y}\left(n\right)& =3\cdot 5^{y-1}+\frac{4}{5}\left[3\cdot 5^y+\frac{2}{5}{T}_{L_{y+1}}\left(n\right)+\frac{1}{5}{T}_{R_{y+1}}\left(n\right)\right]\hfill \\ \hfill & =3\cdot 5^{y-1}+\frac{4}{5}\cdot 3\cdot 5^y+\frac{4}{5}\cdot \frac{3}{5}{T}_{L_{y+1}}\left(n\right)\hfill \\ \hfill & =\cdots \hfill \\ \hfill & =3\cdot 5^{y-1}+4\cdot 5^{y-1}\sum _{k=1}^{n-y}{3}^k+\frac{4}{5}\cdot {\left(\frac{3}{5}\right)}^{n-y}{T}_{L_n}\left(n\right)\hfill \\ \hfill & =3\cdot 5^{y-1}+4\cdot 5^{y-1}\sum _{k=1}^{n-y}{3}^k+\frac{4}{5}\cdot {\left(\frac{3}{5}\right)}^{n-y}\cdot 3\cdot 5^n\hfill \\ \hfill & =18\cdot 5^{y-1}\cdot 3^{n-y}-3\cdot 5^{y-1}.\hfill \end{array}(24)$

Here, $T_{L_n}(n)=3\cdot 5^n$ and the expression for the hitting time of nodes that belong to ${\mathrm{\Omega }}_n^f$ is as follows:

$\begin{array}{ll}\hfill T^f\left(n\right)& =\sum _{i\in {\mathrm{\Omega }}_n^f}{T}_i\left(n\right)\hfill \\ \hfill & =\sum _{y=1}^n{T}_{f_y}\left(n\right)\hfill \\ \hfill & =\sum _{y=1}^n\left(18\cdot 5^{y-1}\cdot 3^{n-y}-3\cdot 5^{y-1}\right)\hfill \\ \hfill & =\frac{33}{4}\cdot 5^n+9\cdot 3^n+\frac{3}{4}.\hfill \end{array}(25)$

Therefore, substituting Eq. 15 and Eq. 25 into Eq. 23, the expression of H_total(n) is as follows:

$H_{\mathit{total}}\left(n\right)=\frac{1}{2}\left[\frac{25}{4}\cdot 5^n\cdot 3^n+\frac{45}{4}\cdot 5^n-\frac{37}{4}\cdot 3^n-\frac{5}{4}\cdot 3^{2n}-\frac{n\left(3^n+1\right)}{2}+1\right].(26)$

Then, the hitting time on HS_n◦N is as follows:

$\bar{H}\left(n\right)=\frac{1}{H_n-1}{H}_{\mathit{total}}\left(n\right)=\frac{1}{H_n-1}\sum _{i=2}^{H_n}{H}_i\left(n\right)=\frac{2}{5\cdot 3^n+2n+3}\left[\frac{25}{4}\cdot 5^n\cdot 3^n+\frac{45}{4}\cdot 5^n-\frac{37}{4}\cdot 3^n-\frac{5}{4}\cdot 3^{2n}-\frac{n\left(3^n+1\right)}{2}+1\right].(27)$

In order to compare it with the exact formula for the hitting time of a random walk on S_n◦N and HS_n◦N, we draw the numerical simulation diagram of $\bar{T}(n)$ and $\bar{H}(n)$ for n = 1, 2, …, 7, as shown in Figure 5. The figure shows that the difference in the hitting time between S_n◦N and HS_n◦N is small when n is large enough. Therefore, the curves of both the networks are nearly merged for large scales, which indicate that the diffusion efficiency of HS_n◦N is consistent with S_n◦N for a large scale.

FIGURE 5

FIGURE 5. Numerical simulation diagram of $\bar{T}(n)$ and $\bar{H}(n)$.

5 Conclusion

In this paper, we study analytically the unbiased random walk on the Sierpinski network (S_n◦N) and the half Sierpinski network (HS_n◦N), which is obtained by vertically cutting S_n◦N along the symmetry axis. After cutting, the global self-similarity of S_n◦N is destroyed, but only the local self-similarity is maintained. We have analytically obtained the closed-form expression of the hitting time for a random walk on the nth generation HS_n◦N, which is $\bar{H}(n)=\frac{2}{5\cdot 3^n+2n+3}[\frac{25}{4}\cdot 5^n\cdot 3^n+\frac{45}{4}\cdot 5^n-\frac{37}{4}\cdot 3^n-\frac{5}{4}\cdot 3^{2n}-\frac{n(3^n+1)}{2}+1]$. However, the hitting time of S_n◦N is found to be $\bar{T}(n)=\frac{25\cdot 5^n\cdot 3^n+12\cdot 5^n-3^n}{2(3^n\cdot 5+1)}$. The hitting time is the quantity that characterizes the diffusion efficiency of the network. The curves of the two networks show that the hitting time of HS_n◦N is practically similar compared with S_n◦N when n is large enough, and our results show that the diffusion efficiency of HS_n◦N has little effect compared with S_n◦N at a large scale. Our work is further helpful in understanding the properties of Sierpinski networks.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Author contributions

All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The handling editor declared a shared affiliation with the authors at the time of review.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

1. Strogatz SH. Exploring complex networks. nature (2001) 410:268–76. doi:10.1038/35065725

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Wang XF, Chen G. Synchronization in scale-free dynamical networks: Robustness and fragility. IEEE Trans Circuits Syst (2002) 49:54–62. doi:10.1109/81.974874

CrossRef Full Text | Google Scholar

3. Guimera R, Amaral LAN. Functional cartography of complex metabolic networks. nature (2005) 433:895–900. doi:10.1038/nature03288

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Grubesic TH, Matisziw TC, Zook MA. Global airline networks and nodal regions. GeoJournal (2008) 71:53–66. doi:10.1007/s10708-008-9117-0

CrossRef Full Text | Google Scholar

5. Deng F-G, Li X-H, Li C-Y, Zhou P, Zhou H-Y. Quantum secure direct communication network with einstein–podolsky–rosen pairs. Phys Lett A (2006) 359:359–65. doi:10.1016/j.physleta.2006.06.054

CrossRef Full Text | Google Scholar

6. Batool A, Pál G, Danku Z, Kun F. Transition from localized to mean field behaviour of cascading failures in the fiber bundle model on complex networks. Chaos, Solitons & Fractals (2022) 159:112190. doi:10.1016/j.chaos.2022.112190

CrossRef Full Text | Google Scholar

7. Cotilla-Sanchez E, Hines PD, Barrows C, Blumsack S, Patel M. Multi-attribute partitioning of power networks based on electrical distance. IEEE Trans Power Syst (2013) 28:4979–87. doi:10.1109/tpwrs.2013.2263886

CrossRef Full Text | Google Scholar

8. Jung J, Jin W, Kang U. Random walk-based ranking in signed social networks: Model and algorithms. Knowl Inf Syst (2020) 62:571–610. doi:10.1007/s10115-019-01364-z

CrossRef Full Text | Google Scholar

9. Bai S, Niu M. Spectral properties of a class of treelike networks generated by motifs with single root node. Mod Phys Lett B (2022) 36:2250038. doi:10.1142/s0217984922500385

CrossRef Full Text | Google Scholar

10. Li B-G, Yu Z-G, Zhou Y. Fractal and multifractal properties of a family of fractal networks. J Stat Mech (2014) 2014:P02020. doi:10.1088/1742-5468/2014/02/p02020

CrossRef Full Text | Google Scholar

11. Zhou Y, Leung Y, Yu Z-G. Relationships of exponents in two-dimensional multifractal detrended fluctuation analysis. Phys Rev E (2013) 87:012921. doi:10.1103/physreve.87.012921

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Chen J, Dai M, Wen Z, Xi L. A class of scale-free networks with fractal structure based on subshift of finite type. Chaos (2014) 24:043133. doi:10.1063/1.4902416

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Wang Q, Xi L, Zhang K. Self-similar fractals: An algorithmic point of view. Sci China Math (2014) 57:755–66. doi:10.1007/s11425-013-4767-x

CrossRef Full Text | Google Scholar

14. Babič M, Mihelič J, Calì M. Complex network characterization using graph theory and fractal geometry: The case study of lung cancer dna sequences. Appl Sci (2020) 10:3037. doi:10.3390/app10093037

CrossRef Full Text | Google Scholar

15. Wang Z, Zhou Y, Bajenaid AS, Chen Y. Design of wireless sensor network using statistical fractal measurements. Fractals (2022) 30:2240092. doi:10.1142/s0218348x22400928

CrossRef Full Text | Google Scholar

16. Klafter J, Blumen A, Zumofen G. Fractal behavior in trapping and reaction: A random walk study. J Stat Phys (1984) 36:561–77. doi:10.1007/bf01012922

CrossRef Full Text | Google Scholar

17. Zhang P, Zhang Y, Huang Y, Xia Y. Experimental study of fracture evolution in enhanced geothermal systems based on fractal theory. Geothermics (2022) 102:102406. doi:10.1016/j.geothermics.2022.102406

CrossRef Full Text | Google Scholar

18. Blumen A, Klafter J, Zumofen G. Trapping and reaction rates on fractals: A random-walk study. Phys Rev B (1983) 28:6112–5. doi:10.1103/physrevb.28.6112

CrossRef Full Text | Google Scholar

19. Rammal R. Random walk statistics on fractal structures. J Stat Phys (1984) 36:547–60. doi:10.1007/bf01012921

CrossRef Full Text | Google Scholar

20. Hughes B, Montroll E, Shlesinger M. Fractal and lacunary stochastic processes. J Stat Phys (1983) 30:273–83. doi:10.1007/bf01012302

CrossRef Full Text | Google Scholar

21. Mandelbrot BB, Mandelbrot BB. The fractal geometry of nature, 1. New York: W. H. Freeman (1982).

Google Scholar

22. Chen J, Dai M, Wen Z, Xi L. Trapping on modular scale-free and small-world networks with multiple hubs. Physica A: Stat Mech its Appl (2014) 393:542–52. doi:10.1016/j.physa.2013.08.060

CrossRef Full Text | Google Scholar

23. Barlow MT, Bass RF. Transition densities for brownian motion on the sierpinski carpet. Probab Theor Relat Fields (1992) 91:307–30. doi:10.1007/bf01192060

CrossRef Full Text | Google Scholar

24. Zhou Z, Wu M. The hausdorff measure of a sierpinski carpet. Sci China Ser A-math (1999) 42:673–80. doi:10.1007/bf02878985

CrossRef Full Text | Google Scholar

25. Kajino N. Heat kernel asymptotics for the measurable riemannian structure on the sierpinski gasket. Potential Anal (2012) 36:67–115. doi:10.1007/s11118-011-9221-5

CrossRef Full Text | Google Scholar

26. Grabner PJ. Functional iterations and stopping times for brownian motion on the sierpiński gasket. Mathematika (1997) 44:374–400. doi:10.1112/s0025579300012699

CrossRef Full Text | Google Scholar

27. Wu S, Zhang Z, Chen G. Random walks on dual sierpinski gaskets. Eur Phys J B (2011) 82:91–6. doi:10.1140/epjb/e2011-20338-0

CrossRef Full Text | Google Scholar

28. Wu S, Zhang Z. Eigenvalue spectrum of transition matrix of dual sierpinski gaskets and its applications. J Phys A: Math Theor (2012) 45:345101. doi:10.1088/1751-8113/45/34/345101

CrossRef Full Text | Google Scholar

29. Guyer R. Diffusion on the sierpiński gaskets: A random walker on a fractally structured object. Phys Rev A (Coll Park) (1984) 29:2751–5. doi:10.1103/physreva.29.2751

CrossRef Full Text | Google Scholar

30. Hattori K. Displacement exponent for loop-erased random walk on the sierpiński gasket. Stochastic Process their Appl (2019) 129:4239–68. doi:10.1016/j.spa.2018.11.021

CrossRef Full Text | Google Scholar

31. Lara PCS, Portugal R, Boettcher S. Quantum walks on sierpinski gaskets. Int J Quan Inform (2013) 11:1350069. doi:10.1142/s021974991350069x

CrossRef Full Text | Google Scholar

32. Bedrosian SD, Sun X. Theory and application of pascal-sierpinski gasket fractals. Circuits Syst Signal Process (1990) 9:147–59. doi:10.1007/bf01236448

CrossRef Full Text | Google Scholar

33. Shima T. On eigenvalue problems for the random walks on the sierpinski pre-gaskets. Jpn J Ind Appl Math (1991) 8:127–41. doi:10.1007/bf03167188

CrossRef Full Text | Google Scholar

34. Chen J, Gao F, Le A, Xi L, Yin S. A small-world and scale-free network generated by sierpinski tetrahedron. Fractals (2016) 24:1650001. doi:10.1142/s0218348x16500018

CrossRef Full Text | Google Scholar

35. Chen J, Le A, Wang Q, Xi L. A small-world and scale-free network generated by sierpinski pentagon. Physica A: Stat Mech its Appl (2016) 449:126–35. doi:10.1016/j.physa.2015.12.089

CrossRef Full Text | Google Scholar

36. Xi L. Differentiable points of sierpinski-like sponges. Adv Math (2020) 361:106936. doi:10.1016/j.aim.2019.106936

CrossRef Full Text | Google Scholar

37. Kozak JJ, Balakrishnan V. Analytic expression for the mean time to absorption for a random walker on the sierpinski gasket. Phys Rev E (2002) 65:021105. doi:10.1103/physreve.65.021105

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Kozak JJ, Balakrishnan V. Exact formula for the mean length of a random walk on the sierpinski tower. Int J Bifurcation Chaos (2002) 12:2379–85. doi:10.1142/s0218127402006138

CrossRef Full Text | Google Scholar

39. Wu B, Zhang Z. The average trapping time on a half sierpinski gasket. Chaos, Solitons & Fractals (2020) 140:110261. doi:10.1016/j.chaos.2020.110261

CrossRef Full Text | Google Scholar

40. Qi Y, Dong Y, Zhang Z, Zhang Z. Hitting times for random walks on sierpiński graphs and hierarchical graphs. Comp J (2020) 63:1385–96. doi:10.1093/comjnl/bxz080

CrossRef Full Text | Google Scholar