Next Article in Journal
Fractional Sequential Coupled Systems of Hilfer and Caputo Integro-Differential Equations with Non-Separated Boundary Conditions
Next Article in Special Issue
Cohen–Grossberg Neural Network Delay Models with Fractional Derivatives with Respect to Another Function—Theoretical Bounds of the Solutions
Previous Article in Journal / Special Issue
Topological Degree via a Degree of Nondensifiability and Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stability Analysis of a Credit Risk Contagion Model with Distributed Delay

1
Department of Mathematics, University of Cape Coast, Cape Coast P.O. Box 046, Ghana
2
Department of Management, Polytechnic University of Marche, 60121 Ancona, Italy
3
Department of Statistics, University of Cape Coast, Cape Coast P.O. Box 046, Ghana
*
Author to whom correspondence should be addressed.
Axioms 2024, 13(7), 483; https://doi.org/10.3390/axioms13070483
Submission received: 21 June 2024 / Revised: 4 July 2024 / Accepted: 15 July 2024 / Published: 18 July 2024

Abstract

:
This research investigates the stability and occurrence of Hopf bifurcation in a credit risk contagion model, which includes distributed delay, using the chain trick method. The model is a generalized version of those previously examined. The model is an expanded version of those previously studied. Comparative analysis showed that unlike earlier models, which only used the nonlinear resistance coefficient to determine the rate of credit risk infection, the credit risk contagion rate is also affected by the weight given to past behaviors of credit risk participants. Therefore, it is recommended to model the transmission of credit risk contagion using dispersed delays.

1. Introduction

