Systemic Risk Contagion in Reconstructed Financial Credit Network within Banking and Firm Sectors on DebtRank Based Model

1is paper is dedicated to building a multilayer financial network within banking sectors and firm sectors (nonbanking) on the balance sheet of two types of agents and to assessing systemic risk contagion in the reconstructed network. Two propagation channels due to interbank credit and counterparty risk via banks’ loans to firms are comprehensively taken into account in systemic risk contagion assessment, which is based on the DebtRankmodel by analyzing the relative loss of each bank’s equity and the vulnerability of the network. 1e computational simulation on how systemic risk contagious process evolves has been conducted, where the possible influential factors of network structure, agent’s initial risk status, external shock ratio, liquidity flow rate, and different layers of the network are considered.1e findings show that the reconstructed network is absolutely vulnerable under the assumed market circumstance without any bailouts and the risk contagion process shows nonlinear behavior. Specifically, when the average degree of the network and the external shock ratio increases, the risk contagion speed becomes relatively high and the resulting negative effects on the network are more intense. Besides, risks originating from the failed firms in bank-firm layer should place more negative effect on the financial system than that only happening in interbank market. Different liquidity rates in financial market could lead to obvious discrepancy of the risk contagion speed and the extent of asset loss. Additionally, the two layers of the network have diverse influences on risk contagious process resulting in totally different banks’ status in each layer.


