Dynamics Evolution of Credit Risk Contagion in the CRT Market

Abstract

This work introduces a nonlinear dynamics model of credit risk contagion in the credit risk transfer (CRT) market, which contains time delay, the contagion rate of credit risk, and nonlinear resistance. The model depicts the dynamics behavior characteristics of evolution of credit risk contagion through numerical simulation. Meanwhile, numerical simulations show that, in the CRT market, the contagion rate of credit risk and the nonlinear resistance among CRT activities participants have some significant effects on the dynamics behaviors of evolution of credit risk contagion. Specifically, on the one hand, we find that the status curve of credit risk contagion that causes some significant changes with the increase in the contagion rate of credit risk, moreover, emerges a series of Hopf bifurcation and chaotic phenomena in the process of credit risk contagion. On the other hand, Hopf bifurcation and chaotic phenomena appear in advance with the increase in the nonlinear resistance coefficient and time-delay. In addition, there are a series of periodic windows in the chaotic interval inside, including Hopf bifurcation, inverse bifurcation, and chaos.

1. Introduction

Over the past few years, with the significant development of nonlinear science, economists have gradually started to use nonlinear theory to study the complex phenomena of social economic system [1–7]. Some far-sighted economists began to apply the nonlinear science research results into economics, which produced the nonlinear economics and the chaos economics. The latest studies of nonlinear theory show that whether interpersonal network, computer network, ecological system, economic system, or disease spread, computer virus spread, forest fire spread, risk spread, complex nonlinear dynamics phenomena, and so forth, can be observed in these social phenomena [8, 9]. The aforementioned phenomena present complex dynamical behavior, involving Hopf bifurcation, inverse bifurcation, chaos, and fractals. Among these behavior types, chaos and bifurcation are complex phenomena that exist in the nonlinear financial system and are important issues in economic and financial dynamics research [10]. Credit risk transfer (CRT) market is a third-party market that connects with the credit markets, the securities market, and the insurance market, in which credit risk contagion has some complex nonlinear characteristics obviously.

At present, participants of the CRT market covered mainly universal banks, commercial banks, securities dealers, insurance companies, investment funds, and parts of nonfinancial institutions. Among them exist close and complicated network relations directly or indirectly, and that constituted a nonlinear giant system. Because the interactions between individuals that have complex nonlinear dynamic properties. Moreover, credit risk contagion is dependent on CRT behaviors of participants of the CRT market and market information dissemination of the relationship network. With the rapid development of the CRT market, the quantity of participants, and the depth and breadth of CRT trading all rapidly increase. This will lead to the increase in the complexity of the CRT market and make the distribution of information and risk of terminal undertaker of credit risk change more complicate. Meanwhile, the rapid development structured products will also increase the complexity. These will make the financial institutions extremely easily cause the superposition or clustering of credit risk in credit risk transfer and cause credit risk contagion. However, credit risk contagion also has complex nonlinearity. It will increase the difficulty of the prediction and control of credit risk in CRT market and bring great challenges to credit risk management departments.

Generally, some participants do not fully understand the potential risk of CRT market or lack of corresponding risk management ability into the CRT market, which will lead to some new risks in the process of CRT behaviors. Moreover, the systemic risk can increase in the CRT market. In the imperfect competition market, CRT behaviors not only did not spread risk, but also added to the system risk and increased the likelihood of the credit risk contagion [11]. The existing literature also showed that the rapid development of the CRT market increased the possibility of credit risk contagion across departments and trade. For example, credit risk transfer in creating contagion between banking and insurance systems and caused contagion, and the spread in systemic risk made everybody worse off. At the same time, credit risk transfer induced insurance companies to hold the same assets as banks [12]. Banks' motive of extensive using CDS (Credit Default Swap) is that improve the diversification of their credit risk. However, this might reduce banks’ stability. The main reasons behind these negative impacts are firstly, that banks are induced to increase their investment in an illiquid, risky credit portfolio and secondly, that these CDS create a possible channel of credit risk contagion [13].

The theory and practice have recognized the serious consequences of the credit default contagion by the US subprime mortgage crisis in 2008. Moreover, a number of studies are also aware of credit risk contagion in the CRT process [11–14]. At present, the study of credit risk contagion mainly focus on the interbank market and credit market. However, the existing credit risk model have not yet discussed and involved nonlinear dynamic problems of the risk contagion process. However, nonlinear dynamic behaviors are obvious in credit risk contagion due to the complex network relationships, the continuous innovation of CRT tools, and the asymmetric information in CRT market. Moreover, network relations of CRT market exist time delay and nonlinear resistance. Therefore, we try to put the nonlinear system theory into the study of the credit risk contagion in CRT market and construct the nonlinear dynamic model of credit risk contagion in CRT market. Then, we conduct numerical simulation to analyze the dynamic behaviors characteristics of evolution of credit risk contagion in CRT market.