The traditional perspective on credit risk involves borrowers failing to make timely payments, resulting in financial losses for the creditors [1]. The industrial and real economy sectors or firms with limited internal funds depend on the financial sector as a crucial source of funding, so their risk and return are significantly influenced by the unpredictability of the financial sector, especially in terms of their profitability and stability. According to Chen et al. [2], as the system’s stability decreases, the infectious scale of credit risk and the wavy frequency of credit risk contagion increase. This leads to chaotic phenomena in the credit risk contagion system, and the increase becomes gradual as the time delay increases. Any credit risk contagion system includes inherent structural instability, meaning that any disruptions within the system can cause a qualitative or topological change in its dynamic properties [2]. Hence, it is crucial to develop a mathematical model with time delays to capture the dynamics of credit risk contagion in order to understand the impact of the internal complexity of the credit risk contagion system and to offer both theoretical and practical guidance for credit risk management [3].
The significance of time delay is widely recognized in mathematical models related to finance, economics, and population dynamics [4,5]. This is because differential equations alone cannot fully capture the complexities of economic and financial phenomena [3]. Time delays commonly occur in real financial markets due to finite information transmission speed, noise, and the transport of matter and energy [6]. Additionally, mathematical models incorporating time delays have been suggested in the relevant literature to account for system behavior at previous time points, such as those observed in economic and biological systems [4,5,7,8,9]. The time delay has had notable impacts on the system, causing a switch in stability from a stable state to an unstable one, as well as leading to fluctuations and Hopf bifurcation [10,11].
Chen et al. [12] modeled the mechanism of time delay and nonlinear resistance in credit risk contagion within the Newman–Watts length scale connection and long-distance connection. They presented its dynamics using the following time-delay differential equation:
N · ( t ) = λ k 1 N ( t ) + λ k 2 N ( t τ ) μ ξ [ λ k 2 N ( t τ ) ] 2 , t 0
N ( t ) = c , τ t 0
where N ( t ) denotes the number of CRT activity participants who are affected by credit risk in the CRT market at time t. The parameters ξ , k 1 , k 2 , λ , μ , c , and τ are all positive, and, respectively, represent the Newman–Watts length scale, the number of instances where the connection distance from a participant infected by credit risk matches the Newman–Watts length scale, the number of instances where the connection distance from an infected participant constitutes a long-distance connection, the effective contagion rate of credit risk in the credit risk transfer (CRT) market, the nonlinear resistance coefficient of the network of relationships among CRT market participants, a real parameter, and the time delay in credit risk contagion within the long-distance connection.
Chen et al. [2] investigated the following stochastic time-delayed system of credit risk contagion, which is influenced by correlated Gaussian white noise:
N · ( t ) = λ 1 N ( t ) + λ 2 N ( t τ ) μ ξ [ λ 2 N ( t τ ) ] 2 + N ( t ) ψ ( t ) + η ( t ) ,
where N ( t ) refers to the number of credit activity participants affected by credit risk in the financial market, while the terms ψ ( t ) and η ( t ) are stochastic and represent noise in the system. The parameters λ 1 and λ 2 , with λ 1 > λ 2 , respectively, represent the effective spread of credit risk among various participants in credit activities, reflecting both direct and indirect business relationships in the financial market. τ represents the time delay for the spread of credit risk among participants through indirect connections. ξ is the Nelength scale and μ is the nonlinear resistance factor within the network of relationships among credit activity participants, indicating the diversity among participants, encompassing their offerings and psychological tendencies. For recent studies on credit risk contagion models, see [13,14,15,16,17,18,19,20,21].
Bianca and Guerrini [22] extended the work of Chen et al. [2] in the form of the following equation:
N · ( t ) = λ 1 N ( t ) + λ 2 N ( t τ ) μ ξ λ 2 2 N ( t τ ) 2 ,
where N ( t ) represents the density of individuals with credit activity in the financial sector who are at risk of credit exposure, λ 1 is the credit risk contagion rate associated with direct business relations, λ 2 is the credit risk contagion rate associated with indirect business relations, μ is alluded to as the Nelength scale, ξ symbolizes the nonlinear resistance factor in the network connecting financial market participants involved in credit activity, and τ 0 is the time delay of credit risk contagion. Bianca and Guerrini [22] used normal form theory and the center manifold argument to establish the presence of the Hopf bifurcation, as opposed to the numerical method employed by Chen et al. in [2]. They also identified the stability, direction, and period of the bifurcating periodic solutions with specific formulas and demonstrated that the examination of the deterministic aspect of model (3) is not limited to the scenario where λ 1 > λ 2 .
In the current literature, time delays are categorized as either fixed or continuously distributed (referred to as distributed delay). Fixed delays apply to contexts where specific time gaps are stipulated for involved parties. In contrast, distributed delays are suitable for scenarios where delays among actors vary. A key challenge lies in the uncertain nature of these time intervals. Distributed delays represent a weighted average incorporating all historical data from time zero to the present, making them more accurate in describing time-delayed systems. As found by Elaiw and Algha in [23], the concept of a distributed delay provides a more comprehensive and realistic way to model processes. Compared to models with a fixed time delay, it is evident that the latter oversimplifies the actual system process [24]. This assertion is bolstered in a study by Hu et al. [25] when a distributed delay was proposed to model time-delayed inherent biology systems. It was discovered that this model was sufficiently general to encompass a variety of characteristics exhibited by the system. Therefore, we generalize Equation (3) by substituting continuously distributed time delays for the discrete delays, as shown below, to assess if the distribution of time delays induces some changes in the rate of credit risk contagion compared to fixed time delay in the literature, which depends only on the nonlinear resistance coefficient. Thus, from the dynamics of the distributed delay model, the time delay of the number of individuals infected by credit risk varies among them, which is an extension of the effect of the discrete delay parameter in the previous models, which were the same for individuals infected by credit risk.
This paper is structured as follows. Section 2 deals with the formulation of the model, while basic properties of the model are derived in Section 3. In Section 4, the discussion revolves around the stability analysis, stability switch, and Hopf bifurcation analysis of the model. In Section 5, detailed numerical solutions and a stability analysis are conducted, comparing our model with the findings of Chen et al. [2]. Finally, Section 6 concludes.

2. Model Formulation

We propose the following credit risk contagion model with distributed delays:
N ˙ ( t ) = λ 1 N ( t ) + λ 2 t N ( r ) f ( t r ) d r μ ξ λ 2 2 t N ( r ) g ( t r ) d r 2 ,
where the delay kernels f ( r ) and g ( r ) are non-negative bounded functions defined on [ 0 , ) , satisfying the normalized conditions:
0 f ( r ) d r = 1 ,
0 g ( r ) d r = 1 .
When
f ( r ) = g ( r ) = δ ( r ) ,
where δ denotes a Dirac delta function, one recovers the discrete delay case (3). Another case of interest is when the kernels are gamma distributions (see Cushing [26]), namely
f ( r ) = g ( r ) = a m r m 1 e a r ( m 1 ) ! ,
with m being a positive integer that defines the form of the weighting function whose average delay is T = m / a . Notice that as T 0 , the distribution functions f , g approach the Dirac distribution. In the literature, researchers focus on the special cases where the delay kernel functions f and g may take the so-called weak kernel form
f ( r ) = g ( r ) = a e a r ( a > 0 ) ,
and strong kernel form
f ( r ) = g ( r ) = a 2 r e a r ( a > 0 ) .
The term weak kernel suggests that the significance of past events diminishes exponentially as one looks further into the past, while the term strong kernel indicates that a specific time in the past holds more importance than any other. In the next sections, we shall consider several special cases given by a combination of the special forms of the kernels.

