Role of Recovery in Evolving Protection against Systemic Risk: A Mechanical Perspective in Network-Agent Dynamics

We propose a model of evolving protection against systemic risk related to recovery. Using the failure potential in network-agent dynamics, we present a process-based simulation that provides insights into alternative interventions and their mechanical uniqueness. The fundamental operating principle of this model is that computation allows greater emphasis on optimizing the recovery within the general regularity of random network dynamics. The rules and processes that are used here could be regarded as useful techniques in systemic risk measurement relative to numerical failure reduction analyses.


Background.
Various contemporary studies have argued that systemic risk and abrupt failure events are related to the highly interconnected systems and networked structures created by agents [1]. Based on a simple set of properties, observations derived from several models indicate that the proportion of protection between nodes in the network can be described as a probability. is is related to how systemic risk should be coped with rather than simply predicting it based on the probability of failure (e.g., a failed bank in a financial system, asymptomatic transmission of a disease, and cognitive bias in decision-making). e risk of propagation is higher than that of independent failure events, extending to interdependent ones, which we refer to as cascading failures among system components [2]. Many scenarios that arise in simulations should be regarded not as indicating uncertainty or mistakes but rather as the consequences of inappropriate settings and interactions. In particular, proper protection against systemic risk with system components can be evolved heuristically through strategy dynamics (social learning and exploration) as a potential means of limiting failure. e present model expands this concept to investigate nonlinear randomness effects due to delayed responses, which may result in sensitivity to small changes that are difficult to prepare or manage [3]. To assess the profound implications that this approach may have for our understanding of dynamic behavior, including protection processes, this study investigates the influence of the necessary heuristics through which a proper response could mediate risk diffusion before a system completely fails [4], as discussed below.

Literature Review.
By investigating how the complexity of networked structures underpins real-world systemic phenomena, various simulation studies have identified implications for individual robustness, the propagation of systemic risk, protection flow, and collective behavior across networks [5]. A distinguishing feature of such phenomena is that they emerge from the complex interactions among individual elements in a system or from their associations with each other [3]. e effect of context-varying mechanical flux on a system's risk is highly complex. e possibility of quantifying such a risk needs to be evaluated in consideration of the distortions and patterns of such effects [4]. e investigations involve a variety of information, e.g., social contacts that are favored as the infection and spreading route, which in turn can be used to infer the characteristics of the underlying networks [6]. As demonstrated by the existence history (e.g., the financial crisis in 2008, the outbreak of COVID-19 in 2019), a reproduction (i.e., failure, bias, or virus) varies greatly from individual to individual, which is generally believed to affect the spreading dynamics significantly [7]. Whether such individual inhomogeneity aggravates the outbreak is a challenging question, and the answer depends on the specific model [8,9]. In particular, the connectivity patterns of individuals are key to understanding [10] how networks are structured and communicated with each other [11]. Other network properties that have been investigated include the concept of evolutionary dynamics, which helps characterize and understand the architecture of artificial systems in relation to the network properties [12]. As most tools for laying out networks are variants of an algorithm, it is difficult to use them to explore how the conditions of a network affect the network's dynamics [13]. e assessment process can be used to make macroscale observations for input performance, while approaches for microscale evaluations to simultaneously obtain more detailed insights must be treated within the structure of the network itself [14]. Several studies have reported such structures in terms of both microscale (e.g., individual incentives and relative gain versus effort) and macroscale (e.g., institutional competition and central intervention) behavior. For example, evolutionary explanations of systemic risk demonstrate how optimal decision makers are constrained when creating biased estimates of their capability and show how individuals alter their strategies in response to perceptions of resource value [15]. Standard evolutionary models in complex environments show that potentially different biases in decision-making expose different experimental groups at different transition probabilities [16]. A recent study found that by employing strong mitigation (i.e., social distancing and isolation of confirmed cases as guided by risk diffusion testing) related to the different response strategies, an outbreak can be suppressed to levels below the normal critical care capacity [10]. Although a triggered cascade can evolve over a certain time scale (i.e., days), it can be mitigated with intervention by the central system [17]. Evidence from many nowcasting and forecasting estimates indicates that in the absence of prevention and control measures other than simply isolating the risk cases, the probability of continued transmission with the projected trajectory remains high (exponential growth of the number of infections) [7,8]. us, there is an urgent need to reduce propagation rates and control the growth of this risk to reduce not only the peak demand on the system but also the total number of eventually affected individuals [18].