The remainder of this paper is organized as follows. In Section 2, the model of credit risk contagion in CRT market and dynamics behavior characteristics of evolution of credit risk contagion are discussed through numerical simulation. In Section 3, we discuss the bifurcation and chaotic behaviors of credit risk contagion. Finally, we conclude the paper in Section 4.

2. Dynamics Evolution of Credit Risk Contagion Based on the Vector Field

With the development of network theory, a number of studies have taken into account the spread and response characters of events in a long-distance connection of network. Newman and Watts [15], Moukarzel [16] have given the dynamic model of constant speed transmission of the events in the network. However, they have not taken into account time-delay and various nonlinear factors. Yang [17] took into account the nonlinear factor and time-delay of events in the long connection and constructed the reasonable dynamic model.

2.1. The Contagion Model of Credit Risk in the CRT Market

We are enlightened by the works [17–19] and propose the dynamic model of credit risk contagion in the CRT market. On the one hand, we assume that the complex network connections among CRT activities participants are Newman-Watts length scale connections and long-distance connections. On the other hand, we take into account the time-delay and nonlinear resistance of long-distances connection between CRT activities participants. In fact, the model is also a nonlinear time-delay differential equation. Therefore, the dynamic model of credit risk contagion is described by the following time-delay differential equation: 𝑑𝑁(𝑡)𝑑𝑡=𝜆𝑘₁−𝑁(𝑡)+𝜆𝑘₂𝑁(𝑡−𝜏)−𝜇𝜉𝜆𝑘₂𝑁(𝑡−𝜏)²𝑁(𝑡)𝑡≥0,=𝑐−𝜏≤𝑡≤0, (1) where 𝑁(𝑡) denotes the number of CRT activities participants that are infected by credit risk in the CRT market, 𝜉 refers to Newman-Watts length scale, 𝑘₁ is the number of instances that the connection distance from the participant infected by credit risk is a Newman-Watts length scale, 𝑘₂ is the number of instances that the connection distance from the participant infected by credit risk is a long-distance connection, 𝜆 is the effective contagion rate of credit risk in the CRT market, 𝜇 is the nonlinear resistance coefficient of the relationship network comprising CRT market participants, 𝑐∈ℝ⁺ is a real parameter, and 𝜏 is the time-delay of credit risk contagion in the long-distance connection. Therefore, the mechanism of the time-delay and the nonlinear resistance of credit risk contagion in Newman-Watts length scale connection and long-distance connection can be described by the time-delay differential equation (1).

According to the general definition, we can derive the balance position and stable point of credit risk contagion when the left side of equation (1) is equal to zero. In fact, this kind of nonlinear dynamics system can be denoted by equation (1), where balance positions may become unstable, periodic solution and the system vibration may emerge, and the phenomenon of Hopf bifurcation and chaos may occur, along with the change in various parameters [20]. Torelli [21], Liu and Spijker [22] have given a numerical Euler method for the solution of delay differential equation as equation (1). We still use the method in this paper. Now, let the stepsize ℎ is such that ℎ=𝜏/𝑚 and 𝜃∈[0,1], where 𝜏 is a time-delay, and 𝑚 is a positive integer. Therefore, according to the one-point collocation rule for delay differential equation (1), we can get 𝑁_𝑛+1=𝑁_𝑛+𝜆𝑘₁ℎ−ℎ(1−𝜃)𝑁_𝑛+𝜃𝑁_𝑛+1+𝜆𝑘₂ℎ(1−𝜃)𝑁_𝑛−𝑚+𝜃𝑁_{𝑛−𝑚+1}−𝜇𝜉𝜆²𝑘²₂ℎ(1−𝜃)𝑁_𝑛−𝑚+𝜃𝑁_{𝑛−𝑚+1}², (2) where 𝑁_𝑛 denotes the approximate value of 𝑁(𝑡) at the point 𝑡_𝑛. Let =(𝑚−𝛿)ℎ+ℎ/2(0≤𝛿<1), then Ω_ℎ={𝑡_𝑛=𝑛ℎ,𝑛∈𝑍}. Thus, we can get 𝑡_𝑛+𝜃ℎ∈[𝑡_𝑛,𝑡_𝑛+1] and 𝑡_𝑛+𝜃ℎ−𝜏∈[𝑡_𝑛−𝑚,𝑡_{𝑛−𝑚+1}]. We apply the 𝜃-collocation method to define the approximate value of 𝑁(𝑡) at the point 𝑡_𝑛+𝜃ℎ and 𝑡_𝑛+𝜃ℎ−𝜏 as follows: 𝑁𝑡_𝑛+𝜃ℎ≈𝜃𝑁𝑡_𝑛+𝑁𝑡_𝑛+1,𝑁𝑡_𝑛+𝜃ℎ−𝜏≈𝜃𝑁𝑡_𝑛−𝑚+𝑁𝑡_{𝑛−𝑚+1}. (3)

