Abstract

The stability of the financial system plays a crucial role in the sustainable economic development. Hence, to identify systemically important banks and firms, we take lending relationships with different loan terms and common asset relationships with different investment cycles into consideration to present a multilayer DebtRank model of the bank-firm system. In the light of simulation research, we can obtain the following results. First, the bank-firm system constructed displays a significant core-periphery structure, which exists in the actual financial system. Then, only very few banks and firms show systemically important characteristics, where “important” subjects hold very high net assets and profits, while "fragile" subjects possess negative net assets and serious losses. Furthermore, the bank-firm multilayer DebtRank model presents a great stability to a certain extent. Overall, the multilayer DebtRank model constructed in this paper has certain theoretical reference value for the supervisory authorities to extract the internal characteristics of systemically important banks and firms and identify them effectively.

1. Introduction

The sustainable and healthy development of the economy cannot do without the stability of financial systems [13]. The global financial markets have presented multiple complex forms because of the expansion and deepening of financial activities in the process of financial globalization. Since the 2008 financial crisis, regulators have begun to pay more attention to systemically important financial institutions. Once the crisis occurs in systemically important financial institutions, it will cause great damage to the whole financial system. Based on this, in order to ensure the stable and healthy development of the financial system, it is necessary to accurately identify the systemically important financial institutions.

In the aspect of identifying systemically important financial institutions, scholars have given diversified measurement methods from different dimensions such as models and indexes. The CoVaR model is the current mainstream research method [4, 5]. In addition, The Shapley index is also widely used to measure the importance of various banking systems and clarify the important factors affecting the systemic importance of the financial system [6]. Similarly, the SRISK index is used to measure the capital shortage degree of financial institutions, which is an important means to reflect systemically important financial institutions [7]. Some scholars also study the systemic importance of financial systems from the perspective of the marginal expected shortfall (MES) [8] and the systemic expected shortfall (SES) [9].

However, the modern financial system is more and more complex, and the relationship between banks and firms is more and more close, which makes the use of financial network methods to study the systemic importance of financial systems more and more urgent [10]. In the field of the financial network research, scholars have achieved fruitful results. For instance, Martinez-Jaramillo et al. [11] measure systemically important banks from the perspective of the strength centrality, the degree centrality, the betweenness centrality, the closeness centrality, the entropy eigenvector centrality, and the PageRank centrality in financial networks. Nicosia et al. [12] propose the concept of the feedback centrality in complex networks based on the PageRank centrality. Referring to the feedback centrality, Rincón and Villalobos [13] use the authority centrality and the hub centrality to measure the importance of financial systems.

For the study of the systemic importance from the perspective of financial networks, in addition to the above network centrality indicators, the DebtRank centrality has also received attention. Considering that the calculation process of the feedback centrality has multiple feedback effects among network nodes, the DebtRank sets different impact states so that the nodes are not affected by the secondary impact, which can well identify the systemically important financial institutions. Battiston et al. [14] propose the DebtRank model to identify the global systemically important financial institutions from 2008 to 2010 and find that 22 institutions are at the core of the financial system during the financial crisis. The DebtRank model can detect the bank-firm default effect by exerting additional shocks on the bank-firm network and obtain more reliable results. For example, small regional banks will also have a greater impact on the economy by connecting with local important firms [15].

The related research of the DebtRank model mainly focuses on the single risk exposure of financial systems, which is difficult to truly describe the complexity and diversity of the business associations between real financial market subjects. In view of this, considering the lending relationships with different loan terms and the common asset relationships with different investment cycles, the multilayer DebtRank model of the bank-firm system is constructed. Compared with the existing research, the contribution of this paper is to construct a multilayer DebtRank model to identify systemically important banks and firms from the complexity and the diversity of bank-firm business associations. In addition, we analyze the internal characteristics of systemically important banks and firms and verify the robustness of the model. This paper is beneficial to deeply mine the internal relationships between the complex association structure and the systemic importance of financial systems, which has certain theoretical reference value for the supervisory authorities to extract the internal characteristics of systemically important banks and firms and identify them effectively.

