Systemic risk and uncertainty quantification

Transcription

1 Systemic risk and uncertainty quantification Josselin Garnier (Université Paris 7) with George Papanicolaou and Tzu-Wei Yang (Stanford University) October 6, 22

2 Modeling Systemic Risk We consider evolving systems with a large number of inter-connected components, each of which can be in a normal state or in a failed state. We want to study the probability of overall failure of the system, that is, its systemic risk. There are three effects that we want to model and that contribute to the behavior of systemic risk: The intrinsic stability of each component The strength of external random perturbations to the system The degree of inter-connectedness or cooperation between components

3 Possible Applications Engineering systems with a large number of interacting parts. Components can fail but the system fails only when a large number of components fail simultaneously. Power distribution systems. Individual components of the system are calibrated to withstand fluctuations in demand by sharing loads, but sharing also increases the probability of an overall failure. Banking systems. Banks cooperate and by spreading the risk of credit shocks between them can operate with less restrictive individual risk policies. However, this increases the risk that they may all fail, that is, the systemic risk. We want to propose a simple model to explain that individual risk does not affect the systemic risk in an obvious way. In fact, it is possible to simultaneously reduce individual risk and increase the systemic risk.

4 A Basic, Bistable Mean-Field Model The (banking) system has N components (banks). The risk variable X j (t) of each component j =,...,N, satisfies the SDE dx j (t) = hv (X j (t))dt +σdw j (t)+θ ( X(t) Xj (t) ) dt. V(x) = x 4 /4 x 2 /2 is a potential function with two stable states ±. With θ = σ =, X j (t) stays in one of the states. We define the normal state and + the failed state. h is the intrinsic stability. The average X(t) := N N i= X i(t) is the risk variable of the system. θ > is the attractive interaction. {W j (t),j =,...,N} are independent Brownian motions serving as the external uncertainties. σ is the noise strength.

5 Why this Model? The three parameters, h, σ and θ control the three effects we want to study: intrinsic stability, random perturbations, and the degree of cooperation, respectively. Why mean field interaction? Because it is the simplest interaction that models cooperative behavior. It can be generalized to include diversity as well as other more complex interactions such as hierarchical ones. Connection with QMU (Quantification of Margins of Uncertainty): In an inter-connected system, individual components can be operating closer to their margin of failure, as they can benefit from the stability of the rest of the system. This, however, reduces the overall margin of uncertainty, that is, increases the systemic risk.

6 Schematic for the Individual Risk and the Systemic Risk - + Status

7 The probability-measure-valued process µ N (t) First idea: to analyze X(t) = N N i= X i(t). Unfortunately, no closed equation for X(t), so we need to generalize this problem into a larger space. We define a probability-measure-valued process µ N (t,dx) := N N j= δ X j (t)(dx). Then X(t) = yµ N (t,dy). (Dawson, 983) µ N (t,dx) converges in law as N to a deterministic process u(t, x)dx, the solution of the nonlinear Fokker-Planck equation: t u = h [ V (x)u ] θ {[ ] 2 2 yu(t,y)dy x u }+ σ2 x x x 2u, recalling that dx j (t) = hv (X j (t))dt +θ( X(t) X j (t))dt +σdw j (t).

8 Existence of two stable equilibria The equilibrium u e ξ (x) = lim t u(t,x) has the form: { uξ e (x) = (x ξ)2 exp Z ξ 2 σ2 2θ } 2h σ 2V(x), with the compatibility condition: ξ = m(ξ) := yue ξ (y)dy. Given θ and h, there exists a critical value σ c such that ξ = m(ξ) has two stable solutions ±ξ b if and only if σ < σ c. Simplification: for small h, we have σc 2 = 2θ/3+O(h) and ξ b = 3σ2 2θ +O(h). Let us say that u ξ e b is the normal state of the system, u+ξ e b is the failed state. We have µ N (t,dx) N u(t,x)dx. If u() = u ξ e b, then u(t) = u ξ e b for all time.

9 Simulation of X: impact of h and N N=, h=.5, σ=, θ= dm()/dξ=.44, 3σ 2 /2θ=.5.5 N=, h=., σ=, θ= dm()/dξ=.86, 3σ 2 /2θ=.5.5 Systemic Risk.5 Systemic Risk t.5 N=2, h=., σ=, θ= dm()/dξ=.86, 3σ 2 /2θ= t.5 N=8, h=., σ=, θ= dm()/dξ=.86, 3σ 2 /2θ=.5 Systemic Risk.5.5 Systemic Risk t t X n+ j = Xj n hv (Xj n ) t +σ Wj n+ +θ( N N k= X n k Xn j ) t.