We apply the midpoint collocation method (one-point collocation with 𝜃=1/2) to equation (1), and can get 𝑁_𝑛+1=𝑁_𝑛+𝜆𝑘₁𝑁ℎ−ℎ^𝑛+𝑁_𝑛+12+𝜆𝑘₂ℎ𝑁_𝑛−𝑚+𝑁_{𝑛−𝑚+1}2−𝜇𝜉𝜆²𝑘²₂ℎ𝑁_𝑛−𝑚+𝑁_{𝑛−𝑚+1}2². (4) Namely, 𝑁𝑡_𝑛+1=𝑁𝑡_𝑛+𝜆𝑘₁𝑁𝑡ℎ−ℎ^𝑛+𝑁𝑡_𝑛+12+𝜆𝑘₂ℎ𝑁𝑡_𝑛−𝑚+𝑁𝑡_{𝑛−𝑚+1}2−𝜇𝜉𝜆²𝑘²₂ℎ𝑁𝑡_𝑛−𝑚+𝑁𝑡_{𝑛−𝑚+1}2². (5) Put equation (3) into equation (5), we can get 𝑁𝑡_𝑛+1≈𝑁𝑡_𝑛+𝜆𝑘₁ℎ−ℎ𝑁𝑡_𝑛+ℎ2+𝜆𝑘₂ℎ𝑁𝑡_𝑛+ℎ2−𝜏−𝜇𝜉𝜆²𝑘²₂ℎ𝑁𝑡_𝑛+ℎ2−𝜏². (6) We put ℎ=𝜏/𝑚 into equation (6), we can get 𝑁𝑡_𝑛+1≈𝑁𝑡_𝑛+𝜆𝑘₁𝜏𝑚−𝜏𝑚𝑁𝑡_𝑛+𝜏+2𝑚𝜆𝑘₂𝜏𝑚𝑁𝑡_𝑛+(1−2𝑚)𝜏−2𝑚𝜇𝜉𝜆²𝑘²₂𝜏𝑚𝑁𝑡_𝑛+(1−2𝑚)𝜏2𝑚². (7)

To understand the effect of nonlinear factors on credit risk contagion further, we have to use equation (7) to conduct numerical simulations under the given parameters 𝜇, 𝜉, 𝑘₁, 𝑘₂, 𝜆, and 𝜏 and the initial condition 𝑁(𝑡)=𝑐  (−𝜏≤𝑡≤0).

2.2. Simulation Analysis of the Dynamics Behavior of Evolution of Credit Risk Contagion in the CRT Market