The remainder of this paper is organized as follows. Section 2 describes the multilayer DebtRank model of the bank-firm system. In Section 3, we discuss the main simulation results, and the conclusion is in Section 4.

2. The Model

In the artificial financial system, we consider two types of subjects: banks and firms. The entire model is divided into three parts, namely, bank-firm balance sheet construction, bank-firm behavior evolution, and bank-firm systemic importance identification.

2.1. Bank-Firm Balance Sheet Construction
2.1.1. Bank Balance Sheet Construction

We assume that an individual bank’s assets include interbank short-term loans , interbank long-term loans , firm short-term loans , firm long-term loans , short-term investments , long-term investments , and liquid assets , and that a bank’s liabilities are composed of interbank short-term borrowings , interbank long-term borrowings , short-term deposits , long-term deposits , and net assets . At time t, there are banks. Drawing on the experiences of the research conducted by Georg [16] and Li et al. [17], the bank balance sheet is presented as follows (Tables 1).

The indicators mentioned above can be measured as follows:where is the set of borrowing banks of the bank ; is the set of borrowing firms of the bank ; is the set of creditor banks of the bank ; and denote interbank short-term loans and long-term loans from the bank of the bank , respectively; and represent bank-firm short-term loans and long-term loans from the bank of the firm , respectively; and denote interbank short-term loans and long-term loans from the bank of the bank , respectively; and represent short-term investments and long-term investments of the bank , respectively; and denote short-term deposits and long-term deposits from depositors of the bank and satisfy the formulas and , respectively, where is a random variable and and are used to measure the short-term and long-term volatility of deposits; and are the maturity of interbank short-term loans and long-term loans, respectively; and are the maturity of bank-firm short-term loans and long-term loans, respectively; and are the maturity of bank short-term investments and long-term investments, respectively; and are the maturity of short-term deposits and long-term deposits, respectively; and , , and represent the ratio of interbank short-term loans, the ratio of bank-firm short-term loans, and the ratio of bank short-term investments, respectively.

2.1.2. Firm Balance Sheet Construction

We assume that an individual firm’s assets include production costs , short-term investments , long-term investments , and liquid assets and that a firm’s liabilities are composed of sale revenues , bank short-term loans , bank long-term loans , and net assets . At time t, there are banks. Drawing on the experiences of the research conducted by Ma et al. [18], the firm balance sheet is presented as follows (Table2).

The indicators mentioned above can be measured as follows:where denotes the labor wage for producing the products; represents the actual labor for producing the products; and denote short-term investments and long-term investments of the firm , respectively; means that the product price per unit output is randomly taken from the interval ; represents the actual output for producing the products and satisfies the formula , where , , and denote the comprehensive technical level, capital output elasticity coefficient, and labor output elasticity coefficient, respectively; denotes the set of creditor banks of the firm ; and represent bank-firm short-term loans and long-term loans from the bank of the firm , respectively; and are the maturity of firm short-term investments and long-term investments, respectively; and are the maturity of bank-firm short-term loans and long-term loans, respectively; and is the ratio of firm short-term investments.

2.2. Bank-Firm Behavior Evolution
2.2.1. Bank-Firm Lending

For simplicity, represents the intraperiod cash of the bank or the firm. At time t, banks and firms with negative liquidity can be classified as potential debt banks and debt firms . In addition, banks with sufficient liquidity can be classified as potential creditor banks . Banks and firms with negative liquidity conduct loan applications to banks with sufficient liquidity.

For the debt bank or the debt firm , each potential borrower randomly selects a certain proportion of potential creditor banks and observes the credit lending rates that they can provide. The short-term lending rate and the long-term lending rate set by the potential creditor banks are shown in equations (3) and (4), respectively [19]:where represents the risk-free interest rate; and denote the sensitivity of the short-term lending rate and the long-term lending rate to the bank-firm risk, respectively; represents the current ratio of creditor bank ; and denotes the debt leverage ratio of the debt bank or the debt firm .