Gap Statement.
e computational modeling technique shows no bridge between the dynamics of agent nodes (with the vertex as a fundamental element) and the emergent properties of failure in recovery [19] (note that "recovery" in the context of financial systemic risk often refers to the fraction of a loan that is recovered after the default of the counterparty; here, it refers to a different quantity called "recovery time delay" similar to the concept of intervention). Most tools for laying out networks are variants of an algorithm and hence cannot easily be used to explore how the conditions affect the dynamics of the network [20] owing to the following factors: (1) Many of them take the form of a theoretical explanatory insight constructed in response to a hypothetical assumption. (2) e type and number of individuals are arbitrary or left undefined. (3) Validation with respect to stylized mechanical parameters cannot explain their potential over-parameterization. (4) ere is an extended transient or burn-in phase that is discarded before analysis [21]. (5) Most importantly, the time units for many of these models may have no clear interpretation. To address these issues, we extend the model to fit an estimation of macro-/microscale variables, such as protections and interventions. e assessment process can be used to make macro-or microscale observations of input performance, while approaches for improving the recovery delay and obtaining more detailed insights should be investigated in the structure of the network itself [22]. is requires the combination of large repositories to construct representations of trajectories that can be analyzed at different scales and from different perspectives [23]. Indeed, the mechanisms and serial algorithms that underpin our understanding of systemic risk in networked agents must be evaluated through various means. Accordingly, we can establish a common ground for the integration of knowledge and methodologies with consistent definitions and reconcile the approaches for studying networks from various fields, which will intuitively enable us to face all the difficulties and pitfalls that are inherent in interdisciplinary work.

Purpose.
is study develops a modeling framework that can account for quantitative measurement in agented networks, allowing us to explore how the recovery time delay affects the risk potential in both macro-and microscale cases. To regard agent dynamics as a random network, this model follows the standard approach in agent-network modeling, where by default, a small event (agent n is hit by a shock at time t) can trigger the initial passage in a risk diffusion process. e mechanism tests the clear implications for different values of the interaction, including interruption, and how protection may be related to a set of interconnectedness with mitigation entities against failure potential, rather than solely focusing on cascading events.

Value.
With the objectives of better risk assessment and effective risk reduction, this model will enable us to not only directly observe the spread of failure in agented network industries but also better understand how protection can be accomplished through intervention. is work is related to the mainstream of research that contributes to the discussion of systemic robustness.

Model (Process-Based Methodology)
e concept of network generation in this study is based on random processes, such as those described in graph theory [24]. Although an arbitrary construction cannot fully capture the local characteristics of individuals observed in real-2 Complexity world networks, everyone in the world is connected to everyone else through a chain of mutual acquaintances or even stronger relationships [25]. To examine the systemic risks that may result from failures originating and cascading on such contact networks, together with how the networked agents may be expected to protect themselves against failure cascades, we consider an agent-based systematic-risk model with evolving protection strategies that have been developed at the International Institute for Applied Systems Analysis (IIASA; Ulf Dieckmann). is model enables agent-based simulations [26], beginning with the simple assumption that pairs of nodes can be randomly connected by an edge with a given connection probability (Table 1). Using a parameter to evaluate the impact of risk on the networked agents, we can estimate the influence of primary risk along the structure as a general failure property [27]. rough scaling for the different evolutionary (Table 2) and nonevolutionary components (Table 3), each step computes a new entity and generates a new proportion in relation to the intervention [3] (Table 4).

Operating Principles.
First, each agent is characterized by two values: capital and strategy. An agent that has lost all of its capital is regarded as failed. Initially, only one agent is given as failed. Each agent is assigned to a node of the network, which is given with a number n of nodes, a connection probability p d , and a resultant adjacency matrix A. Second, at each time step, units of payoff (c) are added to each agent's capital, of which fractions f m and f p are spent on maintenance and protection. e updated capital is expressed as 1 + (1 − f m − f p )c. ird, a failure potential can originate at each node with probability p n and can propagate through each link with probability p l . e system may reach one of the two states after the appearance of a random probability (initially failed node of p n through the link of p l ) and before the second probability occurs. Each failed node with a failure potential becomes a failed node with probability p r according to the equation below. A failed node losses its capital and remains in a failed state for one time step p p � p p,max /(1 + c p,1/2 /(f p c)). We can measure capital in units of c p,1/2 to obtain p p � p p,max /(1 + 1/(f p c)). Fourth, each agent's strategy values, namely, f p0 and f p1 , are updated through social learning and exploration. Social learning represents the process of choosing a random role model with probability p and imitating the role model's strategy. e probability of imitation is p i � 1/(1+ exp(− sΔc)), where Δc is the difference between the role model's capital and the focal agent's capital and s is the selection intensity. en, exploration is performed by altering a randomly chosen strategy with probability p e with a random value chosen from a normal distribution with mean 0 and standard deviation σ e [f(x|μ, σ 2 )]. e protection level for each agent is determined according to the heuristics f p � f p0 + f p1 C, where C is the eigenvector centrality of the graph, which is a measure of the influence of a node in a network. If A denotes the adjacency matrix of the graph, the eigenvector centrality or x must satisfy the equation Ax � λx, where the vector is normalized to 1. We can normalize this vector to a maximum value, which brings the vector components closer to 1; here, f p1 is important for the measurement of f p (note that we did not explicitly demonstrate this because it has been shown in a previous study [1]). e value for f p must be truncated to the interval (0, 1 − f m ). Finally, the failure potential, which lasts for only one time step (r t � 1), is controlled by another variable, which we refer to as the recovery time delay (t r ∈ [0, ∞]). e details of these process-based mechanisms are presented below.