We try to describe the dynamics behavior characteristics of evolution of credit risk contagion and its influencing factors by the nonlinear time-delayed differential equation in this paper. According to the solving process of equation (1), we know that parameters 𝜇, 𝜉, 𝑘₁, 𝑘₂, and 𝜆 and the initial condition 𝑁(𝑡)=𝑐(−𝜏≤𝑡≤0) will affect the stability of the solution of time-delayed differential equations and the trajectory of the process of credit risk contagion. In order to describe the dynamic behaviors and its influencing factors of the process of credit risk contagion in CRT market, we take parameters 𝜆 and 𝜇 as the bifurcation parameter. Then, we conduct numerical simulations to the dynamics system (1) and analyze the dynamics behavior of credit risk contagion in CRT market. Let 𝜏=1, ℎ=0.01, 𝑚=100, 𝛿=0.5, 𝜉=3, 𝑘₁=10, 𝑘₂=25, and the initial condition 𝑁(𝑡)=2(𝑡∈(−𝜏,0)). Figure 1 depicts the effect of the effective contagion rate 𝜆 of credit risk on the trajectory curve of credit risk contagion in the CRT market. We find that the status of credit risk contagion changes gradually from “hyperbolic attenuation” (a piece of the hyperbolic) to “logarithm Gauss attenuation,” and the influence strength and range of credit risk contagion emerge nonlinear velocity increasing with the increase in the effective contagion rate 𝜆 of credit risk in CRT market. However, the influence strength and range attenuate rapidly after a period of time and emerge the fat-tail characteristic. This shows that the effect of the default behaviors of CRT activities participants on other participants weakened gradually after a period of time and the default intensity and default state depend on the company oneself and macroeconomic factors. Figure 2 shows that oscillation amplitude and frequency increase gradually with the increase in the effective contagion rate 𝜆 of credit risk in CRT market. However, the oscillation will weaken after a period of time. Figure 3 shows that the stable state of credit risk contagion will trend to unstable and emerge periodic solution and Hopf bifurcation with the increase in the effective contagion rate 𝜆 of credit risk in CRT market. Namely, the contagion amplitude and range of credit risk will emerge periodic oscillation with the increase in the effective rate of credit risk contagion in CRT market. Moreover, the limit cycle radius of the attractive domain increases gradually, and the shape of the limit cycle becomes increasingly irregular, such that the bifurcation and chaos phenomena occur with the increase in the effective contagion rate 𝜆 of credit risk contagion. In Figure 4, we find that the process of credit risk contagion present different “logarithm Gauss attenuation” feature under the influence of the nonlinear resistance of the relationship network comprising CRT activities participants. In Figure 5, we find that the oscillation of the process of credit risk contagion is not affected with the increase in the nonlinear resistance coefficient 𝜇. However, the number of CRT activities participants 𝑁(𝑡) gradually reduces with the increase in the nonlinear resistance coefficient 𝜇. In Figure 6, we find that the effect of the nonlinear resistance coefficient 𝜇 on the attract factor of balance position of credit risk contagion, and the number of CRT activities participants 𝑁(𝑡) is very sensitive. Namely, the attractive factor of credit risk contagion and the number of CRT activities participants 𝑁(𝑡) will decrease rapidly with the increase in nonlinear resistance coefficient 𝜇.

Figure 1 

The trajectory curve of credit risk contagion where 𝜇=0.03.

Figure 2 

The step response of the process of credit risk contagion where 𝜇=0.03.

Figure 3 

The phase diagram of the process of credit risk contagion where 𝜇=0.03.

Figure 4 

The state trajectory curve of the process of credit risk contagion where 𝜆=0.1.

Figure 5 

The step response of the process of credit risk contagion where 𝜆=0.1.

Figure 6 

The phase diagram of the process of credit risk contagion where 𝜆=0.1.

3. Bifurcation and Chaotic Analysis of Credit Risk Contagion Based on Logistic Mapping

3.1. The Model Analysis of Credit Risk Contagion Based on Logistic Mapping

The model (1) of credit risk contagion used the form of vector field to discuss credit risk contagion in credit risk transfer. However, the previous figures are not intuitive and are difficult to interpret. Thus, analyzing the properties of the dynamic system of credit risk contagion, such as the difference of the trajectory curve of period doubling, may be challenging. However, given the intuition, legibility, and geometrical features of the logistic mapping, we often discretize the nonlinear problem of the continuous vector field to the logistic mapping by using a numerical approximation method to analyze the periodic bifurcation and chaos of nonlinear dynamics system. A number of studies use the Euler [23–25] to analyze bifurcation, periodic solution, and chaotic phenomena of nonlinear time-delayed system. We also adopt the Euler method and take step length ℎ. Therefore, equation (1) can be transformed into the form following form: 𝑁(𝑡−𝜏+ℎ)−𝑁(𝑡−𝜏)=ℎ𝜆𝑘₁−𝑁(𝑡)+𝜆𝑘₂𝑁(𝑡−𝜏)−𝜇𝜉𝜆𝑘₂𝑁(𝑡−𝜏)². (8)

Let ℎ=𝜏, 𝑁(𝑡)=𝑁_𝑛+1, and 𝑁(𝑡−𝜏)=𝑁_𝑛. Thus, equation (8) can be transformed into the form following form: 𝑁_𝑛+1=𝜆𝑘₁𝜏+1+𝜏1+𝜆𝑘₂𝜏𝑁1+𝜏^𝑛−𝜇𝜉𝜏𝜆²𝑘²₂𝑁1+𝜏^𝑛2. (9)