The potential debt banks and debt firms borrow funds from the potential creditor banks based on the optimal partner selection mechanism. If they cannot obtain the sufficient liquidity from the first potential creditor bank, they contact other banks for the remaining funds until their total demand for liquidity is satisfied or all loanable funds are exhausted. As a potential creditor bank, the amount that can be used for fund lending is , where denotes the deposit reserve ratio. If the creditor bank has enough liquidity, it satisfies all its potential borrowers. Otherwise, it allocates its surplus liquidity in sequence according to the rank of the potential borrowers’ net assets from high to low until all its loanable funds are exhausted.

2.2.2. Bank-Firm Investments

At time t, if the bank and the firm have residual liquidity, they will carry out investment activities. The investment amount satisfies the following formulas, respectively:where denotes the net profit of the bank and ; represents the net profit of the firm and .

Suppose there are types of risky assets; the return rate of the kth risky asset is , which is abbreviated as . The return rate of risky assets is affected by common risk factors and their idiosyncratic risk factors. Therefore, the return rate of the kth risky asset satisfies , where the return rate of the common risk factor and the return rate of the idiosyncratic risk factor obey the normal distributions and , respectively. represents the degree and the direction of the kth risky asset affected by the common risk factor and is randomly taken from the interval .

Banks and firms randomly choose a certain proportion of risk assets for equal investments. For simplicity, the selected risky assets are recorded as , and the corresponding return rates are abbreviated as . Then, the short-term returns and the long-term returns of investments by banks and firms are measured as follows:

2.3. Bank-Firm Systemic Importance Identification

In this paper, the systemic importance of the bank-firm system refers to the impact of losses suffered by banks and firms on the whole bank-firm system through lending correlations and common asset correlations. Considering the lending relationships with different loan terms and the common asset relationships with different investment cycles, we use the multilayer pressure diffusion to evaluate the multilayer DebtRank of banks and firms and then identify the systematically important banks and firms.

Drawing on the experience of research conducted by Bardoscia et al. [20], the pressure level of banks and the pressure level of firms are given at time t, respectively. The pressure levels satisfy and , respectively. The default probability of the bank under the pressure level is and satisfies , where represents the default amplification factor. Similarly, the default probability of the firm under the pressure level is and satisfies .

Considering that if there is a loop in the network, the pressure will circulate between nodes. In order to avoid the nodes participating in the pressure diffusion process repeatedly, drawing on the experiences of the research conducted by Battiston et al. [14] and Poledna et al. [21], the nodes in the bank-firm system are given three pressure states, pressured , unpressured , and inactive at time t. For simplicity, and denote the pressure states of the bank z and the firm j at time t, respectively. The whole process of the pressure diffusion is as follows. When the node in the unpressured state is impacted by the external shock, it will be transformed into the pressured state , which will infect other nodes in the next round of the pressure diffusion and transform itself into the inactive state . When all nodes are in the states of unpressured or inactive , the whole process of the pressure diffusion ends. The stress states of banks and firms are calculated as follows:where and . denotes the set of all banks and denotes the set of all firms.

At the initial time , we assume that is the set of impacted nodes of banks and firms. The pressure levels and pressure states of banks and firms are set as follows:

In order to further describe the bank-firm pressure level in different network layers, we set the number of layers of the complex network as and , where , and , respectively, represent the short-term common asset network, the short-term lending network, the long-term common asset network, and the long-term lending network between banks and firms. In addition, we set , and as the short-term common asset scale matrix, the short-term lending scale matrix, the long-term common asset size matrix, and the long-term lending scale matrix between banks and firms, respectively.

At time t, the pressure level of the bank z in the different network layers is expressed as the sum of the pressure level in the previous state and the pressure level caused by all firm nodes in the present state. The calculation formula is as follows:where , and represent the short-term common asset sensitivity coefficient, the short-term lending sensitivity coefficient, the long-term common asset sensitivity coefficient, and the long-term lending sensitivity coefficient of the bank z to the firm j, respectively. They satisfy the calculation formulas as follows:

In addition, denotes the difference of the default probability received by the bank z from the pressured nodes and satisfies .