Process 1. Basic Properties of Created Networks.
Individuals in this model are considered as vertices (nodes), and sets of two elements are drawn as a line (edge) connecting two vertices. Data are stored in the nodes, and the edges represent the connections between them, although they can also store data. e edges between the nodes can describe any connection (adjacency) between them (Figure 1(a)). e nodes can contain any amount of data that is assigned to them, and the edges include the data of the strength of the connection they represent.
Connectivity is another essential property of this structure. A disconnected network has some vertices that cannot be reached by the edges from any other vertex (Figure 1(b)).
A disconnected network might have one vertex connected to no edges at all, or it might have two connected networks that have no connection between them. Similarly, a connected network has no disconnected vertices; thus, a metric called connectivity is used to describe a network as a whole, and it depends on the information being presented, usually identified by (n, d[p])]. Networks also have additional properties, i.e., edges can have a direction, such that a relationship    between two nodes is only one-way and not reciprocal. However, in the present model, we used an undirected network, the edges of which have no direction, because in our case, edges are drawn between two individual nodes who have met; hence, all relationships being represented are reciprocal. us, the relationship created (network) is undirected and begins with edges randomly drawn between one pair of nodes at a time. For example, four nodes may have edges between them, as shown in Figure 2.
By contrast, if we take an adjacency matrix, we may consider the rows and columns of the matrix to be labeled by the vertices (nodes), giving us one, two, three, or four vertices here. We may use any actual labeling; here, we denote our adjacency matrix by A. e definition of our matrix is that an entry in row m and column n is equal to 1 if there is an edge between m and n, and 0 otherwise: Note that the adjacency matrix of a network contains all the information contained in that network. Similarly, note that the presented network is given a random appearance because of the way the computer generates the visuals. If the program was to run a second time, a different picture would be generated; however, regardless of how it is run, the same relationship holds between the vertices (nodes) and the edges (line), resulting in the same degrees: (2) e resulting m × n matrix is obtained as [A � G(n, p[d])], and any node (n) can be randomly linked to any other node. Because the collection of nodes influences the connection probability [p[d] ∈ (0, 1)], the model investigates the distribution of the connection in the network (probability of degree) ( Figure 3). Figure 4 is based on the standard representation in network theory. Note that the vertices have degree d . is d-regular network has degree d for all vertices, and the resulting matrix can be compared to another well-known random network called the Erdös-Rényi network [28], obtained from [A � G(n, p)], shown in Figure 5, where any node (n) can be randomly exposed to any other nodes, creating random connections.
Technically, with respect to the collection of nodes influenced by the connection probability [p ∈ (0, 1)], it can be observed that both cases with a higher degree (on the right-hand side of the plots) show a higher probability of being connected to an agent. However, for simplicity, the model presents a random regular structure owing to the connection degree [d ∈ (0, n − 1)] because the network has all vertices of a certain degree (the same connection probability). In general, various network types (e.g., institutions, firms, banks, food distributors, and supply chains) have ambiguous effects in relation to individual ties because they allow for different ways of diversifying risk; this is particularly true as they influence each other in the same procedure of originating strategies and distributing security [18]. is is considered in the literature to provide strong empirical evidence for a technical approach for comparison with another network (note the lattice and calculation of eigenvector_centrality (Table 5)) to obtain a simple numerical solution ( Figure 6).
In the following section, analyses are conducted for the case of a random regular network such that agents within the strategy have the same level of diversification [17]. 4 Complexity (p j). Each node can be in one of the two states: not failed or failed. All nodes are initially without failure.
Next, the probability of failure can be determined by the number of links from the node of the specification scaled by R/S. If we retain the individual characteristics as constant (K), R/S is equal to the risk (failure probability: p ∈ (0, 1)) as a function of the connectivity created by the connections (λx), and R is set to K/S. Nodes at lower (higher) links can be expected to have a lower (higher) connectivity to their risk.

Complexity
In other words, if we remove nodes from the network, the bias is reduced where the links are lower, even if they retain their individual characteristics throughout the process.

Process 3. Impose Protection against Systemic Risk.
Along with the basic intuitions mentioned above, protection dynamics is also implemented. First, we divide the program into the subdynamics (payoff, failure, and strategy). en, the result of each subdynamics is saved. ese subdynamics are trivial problems that add complexity to the dynamics. To implement each of these, we use simple equations. ese equations combine the previously computed variables with the newly added or computed variables, for example, In the example given in Figure 9, we use the values already stored in the table to compute new variables. is technique is often referred to as memorization ( Figure 9).

Payoff
Dynamics. An agent is associated with each node and is characterized by its capital and strategy as follows: for each time step, each agent receives one unit of payoff, which is added to that agent's capital c, of which fractions f m and f p are spent on maintenance and protection, respectively; thus, the capital value is updated We used an elementwise computation with arrays for vectorization ( v → ) instead of a loop: e previous initial random network property is as follows: us, the applied output given the payoff dynamics becomes where the vector ( v → � payoff_dynamics) components are equal to matrix B. is product is equal to B v → ( Figure 10).