Therefore, there exists the logistic mapping 𝑓 as follow: 𝑓∶𝑁_𝑛↦𝑁_𝑛+1. (10)

According to the definition of the fixed point of the logistic mapping, we know that the fixed point of the logistic mapping 𝑓 should meet 𝑁_𝑛+1=𝑁_𝑛=𝑁^∗. Therefore, we can get the analytic equation of the fixed point of the logistic mapping as follow: 𝜇𝜉𝜏𝜆²𝑘²₂𝑁1+𝜏^∗2−𝜆𝑘₂𝜏−𝜏𝑁1+𝜏^∗−𝜆𝑘₁𝜏1+𝜏=0. (11)

Therefore, we can get the fixed point of the logistic mapping 𝑓 by equation (11) as follow: 𝑁^∗₁=𝜆𝑘₂𝜏−𝜏+𝜆𝑘₂𝜏−𝜏²+4𝜇𝜉𝑘₁𝜏²𝜆³𝑘²₂2𝜇𝜉𝜏𝜆²𝑘²₂,𝑁^∗₂=𝜆𝑘₂𝜏−𝜏−𝜆𝑘₂𝜏−𝜏²+4𝜇𝜉𝑘₁𝜏²𝜆³𝑘²₂2𝜇𝜉𝜏𝜆²𝑘²₂. (12)

Obviously, 𝑁^∗₂<0 is unrealistic. Therefore, the fixed point 𝑁^∗₁ is sole fixed point of the logistic mapping 𝑓. According to the definition of the logistic mapping and the Lyapunov movement stability, we know that the movement stability of the fixed point depends on the characteristic root of the derived operator of the logistic mapping, which is Floquet multiplier [26, 27]. Therefore, the Floquet multiplier will determine the stability of the fixed point 𝑁^∗₁. Namely, 𝐷𝑔_𝑁∗=1−𝜆𝑘₂𝜏−𝜏²+4𝜇𝜉𝑘₁𝜏²𝜆³𝑘²₂.1+𝜏 (13)

According to the nonlinear system theory [27, 28], if |𝐷𝑔|_𝑁∗>1, then the fixed point 𝑁^∗ will become unstable; if |𝐷𝑔|_𝑁∗<1, then the fixed point 𝑁^∗ is asymptotically stable; if |𝐷𝑔|_𝑁∗=1, then the fixed point 𝑁^∗ is criticality stable. So, for the fixed point 𝑁^∗ of the mapping 𝑓, the fixed point 𝑁^∗ is asymptotically stable when 𝜇<(4(1+𝜏)²−(𝜆𝑘₁𝜏−𝜏)²)/4𝜉𝑘₁𝑘²₂𝜏²𝜆³, is criticality stability when 𝜇=(4(1+𝜏)²−(𝜆𝑘₁𝜏−𝜏)²)/4𝜉𝑘₁𝑘²₂𝜏²𝜆³, or is unstable when 𝜇>(4(1+𝜏)²−(𝜆𝑘₁𝜏−𝜏)²)/4𝜉𝑘₁𝑘²₂𝜏²𝜆³.

According to the nonlinear dynamic related theory [26–28], if there exists a series of period-doubling bifurcation phenomena, then a series of period-doubling bifurcation leads to chaos. In recent years, much works used topological horseshoes embedded method to study chaos rigorously [28–34]. By this method, one can not only prove the existence of chaos, but also reveal the mechanism of chaotic phenomena by showing the structure of chaotic attractors [31–34]. Beyond that, some works used the Lyapunov exponents [35] and set oriented numerical methods [36, 37] to prove the existence of chaos. Li and Yorke [38] gave a definition of chaos that the existence of a point of period 3 implies the existence of chaos. Therefore, according to this definition, we use numerical simulation to discuss the fixed point and its stability, bifurcation, and chaos of the mapping from the intuitive.

3.2. Numerical Simulation Analysis