Introduction
e financial crisis of 2007-2008, originally sparked by the collapse of relatively small investment banks, has caused magnified losses for numerous financial and nonfinancial institutions eventually and ultimately led to a worldwide economic slump.
e Bank for International Settlements (BIS) and the International Monetary Fund (IMF) believed that financial systemic risk would lead to functional disorder of the financial system and actually place serious impact on real economy, due to the interrelationships between them [1]. erefore, the emergence of systemic risk in financial networks has been receiving increasing attention in the literature [2,3] and regulators (Yellen [4]). Generally, systemic risk in financial markets can be seen as a specific risk that a significant fraction of the financial system can no longer perform its function as a credit provider and collapses. Briefly speaking, systemic risk is the notion of contagion or impact that starts from the failure of a financial institution and propagates through the financial system, potentially to the real economy [5,6]. Network places paramount effect on systemic risk and financial stability because shocks can be amplified and transmitted across various channels both directly and indirectly, including credit lending [7][8][9][10], common funding [11,12], balance sheet interlock [13,14]), exposures to common assets [3,[15][16][17], derivative exposures [18], and deleveraging across financial institutions [19,20]. In view of this, network-based systemic risk is extremely harmful and can be potentially applied to evaluate the financial risk, which has already been testified by the academic circles and regulation authorities.
Until now, lots of researches have been done on whether the network structure can amplify or mitigate risk within the financial systems, at both regulatory level [21] and institution level [19,22]. ese achievements have accelerated the research on systemic risk and financial networks for different purposes, especially on the measurement of systemic risk from the network perspective [14,18,[23][24][25]. However, other than the most works on systemic risk in single layer financial network, some works have already extended the notion of systemic risk to be explored in multilayer networks, which is much more complicated than the simple accumulation of the risks in homogeneous single layer network [26]. Concerning the multilayer network-based financial risk analysis, recent literatures have appeared and spurred the future exploration [26][27][28][29][30][31][32][33]. León et al. [29] and Bargigli et al. [30] analyzed the interactions of financial institutions on different markets in Colombia and Italy, respectively. Bluhm and Krahnen [28] designed an agent model of multilayer interbank network with the channels of common asset exposure, direct lending exposures, and fire sales. Montagna and Kok [27] studied the individual contagion layer to systemic risk through the multilayer networks with three layers: long-term direct lending exposures, shortterm direct lending exposures, and common asset exposures. Poledna et al. [26] showed that the risk was greatly underestimated only in a single layer network and quantified the contribution of the four-tier network of Mexican banking system on systemic risk. Peralta and Crisóstomo [31] divided the bank network structure into two layers based on the type of mortgage-free transaction and showed that the speed of risk transmission under the multilayer network structure was more rapid. Korniyenko et al. [32] constructed a multilayer network based on the division of different financial instruments and analyzed the transmission of risks in global financial networks. Li et al. [33] constructed a multiplex network model of banks with arbitrary structural characteristics by long-term and shortterm interbank lending.
In real markets, especially in some emerging markets, banks are the main creditors to the real economy in financial system, but part of the large firms (nonbanking) also has these similar financial functions to the banks, which means that financial crisis could also originate from the real economy and spread to the financial system even the whole economy in back and forth way. erefore, the topology of financial networks can be adopted to identify and quantify systemic collapse systematically, and the networks should not only comprise financial institutions, but also include certain industrial firms [34][35][36]. To our knowledge, only a few works empirically and theoretically study the interrelationship between the financial economy and the real economy [37][38][39][40]. Moreover, De Masi et al. [41] and Miranda and Tabak [42] studied credit networks in Italy and Brazil empirically, while Lux [43] developed a theoretical model of bipartite credit networks. Marotta et al. [40] investigated the evolution of the network structure in Japan and then identified the communities in the composed financial network between banks and firms. Except for the concentration on the structure of bank-firm networks, systemic risk assessment in them has also been put forward. Miranda and Tabak [42] and Aoyama [39] firstly attempted to study systemic risk in credit networks of Japan and Brazil empirically. In these works, several measures for systemic risk have been proposed that are mostly based on the traditional statistics of losses, accompanied by a potential distress during the periods of risk contagion [44][45][46]. However, none of these measures take cascading defaults into account. As a complement solution, network-based theory and measures, such as network centrality approach and its derivations, have been great contributions to understand this issue. erefore, considering the interrelationship between the financial department and the real economy sector, financial risk not only happens when the institutions directly are shocked, but also has adverse impact on the deterioration of credit quality [47] in this complicated financial system. We extend the existing literature by reconstructing and analyzing a synthetic financial network that not only includes interbank liabilities but also covers the lending relationship from banks to firms. e theoretical reconstruction of financial network within banking and firm sectors allows identifying the assessment systemic risk contagion by employing DebtRank algorithm as proposed in Battiston et al. [14] and Bardoscia et al. [48], whose derivation has been approved useful in several works [26,49].
is paper is organized as follows: the two-layer financial credit network model between banks and firms is reconstructed theoretically in Section 2, and the systemic risk contagion mechanism in the reconstructed financial network is presented in Section 3. Section 4 builds a systemic risk contagion model on DebtRank model. In Section 5, simulation analysis is conducted. Section 6 draws the conclusions.

Reconstruction of Financial Network within Banking and Firm Sectors
In order to capture the systemic risk contagion effects in financial system precisely, joint exposures via different channels to the same set of correlated agents should be included. With regard to this, a consistent and extremely close to realistic network structure for each type of agents within different but combined channels should be designed, especially for the financial bipartite network of banks and firms. Some basic and reasonable stylized realities can be inferred from the analysis of comprehensive empirical studies on Italian data [41], Japanese data [37][38][39], Brazil data [42], and Austrian data [50], which have been simply described and explored in some related articles [43,50]. In most of the empirically analyzed articles on bipartite networks, the distribution of degrees is always wider for banks than for the correlated firms [38,50], and the number of links of the networks also depends on the size of each agent in the network system [37,41]. Based on the above literatures, we aim to reconstruct a multilayer financial credit network within banking and firm sector theoretically by adopting the hypothesis and the algorithm proposed in the directional articles of Montagna and Lux [51], De Masi and Gallegati [38], and Lux [43]. A two-layer financial directed network is to be established as G � (V, E) in this article, where V and E represent the nodes and links, respectively, consisting of interbank credit layer and bank-firm credit layer. Basically, the schematic diagram of the reconstructed network can be depicted in Figure 1.
Above all, we set a fix number of banks and firms in our reconstructed network system, and it is certain that number of the firms in real economy exceeds that of banks by two or more orders of magnitude. Meanwhile, we assume that the balance sheet size of the banks and firms both follows Pareto distribution which is well-known insights in reality [43,51]. erefore, interbank credit layer network and bank-firm credit layer network are reconstructed specifically in the following sectors. For all, we denote the number of banks and firms in the network by N b and N f , respectively.

Interbank Credit Network.
In interbank network, the edges are directed and weighted, and the direction of a link from node i to node j means the total amount of money that bank i lends to bank j. Following Montagna and Lux [51], we adopt that simplified scheme of the agents' balance sheet. e assets side A i of each bank can be partitioned into two parts: interbank loans l i and external assets x i (loans to firms), so that the total assets A i can be given as A i � l i + x i . In the same way, the liabilities L i of each bank can be divided into three parts: interbank borrowing b i , deposit d i , and equity e i , so the liabilities side can be obtained as erefore, the links between the banks can be generated by the probability generating functions P(A i , A j ) according to the relative sizes A i and A j of balance sheets [43], which is in accordance with the fact that relatively small and medium sized financial entities usually lend money to much bigger banks so as to reallocate liquidity to the market more efficiently [10,11,52]. erefore, the probability of the interbank credit link between each node i and node j can be drawn as where α, β, d are parameters and A max is the balance sheet size of the largest bank. With the probability generation function, en, each entry of the interbank network adjacency matrix A can be drawn with the probabilities p ij and 1 − p ij so as to make sure that the interbank credit link can be created. Meanwhile, the selected probability generation function, which can be used to capture the right-skew of the degree distribution and the disassortative feature of link formation, has been explored and testified in recent researches [10,53].
Naturally, it is assumed that financial entities should always have more intense links to the banks with high balance sheet size; the volume of credit l ij between each bank i and j can be calculated as l ij � (l i p ij A j / k∈Ω i p ik A k ), where Ω i denotes the set of nodes with a ij � 1. us, the interbank credit network can be constructed theoretically.

Bank-Firm Credit Network.
On bank-firm credit network, let the average number of links of firms be λ f , so the average number of links for each bank can be λ b � λ f (N f /N b ), but does not follow uniformly and shows the means of linking probabilities across the banks and firms varying by their balance sheet sizes. erefore, the procedure of reconstructing the bank-firm financial network should be conducted on banks and firms separately.
For banks, each one's total assets ) that the balance sheet sizes of banks could be got from a truncated Pareto distribution randomly, where L and H are the lower and the upper limit, respectively, and α is the Pareto index. Besides, we also assume that the degree of each bank (the links number of each bank) is distributed proportionally to its balance sheets sizes, so that bank i, i ∈ (1, N b ) has an expected degree λ i � λ b A i . erefore, the average expected degree across the banking sectors can be calculated as representing the average degree λ b divided by the mean balance sheet size A i within the financial system.
For firms, the distribution of loan sizes and degrees can also be consistent with those obtained for banks and the typical facts represented above. erefore, the mean size of firms' loan can be obtained as where θ is the proportion of the assets as loans to firms in bank's balance sheets. Supposing that the distribution of firm size also follows a Pareto distribution with the same parameter α as for banks, the truncated Pareto distribution for firms can be therefore obtained by adopting the minimum threshold l and maximum threshold h in the same way that are l � θL(N b /N f ) and h � θH(N b /N f ). For the number of firms being absolutely larger than that of banks, it is guaranteed that the loans by firms should appear on the left side of the size distribution of banks, with the distribution of aggregate loans in banks' balance sheet. us, the realization for the ensemble of firms can be obtained just like the way for banks, and the number of links of each firm j � 1, 2, . . . , N f can be drawn with the above parameters λ j � λ f B i , so that the average number of that can be obtained as λ f � λ f /B i , where λ f and B i are the average degree and the average loan size across firms.
Overall, in addition to the total loans from the banks to firms, the obtained numbers of the degrees for all banks and firms should be roughly accordingly set up to the initial expectation. In order to make it more consistent with the reality, especially for the difference of both numbers exactly, the approach of static model for network reconstruction is applied [54], which is operated by taking the minimum of the aggregate links and adding connections one after another. From this, the steps of network reconstruction can be conducted as follows: first, each bank or firm will be assigned a weight according to the realization of its degree. en, we pick one node from the weighted ensemble of banks and firms, respectively, and create a link between them. Hereafter, the previous used link from the distribution of degrees of banks and firms will be excluded from the set of possible links each time, and the same operation of building the connection will be preceded until the available links from either banks' or firms' side run out. At last, the total loan amount of each bank is allocated to its borrowers (firms) proportionally to their loan sizes by Pareto distribution of firms. If so, this could guarantee the possible increase of degrees with loan size within firms to large extent, avoiding the potential errors that the exact numbers would not be absolutely the sum pertinent loan position of banks that have been generated independently.

Systemic Risk Contagion Mechanisms in Reconstructed Network
On the existing literatures, we go through the whole banking and firm sectors and analyze how default or distress is evolving in the two-layer financial credit network. Normally, it is assumed that external shock on the financial network system should damage the equity of the shocked firm(s) or bank(s). Moreover, it is expected that the shocked firms or banks should linearly transfer the loss of equity to the creditor banks. In the network, the risk exposure is supposed to be short-term loan, and the credit from banks to firms is assumed to be a highly liquid asset for creditors. Normally, if banks are faced with liquidity risk, credit risk, and the parallel influences, they could adopt several solutions to avoid these possible risks. Meanwhile, if a bank suffered from liquidity crisis, regardless of the solutions to be taken by the bank, negative effect would be placed on for different reasons, which may be resulting from high interbank interest's rate, asset devaluation for fire sale, etc. In general, the diagram of systemic risk contagion mechanism in financial credit network within banks and firms can be portrayed as shown in Figure 2.
According to the diagram of the systemic risk contagion mechanism, the institutions face losses in all classes of assets, both interbank market and the bank-firm financial market, and the framework of dynamics of the relative loss in equity for each institution as well as the whole system is concerned. Systemic risk contagious process and how it places impact on the stability of the financial network system can be depicted as the following two rounds.
(1) e first round of shock is on the external part of the assets. Initially, it is supposed that one or more nodes in the financial network that may be firms or (and) banks suffered from external shock and we can let m be the shock ratio, representing the loss extent of the shocked agents' external assets.
For banking system, after the first round of shock, shocked banks may lose part of the assets amount to mV OI (0), where m represents the extent of the damage caused by the first round shock. However, the liabilities will be unchanged for the book value regardless of how the shocks are going. In the same way, the equity of the shocked firms should be also damaged and suffer from the loss of Due to the external shock, even if the shocked bank(s) should not default, the panic of default or the probability of default still spreads within the network system, placing especially negative impact on the creditor banks or other correlated agents. Meanwhile, creditor banks also face the dilemma of repayment for their own debts and have to sell assets. However, it is difficult for banks to sell assets without realizing losses under market liquidity deteriorate; that is, their assets should be traded at a certain discount rate. erefore, in this paper, we let q be the asset flow rate reflecting the market liquidity, so that the banks would sell the assets with the book value ≤ 0 a bi (t))/q) to make up for the funding shortage. (2) e second round of shock is on the assets within the bank-firm credit network. In this round, shocks on the asset side of the balance sheet of the firm(s) or bank(s) propagate along the network; even such shocks are not large enough to trigger the default of these agents. Banks would make the evaluation of the shocked counterparties' ability of paying the loans. is is motivated by the fact that as one agent's equity decreases, so does its distance to default [55], and it will be less likely to repay its obligations in case of future distress. erefore, banks should have incentives to draw back credit from the correlated institutions or make liquidity hoarding to them in case of the default. If so, banks can make self-rescue when their assets are lower than a certain value, which could be the benchmark of the bank's compulsory leverage. Meanwhile, it is also presumed that the banks should make the decision of drawing back the credit from the firms randomly, whoever has credit-lending relationship with them. As far as banks are concerned, drawing back credit may avoid potential loss of assets, but it could place negative impact on the counterparty agents. It is assumed that being drawn back credit firms or banks should go bankrupt potentially and their liquidation of assets should meet the repayment of their loans to the banks. So, the relative changes of the shocked firms' equity and the survival ones within the network can be presented as eoretically, systemic risk contagion could last until all banks go bankrupt, and E i (t) ≤ 0 will be the exact bankrupt benchmark for each agent. As a matter of fact, except for credit risk and liquidity risk mentioned above, there are still lots of channels or factors resulting in risk spreading in financial networks. Specifically, the synchronization of these possible factors should accelerate or strengthen the risk contagion, which means that the impact of the outbreak of financial systemic risk in reality is likely to be worse than the analysis in this paper with certain constraints.