Similarly, at time t, the pressure level of the firm j in the different network layers is expressed as the sum of the pressure level in the previous state and the pressure level caused by all bank nodes in the present state. The calculation formula is as follows:where , and represent the short-term common asset sensitivity coefficient, the short-term lending sensitivity coefficient, the long-term common asset sensitivity coefficient, and the long-term lending sensitivity coefficient of the firm j to the bank z, respectively. They satisfy the calculation formulas as follows:

In addition, denotes the difference of the default probability received by the firm j from the pressured nodes and satisfies .

Assume that the pressure diffusion stops at . The DebtRank of the impacted node in the different network layers is expressed as the difference between the total pressure level of the four network layers at time and the pressure level of the impacted node at the initial time . The corresponding calculation formula is as follows:where represent the short-term investment weight, the short-term lending weight, the long-term investment weight, and the long-term lending weight held by the bank node z in all banks at time t, respectively. represent the short-term investment weight, the short-term lending weight, the long-term investment weight, and the long-term lending weight held by the firm node j in all firms at time t, respectively. In addition, and can be calculated as follows:

Based on the above analysis, the multilayer DebtRank of the impacted node i in the bank-firm multilayer network can be expressed as follows:

3. Simulations

3.1. Model Parameters

Drawing on the experiences of research led by Li et al. [17] and Ma et al. [18], the simulation parameter values are set in this section. The simulation is performed with 50 banks and 100 firms over a time span of 200. At the initial time of the simulation, 50 banks are divided into two categories. The first category includes 10 banks with , , , and . The second category has 40 banks with , , , and . In addition, 100 firms are divided into two categories. The first category includes 20 firms with and . The second category has 80 firms with and . Table 3 summarizes the other benchmark parameters of the model.

3.2. Network Structures of Bank-Firm System

The ban-firm system constructed in this paper contains four network layers, including the short-term common asset network, the short-term lending network, the long-term common asset network, and the long-term lending network. The lending network describes the lending relationship between banks and firms, and the common asset network represents the common holding relationship of external assets between banks and firms. Figure 1 illustrates the complex network structure diagram of the bank-firm system at t = 200. As is shown, the network is represented by nodes and edges, with the left nodes representing banks and the right nodes representing firms. The connected sides indicate the common asset relationships and the lending relationships between banks and firms. As can be seen from Figure 1, the bank-firm network is closely connected, which shows that the financial business between banks and firms is closely related. The characteristics of bank-firm network structures need to be further explored.

The node degree is a simple and important concept to describe the characteristics of bank-firm network structures. The larger the node degree, the more active the node in the bank-firm network. Drawing on the experience of research conducted by Boccaletti et al. [22], we use the multilayer network node degree to analyze the bank-firm network with four network layers. Figure 2 illustrates the node degree distributions of the bank-firm multilayer networks at different times (t = 80, t = 120, t = 160, and t = 200), where the horizontal and vertical coordinates represent the bank-firm number and the multilayer network node degree, respectively. Numbers 1–50 are banks and numbers 51–150 are firms.

It can be seen from Figure 2 that the operation results of the bank-firm system at different times are similar. In other words, most banks possess larger node degrees and are at the core of the bank-firm network while a few banks and all firms have smaller node degrees and are at the periphery of the bank-firm network. This indicates that the bank-firm system constructed in this paper has always shown a significant core-periphery structure at different times, which is also verified in the empirical study of Bargigli et al. [23]. This also means that most banks in the whole bank-firm system are highly active while a few banks and all firms have relatively low activities.

In addition, Figure 3 illustrates the complementary cumulative distribution function (CCDF) of total assets of the bank-firm system at different times (t = 80, t = 120, t = 160 and t = 200). Through the fitting analysis, we can find that the distribution of total assets of the bank-firm system in Figures 3(a) and 3(d) can be described as lognormal distributions with Pareto tails, where the tail exponents are 6.6040, 8.4225, 4.4721, and 6.7831, respectively. Similar conclusions have appeared in the empirical study of Chinese financial market by Xu et al. [24]. This suggests that most banks and firms have large-scale assets while only a few banks and all firms possess a small amount of wealth.

