Extreme Events and Event Size Fluctuations in Resetting Random Walks on Networks

1. Introduction

Extreme events (EE) are typically witnessed in many geophyscial phenomena such as the earthquakes, cyclones, drought and so on [1,2,3]. These represent extreme events recorded in univariate time series, and a large amount of theoretical and practical results have been obtained regarding the statistical and dynamic analysis of these extreme events [4,5]. However, as early as a century ago, extreme events in univariate time series arised in the research of classical extreme value theory, the extreme events in complex networks have only attracted attention in recent years [6,7,8]. Extreme events take place on complex networks is a fairly common place experience, such as traffic jam in roads and other transportation networks, web servers not responding due to the heavy load of http requests, floods in the network of rivers, power black outs due to tripping of power grids, etc [9,10]. Though such events typically have low probability of occurrence, they are regarded as significant owing to the social and financial losses suffered on account of most of these extreme events. Generally, these physical phenomena can be summed up as an emergent phenomena arising due to flux on networks and could be regarded as extreme events arising primarily because of the limited handling capacity of the node [11,12]. However, extreme event, defined as exceedences above the threshold, is not only affected by the capacity of the node, but also is attributable to natural fluctuations passing through a node. When the fluctuation in the flux passing through a node exceeds a prescribed threshold, extreme events will occur accordingly [6,13,14].

The analysis model we introduce in this research is the random walk on complex networks [15]. Random walks are ubiquitous in nature and have applications in a broad range of fields [16,17,18]. Recent research has shown that complex network exploration using random walks can be defined in terms of moving to nearest-neighbor nodes or as long-range hops between distant nodes [31,32]. Random walks are also interesting as a mechanism of transport and search within networks [19,30]. Recently, stochastic resetting has emerged as a significant concept in understanding a variety of random walk processes on complex networks [20,21,22,23,37,38]. This mechanism, which involves reinitializing a random walk to a particular node, has been explored to control and optimize diffusion-based search dynamical processes [24,33]. In addition to these developments, a broad theoretical foundation for stochastic resetting has been established in recent years. In particular, Pal and collaborators have made significant contributions to the understanding of resetting dynamics, first-passage under resetting, and optimal resetting strategies [34,35,36]. The effect of stochastic resetting on network-based processes, especially random walks, has implications for a range of applications, from search algorithms to biological transport phenomena [27,28]. Recent contributions relating to optimal reset protocols further emphasize the importance of tuning the reset conditions to achieve specific goals, such as minimizing the search time or optimizing resource allocation on the network [25,26]. By optimizing the reset protocol, the efficiency of random walks on complex networks can be expedited.

In this work, we model the transport process as resetting random walks on networks and delve into the study of the probability for EE and event size fluctuations. Relating to how the connectivity of the network affects the probability for the extreme event occurrence, by modeling the transport as the standard random walk, the V. Kishore has displayed in Ref [6] that the the probability for the occurrence of EE which depends only on the stationary distribution of a node arising due to inherent fluctuations. In this research, the threshold was choosen to be proportional to typical fluctuation size on i-th node and the extreme events are identified after taking care of the natural variability of the flux passing through the given node. Furthermore, it was shown that the probability for EE is higher for small degree nodes than for hubs or large degree nodes. This hints that the large flux does not translate into extreme events and even though hubs or large degree nodes draw large flux they are less liable to extreme events on a connected network. By leveraging the stochastic resetting, there is a surprising result that the stochastic resetting can reduce the probability of EE occuring and event size fluctuations in our research. And our work can be generalized to study the impact of stochastic resetting on complex networks. Specifically, we investigate the impact of stochastic resetting on the probability for EE and event size fluctuations within the Barabási-Albert network [29], which exhibits scale-free characteristics.

The remainder of this paper is organized as follows. Section 2 introduces random walks on networks with stochastic resetting. Section 3 describes the probability for extreme events with stochastic resetting. Section 4 elaborates fluctuations in event size with stochastic resetting. Section 5 summarizes the conclusions derived from this research.

2. Random Walk on Networks with Stochastic Resetting

A connected, undirected network with N nodes and E edges is considered. Its connectivity structure is determined by an adjacency matrix A whose element $A_{ij}=1$ if nodes i and j are connected by an edge, and is 0 otherwise.