Model of Systemic Risk Contagion in
Financial Network e above reconstructed resulting bipartite network G � (B, F) consists of two sets of nodes: banks B and firms F. In network, the links connect banks not only with other banks (interbank liabilities), but also with the firms (liabilities of firms or credit from the banks). e weighted links between each agent in all can be defined as liability matrix where each entry L ij indicates the liability that node i has to node j, excluding the null entries between firms and the entries from firms to banks. erefore, the reconstructed multilayer financial network can be partitioned into the following two types: interbank network and bank-firm network.
Let A ji,c � a ji,c N * N(M) be the adjacency matrix of the financial network; the credit relationship between each agent can be represented by these entries, where a ji,c indicates the credit level of node i to node j in layer c and c indicates the network layer with different value. Specifically, it is assumed that there should be no circle credit relationship within each node's own, which means that a ji,c � 0 is a default setting.
For simplicity, in asset side of a bank i's balance sheet, the asset V i is subdivided into the following three parts: interbank assets V BI � N b�1 a bi , bank-firm assets V FI � M f�1 a fi , and other assets V OI , while in liability side of the bank i's balance sheet it is also divided into three parts: interbank liabilities L BI , other liabilities L OI , and equity E i , as is shown in Table 1. For interbank asset in each bank's balance sheet, there is a corresponding liability in another bank's balance sheet; that is, L ij � a ji , thus, interbank liabilities is L BI � N b�1 a ib and the rest liabilities can be cal- equity E i is used as a buffer of loss, and the relationship V BI + V FI + V OI � L BI + L OI + E i is constant in common sense. Similarly, the simplified firm balance sheet can also be defined as that in Table 2.