3. Existence and Uniqueness of Solution

We use the contraction mapping theorem (see, e.g., [27,28,29]) to determine the existence and uniqueness of the solution of Equation (4). Set T N = N and rewrite Equation (4) as
λ 1 N = N · λ 2 t N g ( t r ) d r + μ ξ λ 2 2 t N g ( t r ) d r 2 .
Since g ( t r ) is continuous, we have that g r H . Given any N 1 and N 2 ,
T N 1 = 1 λ 1 N · 1 λ 2 t N 1 g ( t r ) d r + μ ξ λ 2 2 t N 1 g ( t r ) d r 2
and
T N 2 = 1 λ 1 N · 2 λ 2 t N 2 g ( t r ) d r + μ ξ λ 2 2 t N 2 g ( t r ) d r 2 ,
implying
T N 1 T N 2 = 1 λ 1 λ 2 t ( N 1 N 2 ) g ( t r ) d r + μ ξ λ 2 2 t ( N 1 2 N 2 2 ) g 2 ( t r ) d r = 1 λ 1 λ 2 t ( N 1 N 2 ) g ( t r ) d r + μ ξ λ 2 2 t ( N 1 N 2 ) ( N 1 + N 2 ) g 2 ( t r ) d r = 1 λ 1 ( N 1 N 2 ) t [ λ 2 g ( t r ) + μ ξ λ 2 2 ( N 1 + N 2 ) g 2 ( t r ) ] d r ,
consequently, we obtain
T N 1 T N 2 = 1 λ 1 | N 1 N 2 | t λ 2 g ( t r ) + μ ξ λ 2 2 ( N 1 + N 2 ) g 2 ( t r ) d r 1 λ 1 t λ 2 g ( t r ) + μ ξ λ 2 2 ( N 1 + N 2 ) g 2 ( t r ) d r β λ 1 λ 2 g ( t r ) + μ ξ λ 2 2 ( N 1 + N 2 ) g 2 ( t r ) | N 1 N 2 | η | N 1 N 2 | ,
where
η = β λ 1 λ 2 g ( t r ) + μ ξ λ 2 2 ( N 1 + N 2 ) g 2 ( t r ) .
For contraction,
β λ 1 λ 2 g ( t r ) + μ ξ λ 2 2 ( N 1 + N 2 ) g 2 ( t r ) 1 .
Since (5) holds, according to the contraction mapping theorem, it can be concluded that there is a single fixed point that serves as the only solution to the original integro-differential Equation (4). As a result, in cases where the parameters are chosen such that the kernel is small, ensuring contraction, the existence and uniqueness of the solution to Equation (4) are guaranteed.

4. Stability and Hopf Bifurcation Analysis

Since steady states are unaffected by the delay, Equation (4) admits the following unique positive equilibrium:
N = λ 1 + λ 2 μ ξ λ 2 2 .
In the following, we examine the stability and Hopf bifurcation of our model in certain specific scenarios. Our method utilizes the well-established linear chain trick (refer to [30]), which enables a reduction in a system with distributed delays to a system of ordinary differential equations.

4.1. Case f , g are Weak Kernels