3.3. Systemic Importance of Bank-Firm System

The above research shows that the bank-firm system constructed in this paper shows the core-periphery structure, and most banks are active at the core of the network, while a few banks and firms are relatively less active at the periphery of the network. It is also very important to identify the systemically important banks and firms. Based on this, we use the multilayer DebtRank to analyze the problem. In order to explore whether there are differences in the systemic importance between banks and firms at different times during the operation of the ban-firm system, Figure 4 calculates the bank-firm multilayer DebtRank at different times (t = 80, t = 120, t = 160, and t = 200), where the horizontal and vertical coordinates represent the bank-firm number and the multilayer DebtRank, respectively. Numbers 1–50 are banks and numbers 51–150 are firms.

Figure 4 presents the operation results of the bank-firm system at different times are similar. In other words, most nodes of banks and firms have a multilayer DebtRank between 0 and 0.3, showing weak systemic importance. However, only very few nodes of banks and firms (the bank of No. 7 and the firm of No. 137) have a multilayer DebtRank that exceeds 0.5, showing strong systemic importance. This indicates that only very small numbers of banks and firms show the systemic importance in the bank-firm system composed of 50 banks and 100 firms. Therefore, we should strengthen the monitoring of the financial situations and the operation conditions of these systemically important banks and firms and identify their “importance” and “vulnerability” significantly. This also reflects the problem of “too connected to fail” in the real bank-firm system to prevent the financial systemic risk effectively.

The aforementioned research results show that only very small numbers of banks and firms show the systemic importance in the bank-firm system while the internal characteristics of systemically important banks and firms need to be further explored. Based on this, considering that the financial situations and the operation conditions of banks and firms may have significant impacts on the systemic importance of the bank-firm system [21], we select two indicators of net assets and profits for in-depth analysis.

In the actual operation process of the bank-firm system, due to the differences in asset scales and operation conditions between banks and firms, there may be large differences in the net assets between banks and firms. Whether this will affect the systemic importance of banks and firms remains to be further studied. Based on this, Figure 5 presents the distribution of the net assets and their multilayer DebtRank of the bank-firm system at t = 200. The abscissa represents the net assets of banks and firms, and the ordinate represents the multilayer DebtRank of banks and firms. In Figure 5, all banks are represented by blue diamonds, and all firms are represented by red squares.

It can be seen from Figure 5 that the net assets of most banks and firms are in the interval (0, 200), and the net assets of most banks are higher than those of the firms. The multilayer DebtRank of most banks and firms are between 0 and 0.3, while only the bank of No. 7 and the firm of No. 137 show a strong systemic importance. The bank of No. 7 owns the highest net asset and the highest multilayer DebtRank, which shows an obvious “importance.” However, the firm of No. 137 has the negative net asset and its multilayer DebtRank remains a high level, showing a significant “vulnerability.” It can be concluded that for the systemically important banks and firms, “important” subjects possess high net assets, while “vulnerable” subjects hold negative net assets.

Banks and firms in good operating conditions may gradually evolve into systemically important banks and firms through the accumulation of profits while poorly managed banks and firms may gradually become systemically vulnerable banks and firms due to frequent losses. Whether this will affect the systemic importance of banks and firms remains to be further studied. Based on this, Figure 6 presents the distribution of the profits and their multilayer DebtRank of the bank-firm system at t = 200. The abscissa represents the profits of banks and firms, and the ordinate represents the multilayer DebtRank of banks and firms. In Figure 6, all banks are represented by blue circles, and all firms are represented by red triangles.

It can be seen from Figure 6 that the profits of most banks and firms are in the interval (−2, 5), and the multilayer DebtRank of most banks and firms are between 0 and 0.3. We can see that only the bank of No. 7 and the firm of No. 137 show a strong systemic importance. The bank of No. 7 owns the highest profit and the highest multilayer DebtRank, which shows an obvious “importance.” However, the firm of No. 137 has the negative profit and its multilayer DebtRank remains a high level, showing a significant “vulnerability.” This indicates that, for the systemically important banks and firms, “important” subjects hold high profits while “vulnerable” subjects have serious losses.