General DebtRank Model.
In order to measure the systemic risk and its dynamics of contagious process in financial network system, the DebtRank model, originally proposed by Battiston et al. [14] and Bardoscia et al. [48], is applied to be a special tool for measurement in this article. e original DebtRank model, transformed from PageRank algorithm [56], can be used to describe the dynamic process of stress propagation within interbank market successfully in several literatures. Briefly speaking, the DebtRank model could be adopted to quantify the extent of the financial distress that a particular node should face under the external shocks and the corresponding risk contagion within different channels. erefore, a transition matrix called interbank market's vulnerability matrix V ∈ B × B can be created in the model, in which B is the set of banks, and the vulnerability matrix can be defined as V ij � (A ij /e i ), where ∀i, j ∈ B and V ij ∈ [ 0, ∞ ).
e entry A ij represents the unsecured exposure of creditor bank i to debtor bank j, while e i > 0 denotes the capital buffer of bank i. us, when V ij ≥ 1, bank j leads bank i to default as well and the intermediate value of V ij ∈ ( 0, 1 ] leads bank i into distress but not into default. In the general DebtRank model, two indicators are defined to measure the risk status of each node at each time step. One indicator is H i (t) ∈ [0, 1], which represents the  Discrete Dynamics in Nature and Society 5 vulnerability of node i, and the other one is C i (t) � U, D, I { }, a discrete variable, representing the status of each shocked node i. By the way, the entries in the set of C i (t) represent the status of "Undistressed," "Distressed," and "Inactive," respectively.
e specific modal indicators in the general DebtRank model can be seen as follows: e above DebtRank model always converges for H i (t) within the value by H i (t) � 1, due to the operator min, and shows nondecreasing monotonically, as a result of the nonnegative value in vulnerability matrix V ij � A ij /e i .
Consequently, for a sufficiently large number of iteration steps T < ∞, the system converges to a stable status. e resulting DebtRank index DR is simply the additional stress that the interbank market suffers apart from the initial exogenous shock b(0) that can be defined as follows: where φ i represents the economic importance of agent i, which can be simplified and designed to be the share of liabilities that agent i has in the financial system. In this formula, DR(b(0)) is conditional on the initial shock b(0).