On this network structure, diffusion of W independent walkers is considered who perform random walks in the sense explained below. We consider a random walker starts from node i at $t=0$ in discrete time.

The random walker performs at $t=1,2,3,\dots$ in two types of steps: a hop to one of the neighbors of the node presently occupied with probability $1-\gamma$ or a resetting to a given node r with probability $\gamma$.

Without resetting ($\gamma =0$), the probability of the random walker hops from i to j is $b_{i\to j}={\displaystyle \frac{A_{ij}}{k_i}}(i,j=1,2,\dots ,N)$ which is the element of the transition matrix B, depicting the random walk without resetting and possessing the characteristic of normalization (each row of B adds up to 1). Thus, the time evolution of probability is given by the master equation [20,21]:

$$P_{ij}(t+1;r,\gamma )=(1-\gamma )\sum _{l=1}^N{P}_{il}(t;r,\gamma )b_{l\to j}+\gamma {\delta }_{rj},$$

(1)

where $P_{ij}(t+1;r,\gamma )$ denotes the probability of the random walker starting from node i and reaching node j at time t given the resetting node r and the resetting rate $\gamma$. The first term of the equation represents the walker hoping to the nearest nodes from current position, and the second term indicates resetting to the given node r(Where ${\delta }_{rj}$ is the Kronecker delta). Let define $\Pi (r,\gamma )$ as the transition matrix with elements ${\pi }_{l\to j}(r,\gamma )=(1-\gamma )b_{l\to j}+\gamma {\delta }_{rj}$. Thus, the master equation can be rewritten as

$$P_{ij}(t+1;r,\gamma )=\sum _{l=1}^N{P}_{il}(t;r,\gamma ){\pi }_{l\to j}(r,\gamma ).$$

(2)

The form of Equation (2) satisfies the property of a Markov chain, which means ${\sum }_{j=1}^N{\pi }_{i\to j}(r,\gamma )=1.$ Like W, the matrix $\Pi (r,\gamma )$ is also a stochastic matrix: knowing its eigenvalues and eigenvectors expedites us to calculate the occupation probability at any time, even including the stationary distribution.

Owing to the fact that the matrix $\Pi (r,\gamma )$ is related to B, we can obtain the relationship between the two stochastic matrices. Let us denote the eigenvalues of the matrix B as ${\lambda }_i$, and its left and right eigenvectors are denoted as $\langle {\overline{\phi }}_i|$ and $|{\phi }_i\rangle$ respectively, for $i=1,2,\dots ,N$.

Similarly, the eigenvalues of $\Pi (r,\gamma )$ are denoted as ${\zeta }_i(r,\gamma )$ and its eigenvectors are denoted as $\langle {\overline{\psi }}_i(r,\gamma )|$ and $|{\psi }_i(r,\gamma )\rangle$. Based on the master equation, we import the identity

$$\Pi (r,\gamma )=(1-\gamma )B+\gamma \Psi \left(r\right),$$

(3)

where the elements of the matrix $\Psi \left(r\right)$ are ${\Psi }_{lm}\left(r\right)={\delta }_{mr}$, so the matrix has entries 1 in the $r^{\mathrm{th}}$-column and null entries everywhere else. Then the ${\zeta }_i(r,\gamma )$ can be described as [20]