3.4. Robustness Analysis

In order to ensure the robustness of the simulation results, it is necessary to further verify the robustness of the bank-firm multilayer DebtRank. Considering that the parameter represents the bank-firm default amplification factor, it depicts the correlation between the default probability of banks and firms and the pressure level, which plays an important role in the pressure diffusion of the bank-firm system. Particularly, when , the default probability of banks and firms is linearly related to the pressure level. When , there is a nonlinear correlation between the default probability of banks and firms and the pressure level. In order to further verify the robustness of the bank-firm evolution process, Figure 7 presents the bank-firm multilayer DebtRank with different values (, , and ) at different times (t = 80, t = 120, t = 160, and t = 200).

It can be seen from Figure 7(a) that the multilayer DebtRank of banks and firms are similar under different values (, , and ) at t = 80. For the systemically important bank of No. 7, the multilayer DebtRank is 0.593, 0.590, and 0.585 under , , and , respectively. For the systemically vulnerability firm of No. 137, the multilayer DebtRank is 0.604, 0.602, and 0.596 in the three cases, respectively. For all nodes of banks and firms, with the gradual increase of , the multilayer DebtRank decreases slightly. Under different values, the distribution of the bank-firm multilayer DebtRank is approximately the same. This shows that the increase of the bank-firm default amplification factor can only reduce the multilayer DebtRank very limitedly, but it can maintain a good robustness to the overall systemic importance of banks and firms. In addition, it can be seen from Figures 7(a)7(d) that the overall systemic importance of the bank-firm system keeps a very good stability at different times (t = 80, t = 120, t = 160, and t = 200). In other words, the multilayer DebtRank constructed in this paper shows a good robustness to a certain extent.

4. Conclusions

The multilayer DebtRank is beneficial to deeply mine the internal relationships between the complex association structure and the systemic importance of financial systems. Considering the lending relationships with different loan terms and the common asset relationships with different investment cycles, the multilayer DebtRank model of the bank-firm system is constructed. The entire model is divided into three parts, namely, bank-firm balance sheet construction, bank-firm behavior evolution, and bank-firm systemic importance identification. This paper focuses on the identification of systemically important banks and firms, the internal characteristics of systemically important banks and firms, and the robustness test of the model.

The simulation results are as follows. Firstly, the bank-firm system constructed shows a significant core-periphery structure at different times, which means that most banks in the whole bank-firm system are highly active, while a few banks and all firms have relatively low activities. Secondly, we find that the distributions of total assets of the bank-firm system are lognormal distributions with Pareto tails, which suggests that most of the banks and firms have large-scale assets, while only a few banks and all firms possess a small amount of wealth. Thirdly, we discover that only very small numbers of banks and firms display the systemic importance in the bank-firm system, where “important” subjects hold very high net assets and profits, while “fragile” subjects possess negative net assets and serious losses. This suggests that we should pay more attention to the systemically important banks and firms that are “too connected to fail.” Finally, we verify that the overall systemic importance of the bank-firm system keeps a very good stability at different times.

It is urgent for the government departments to better understand the number and scale distribution of the current Chinese banks and firms, balance the development mode of related business, reasonably optimize the loan mechanism, and establish a multidirectional risk monitoring system. This study reveals the microbasis of complex business associations between banks and firms, which has certain theoretical reference value for the regulatory authorities to extract the internal characteristics of systemically important banks and firms and prevent financial systemic risks.

Conflicts of Interest

The authors declare no conflicts of interest.

Authors’ Contributions

Qianting Ma and Kun Yang contributed equally to this work. They are cofirst authors.

Acknowledgments

The authors are grateful to the financial support from the Postgraduate Research and Practice Innovation Program of Jiangsu Province (KYCX20_0165) and the Scientific Research Foundation of the Graduate School of Southeast University (YBPY1971).