Failure
Dynamics. e failure potential can originate at each node with probability p n ∈ [0, 1], and it also propagates along each link with probability p l ∈ [0, 1] at each time step. e failure potential turns into failure with probability 1 − p p depending on the agent's investment in protection: a possible choice is [p p � p p,max /(1+ c p,1/2 /(f p c))], where protection (p p ) is equal to the applied (saturation) function. Here, p p,max is a designated protection maximum, c p,1/2 denotes an allocated reference point, and f p c represents an evolutionary protection level multiplied by the updated capital: e applied output given the failure dynamics is as follows: where the vector ( u → � failure_dynamics) components are equal to matrix B. is product is equal to B u → . In this section, a prewritten function (regular) is used to create a short iterate 1 � D array for vectorization ( u → ) instead of using the adjacency matrix directly.
is substitution shortens the loop. e failure lasts for one time step (default) and results in the loss of an agent's capital ( Figure 11).  10 Complexity

Strategy Dynamics. Each agent chooses its protection level according to the heuristics
For initialization of the strategy values, two arrays are added for vectorization ((f p0 � w i ), (f p1 C � Cw i )): where w → i1 is a vectorization as the designated strategy of f p0 and w → ii is a vectorization as the designated strategy of f p1 , multiplied by the eigenvector centrality from the random graph (C), which is a measure of the centrality of the agent's node, normalized to the interval (0, 1): e eigenvector centrality for node I is Ax � λx, where A is the matrix of the network with eigenvalue λ. e strategy values f p0 and f p1 evolve through social learning and strategy exploration as follows: at each time step, each agent randomly chooses another agent as a role model with probability p r ∈ [0, 1] and imitates that agent's strategy values with the following probability: where p I is the probability of acceptance of the role model for imitation, π f is the capital of the focal individual, π r is the capital of the role individual, e denotes the exponential, and ω is the intensity of the selection (ω < 1 � weak selection; ω ∞ � strong selection). e focal individual imitates the strategy of the nearby role individual, comparing its new capital (large Δπ � large capital difference; small Δπ � small capital difference); then, the focal individual chooses to imitate the strategy of the role individual. In relation to this imitation, a temporary matrix is employed to avoid changing and using one matrix in the loop ( Figure 12).
Finally, at each time step, each agent with probability p e ∈ [0, 1] randomly chooses one of its two strategy values and alters it with a normally distributed increment with mean 0 and the following standard deviation ( Figure 13):

Process 4. Reset Failure (and/or Protection) Potential.
A failure lasts for one time step and results in the loss of an agent's capital, according to the reset failure potential and/or failure, as follows: where random(n) < rec � 1 denotes randomly chosen individuals with a certain probability [random(n)], and the failure potential of any individual is (1 − p p ), which is chosen according to a certain probability that approaches 0: By default, this recovery rate is implemented by resetting the failure potential after every r t � 1 time step. At the same time, to control this intervention, we allowed the number of time steps to be controlled by another parameter (t r ∈ [0, ∞]), which represents the recovery time delay (Figure 14).

Mechanical Insight of the Programming.
e broader objective of this step-by-step procedure is to show the computerized process underlying the fundamental mechanisms that are used. e Tables 1-5 and matrix mentioned above account for the details of the rules and procedures (see Appendix 6 for details on the coding). In relation to the technical insights of the implementation, we present the details in Figure 15 that may interest program developers ( Figure 15). e versions on the left-and right-hand sides in Figure 15 represent vectorization using a loop and an array, respectively. Both versions are used in the development of our programming, as can be seen in the code book, which presents simple calculations. Specifically, following the stepby-step procedures of the operation, the version on the lefthand side can be recommended in general; however, it requires more time to simulate many individuals across time steps (see Output in Figure 15). us, in some dynamics implementations of models, we used the version on the right-hand side owing to its efficiency.

Performance (Results)
Given the set of features in relation to the dynamics presented according to the model description, the results of the simulation show the fundamental characteristics of risk diffusion in a randomly networked system and present a framework that enables us to examine the assumptions efficiently while imposing realistic protections against failure. In particular, the simulation characterizes the role of the recovery time delay from the observations of applied dynamics over time to indicate how failure spreads.
First, individuals in the model are considered as vertices, and a set of two elements is drawn as edges connecting two vertices in relation to the information given in the graph.
is representation involves two parameters: the number of nodes (n) and the probability that an edge is present (d). In network analysis, indicators of centrality identify the most important vertices within a graph, and their applications include the identification of the most influential node(s) in a network. e eigenvector centrality of a node is defined as Ax � λx, where A represents the adjacency matrix of the network with eigenvalue λ. e principal has an entry for each of the n vertices. e larger the entry for each vertex, the higher it's the ranking with respect to the eigenvector centrality ( Figure 16). Figure 16(a) represents the influences of failures, computed as the frequency of failure nodes in relation to their initial and final states, with the computed relative frequencies f j � h j /N, h j � k�1,...,N;X (k) �j 1, where h j is the frequency of the failure state (j) and the percentage of occurrences of that failure outcome in the statistical ensemble corresponds to the relative frequency f j . Figure 16 Figure 15: Numerical results using a loop (left side) and an array (right side).