Risk Measurement Model Based on General DebtRank
Model. e basic assumption of the general DebtRank model is that the loss or distress should be linearly transmitted from the debtor to the creditor and the devaluation of the debtor's equity leads to the devaluation of the creditor's assets on the same scale. Based on this, systemic risk contagion may result in potential loss of creditor assets, and the change of creditor's asset is related to that of debtor's equity. erefore, in this paper, changes of each bank's asset within interbank network and bank-firm market are described in the following formulas, respectively. Specifically, the credit amount of banks' drawing from the firms is considered to be the changes of assets in bank-firm market.
Different from the original DebtRank model, where only the credit risk was considered, both liquidity risk and credit risk are taken into account in our optimized model, and the loss caused by the liquidity risk can be calculated as erefore, based on the risk contagion process through the balance sheet in both asset side and liability side, the dynamic changes of equity in financial network can be microscopically drawn as where . e last two terms show the assumption that the debt is unchanged according to the face value.
According to the above analysis, the ration of risk exposure to equity is defined as the risk leverage matrix [48]. Concerning the deterioration of liquidity risk, the risk exposure should be magnified by (1/q) times potentially, which is mostly significant in short-term loans with high liquidity. According to formulas (3) and (4), the recursion process of the risk contagion can be drawn as in the following formula, on the condition that the risk exposure should not change in the first step; that is, a ji,c (1) � a ji,c (0): us, the maximum eigenvalue |λ max | of the interbank leverage matrix can be calculated to characterize whether the risk of the external shock could be augmented in the risk where In original DebtRank algorithm, the node propagates the shock only once after receiving the shock. However, in this article, shocks could propagate until the agents default. Correspondingly, we optimize the DebtRank model so as to fit the new designed mechanism; therefore, the indicator C i (t) of the risk measurement can be adjusted into the one in the following formula: It is assumed that there be X firms and Y banks in initial distress state in this article. After a finite time step T, all banks in the multilayer financial network are in status U or status I, and eventually the dynamics could tend to be stable in a certain scale.