10 Simulation of X: impact of θ and σ.5 N=, h=., σ=, θ= dm()/dξ=.86, 3σ 2 /2θ=.5.5 N=, h=., σ=, θ=5 dm()/dξ=.4, 3σ 2 /2θ=.3 Systemic Risk.5.5 Systemic Risk t.5 N=, h=., σ=, θ=6 dm()/dξ=.26, 3σ 2 /2θ= t 2 N=, h=., σ=4, θ=6 dm()/dξ=.9572, 3σ 2 /2θ=4.5 Systemic Risk.5.5 Systemic Risk t t X n+ j = Xj n hv (Xj n ) t +σ Wj n+ +θ( N N k= X n k Xn j ) t.

11 Probability of system failures Assume that σ < σ c so the limit of µ N (t) has two stable equilibria u e ξ b, the system s normal/failed states. For large N, let the empirical density µ N () u e ξ b. Then we expect that µ N (t) u e ξ b for all t >. However, as long as N is finite, the system collapse: µ N () u e ξ b, µ N (T) u e +ξ b happens in the time interval [,T] with small but nonzero probability. We use large deviations to compute this small probability.

12 Large deviation principle (Dawson & Gärtner, 987) Given a rare event A in the suitable space (smaller than C([,T],M (R))) Detail, then µ N (t,dx) = N N i= δ X j (t)(dx) satisfies the large deviation principle with the rate function I h : Detail ( ) P(µ N A) N exp N inf I h(φ). φ A The rate function I h has the following variational form: I h (φ) = T 2σ 2 sup f: φ,f 2 L h φ = 2 {[ σ2 2 x 2φ+θ x x µ,f = φ t L h φ,f 2 x φ,fx 2 dt, ] } yφ(t,dy) φ +h [ V (x)φ ], x µ (dx)f(x).

13 How to interpret the rate function? I h (φ) = T 2σ 2 sup f: φ,fx 2 φ t L h φ,f 2 / φ,fx 2 dt. Recall that the limit u(t,x)dx µ e ξ b (x)dx of µ N (t,dx) satisfies u t L h u =, and I h(u) =. So P(µ N A) N exp( Ninf φ A I h (u)) = if u A. Given a smooth probability-measure valued process φ u, let g satisfy φ t L h φ = (φg) x. () g is the external force/randomness making φ deviate from u. Note that by (), sup f: φ,f 2 φ t L h φ,f 2 y φ,fx 2 φ,f x = sup g 2 f: φ,fx 2 φ,fx 2 = φ,g 2. Then I h (φ) = φ,g2 dt, which is the L 2 -norm of g (the randomness) and indicates the price to pay to have µ N φ. This expression () exists only if (very unfortunately) φ is a density function. 2σ 2 T

14 Probability of system failures and small h analysis P(µ N A) N exp( Ninf φ A I h (φ)). The rare event A of system failures is the set of all possible paths starting from u e ξ b (the normal state) to u e +ξ b (the failed state): A = { φ : φ() = u e ξ b, φ(t) = u e +ξ b }. The major work is to compute inf φ A I h (φ), which is a nonlinear and infinitely dimensional problem. Here we assume that the intrinsic stability, h, is small so that we can make this problem tractable. Two good reasons to consider a small intrinsic stability h:. P(µ N A) is extremely small for large h, which is not the regime we have in mind. 2. The case h = is analytically solvable and the cases h small are perturbations of it.