Complexity 13
failures for 100 realizations (black lines) until the number of failed nodes is stabilized. is implies that the greater the propagation from a higher rank (with respect to degree or centrality), the greater the area drawn between the failed nodes (red) and the nonfailed nodes (blue). is result intuitively implies that, primarily, higher-degree agents have greater exposure to cascading failure risk than lowerdegree agents; this increases the potential for cascading failure [29].

Protection against Failure (Imposing Realistic Dynamics).
In line with our proposition, protection dynamics was applied in reference to the risk of failure. is model allows an agent to make a large investment in protection. Note that we assume that the failure potential from the dynamics becomes failure with probability 1 − p p depending on the agent's investment in protection; p p � p p,max /(1 + c p,1/2 / (f p c)) when a failure lasts for one time step (Figure 17). e probabilities of these scenarios can be proven as follows.
First, scenario A's failure potential 1 − p p is 0.527 because the investment in protection p p is 0.473, based on the possible choice p p : p p,max /(1 + c p,1/2 /(f p c)) � 1/(1 + 1/(0.9 * 1)), for certain parameter values (p p,max � 1, c p,1/2 � 1, p r , p e � 0.9, C � 1). Second, scenario B's failure potential 1 − p p becomes 0.953 because the investment in protection p p is 0.047 based on the possible choice p p : p p,max /(1 + c p,1/2 / (f p c)) � 0.1/(1 + 1/(0.9 * 1)), for certain parameter values (p p,max � 0.1, c p,1/2 � 1, p r , p e � 0.9, C � 1). As can be seen from the set on the left-hand side of the plot, this scenario has a better result than the set on the righthand side, corresponding to the coexistence of failure and the absence of the 20 eurist with time. e results imply that contagion and systemic risk are likely to be intensified by the parameter setting (i.e., p p,max , c p,1/2 , and f p ), resulting in a significant failure cost [30]. Note that the failure potential 1 − p p also changes in relation to the choice of strategy values f p0 and f p1 based on social learning (p r ∈ [0, 1]) and exploration (p e ∈ [0, 1]) because these probabilities can lead to different levels of protection in f p [1].

Role of Recovery (Imposing Recovery Time Delay).
Here, we present a specific characterization to determine whether an alternative intervention could affect the failure trends observed in coexistence scenario A (p p,max � 1, c p,1/2 � 1, p r , p e � 0.9, C � 1). e model shown in Figure 18 implies that contagion and systemic risk are also likely to be intensified with the recovery time delay, resulting in a significant failure cost. We observe that immediate or delayed intervention is associated with microscale propagation criteria for each node (Figure 18).
In Figure 18, the individuals can be distinguished clearly. Obviously, nodes that have immediate recovery (plots on the left-hand side) appear to have the potential to protect against the propagation of failure; on the contrary, nodes with a malfunction (plots on the right-hand side) do not have such potential. e simulation consisted of repeated trials, and each trial had two possible outcomes: failure (� red) and absence of failure (� blue). e probability of failure is the same for every trial, as in flipping a coin n times, based on the binomial variable that we defined in the model's basic structure. e probability of failure in each trial is given by P(X) � n!/r!(n − r)!p r (1 − p) n− r � C(n, r)p r (1 − p) n− r , where n is the total number of trials, r is the total number of failure events, and p is the probability of failure in a single trial.
In the plots on the left-hand side in Figure 18, the probability of absence of failure is approximately 0.6, and the probability of failure is 0.4 (the case probability for Scenario A is calculated by the number of failures from the time steps: f j � h j /N, h j � k�1,...,N;X (k) �j 1). We assume a random variable X, which is equal to the number of failures after a certain number of time steps. One of the first conditions of this result is that it consists of a finite number of independent trials. is means that the probability of obtaining failure or absence of failure in each trial is independent of whether we obtained failure in a previous trial. us, in the case of the plots on the left-hand side with recovery at every time step (immediate intervention), the results of the simulation show independent trials. Another condition is that each trial clearly has one of the two discrete outcomes in which the variable X should be clearly classified as showing either failure or absence of failure with a given fixed number of trials. en, the final condition of the probability of failure and absence of failure for each trail is constant, as we have already obtained measurements for each trial from the cases of Scenarios A and B.
Proof of the role of the recovery rate. We now examine the frequency of failed agents observed on the left-hand side in Figure 18. Suppose that at a given time, there are failed nodes (N f ) and nonfailed nodes (N − N f ). e number of failed nodes at t + 1 is given by where p p denotes the protection probability (arbitrarily designated as 0.5 only for this numerical calculation); conversely, the number of failed nodes is Failure is propagated through an existing link with probability p l ; at the same time, a link exists with a certain probability through the created network (p RR , RR � random regular graph). erefore, if we impose the condition that the failure potential can propagate through each link, and the propagation probability p l is less than 1 (when the original probability p n � constant), this equation can be modified simply by replacing p p with For example, when we arbitrarily consider the protection p p � 0.5 and failure propagation (p l ) through a random network (p RR ) to be sufficiently high (p l p RR � 0.9), the protection influenced by the failure propagation p p ′ becomes small because 1 − (1 − 0.5)(1 − (1 − 0.9)) N f � 0.09. On the contrary, if a weak failure propagation occurs, such as p l p RR � 0.1, the protection influenced by the failure propagation p p ′ increases because However, in the plots on the right-hand side in Figure 18, the probability is no longer the same but changes from trial to trial. We have the variable X(p l ), which is equal to the number of failures from a designated population. It seems to represent the same operation because, with it, each trial can be classified as either a failure or an absence of failure over a fixed number of trials (t � 16). At the same time, there is the probability that the variable p l p RR is not constant for each trial owing to the recovery time delay, which does not consist of independent trials. e probability of failure or absence of failure in the first trial would be equal to the total number of individuals between the two simulation cases: Small Large rec_t p p = p p,max /(1 + c p,1/2 /(f p c)), p p,max = 1, c p,1/2 = 1, p r = 0.9, p e = 0.9 Scenario A: large p p,max , large f p , and large c p,1/2 → enough protection → coexistence p p = p p,max /(1 + c p,1/2 /(f p c)), p p,max = 1, c p,1/2 = 1, p r = 0.9, p e = 0.9 Scenario C: small p p,max , large f p , and large c p,1/2 → low protection → all failure ). e graph on the right-hand side represents the dynamics in a random regular network, such that the node number � random label for each node and the node color � state (failure � red, absence of failure � blue) (see Appendix 3 for the results of the entire simulation).
However, the probability of the second trial (and the following trial) would not be the same because the simulation of the right-hand side case depends on what happened in the first (previous) trial: In other words, each trial is carried out without replacement, which results in an exponentially large difference between the two cases. us, this result does not meet the independent condition: the probability present in the following trial depends on what happened in the previous one. Because replacement does not take place, the probability of failure for each trail is also not constant, unlike the simulation on the left-hand side, where the probability of failure is constant: Inspired by the plausible scenarios included in the recovery delay for cases A and C in Figure 18, we focus on the comparison of the parameter p p,max (when recovery delay constant � 1) with more time steps (t�100). is is because of the applied function p p [� p p,max /(1 + c p,1/2 /(f p c))] for this application, which will primarily be decided by p p,max when we consider c p,1/2 /(f p c) and the recovery time as constants ( Figure 19).
is refers to the individual for control by interventions to be protected against the risk of failure when the failure propagation mechanism influences agents' decisions as they are generated. Figure 19 shows a scatter plot constructed using p p,max (Figure 19(a)) with a clear negative correlation between failure and capital. Although it is constructed by the recovery time delay (Figure 19(b)), strong power law relationships are still obeyed in cascading failure, even if the capital values remain strong (see the averaged values for the failure and capital in the lower part of Figure 10). In the following results, this trend (failure and capital) is determined by two parameters (p p,max , rec t); however, their influences do not have the same weight. Here, we should note that intervention applied to recovery can persist for a short period of time, which we could cite as a bias or rationality.