Let
f ( r ) = g ( r ) = a e a r , a > 0 .
Equation (4) writes as
N ˙ ( t ) = λ 1 N ( t ) + λ 2 t N ( r ) a e a ( t r ) d r μ ξ λ 2 2 t N ( r ) a e a ( t r ) d r 2 .
Introducing the variable
u ( t ) = t N ( r ) a e a ( t r ) d r ,
by the linear chain technique, Equation (7) gets converted to the equivalent system below:
N · ( t ) = λ 1 N ( t ) + λ 2 u ( t ) μ ξ λ 2 2 u ( t ) 2 , u ˙ ( t ) = a N ( t ) u ( t ) .
Linearizing (8) at its equilibrium point ( N , u ) , where u = N , with N as in (6), we obtain the Jacobian matrix:
J = λ 1 λ 2 2 μ ξ λ 2 2 N a a .
Let ρ denote a characteristic root of (9). The related characteristic equation is of second order in ρ ,
ρ 2 + a λ 1 ρ + a 2 N μ ξ λ 2 2 λ 1 λ 2 = 0 ,
resulting in
ρ 2 + a λ 1 ρ + a λ 1 + λ 2 = 0 .
The determinant of J equals the product of the eigenvalues and is given by
ρ 1 ρ 2 = a λ 1 ,
while the trace of J is equal to the sum of the eigenvalues
ρ 1 + ρ 2 = a λ 1 + λ 2 > 0 .
The model’s local asymptotic stability is assured when both roots have negative real parts, leading to the following outcome.
Theorem 1. 
For a > 0 , system (8) is locally asymptotically stable if a > λ 1 , unstable if a < λ 1 .
We now use the Hopf bifurcation theorem to demonstrate that a limit cycle could start at a = λ 1 . This theorem says that we can prove the existence of a cyclic solution if the characteristic equation has a pair of purely imaginary roots and the real parts of these roots switch signs with a bifurcation parameter. Suppose ρ = i ω , ω > 0 , is a root of (10). Substituting it into (10), we have
ω 2 + a λ 1 ( i ω ) + a λ 1 λ 2 = 0 .
Separating the real and imaginary parts, we obtain
ω 2 + a λ 1 + λ 2 = 0 ,
a λ 1 ω = 0 ,
which yield
a = λ 1 = a ,
ω = λ 1 λ 1 + λ 2 = ω .
Accordingly, we conclude that Equation (10) has a pair of purely imaginary roots ρ = ± i ω when a = a .
Next, we select a as a bifurcation parameter and consider the eigenvalue as a function of a, ρ = ρ ( a ) . By implicitly differentiating Equation (10) with respect to a, we obtain
d ρ d a = ρ + λ 1 + λ 2 2 ρ + a λ 1 ,
then,
R e d ρ d a ρ = i ω = R e i ω + λ 1 + λ 2 2 i ω = 1 2 .
so that
s i g n R e d ρ d a ρ = i ω < 0 .
As the value of a decreases, the roots of the characteristic equation will transition from the right side to the left side of the imaginary axis. Based on the prior examination, we can draw the following conclusions.
Theorem 2. 
The equilibrium point of (7) bifurcates to a limit cycle through a Hopf bifurcation when a = λ 1 .

4.2. Case f , g are Strong Kernels