Let 𝜉=3, 𝜏=1, 𝑘₁=10, 𝑘₂=25, and the initial condition 𝑁(𝑡)=2(𝑡∈(−𝜏,0)). We use equation (8) to conduct numerical simulations. The Figure 3 reflects the Hopf bifurcation process and its variation characteristics of credit risk contagion with parameter 𝜆 and 𝜇. Figures 7(a) and 7(b) reflect the Hopf bifurcation and chaos characteristics of credit risk contagion with the increase in the effective contagion rate 𝜆 of credit risk. Figure 7(a) reflects that the process of credit risk contagion exists the only stable constant state when parameter 𝜆 is kept at a proper level. Moreover, the process of credit risk contagion emerge different types of period bifurcation and periodic oscillation with the increase in the effective contagion rate 𝜆 of credit risk in CRT market. According to the definition of Li-Yorke [29], the process of credit risk contagion can occur chaos phenomenon when the effective contagion rate 𝜆 reaches to a proper value. Figure 7(b) reflects a series of similar characteristics with Figure 7(a). However, we also find that the Hopf bifurcation and chaotic phenomena of credit risk contagion emerge in advance with the increase in the nonlinear resistance coefficient 𝜇. In Figure 7(c), we find that the Hopf bifurcation and chaotic phenomena of credit risk contagion emerge in advance with the increase in time-delay 𝜏. In Figure 7(d), we find that the process of credit risk contagion exists the only stable constant state when parameter 𝜆 is kept at a proper level. Moreover, the process of credit risk contagion emerges different types of period bifurcation and periodic oscillation with the increase in the nonlinear resistance coefficient 𝜇 among CRT activities participants. According to the definition of Li-Yorke [29], the process of credit risk contagion can occur chaos phenomenon when the nonlinear resistance coefficient 𝜇 reaches to a proper value.

(a)

(b)

(c)

(d)

Figure 7 

(a) The bifurcation diagram of the process of credit risk contagion with 𝜆 when 𝜇=0.01; (b) the bifurcation diagram of the process of credit risk contagion with 𝜆 when 𝜇=0.015; (c) the bifurcation diagram of the process of credit risk contagion with 𝜆 when 𝜇=0.01, 𝜏=1.5; (d) the bifurcation diagram of the process of credit risk contagion with 𝜇 when 𝜆=0.1.

According to numerical simulation and comparative analysis, we find that the process of credit risk contagion can emerge three states, including the stable constant state, Hopf bifurcation, and chaos with the increase in parameter 𝜆 and 𝜇. However, these cannot more directly depict the nonlinear dynamic behavior characteristics after occurring chaotic phenomena. Therefore, we further discuss the effect of these parameters on the chaotic state and the period window of the process of credit risk contagion. In Figures 8(a) and 8(b), we find that Hopf bifurcation, pour bifurcation, and chaos mixed emerge in chaotic interval internal period window. Moreover, Hopf bifurcation, pour bifurcation, and chaos phenomena emerge in advance in chaos interval inside with the increase in nonlinear resistance coefficient 𝜇. Figure 8(c) shows that chaos states are significant in the process of credit risk contagion with the increase in nonlinear resistance coefficient 𝜇. However, Hopf bifurcation and pour bifurcation features become relatively obscure comparing to the Figures 8(a) and 8(b).

(a)

(b)

(c)

Figure 8 

(a) The bifurcation diagram in the chaos area when 𝜇=0.01; (b) the bifurcation diagram in the chaos area when 𝜇=0.015; (c) the bifurcation diagram in the chaos area when 𝜆=0.1.

4. Conclusion

In this paper, we constructed a nonlinear dynamic model of credit risk contagion based on literatures [17–19]. Moreover, the dynamical properties of the nonlinear dynamics system of credit risk contagion were investigated. We found that the effective rate of credit risk contagion and nonlinear resistance between CRT market participants have significant effect on dynamics behavior of credit risk contagion. Moreover, we found a series of complex Hopf bifurcation, inverse bifurcation, and chaos phenomena in the nonlinear dynamics system of credit risk contagion through a numerical simulation. At the same time, there are a series of period window in chaos interval inside, and that emerge intertwined state including Hopf bifurcation, pour bifurcation, and chaos. The study of dynamics behavior of evolution of credit risk contagion can help us to understand the effect of the interaction between the internal nonlinear factors and external disturbance of credit risk contagion, which has important theoretical and practical value.

There is still much work that is worth further research. For example, in the real world, a variety of noises usually influence the process of credit risk contagion and its dynamics behaviors, such as Gaussian noise, random noises, and so forth. For the kind of credit risk contagion with both time-delay and noises, we leave it for the future work.

Acknowledgments

The authors wish to express their gratitude to the referees for their invaluable comments. This work was supported by the National Natural Science Foundation of China Grant (nos. 71071034, 71173103, and 71201023), the Humanities and Social Science Youth Foundation of the Ministry of Education of China (no. 12YJC630101), the Funding of Jiangsu Innovation Program for Graduate Education (no. CXZZ12-0131), and the Scientific Research Foundation of Graduate School of Southeast University (no. YBJJ1238).