Simulation Parameters.
Based on the above theoretical analysis, a series of computational simulations on systemic risk contagion and its impact on the stability of the financial network are conducted. Although various simulations with different settings will run, we confine some baseline set of the following parameters: the number of firms and banks should be set as N f � 1000 and N b � 50, respectively. e bank size distribution is assumed to be followed by a truncated Pareto distribution with shape parameter α � 1.5, and the boundaries of the nodes' links are set as L � 5 and H � 200. According to the empirical results in existing literatures, let the average number of links be λ b � 100 and λ f � 10 for banks and firms, respectively. Meanwhile, we make the assumption that each agent should possess at least one connection to the other ones and that each bank should be involved in the interbank credit market. Meanwhile, the parameter θ (initially set up as 28%) of the amount of external assets will change systematically in the following experiment so as to explore how the different scenarios influence the vulnerability of the financial system. Additionally, the parameters of link-creating in interbank credit market should be set as d � 0.5, α 1 � 0.25, and α 2 � 1 that can be employed to generate the disassortative features found in the actual market. Other than the above setting, we presume that the initial proportion of equity of total assets of the bank should be φ 2 � 12% and the ratio of the credit drawn from the firm by the banks should be φ 3 � 8%.
Besides the above parameters, eight types of scenarios about the financial correlation condition within banking and firm sectors are assumed and conducted in simulation analysis, so that how the difference in risk contagion among each scenario can be compared and explored. Initially, all banks and firms are considered to be in the status of "Undistress." Concerning the eight types of scenarios, except from the network structure and the behavior of the agents, market liquidity and shock ratio are also the specific factors to be taken into account in the simulation. erefore, these eight types of scenarios are divided into three categories in this article that have different value definition for each focused factor. e first category consists of four scenarios (Scenario 1-Scenario 4), which are conducted in the simulation to compare how the nodes' different initial distress states influence the risk contagion. e second category consists of two scenarios (Scenario 5 and Scenario 6), which are conducted in the simulation to compare how the different external shock ratios on assets influence the risk contagion process. e third category consists of two scenarios (Scenario 7 and Scenario 8), which are conducted in the simulation to compare how the different liquidity flow rates influence the risk contagion process. By the way, all the value setup of the related parameters is shown in Table 3.
Overall, in the whole eight scenarios, ten or more time steps are allowed in the simulation, so that we can observe the changeable state of banks' equity E i and risk measure indicators H i (t) and C i (t) until the time steps stop or the financial network system comes to a relative stable state. erefore, the evolvement of systemic risk contagion process influenced by the factors of network structure, the nodes' initial risk status, the external shock ratio, and the asset flow rate in the network is comprehensively analyzed in the simulation.

Simulation
Results. Based on the above assumption and definition, the following figures derived from the simulation experiment can absolutely prove the above theoretical viewpoints and give some additional meaningful information that all banks in the multilayered network would be affected in different scales in risk contagion process. erefore, the simulation analysis is on how the systemic risk contagion is comprehensively conducted from the factors of the network structure, initial distress status of Discrete Dynamics in Nature and Society agents, the external shock ratio, the liquidity rate, and the layers of the network, by which the risk contagion process can be compared and analyzed.

Impact on Risk Contagion from Network Structure.
How the linkage between each agent directly impacting the systemic contagious process is extremely important to both the speed of contagion and the extent of agent's failure. erefore, in order to measure how the network structure (there are many indicators that can be applied to reflect the network structure, such as node density, degree distribution, diameter, and clustering coefficient; in this paper, two indicators of the average degree and the power law exponent of degree distribution are used) influences the risk contagion process, the related variables in simulation experiment could be assigned specific value. Specifically, the external shock ratio and the liquidity rate are defined as 0.01% and 95%, respectively. Meanwhile, the average degree of interbank layer network and bank-firm layer network is defined as Avd Interbank � 5.8 and Avd bank− firm � 4.5, and the power law exponent of the degree distribution is defined as c � 2.6 in Figure 3(a). On the contrary, the average degree of interbank layer network and bank-firm layer network is defined as Avd Interbank � 8.9 and Avd bank− firm � 14.9, and the power law exponent of the degree distribution is defined as c � 2.9 in Figure 3(b). On this assumption and definition, Figure 3 shows the losses of banks' equity over the time influenced by the different network structure of interbank credit network and bank-firm credit network, when five banks and five firms in the system are shocked to be in the status of "distress." From Figures 3(a) and 3(b), it can be seen that almost all of the banks in reconstructed financial multiplex network go bankrupt within the fourth and tenth time step, which indicates that the speed of contagion is very high and the influential scope is relatively wide in both types of network structure. In Figure 3(b), the banks' equity drops suddenly from the start of the simulation that is different and a little bit slow in Figure 3(a). Comparing the banks' relative loss in Figure 3(b) with that in Figure 3(a), we can see that, in the two scenarios, as the average degree of the network increases, the speed and strength of risk contagion are significant. erefore, for the more linkages within bank-firm credit network, distress or default from the bank-firm financial network could cause relatively more severe risk to the entire system, especially the postderived risk combined with the interbank market, which is extreme disaster and may be more intensive if the assumed credit relation is decomposed by actual complex financial correlation.