Let
f ( r ) = g ( r ) = a 2 r e a r , a > 0 .
Equation (4) becomes
N ˙ ( t ) = λ 1 N ( t ) + λ 2 t N ( r ) a 2 ( t r ) e a t r d r μ ξ λ 2 2 t N ( r ) a 2 ( t r ) e a t r d r 2 .
To apply the linear chain trick on the integro-differential Equation (11), we define
u ( t ) = t N ( r ) a 2 ( t s ) e a t r d s , v ( t ) = t N ( r ) a e a t r d s ,
so that (11) transforms into
N · ( t ) = λ 1 N ( t ) + λ 2 u ( t ) μ ξ λ 2 2 u ( t ) 2 , u ˙ ( t ) = a v ( t ) u ( t ) , v ˙ ( t ) = a N ( t ) v ( t ) .
The equilibrium point of (12) takes the form ( N , u , v ) , with u = v = N , N defined in (6). To determine its stability, the characteristic equation of the linearized system is taken into consideration, leading to the following polynomial equation in ρ :
ρ 3 + 2 a λ 1 ρ 2 + a a 2 λ 1 ρ + a 2 λ 1 + λ 2 = 0 .
The Routh–Hurwitz criterion implies that the real parts of the solutions of (13) are negative if and only if the following conditions hold:
2 a λ 1 > 0 ,
a 2 λ 1 + λ 2 > 0 ,
2 a λ 1 a a 2 λ 1 a 2 λ 1 + λ 2 > 0 .
These can be reduced to
a > λ 1 2 = a 0
and
2 a 2 6 λ 1 + λ 2 a + 2 λ 1 2 > 0 .
The discriminant of the left-hand side of (14) is given by
Δ = 10 λ 1 + λ 2 2 λ 1 + λ 2 > 0 .
Thus, the inequality (14) is solved by
a < 6 λ 1 + λ 2 10 λ 1 + λ 2 2 λ 1 + λ 2 4 = a 1 ,
a > 6 λ 1 + λ 2 + 10 λ 1 + λ 2 2 λ 1 + λ 2 4 = a 2
Noticing that 0 < a 1 < a 0 < a 2 , we arrive at the following result.
Proposition 1. 
Let a 2 be defined as in (15). The equilibrium point of (12) is locally asymptotically stable if a > a 2 , unstable if a < a 2 .
We will now go back to (13) and demonstrate via the Hopf bifurcation theorem that a limit cycle might begin at a = a 2 . At the value a = a 2 , where one has
2 a 2 λ 1 a 2 a 2 2 λ 1 = a 2 2 λ 1 + λ 2 ,
the characteristic Equation (13) can be rewritten as
ρ 3 + 2 a 2 λ 1 ρ 2 + a 2 a 2 2 λ 1 ρ + 2 a 2 λ 1 a 2 a 2 2 λ 1 = 0 ,
which factors as
ρ + 2 a 2 λ 1 ρ 2 + a 2 a 2 2 λ 1 = 0 .
Thus, one of (16) is real and negative, i.e., ρ = 2 a 2 λ 1 , and the other two are purely imaginary:
ρ = ± i a 2 a 2 2 λ 1 = ± i ω ,
where ω > 0 . By choosing a as the bifurcation parameter, taking the derivative of (13) with respect to a, we derive
d ρ d a = 2 λ 2 + 2 a λ 1 λ + 2 a λ 1 + λ 2 3 λ 2 + 2 2 a λ 1 λ + a a 2 λ 1 .
A direct calculation implies
R e d ρ d a λ = i ω = R e 2 ω 2 + 2 a λ 1 + λ 2 + 2 a λ 1 i ω 3 ω 2 + a a 2 λ 1 + 2 2 a λ 1 i ω = λ 1 2 + 3 a 2 2 6 λ 1 a 2 λ 2 a 2 2 a 2 λ 1 2 + a 2 2 a 2 2 λ 1 = 10 λ 1 + λ 2 2 λ 1 + λ 2 + 10 λ 1 + λ 2 2 λ 1 + λ 2 6 λ 1 + λ 2 8 2 a 2 λ 1 2 + a 2 2 a 2 2 λ 1 ,
namely
s i g n R e d ρ d a ρ = i ω < 0 .
As a result, the imaginary axis is crossed by the roots at a = a 2 from left to right as a decreases. In summary, the outcome is as follows.
Theorem 3. 
Let a 2 be as in (15). The equilibrium point of (12) undergoes a Hopf bifurcation at a = a 2 .

5. Numerical Simulations

The first part of the numerical analysis of Chen et al. [2] was restricted to the deterministic part of the model (3) using the following parameters: time delay τ = 4 , Nelength scale ξ = 3 , spread of credit risk among different participants engaged in credit activities is closely linked to their direct business relationships λ 1 = 0.12 , transmission rate of credit risk among various participants in credit activities is associated with their indirect business relationships λ 2 = 0.07 , and the initial number of credit activity participants who are infected by credit risk in the financial market N 0 . We keep τ = 4 and vary the nonlinear resistance coefficient μ and the weight γ on the previous activities of the credit risk participant to assess the dynamics of the credit risk contagion in the financial market on the assumption that λ 1 = 0.12 > λ 2 = 0.07 , as we compare model (4) to model (3) (see [2]).

5.1. Time Process Diagrams of Credit Risk in Financial Market on λ 1 > λ 2

The actual number of credit activity participants infected by credit risk in the financial market can be determined by multiplying the simulated value from model (4) by 10 from the N ( t ) axis compared to the values in Chen et al. [2].
From Figure 1, it is observed that at a time delay of τ = 4.0 and nonlinear resistance coefficient of μ = 0.15 , the range of credit risk and the infectious scale are diminished compared to the subsequent simulation analysis of the diagrams in the study.
It is observed in Figure 2 that as the nonlinear resistance coefficient is reviewed downward to μ = 0.08 from μ = 0.15 , at a time delay of τ = 4.0 , the infectious scale and range of credit risk are gradually increased.
In Figure 3, the nonlinear resistance coefficient is further reviewed downward to μ = 0.05 with all other parameters held at the same level compared to Figure 2 and Figure 3, and the rate of credit contagion increased to 26.25 from 8.75 in Figure 1.

5.2. Time Process Diagrams of Credit Risk in Financial Market on λ 1 < λ 2