Generalize with Stationarity.
Let us assume that the system has evolved and reaches stationarity, which implies that all variables have nearly constant fluctuation around their mean values. In stationarity, the state of the system can be considered independently of the initial conditions. It is useful to obtain a relationship among the variables in stationarity. First, we denote the number and fraction of failed agents by N f and f, respectively. For stationarity, the average number of failed agents remains constant. Consider the two successive steps t and t + 1. Suppose that at time t, there are N f and N − N f failed and nonfailed nodes, respectively. erefore, the number of nonfailed nodes at t + 1 is N f + p p (N − N f ), and the number of failed nodes is (1 − p p )(N − N f ). Stationarity means that the number of failed (or nonfailed) nodes remain constant (on average), i.e., Let us define q p � 1 − pp. e above formula for f can be written as f � q p /(1 + q p ), which is simple. Here, p p is the average probability of protection among the nonfailed nodes. Note that the protection probability for the failed nodes is zero because they have no capital. In other words, p p � p p (i)/N − Nf, where p p (i) represents the protection probability for node i. In the second equality, f p and c are defined similarly.
is relationship shows that the stationary value for f is bounded in the interval [0, 1] and is a decreasing function of p p . e limiting values seem promising. For p p � 1, i.e., full protection, there is no failed agent. For p p � 0, i.e., no protection, there is a failed agent. Note that at the end of each time step, all the failed agents recover (as an initial setting). If the propagation probability p l p RR is less than 1, this equation should be modified by replacing p p as follows: Again, by defining q p ′ � 1 − p p ′ , we implicitly obtain f. To check the validity of these formulae with a simulation, we increase the total simulation time (t � 1,000); for smaller p l p RR < 1, a longer time is required to reach stationarity ( Figure 20).
Proof of the stationary state. we calculate the average value of the capital in the stationary state to prove the long-term role of the recovery time delay. First, consider the case p l p RR � 1. It can be easily seen that the average value of Complexity the capital in the stationary state satisfies the following equation: where c is the average value of the capital among nonfailed individuals. Combining the above equations, we can obtain c and f p as a function of f in the stationary state. However, it is easier to obtain c as a function of p p , such that the protection probability should be replaced by erefore, we obtain the following (note that in the simulation, we use a truncation of f p to constrain it to the interval [0, 1 − f m ] [see code book]): According to the trajectories (upper part of Figure 20), the outcomes such as failure (or capital) reach a stationary state even if the trajectories from the strategies continue to evolve and exhibit different behaviors (see the inset of the upper plots). erefore, the presence of the detailed dynamic features of the system (capital and failure) indicates the possibility of stationarity in the coevolutionary process. is observation can enable us to gain a sense of the qualitatively different nature of propagation (compared to p p,max ) within systems, such that failure is not limited to protection or strategy but is due to regulation over time.
In other words, there is another rate that continues to increase the bias as the length of the delay causes the cascading failure, even where the individuals still have potential ( Figure 21). As indicated by the protection (Figure 21; bottom, marker of green, nodes that are controlled by p p,max � 0.1 (left-hand side) seem to lose their ability to protect, owing to the propagation of failure. On the contrary, nodes that are controlled by t r � 2 (right-hand plots) seem to maintain their potential, even if they all feature cascading failure. We noted that the different levels of potential (between the two parameters [p p,max , t r]) came from their capital, as marked in blue in Figure 20 (see bottom plots of the figure). us, to reduce the potential ramifications (p p ′ ) of such additional losses to others, an immediate recovery intervention may not only be preferred to the potential damage (1 − p p ′ ) from individual failures but also guarantee that it is strategically possible for even large insolvent individuals to recover losses to uninsured connectors. Before the potential for failure can advance the propagation value, it must be identified, and the recovery value of the individual protection potential (i.e., capital) must be estimated.

Discussion
We presented a simple general model to quantify the protection for mitigating systemic risk, with the recovery time delay as the result. Using a simple set of properties, the observations from this model indicate how probability describes the proportion of protection that can be characterized by how systemic risk should be coped with rather than being predicted by the probability of failure [31].

Summary of the Model.
We proposed the following process-based steps (1-4) to create a model. First, networkagent properties were established: when the basic data structure was constructed, the functionalities of the mechanism began with a specific undirected relationship between agents through agent-based simulations, in a coexisting macroscale structure with individuals interconnected at the microscale level in the network. Next, the primary risk influence was established using a parameter to evaluate the impact of risk on the networked agents; the influence of primary risk was estimated along the structure as a general failure property. en, it was necessary to implement protection against systemic risk by embedding protection dynamics that emphasizes the roles of payoff, failure, and strategy dynamics. Finally, recovery was considered through scaling for the different evolutionary and nonevolutionary components; this is crucial to the way this system functions, where each step computes a new entity and e graph on the right-hand side represents their dynamics in a random regular network; node number � random label of each node marked with a yellow background; lines � connections embedded by eigenvector centrality as the thickness; node color � states (failure � red; absence of failure � blue).
generates a new proportion in relation to the intervention. e implications of this model are as follows.

eoretical Implications.
A suite of plausible dynamics and decentralized bottom-up mechanisms was constructed by establishing appropriate rules for the interaction, within which the system components can self-organize, including mechanisms for ensuring rule compliance (vectorized microscale implementation).
is evolutionary heuristic promotes balance with respect to the interactions [32]. For example, when investment in protection is weak (low investment), a pattern of strong systemic risk emerges in the form of agent failure in networked conditions. By contrast, when investment in protection is strong (high investment), a pattern of protection emerges, with little diversification against all challenges (Figure 17). e simulation also reflected a clear correlation between a set of parameter values (between those for capital and failure and between those for the strategies of f p0 and f p1 ). us, the strategy of social learning could be another crucial factor in resource provision that violates expectations and leads to novel trends with high impact. ese results shed light on the propagation modes. e observed contagion and persistence patterns can be regarded not as a direct causal link but in relation to an accumulated rational driven by interconnectedness [33].
More importantly, in this study, the failure potential also reflects the time delay after the official failure of an individual [34]. Newly damaged individuals may affect others or may recover from the next time step onwards (recovery occurs with the rate t r ). e potential needed for recovery is related to the number of healthy individuals n(t), and if individuals recover at time t r � 1, the protection potential changes as follows: Once the capital is exhausted, the systemic risk increases rapidly.
is is based on suitable real-time intervention related to short-term anticipation of risk flows [35]. However, because regulation occurs over time, another rate itself increases the failure potential even if the capital is still exhausted (note that the recovery rate mainly influences failure and does not have the same weight for capital). From our observations, we assume that the protection (or robustness) of an agent π i decreases (or increases) by an amount proportional to the relative exposure of risk sharing η ij . Consistent with the network of dynamics described above, the dynamics, including the recovery time delay (in time series), can be formalized as follows: For convenience, with respect to the results, k substitutes the protection potential, defined by c(t + 1) � c(t)+ n(t) − [1 − pp(t)], such that the function a determines the extent of the loss caused by the default χ j (τ) when we fix risk sharing η ij . is indicates that the agent has defaulted at a previous time, denoted by the time variable τ. In the simulations, there was no delay during the development of the cascade (τ ≈ 0, τ � t 2 − t 1 ). However, this dynamic differs across the default of the recovery [36]. If the system delays the intervention, the cascade tends to be much larger (following the power law) because the defaults of the neighbors of a given agent are statistically dependent on j, such that their default is more likely to cause the default of others: In Table 6, such a range [f ′ (x)] will be determined by failure propagation through each link in the random regular network (p l p RR ), which was determined by the recovery time delay. We simply proved that this numerical trend of the protection probability defined by p p ′ � 1 − (1 − p p )(1 − (1 − p l p ER ) N f ) could not be the same as that of p p max [when p l p ER � constant ( Figure 20)]. us, the failure influences in these simulations are not identical, as p p max changes according to these parameters, including where increases or decreases occur following a nonlinear curvature [37].

Practical Implications.
We propose a modified version of the protection model against systemic risk in terms of its failure recovery mechanics, as demonstrated above. e mechanisms here are based on a few simple basic rules: a node fails at a time with a probability of failure if its failure potential is greater than or equal to that of its nearest neighbors [1], and it fails owing to the interconnected potential, although it spontaneously recovers at an external recovery or according to an intervention probability [3]. e consequences of the damage are crucial for systemic risk and controllability; however, they have not been explored systematically thus far [38]. We show how the process of embedding and the related recovery times impact the dynamics of failure processes in a network (see the section on the proof of the role of the recovery time delay). erefore, we propose that the extent of the intervention regime in the system can be a source of evolutionary (in)stability, wherein the immediate recovery of components can mitigate damage and propagate failure in dynamic networks [17]. ese decentralized management principles could be applied to logistical and production systems, or even to administrative processes and governance [18].
Governments are often reluctant to resolve insolvency in institutions [e.g., banks, firms, supply chains, and infected virus individual (i.e., Ebola, SARS, and COVID-19)], and they permit them to continue operating despite their negative effects. e length of the attendant delay causes a cascading failure, regardless of each individual's capital (i.e., immunity), which then leads to reductions in network welfare. ese delays may propagate the damage to many other individuals in the network, increasing their fragility and probability of failure [39]. Evidence from governing systems suggests that if individuals' troubles are estimated before their capital becomes negative, institutions could, based on their risk potential, weed out the inefficient or unfortunate individuals to avoid more serious adverse effects [40]. is indicates the importance of resolving challenges to individuals as quickly as possible and developing more rapid responses to certify protected individuals by providing immediate intervention. On the contrary, prompt corrective action can increase the willingness to supply protection to reduce the chances of systemic risk [41]. Given the evolutionary mechanisms shown in this simulation model, we observed that the evolutionary response often obtains a critical value for a plausible protection potential. We demonstrated that, although structures have high potential for individual protection and a strategy to maintain their capital, the function of the interconnected recovery delay makes them weak amplifiers, where unprofitable intervention shows a high bias.

Concluding Remarks
Cascading failure was used to assess how rules and processes can collectively add up to a series of interconnected unsuitable outcomes. In the simulations, although every object was considered to be doing well individually, systemic risks could emerge, as the whole was vast. In other words, even if every single object is considered to be highly fitted and behaving properly (i.e., capital), there is no standard solution for proper management for everyone's benefit. erefore, the elimination of resources needed for recovery requires a strong effort to halt cascades as they begin, when the damage is still small, and the problem may not yet be perceived as threatening. ese issues can and must be treated with a proper (re)design of the system and adoption of management principles, as shown by the simulation results [42].
Many disasters in the real-world result from incorrect thinking and inappropriate system design. To address such risks more completely, a better understanding of proper intervention and resilience is crucial. However, there remains a lack of an effective method for calculating networked risk. e proposed model can facilitate realistic calculation of the interdependence and propagation of risks in a network and how they can be absorbed, and it can mediate such that both the system components and the plausible systemic intervention and outcomes will work well.
We derived a unifying framework for the interplay of observations that embed realistic dynamics such as payoff, failure, strategy, and recovery in a random regular network. e mechanics showed that the complexity should be described by the essential features of the model's processes that capture evolutionary damage spread owing to the coexistence of crucial hypotheses in the system. us, it may be possible to create an account where protection and failure are not static quantities, and propagation is likely to be understood in terms of how frequently the system is in a condition that leads to a large diffusion. is availability of mechanisms raises additional expectations that predictability and controllability are a simple matter of proper system design and operation [43]. More intuitively, this study can provide a better understanding of recovery, where real-time management can overcome instabilities caused by delays in feedback or lack of information.

Data Availability
e data used to support this study can be obtained from the corresponding author upon request.

Conflicts of Interest
e author declares that there are no conflicts of interest regarding the publication of this article.

Supplementary Materials
Appendix 1: the results of the entire simulation for protection dynamics against systemic risk (left-hand side of Figure 17). Appendix 2: the results of the entire simulation for protection dynamics against systemic risk (right-hand side of Figure 17). Appendix 3: the results of the entire simulation for protection dynamics against systemic risk (right-hand side of Figure 18). Appendix 4: stationarity controlled by pp,max: (nodes � 10 * 10, connection d � 9, time step � 1,000). Appendix 5: stationarity controlled by rec_t: (nodes � 10 * 10, connection d � 9, time steps � 1,000). Appendix 6: code book. (Supplementary Materials)