Impact on Risk Contagion from Nodes' Initial Distress.
In financial network, the initial status of the agents should lead the whole financial system to different endings no matter what kind of intervention is and how it will be carried out. Consequently, let the average degree of interbank network and the bank-firm network be Avd Interbank � 5.8 and Avd bank− firm � 14.9, respectively. Meanwhile, four types of scenarios with different quantity and category of shocked agents will be selected to explore how the nodes' initial status influences the risk contagion, whose corresponding simulation results can be seen in Figure 4. As shown in Figure 4, most banks in the reconstructed two-layer financial network almost have significant loss between the fourth and the fifth time steps, which shows that the speed of risk contagion is absolutely high and the influence is extremely wide no matter which scenario is tested. Figures 4(a) and 4(b) show the relative loss of the banks with external shocks randomly on one bank and one firm, respectively, while it can be seen that the losses of banks' equity come out in a little bit faster speed in Figure 4(b) with almost no obvious difference. Meanwhile, the significant different simulation results can be obtained between Figures 4(d) and 4(c), where the scenario with five banks and one firm being shocked simultaneously is assumed in Figure 4(c), and the scenario with one bank and five firms   Discrete Dynamics in Nature and Society being shocked simultaneously is shown in Figure 4(d). In Figure 4(c), the decrease of the banks' equity is slower than that in Figure 4(d), and part of the banks survived as the time goes at last. Besides, we can find that the speed of systemic risk contagion in Figure 4(d) is nearly the fastest among all the four scenarios, which implies that some firms in the financial system also play crucial roles in risk contagion. erefore, from the simulation on all four types of scenarios about the systemic risk contagion, the conclusion is that the risks resulting from the failed firms should place more negative effect on the financial systems than that only happening in interbank market. According to that, risks from the real economy in the measurement of systemic risk contagion in financial system should be also taken into account for the complex financial correlation in the economy especially the potential feedback mechanism in the contagious process, which may cause unimaginable disaster.

Impact on the Systemic Risk Contagion from External Shock Ratio.
A sudden incidence on economy may cause a huge loss even if the whole system collapses. us, how the accidental events affect the systemic risk contagion should be analyzed theoretically for the practical regulation. Based on the fifth and the sixth scenarios, average degree of the interbank network and the bank-firm network is also defined as Avd Interbank � 5.7 and Avd bank− firm � 14.9, respectively, and the other parameters follow the default definition. Figure 5 presents the changes of proportion of the banks' status during the systemic risk contagion process under different external shock ratios, which are defined as 0.1% in Figure 5(a) and 10% in Figure 5(b), respectively.
From Figure 5, we can see that almost most of the banks turn to the state of "I," which means that these banks go bankrupt, but the evolving process shows difference between the two scenarios. In both Figures 5(a) and 5(b), there will be no banks being totally healthy within the first three time steps, and meanwhile almost all of the banks come to the state of "D," no matter what the external shock ratio is. As time goes, the proportion of the banks with the state of "I" goes up quickly and the proportion of the banks with the state of "D" comes down accordingly. However, within the twenty time steps, the proportion of banks with the state of "D" is almost zero in Figure 5(b), while that proportion in Figure 5(a) comes to 30%. Besides, the changes of proportion of the banks' state are rapid in Figure 5(b), which conveys the implication that relatively high external shock ratio results in more harmful risk contagion on the condition of the simulation. erefore, some suggestion on financial risk management for the regulation authorities can be drawn that proper intervention should be implemented in the specific period of systemic risk contagion process. However, the intervention from the outside of the financial system should not be simply considered as rescuing some banks or firms unconditionally.

Impact on Systemic Risk Contagion from Liquidity
Rate. For liquidity rate could be an influencing factor in credit network-based risk contagion, how it works should be evaluated in the simulation more in detail. According to the seventh and the eighth scenarios presumed above and all other parameters being defined as the default setup, Figure 6 shows the losses of the proportion of banks' equity over time under different liquidity rates.
In Figures 6(a) and 6(b), most of the banks go bankrupt within the fourth and fifteenth time step, and the changes of banks' equity follow almost similar track, while the liquidity ratios are defined as 50% and 95%, respectively, for the comparative simulation analysis. However, there is also obvious difference between the two evolvement processes. Compared with the evolving tract in Figure 6(a), the time span of the banks' bankruptcy lasts a little bit longer, and the time span for banks' failure and the banks' equity decreases to a relative state which are presented in more concentration in Figure 6(b). Specifically, the survival of banks' equity in Figure 6(b) is also larger than that in Figure 6(a). From the simulation results, we can find that liquidity rate indeed places impact on the speed of risk contagion and the degree of asset loss in that process.
us, some advice may be provided for risk management that liquidity control should be a paramount means to prevent or manipulate risk contagion, implemented dynamically not only in the risk contagious process, but also in the specific polices or regulations for firm sectors and financial sectors in various industries.