In this section, we maintain the parameter values used in Figure 1, Figure 2 and Figure 3 on the condition that λ 1 = 0.07 < λ 2 = 0.12 and vary the weight γ on the previous activities of the credit risk participant to assess its effects on the credit risk contagion in the financial market as opposed to the assumption λ 1 > λ 2 of Chen et al. [2].
It is observed in Figure 4 that at a time delay of τ = 4.0 , on the condition that λ 1 < λ 2 , the weight on the previous activities of the credit risk participant was reduced from γ = 0.290 to γ = 0.149 to keep the infectious scale at 8.75 compared to Figure 1. This implies that the generalized model (4) is not adherent to the condition λ 1 > λ 2 as claimed in Chen et al. [2].
It is observed in Figure 5 that at a nonlinear resistance coefficient of μ = 0.08 and a time delay of τ = 4.0 , little initial oscillations were observed as the infectious scale converged to 16.25 on the condition that λ 1 = 0.07 < λ 2 = 0.12 , as the weight on past credit risk activities is reduced by fifty one ( 51 % ) .
It is indicated in Figure 6 that at a time delay of τ = 4.0 , on condition that λ 1 = 0.07 < λ 2 = 0.12 , the same rate of credit risk contagion of the participants (refer to Figure 3) can be obtained at a few wavy frequencies if the initial weight on previous credit risk activity is reduced by fifty four percent ( 54 % ) .

5.3. Real-World Application of the Model

The practical usefulness of model (4) lies in its ability to help financial creditors control the spread of credit risk contagion. This can be achieved by reviewing the weight assigned to past activities of various credit activity participants, regardless of whether they are involved in direct business relationships ( λ 1 ) or indirect business relationships ( λ 2 ).
It is obvious in Figure 7 that when the weight on past activities of the credit risk population is increased from 0.290 to 0.390 , the infectious scale of credit risk is reduced to 7.37 from 8.75 compared to Figure 1, even under the restricted condition of λ 1 > λ 2 . In all cases, the oscillations of the infectious scale decreased drastically compared to the time process diagrams in Chen et al. [2].

6. Conclusions

This study examines the stability of a credit risk contagion model that incorporates distributed delay. The model is a generalized version of the models that were previously examined. The comparative research revealed that, in contrast to the previous models that solely used the nonlinear resistance coefficient to calculate the rate of credit risk infection, the weight placed on the past behaviors of credit risk participants also influences the credit risk contagion rate. This study establishes that existing models perform best when credit risk individuals with direct business relationships have a higher contagion rate than those with indirect business relationships. Nonetheless, regardless of whether the credit risk individuals’ contagion rate is higher or lower in direct or indirect business relationships, the distributed delay model can quantify the infection rate. Our model (4) is practically useful as the financial creditors can check the spread of the contagion rate of credit risk by reviewing the weight they put on past activities among different credit activity participants exposed to risk. By doing this, the financial institution can lessen losses brought on by the borrowers’ failure to make the agreed-upon payments. As a result, modeling the transmission of credit risk contagion with dispersed delays is advised.

Author Contributions

This article is the result of the joint efforts of M.A., L.G., A.L.S., S.E.A. and H.A. Conceptualization, L.G. and M.A.; methodology, L.G. and A.L.S.; software, M.A.; validation, S.E.A., M.A. and H.A.; formal analysis, M.A. and L.G.; investigation, S.E.A. and A.L.S.; writing—original draft preparation, M.A. and L.G.; writing—review and editing, H.A. and S.E.A.; visualization, L.G. and S.E.A.; supervision, H.A. and L.G.; project administration, A.L.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The authors declare that all data supporting the findings of this study are from previously reported studies which have been cited.

Acknowledgments