$${\zeta }_i(r,\gamma )=\left\{\begin{array}{cc}1\hfill & \mathrm{for}i=1\hfill \\ (1-\gamma ){\lambda }_i\hfill & \mathrm{for}i=2,3,\dots ,N.\hfill \end{array}\right.$$

(4)

Furthermore, the left and right eigenvectors of $\Pi (r,\gamma )$ are described as

$$\langle {\overline{\psi }}_1(r,\gamma )|=\langle {\overline{\phi }}_1|+\sum _{i=2}^N{\displaystyle \frac{\gamma }{1-(1-\gamma ){\lambda }_i}}{\displaystyle \frac{\langle r|{\phi }_i\rangle }{\langle r|{\phi }_i\rangle }}\langle {\overline{\phi }}_i|,$$

(5)

where $\langle {\overline{\psi }}_i(r,\gamma )|=\langle {\overline{\phi }}_i|$ for $i=2,3,\dots ,N$. Similarly, the right eigenvectors are given by

$$|{\psi }_i(r,\gamma )\rangle =|{\phi }_i\rangle +{\displaystyle \frac{\gamma }{1-(1-\gamma ){\lambda }_i}}{\displaystyle \frac{\langle r|{\phi }_i\rangle }{\langle r|{\phi }_1\rangle }}|{\phi }_1\rangle ,$$

(6)

for $i=2,3,\dots ,N$, and $|{\psi }_1(r,\gamma )\rangle =|{\phi }_1\rangle$. $|r\rangle$ is defined as the vector whose elements equal to 0 except the r-th one, which equals to 1. With the left and right eigenvectors at hand, the matrix $\Pi (r,\gamma )$ can be characterized by the spectral representation

$$\Pi (r,\gamma )=\sum _{i=1}^N{\zeta }_i(r,\gamma )|{\psi }_i(r,\gamma )\rangle \langle {\overline{\psi }}_i(r,\gamma )|.$$

(7)

Within the discrete step, the master equation of the occupation probability can be rewritten as

$$P_{ij}(t;r,\gamma )={\langle i|\Pi (r,\gamma )}^t|j\rangle ,$$

(8)

Using the spectral representation, this can be restated as

$$P_{ij}(t;r,\gamma )=\langle i|{\psi }_1(r,\gamma )\rangle \langle {\overline{\psi }}_1(r,\gamma )|j\rangle +\sum _{l=2}^N{\left[(1-\gamma ){\lambda }_l\right]}^t\langle i|{\psi }_l(r,\gamma )\rangle \langle {\overline{\psi }}_l(r,\gamma )|j\rangle ,$$

(9)

where the stationary distribution, which we define as $p_j(r,\gamma )$, is equal to the first term of Equation (9). Utilizing the expression of the left and right eigenvectors of $\Pi (r,\gamma )$, the stationary distribution is given by

$$p_j(r,\gamma )={\displaystyle \frac{k_j}{{\sum }_{m=1}^N{k}_m}}+\gamma \sum _{l=2}^N{\displaystyle \frac{\langle r|{\phi }_l\rangle \langle {\overline{\phi }}_l|j\rangle }{1-(1-\gamma ){\lambda }_l}},$$

(10)

where $\langle i|{\phi }_1\rangle \langle {\overline{\phi }}_1|j\rangle ={\displaystyle \frac{k_j}{{\sum }_{m=1}^N{k}_m}}$ is the stationary distribution of the random walk without resetting. The second term of Equation (10) is the consequence of the resetting dynamical process, which disrupts the equilibrium of the general distribution.

3. Probability for Extreme Events with Stochastic Resetting

In this section, we focus the research emphasis on the probability for EE with stochastic resetting. First of all, we consider the distribution of the number of random walkers at a given node. We assume the dynamical process that there are a total of W random walkers wandering in the network with stochastic resetting. We define the probability of finding w random walkers at node i (the flux passing through a node) where the resetting node is r and the resetting rate is $\gamma$ as $f_i(w,r,\gamma )$. It is obvious that the walkers are relatively independent when performing the random walk process which leads to the variable w follows the binomial distribution. Thus, the probability $f_i(w,r,\gamma )$ is obtained by

$$f_i(w;r,\gamma )=\left({\displaystyle \genfrac{}{}{0pt}{}{W}{w}}\right)p_i{(r,\gamma )}^w{[1-p_i(r,\gamma )]}^{W-w}.$$

(11)

We can also derive the mean and standard deviation of finding w walkers at a given node is

$$\begin{array}{ccc}\hfill \langle f_i(r,\gamma )\rangle & =& W{p}_i(r,\gamma ),\hfill \\ \hfill {\sigma }_i(r,\gamma )& =& \sqrt{W{p}_i(r,\gamma )[1-p_i(r,\gamma )]}.\hfill \end{array}$$

(12)

We can observe that the relationship between the mean and standard deviation can be approximately expressed as ${\sigma }_i{(r,\gamma )}^2\approx \langle f_i(r,\gamma )\rangle$ owing to the value of $p_i(r,\gamma )$ is much smaller than 1. This relation can be thought of as a generalization of a similar relation for the standard random walks reported in Ref. [14]. Moreover, the expression of Equations (11) and (12) depend only on the stationary distribution $p_i(r,\gamma )$ which portrays a node, and are independent of the network topology, making them applicable to different types of networks, such as scale-free networks, random networks, or fractal networks. Figure 1 plots the binomial distribution of random walkers on the non-resetting node (k=4) and the resetting node (k=47) which is the central node (the node with the highest degree) in the network. Different colors represent different resetting rates. These theoretical analysis are based on the BA scale-free network with 200 nodes, 780 edges, and 1560 random walkers.

As is vividly depicted in Figure 1a,b, we can identify that the resetting strategy has distinct impacts on the central node and low-degree node. Figure 1a shows that resetting can reduce the probability of a large number of random walkers appearing at the low-degree node (k=4) to a certain extent, while it can significantly increase the probability of finding a large amount of random walkers at the central node on the contrary. This research result has laid a certain foundation for our better understanding of the probability for EE with stochastic resetting and the fluctuations in event size with stochastic resetting in the subsequent research.

We determine that the the probabilities for the occurrence of EE are very small, which are reflected at the tails of the binomial distribution. We apply this principle to each node in the BA scale-free network with the resetting node r and the resetting rate $\gamma$. As is mentioned above, at any given time, the number of random walkers w passing through a node in the network follows the binomial distribution given by Equation (11), an extreme event is considered to occur at a node if more than q random walkers pass through it, i.e., $w\ge q$. Therefore, the probability of an extreme event can be expressed as

$$F_i(r,\gamma )=\sum _{w=q_i}^W\left({\displaystyle \genfrac{}{}{0pt}{}{W}{w}}\right)p_i{(r,\gamma )}^w{(1-p_i(r,\gamma ))}^{W-w}.$$

(13)

First, let examine the scenario where the threshold $q_i$ is constant. If we set $q_i=0$, then according to Equation (13), the probability $F_i(r,\gamma )$ equals to 1 for all nodes, meaning that every node would experience extreme events all the time. Conversely, if we set $q_i=W$, the probability $F_i(r,\gamma )=p_i{(r,\gamma )}^W$, which approaches to zero since $p_i(r,\gamma )<<1$, implying that extreme events would virtually never occur. Both of these extreme cases are not physically meaningful. Similarly, setting a fixed threshold like $q_i=q_0$ (where $q_0$ is a constant) would result in some nodes experiencing extreme events almost all the time, while others rarely. Such a result is not what we intend to investigate. These considerations suggest that a more meaningful and realistic approach is to define the threshold based on the natural variability of the flux passing through a node. Therefore, we define the threshold for extreme events as [14]

$$q_i=\langle f_i(r,\gamma )\rangle +m{\sigma }_i(r,\gamma ).$$

(14)

We consider $m\ge 0$, and the standard deviation ${\sigma }_i(r,\gamma )$ and the mean $\langle f_i(r,\gamma )\rangle$ are given by Equation (12). Based on Equation (13), we can intuitively analyze that the probability of EE ($F_i(r,\gamma )$) is directly determined by the stationary distribution $p_i(r,\gamma )$. Figure 2 vividly presents the impact of different values of $p_i(r,\gamma )$ on the probability of EE occurring. Both Figure 2a,b show that the the nodes with large values of $p_i(r,\gamma )$ are less likely to experience extreme events and resetting can reduce $F_i(r,\gamma )$ occurring to a certain extent by comparing Figure 2a,b. Meanwhile, the extreme events with large values of parameter m are less likely to occur analyzed from the Figure 2. By substituting Equation (14) into Equation (13), we find that the probability of EE mainly depends on the resetting rate $\gamma$, the parameter m and the total number of random walkers W. However, once given $p_i(r,\gamma )$, the extreme events are a direct consequence of the shape of the binomial distribution independently of the walker evolution. Furthermore, in particular in the limit of large W, the central limit theorem holds and the binomial distribution of w converge to the Gaussian distribution that fully determines the value of the probability of EE [39], independently of $p_i(r,\gamma )$ which can be confirmed from Figure 1. Thus the probability of EE $F_i(r,\gamma )$ can be expressed as an integral as ${\int }_{q_i}^{\infty }{f}_i(w;r,\gamma )dw$. Figure 3 presents $F_i(r,\gamma )$ for different resetting nodes as the resetting rate $\gamma$. Figure 3a–c are based on simulations conducted on the BA scale-free network consisting of 200 nodes, 780 edges, and 1560 random walkers. The simulation results show excellent agreement with the analytical distribution obtained from Equation (13). The simulation results shown in Figure 3 are obtained by averaging over 50 independent runs.

By examining Figure 3a, we can intuitively obtain that the probability of EE on the low-degree resetting node in the BA network decreases significantly as the resetting rate $\gamma$ increases. When the resetting rate $\gamma$ reaches 1, the probability drops to zero, indicating that this resetting nodes no longer experience any extreme events. At the same time, as the parameter m increases, the threshold $q_i$ also increases, leading to a lower probability of EE for the node $(k=4)$. This is consistent with the analytical results derived from Equation (13). The same conclusion applies to Figure 3b,c as well. Moreover, even if a large amount of random walkers passing through the central node ($k=47$), the probability for the occurrence of extreme events is lower than that in nodes with a lower degree than it by comparing Figure 3a–c. Without loss of generality, we draw a conclusion that the probabilities of EE across all nodes in the BA scale-free network decrease as the resetting rate $\gamma$ increases. Furthermore, when $\gamma =1$, no extreme events occur in the network at all, as all random walkers are attracted by the resetting node.

To further validate the conclusions drawn from Figure 3 and to investigate the effect of stochastic resetting on extreme events, we compared the probability of EE under standard random walks and random walks with resetting. Figure 4 presents the probability of EE across all nodes in the BA scale-free network consisting of 100 nodes, 562 edges, and the results are performed by 1124 random walkers. The left panel corresponds to the case without resetting node, while the right panel shows the case where the central node serves as the resetting node with a resetting rate of $\gamma =0.1$. The color change of the energy bar represents the change of the probability for the occurrence of EE: warmer colors indicate a higher probability of extreme events occurring at a node, while cooler colors indicate a lower probability. It is clearly observed that in the presence of a resetting mechanism, the probability of extreme events is significantly reduced across the entire BA network. This observation is in strong agreement with the conclusion drawn from the Figure 3. Actually, we can clearly derive that the stochastic resetting obviously reduce the average occupation across all the nodes in the BA network by comparing the two panels in Figure 5 which corroborates with the results in Figure 4. Besides, we can also intuitively analyze that the probability of EE occurring in low degree nodes is higher than that in central nodes and nodes adjacent to the center. Further, notice that Equation (13) does not depend on the large scale structure of the topology, and therefore it is conducive to delve into other complex networks through theoretical and simulation analysis.

4. Fluctuations in Event Size with Stochastic Resetting

In this section, we focus on the event size fluctuations in the flux passing through a node with stochastic resetting. To explore how event size fluctuations depend on node connectivity in the presence of stochastic resetting, we analyze the distribution of event sizes with nodes arranged in descending order of their degrees and nodes with the same degree in ascending order of their distance to the resetting node. The visualization of event-size distributions in this order highlights how highly connected nodes (hubs) respond to fluctuations in walker flux compared to peripheral nodes.

Figure 6 shows the distribution of event size m, measured in units of the standard deviation $\sigma$, as a function of node degree and their distance to the resetting node. This theoretical analysis is based on the BA scale-free network used in Figure 4. The color map displays the probability of an event of size m occurring on nodes with different degrees or same degrees with different distances to the resetting node. In our research, an event of size m is defined by the condition $m\sigma \le w-\langle f_i(r,\gamma )\rangle \le (m+1)\sigma$, where w is the number of walkers discovered on the node. Then, the probability for the occurrence of an event of size m can be described as

$$\begin{array}{ccc}\hfill P_m& =& F_m-F_{m+1}\hfill \\ & =& \sum _{w=q_m}^{w=q_{m+1}}\left({\displaystyle \genfrac{}{}{0pt}{}{W}{w}}\right)p_i{(r,\gamma )}^w{[1-p_i(r,\gamma )]}^{W-w}.\hfill \end{array}$$

(15)

To visualize this theoretical result, we plot the distribution of event sizes as the node degrees and their distances to the resetting node with different resetting rate $\gamma$. In Figure 6a, we observe a pronounced asymmetry in the fluctuation landscape:the nodes on the left exhibit a narrow range of event sizes, typically within $\left|m\right|\le 3$, and the probabilities for larger event sizes $(m>5)$ rapidly decay below ${10}^{-6}$. This reflects the fact that hubs receive a relatively stable and high average flux of walkers, making large deviations statistically rare which leads to the probability for the occurrence of events of large sizes are relatively small .This is consistent with the results shown in Figure 1. the node in the middle begin to show moderate broadening in event size distribution. Occasional events up to $m\approx 6$ are observed, with probabilities $P_m\approx {10}^{-6}$, indicating increased flux volatility compared to the hubs. the nodes on the right display the most dramatic fluctuations. Some of these nodes experience extreme events as large as $m=7$ or more, with non-negligible probabilities (e.g., $P_6>{10}^{-5}$ for certain nodes). This suggests that, although these nodes are rarely visited in the dynamical process performed by random walkers, the few walker arrivals they do receive are sufficient to produce disproportionately large statistical deviations which implies that even occasional visits by a few random walkers contribute to extremely events of large sizes.

Compared to Figure 6a,b demonstrates that event size fluctuations under resetting random walks become more concentrated. Notably, the probability of unit event size $(m=1)$ increases significantly, indicating that stochastic resetting effectively reduces size variability. Furthermore, we observe a marked decrease in extreme events $(m>4)$, consistent with the trends shown in Figure 3. These results collectively demonstrate that resetting suppresses both overall fluctuations and extreme event probabilities.

Interestingly, these effects exhibit degree and distance dependence: while small-degree and peripheral nodes show clear changes in their event size distributions, hubs display less pronounced variations. This suggests that stochastic resetting differentially affects nodes based on their connectivity in BA scale-free networks. Our findings confirm earlier analytical predictions that those nodes with small degrees and far from the resetting node are particularly vulnerable to extreme fluctuations. Furthermore, as for which factor has a greater impact on the event size fluctuations-the size of the degree and the distance to the resetting node, the research of another sorting method:sorting in ascending order of the distance to the resetting node, with nodes at the same distance sorted in the descending order of their degrees is considered. We found that the result of this sorting method are not significantly different from the former case, meaning that the two factors have little difference in their impact on the event size fluctuations.

The resetting mechanism produces two key effects. First, it diminishes the dominance of central nodes. Second, it makes fluctuations in peripheral regions more discernible while simultaneously mitigating traffic variations network-wide. The specific order visualization in our research reveals these dynamics more clearly, showing a rightward shift in the probability distribution of large events. This spatial redistribution indicates that extreme fluctuations become increasingly concentrated in the network’s periphery, highlighting the systemic vulnerability of low-degree and peripheral nodes under stochastic perturbations.

This analytical approach serves dual purposes: it not only validates previous theoretical predictions but also offers a powerful diagnostic tool for pinpointing fluctuation hotspots and evaluating network resilience in heterogeneous systems.

5. Conclusions

This work investigates extreme events occurring on network nodes under the mechanism of stochastic resetting in random walks. We explore how resetting to a central node affects both the probability of extreme events and the fluctuations in event sizes across the network.

Our key finding is that stochastic resetting serves as an effective control mechanism for extreme events in complex networks. The resetting process significantly reduces the probability of EE across all nodes, with complete suppression occurring when the resetting rate $\gamma =1$. This reduction is particularly pronounced for low-degree nodes and nodes distant from the resetting node, which are naturally more vulnerable to extreme fluctuations in standard random walks.

We have demonstrated that resetting not only lowers extreme event probabilities but also concentrates the distribution of event sizes. The fluctuations become more tightly distributed around the mean, with a marked decrease in extreme deviations. This effect arises because resetting redistributes walkers across the network, preventing the accumulation of large fluctuations at any single node.

An important insight from our work is the differential impact of resetting based on both node degree and distance to the resetting node. While all nodes benefit from reduced extreme events, the effect is more substantial for peripheral nodes—those with low degree and located far from the resetting node—compared to hubs or nearby nodes. This suggests that resetting could be particularly valuable for protecting vulnerable nodes in infrastructure networks, where extreme events (such as traffic congestion or overloads) are most likely to occur at poorly connected or remote nodes.

Our analytical framework, validated by numerical simulations, provides a foundation for designing optimal resetting strategies in real-world networks. The results suggest that by carefully tuning the resetting rate $\gamma$, network operators can achieve a desired balance between search efficiency (often improved by resetting) and extreme event mitigation. This could have practical applications in designing robust communication networks, transportation systems, and power grids where extreme fluctuations need to be controlled.

Future work could extend these findings to other resetting protocols, including multiple resetting nodes [40,41] or resetting biased random walks [21]. The methods developed here may also be applicable to studying extreme events in other stochastic processes with resetting mechanisms.