15 The case h = In this case, inf φ A I (φ) can be solved without approximations. Result: For h =, inf φ A I (φ) has the unique minimizer p e (t,x)dx, a path of Gaussian measures: { } p e (t,x) = exp (x ae (t)) 2, 2π σ2 2 σ2 2θ 2θ where a e (t) = 2ξ b t/t ξ b and ξ b = 3σ2 2θ. Therefore, inf I (φ) = I (p e ) = 2ξ2 b φ A σ 2 T = 2 σ 2 T ) ( 3σ2. 2θ Recall that I h (φ) = sup f: φ,fx 2 φ t L h φ,f 2 / φ,fx 2 dt } and A = {φ : φ() = u e, φ(t) = u e. ξb +ξb 2σ 2 T

16 The case h small When h is small, good candidates for the transition path of empirical densities are Gaussian with small perturbations: { [ (x a(t)) 2 φ = p +hq : p(t,x) = 2πb2 (t) exp 2b 2 (t) ],a() = ξ b,a(t) = ξ b } For h small, the large deviation problem is solvable approximately, and the transition probability is P(µ N A) ( [ exp N ) 2 ( 3σ2 + 6h σ 2 T 2θ σ 2 ( σ 2 Here are some comments of this result: θ ) 2 ) ]) ( σ2 +O(h 2 ). θ A large system is more stable than a small system. In the long run (T large), a transition will happen. Increase of the intrinsic stabilization parameter h reduces systemic risk. Mean transition times are simply related to transition probabilities in this approximation (Williams 82)

17 What about the individual risk? Recall that the risk variable X j (t) of component j satisfies dx j = h(x 3 j X j )dt +θ( X X j )dt +σdw j. Assume X j () = (the normal state) and linearize X j around : X j (t) = +δx j (t), X(t) = +δ X(t) and δ X(t) = N N j= δx j(t). δx j (t) and δ X(t) satisfy linear SDEs: dδx j = (θ+2h)δx j dt+θδ Xdt+σdW j, dδ X = 2hδ Xdt+ σ N δx j (t) is a Gaussian process with the stationary distribution σ N(, 2 2(2h+θ) ) as N. Qualitatively speaking, the individual risk is external risk(σ 2 ) intrinsic stability(h)+risk diversification(θ), and note that σ 2 2(2h+θ) σ2 2θ for small h. N dw j. j=

18 Why and when is risk diversification undesirable? For h small: Systemic Risk exp ( 2N σ 2 T )) ( 3σ2, Individual Risk σ2 2θ 2θ. If a banker wants to take more external risk (to earn more profits) or the economy is more uncertain (increase σ 2 ), in order to keep the individual risk low, he can diversify the risk to the other banks (increase θ). Let us assume that everyone simultaneously increase σ 2 and θ but carefully keep the ratio σ2 2θ at a low level. The individual risk is kept low. The systemic risk increases, although this cannot be observed by the banks (until the catastrophic transition happens).

19 Modeling of Diversity in Cooperative Behavior t u t u K The cooperative behavior of components can be different across groups: dx j (t) = hv (X j (t))dt +σdw j (t)+θ j ( X (t) X j (t) ) dt. The components are partitioned into K groups. In group k, the components have cooperative parameter Θ k. In the limit N the empirical densities of each group converge to the solution of the joint Fokker-Planck equations: = h [ V ] 2 (x)u + x 2 σ2 x 2u Θ x. = h [ V ] 2 (x)u K + x 2 σ2 x 2u K Θ K x where ρ k N is the size of group k. {[ y {[ y } K ρ k u k (t,y)dy x ]u k= } K ρ k u k (t,y)dy x ]u K k=

20 Impact of Component Diversity on the Systemic Risk Why is the diversity interesting? The model is more realistic and more widely applicable. Diversity significantly affects the system stability by reducing it. Impact from the diversity: Analytical and numerical studies show that even with the same parameters and with {θ j } whose average equals θ the system still changes significantly.

21 Simulation 2 - Impact of Diversity, Change of Θ k and ρ k N=, h=., σ=, Θ=[;;] ρ=[.33;.34;.33], stdev(θ)/mean(θ)= dm()/dξ=.86, Σ (ρ /Θ )(3σ 2 /2Θ )=.85 i i i i N=, h=., σ=, Θ=[6;;4] ρ=[.33;.34;.33], stdev(θ)/mean(θ)=.3266 dm()/dξ=.92, Σ (ρ /Θ )(3σ 2 /2Θ )=.996 i i i i N=, h=., σ=, Θ=[2;;8] ρ=[.33;.34;.33], stdev(θ)/mean(θ)=.6532 dm()/dξ=.89, Σ (ρ /Θ )(3σ 2 /2Θ )= i i i i Systemic Risk.5.5 Systemic Risk.5.5 Systemic Risk t t t N=, h=., σ=, Θ=[5;;5] ρ=[.;.8;.], stdev(θ)/mean(θ)= dm()/dξ=.89, Σ i (ρ i /Θ i )(3σ 2 /2Θ i )=.88 N=, h=., σ=, Θ=[5;;5] ρ=[.3;.4;.3], stdev(θ)/mean(θ)= dm()/dξ=.95, Σ i (ρ i /Θ i )(3σ 2 /2Θ i )=.94 N=, h=., σ=, Θ=[5;;5] ρ=[.45;.;.45], stdev(θ)/mean(θ)= dm()/dξ=.99, Σ i (ρ i /Θ i )(3σ 2 /2Θ i )= Systemic Risk.5.5 Systemic Risk.5.5 Systemic Risk t t t

22 Analysis in the Diversity Case The critical value of the system fluctuation is lower: of the homogeneous case is θ, and The critical value σ homo c is σ div c = [( K k= ρ k Θ k ) / ( K k= 3ρ k 2Θ 2 k )] /2. K If θ = k= ρ kθ k, then σc homo σc div. System with diversity have larger transition probabilities: When h is zero, the average of θj is θ, and the diversity of θ j is small, then the system has a higher transition probability. Mathematically, if h =, Θk = θ(+δα k ) with δ, and K k= ρ kα k =, then ( ) P µ div N A { [ exp N ) ( K ) 2 ( 3σ2 (3σ 2δ 2 ρ σ 2 k α 2 2 k T 2θ 2θ + T k= T 2 3 ( e θs) 2 ds ) ]}.

23 A Hierarchical Model of Systemic Risk Here we consider a hierarchical model with a central agent: dx = σ dw h V (X )dt θ ( X N dx j = σdw j hv (X j )dt θ ( X j X ) dt, N ) X j dt j= j =,...,N X models the central stable agent. It is intrinsically stable (h > ), and not subjected to external fluctuations (σ = ). It interacts with the other agents through a mean field interaction. X j, j =,...,N model individual agents that are subjected to external fluctuations. They are (h > ) or are not (h = ) intrinsically stable. They interact with the stable agent X.

24 A Hierarchical Model of Systemic Risk - Analysis Nonlinear Fokker-Planck equations: In the limit N the pair ( X (t), N N j= δ X j (t)(dx) ) converges to (x (t),p(t,x)dx) solution of the nonlinear Fokker-Planck equation t p = σ2 [ 2 2 xxp + x hv (x)+θ(x x (t))p ], with dx dt = h V (x ( ) θ x xp(t,x)dx ). Existence of two equilibirum states when σ is below a critical level: x x e, p(x) ( Z(x e ) exp 2hV(x)+θ(x x e) 2 ) σ 2 with the compatibility equation xp(x)dx = x e + h θ V (x e). Large deviations principle to compute the probability of transition.

25 A Hierarchical Model of Systemic Risk - Results Expansions for small h (for the optimal paths and for the probability of transition). Dimension reduction: min-max problem for the rate function can be reduced to a the resolution of a fourth-order ODE for the optimal x (t) with boundary conditions. For the optimal path the mean of the individual agents x(t) is ahead of x (t): the individual agents drive the transition. Stability increases with θ and decreases with θ.

26 On-going and future work Compute the transition probability of the diversity case for small h. Study models with diversities everywhere: The most general case (when the limit exists) is that V, σ and θ in each group can be different, i.e. if X j is in the group k, then the SDE is dx j (t) = hv k (X j (t))dt+σ k dw j (t)+θ k ( X (t) Xj (t) ) dt. Show that it is unconditionally true that systems with diversity are less stable than homogeneous ones? Using the analysis as a guide, we design importance sampling algorithms for computing efficiently (very) small systemic failure probabilities.

27 Topological spaces for the mean field model M (R) is the space of probability measures on R with the Prohorov metric ρ, associated with the weak convergence. C([,T],M (R)) is the space of continuous functions from [,T] to M (R) with the metric sup t T ρ(µ (t),µ 2 (t)). M (R) = {µ M (R), ϕ(y)µ(dy) < }, where ϕ = +y 2 +λy 4, λ h/2, serves as a Lyapunov function. M (R) is endowed with the inductive topology: µ n µ in M (R) if and only if µ n µ in M (R) and sup n ϕ(y)µ(dy) <. C([,T],M (R)) is the space of continuous functions from [,T] to M (R) endowed with the topology: φ n ( ) φ( ) in C([,T],M (R)) if and only if φ n ( ) φ( ) in C([,T],M (R)) and sup t T sup n ϕ(y)φn (t,dy) <. Go Back

28 Large Deviations Principle Exact statement: inf I h (φ) liminf φ A N limsup N N logp(µ N A) N logp(µ N A) inf I h (φ) φ A Go Back