The authors are highly grateful to the editor and the anonymous reviewers for their valuable suggestions and comments that have significantly enhanced this article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Kai, X.; Mo, J.; Qian, Q.; Zhang, F.; Xie, X.; Zhou, Z. Associated credit risk contagion with incubatory period: A network-based perspective. Complexity 2020, 2020, 5642730. [Google Scholar]
  2. Chen, T.; Li, X.; He, J. Complex dynamics of credit risk contagion with time-delay and correlated noises. Abstr. Appl. Anal. 2014, 2014, 456764. [Google Scholar] [CrossRef]
  3. Gakkhar, S.; Singh, A. Complex dynamics in a prey predator system with multiple delays. Commun. Nonlinear Sci. Numer. Simul. 2012, 17, 914–929. [Google Scholar] [CrossRef]
  4. Anokye, M.; Donkoh, E.K.; Ofori, M.F.; Adom-Konadu, A. Full delay logistic population model with sustainable harvesting. J. Appl. Math. 2022, 2022, 6777827. [Google Scholar] [CrossRef]
  5. Anokye, M.; Barnes, B.; Boateng, F.O.; Adom-Konadu, A.; Amoah-Mensahand, J. Price dynamics of a delay differential cobweb model. Discret. Dyn. Nat. Soc. 2023, 2023, 1296562. [Google Scholar] [CrossRef]
  6. Wu, D.; Zhu, S. Stochastic resonance in FitzHugh-Nagumo system with time-delayed feedback. Phys. Lett. A 2008, 32, 5299–5304. [Google Scholar] [CrossRef]
  7. Anokye, M.; Adom-Konadu, A.; Sackitey, A.L.; Ofori, M.F. Delay logistic model with harvesting on weighted average population. Adv. Dyn. Syst. Appl. 2022, 17, 467–481. [Google Scholar]
  8. Wangersky, W.J.; Cunningham, P.J. Time lag in population models. Cold Spring Harbor Symp. Quant. Biol. 1957, 22, 329–338. [Google Scholar] [CrossRef]
  9. Erneux, T. Applied Delay Differential Equations; Springer: Berlin/Heidelberg, Germany, 2009. [Google Scholar]
  10. Ma, Y. Global Hopf bifurcation in the Leslie-Gower predator-prey model with two delays. Nonlinear Anal. Real World Appl. 2012, 13, 370–375. [Google Scholar] [CrossRef]
  11. Song, Y.; Wei, J.; Han, M. Local and global Hopf bifurcation in a delayed hematopoiesis model. Int. J. Bifurc. Chaos Appl. Sci. Eng. 2004, 14, 3909–3919. [Google Scholar] [CrossRef]
  12. Chen, T.; He, J.; Yin, Q. Dynamics evolution of credit risk contagion in the CRT market. Discret. Dyn. Nat. Soc. 2013, 2013, 206201. [Google Scholar] [CrossRef]
  13. Fanelli, V.; Maddalena, L. A nonlinear dynamic model for credit risk contagion. Math. Comput. Simul. 2020, 174, 45–58. [Google Scholar] [CrossRef]
  14. Ma, J.; Liu, Y.; Zhao, L.; Liang, W. Research on the mechanism and application of spatial credit risk contagion based on complex network model. Manag. Decis. Econ. 2024, 45, 1180–1193. [Google Scholar] [CrossRef]
  15. Chen, N.; Fan, H. Credit risk contagion and optimal dual control-An SIS/R model. Math. Comput. Simul. 2024, 210, 448–472. [Google Scholar] [CrossRef]
  16. Irakoze, I.; Nahayo, F.; Ikpe, D.; Gyamerah, S.A.; Viens, F. Mathematical modeling and stability analysis of systemic risk in the banking ecosystem. J. Appl. Math. 2023, 2023, 5628621. [Google Scholar] [CrossRef]
  17. Aliano, M.; Canana, L.; Cestari, G.; Ragni, S. A dynamical model with time delay for risk contagion. Mathematics 2023, 11, 425. [Google Scholar] [CrossRef]
  18. Chen, T.; Yang, Q.; Wang, Y.; Wang, S. Double-layer network model of bank-enterprise counterparty credit risk contagion. Complexity 2020, 2020, 3690848. [Google Scholar] [CrossRef]
  19. Ma, G.; Ding, J.; Lv, Y. The credit risk contagion mechanism of financial guarantee network: An Application of the SEIR-epidemic model. Complexity 2022, 2022, 7669259. [Google Scholar] [CrossRef]
  20. Sui, X.; Wen, H.; Gao, J.; Lu, S. Research on contagion and the influencing factors of personal credit risk based on a complex network. Discret. Dyn. Nat. Soc. 2022, 2022, 4730479. [Google Scholar] [CrossRef]
  21. Ebrahimi Dehshalie, M.; Kabiri, M.; Ebrahimi Dehshali, M. Stability analysis and fixed-time control of credit risk contagion. Math. Comput. Simul. 2021, 190, 131–139. [Google Scholar] [CrossRef]
  22. Bianca, C.; Guerrini, L. Hopf bifurcations in a delayed microscopic model of credit risk contagion. Appl. Math. Inf. 2015, 9, 1493–1497. [Google Scholar]
  23. Elaiw, A.M.; Al Agha, A.D. A reaction-diffusion model for oncolytic M1 virotherapy with distributed delays. Eur. Phys. J. Plus 2020, 135, 1–40. [Google Scholar]
  24. Sargood, A.; Gaffney, E.A.; Krause, A.L. Fixed and distributed gene expression time delays in reaction-diffusion systems. Bull. Math. Biol. 2022, 84, 1–29. [Google Scholar] [CrossRef] [PubMed]
  25. Hu, S.; Dunlavey, M.; Guzy, S.; Teuscher, N. A distributed delay approach for modeling delayed outcomes in pharmacokinetics and pharmacodynamics studies. J. Pharmacokinet. Pharmacodyn. 2018, 45, 285–308. [Google Scholar] [CrossRef] [PubMed]
  26. Cushing, J.M. Integrodifferential Equations and Delay Models in Population Dynamics; Springer: Berlin/Heidelberg, Germany, 1977. [Google Scholar]
  27. Hernandez, E.; Wu, J. Existence, uniqueness and qualitative properties of global solutions of abstract differential equations with state-dependent delay. Proc. Edinb. Math. Soc. 2019, 62, 771–788. [Google Scholar] [CrossRef]
  28. Podlubny, I. Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications; Elsevier: Amsterdam, The Netgerlands, 1998. [Google Scholar]
  29. Ahmadian, A.; Ismail, F.; Senu, N.; Salahshour, S.; Suleiman, M. Toward the existence and uniqueness of solutions for fractional integro-differential equations under uncertainty. AIP Conf. Proc. 2016, 1739, 020034. [Google Scholar]
  30. MacDonald, N. Time Lags in Biological Models; Springer: New York, NY, USA, 1978. [Google Scholar]