Impact on Risk Contagion from Network's Different
Layers. Based on the above definition about the parameter setting, this section aims to find out how the layers of the reconstructed financial network, both interbank network layer and bank-firm network layer, affect the systemic risk contagion process. Just like the above sectors, let the average degree of the interbank network and the bank-firm network be Avd Interbank � 5.8 and Avd bank− firm � 14.9, respectively. Meanwhile, it is assumed that liquidity rate be defined as 95% and five banks and five firms be shocked by external shock. Further, we conduct the simulation analysis to measure the systemic risk by the variable H within contagious process, which means the relative loss of the agents or the system along the time steps shown in Figure 7.    Discrete Dynamics in Nature and Society of the reconstructed multiplex networks, the bank-firm layer is much more likely to cause banks' bankruptcy eventually. In Figure 7(a), only a few parts of banks will absolutely go bankrupt, and most of the banks just come to be the stable state of distressed or infectious ones as the time steps go. However, in Figure 7(b), when it almost comes to the third time step, all of the banks turn to be bankrupt suddenly. To some extent, the simulation results can also propose an obvious viewpoint that bank-firm credit market plays a significant role in the process of risk contagion, so that the regulatory authorities should take into account the complex correlations between the industrial sectors and the financial sectors systematically when establishing and implementing related polices, and they cannot ignore the impact on the financial system from real economy at all.

Conclusion
e theory of multilayer network provides a comprehensive visualization and solution to explore complex financial markets. In this article, a two-layer financial credit network is reconstructed between banking and firm sectors on the structure of each agent's balance sheet generally, consisting of interbank lending layer and bank-firm lending layer. Based on the reconstructed financial network, systemic risk is considered as a status of loss or distress of the financial system including credit risk and liquidity risk in this article, and its contagion mechanism has been systemically analyzed.
ereafter, systemic risk contagion model is constructed on the DebtRank model, originally proposed by Battiston et al. [14] and Bardoscia et al. [48]; and the simulation analysis is conducted to investigate how the key factors influence risk contagion process due to external and internal shocks.
e main conclusions are drawn as follows: (1) in reconstructed financial network, network structure is an important factor influencing the risk contagious process, no matter how different the average degree of each layer network is. For the scale-free characteristic of the designed network, nearly all of the banks should be bankrupt over time in the simulation without any external intervention on the basic definition of the risk contagion circumstance, but the changes of banks' equity show different evolving process. Overall, the higher the average degree of the network is, the more influential and faster the risk contagion is. (2) In the two-layered financial networks, risks that originated from the failed firms in bank-firm layer should place more negative effect on the financial systems than that only happening in interbank market, especially the firm sectors with highly financial linkage with the financial departments. Meanwhile, how the external shock ratio placing impact on the risk contagion can also be obtained that the higher the initial shock ratio, the lower the peak proportion value of the distressed banks. (3) As a major cause of systemic risk, different liquidity rate in market could lead to discrepancy of the risk contagion speed and the degree of asset loss during the process, which implies that liquidity control should be a paramount and effective means to risk manipulating and the methods should be developed and implemented for the industries collaboratively. (4) In the same assumed market condition, track of H values' changes over time presents obvious difference in interbank layer and bank-firm layer, respectively, where nearly all of the banks go bankrupt eventually in bank-firm layer and most of the banks are not in bad condition in interbank layer. erefore, these simulation results give clear evidence that a series of diverse and associated policies for risk management should be established and implemented in the market.
However, the reconstructed financial credit network within banking and firm sectors in this article mainly focuses on the credit relationship between banks and from banks to firms, while the financial correlation between firms and from firms to banks is not considered comprehensively that would result in underestimating the systemic risk. Besides, there are more financial correlations within different types of agents, directly or indirectly, such as common exposure, derivatives holding, and payment transaction. What is more, concerning the systemic risk model and its simulation analysis, attention has been paid to how banks go in risk contagious process, while it is also important to care about the firms' status and the feedback mechanism between both sectors.
erefore, all these gaps should be fulfilled in future research.

Data Availability
All data used in this article are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.