Figure 1. Dynamics of credit risk contagion under nonlinear resistance coefficient of μ = 0.15 .
Figure 1. Dynamics of credit risk contagion under nonlinear resistance coefficient of μ = 0.15 .
Axioms 13 00483 g001
Figure 2. Dynamics of credit risk contagion under nonlinear resistance coefficient of μ = 0.08 .
Figure 2. Dynamics of credit risk contagion under nonlinear resistance coefficient of μ = 0.08 .
Axioms 13 00483 g002
Figure 3. Nonlinear resistance coefficient influence on credit risk contagion at μ = 0.05 .
Figure 3. Nonlinear resistance coefficient influence on credit risk contagion at μ = 0.05 .
Axioms 13 00483 g003
Figure 4. Effects of past credit risk activities on credit risk contagion at γ = 0.149 .
Figure 4. Effects of past credit risk activities on credit risk contagion at γ = 0.149 .
Axioms 13 00483 g004
Figure 5. Effects of past credit risk activities on credit risk contagion at γ = 0.079 .
Figure 5. Effects of past credit risk activities on credit risk contagion at γ = 0.079 .
Axioms 13 00483 g005
Figure 6. Effects of past credit risk activities on credit risk contagion at γ = 0.052 .
Figure 6. Effects of past credit risk activities on credit risk contagion at γ = 0.052 .
Axioms 13 00483 g006
Figure 7. Dynamics of higher weight on contagion rate of credit risk at γ = 0.390 .
Figure 7. Dynamics of higher weight on contagion rate of credit risk at γ = 0.390 .
Axioms 13 00483 g007
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Anokye, M.; Guerrini, L.; Sackitey, A.L.; Assabil, S.E.; Amankwah, H. Stability Analysis of a Credit Risk Contagion Model with Distributed Delay. Axioms 2024, 13, 483. https://doi.org/10.3390/axioms13070483

AMA Style

Anokye M, Guerrini L, Sackitey AL, Assabil SE, Amankwah H. Stability Analysis of a Credit Risk Contagion Model with Distributed Delay. Axioms. 2024; 13(7):483. https://doi.org/10.3390/axioms13070483

Chicago/Turabian Style

Anokye, Martin, Luca Guerrini, Albert L. Sackitey, Samuel E. Assabil, and Henry Amankwah. 2024. "Stability Analysis of a Credit Risk Contagion Model with Distributed Delay" Axioms 13, no. 7: 483. https://doi.org/10.3390/axioms13070483

APA Style

Anokye, M., Guerrini, L., Sackitey, A. L., Assabil, S. E., & Amankwah, H. (2024). Stability Analysis of a Credit Risk Contagion Model with Distributed Delay. Axioms, 13(7), 483. https://doi.org/10.3390/axioms13070483

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop