Skip to main content
Log in

An SPDE model for systemic risk with endogenous contagion

  • Published:
Finance and Stochastics Aims and scope Submit manuscript

Abstract

We propose a dynamic mean-field model for ‘systemic risk’ in large financial systems, derived from a system of interacting diffusions on the positive half-line with an absorbing boundary at the origin. These diffusions represent the distances-to-default of financial institutions, and absorption at zero corresponds to default. As a way of modelling correlated exposures and herd behaviour, we consider a common source of noise and a form of mean-reversion in the drift. Moreover, we introduce an endogenous contagion mechanism whereby the default of one institution causes a drop in the distances-to-default of the other institutions. In this way, we aim to capture key ‘system-wide’ effects on risk. The resulting mean-field limit is characterised uniquely by a nonlinear SPDE on the half-line with a Dirichlet boundary condition. The density of this SPDE gives the conditional law of a non-standard ‘conditional’ McKean–Vlasov diffusion, for which we provide a novel upper Dirichlet heat kernel type estimate. Depending on the realisations of the common noise and the rate of mean-reversion, the SPDE can exhibit rapid accelerations in the loss of mass at the boundary. In other words, the contagion mechanism can give rise to periods of significant systemic default clustering.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Institutional subscriptions

Fig. 2.1
Fig. 2.2

Similar content being viewed by others

References

  1. Acharya, V.V., Pedersen, L.H., Philippon, T., Richardson, M.: Measuring systemic risk. Rev. Financ. Stud. 30, 2–47 (2017)

    Article  Google Scholar 

  2. Allen, F., Gale, D.: Financial contagion. J. Polit. Econ. 108, 1–33 (2000)

    Article  Google Scholar 

  3. Azizpour, S., Giesecke, K., Schwenkler, G.: Exploring the sources of default clustering. J. Financ. Econ. 129, 154–183 (2018)

    Article  Google Scholar 

  4. Banerjee, T., Bernstein, A., Feinstein, Z.: Dynamic clearing and contagion in financial networks. Preprint (2018). Available online at: arXiv:1801.02091

  5. Benoit, S., Colliard, J.-E., Hurlin, C., Pérignon, C.: Where the risks lie: a survey on systemic risk. Rev. Finance 21, 109–152 (2017)

    Article  MATH  Google Scholar 

  6. Berg, T., Gider, J.: What explains the difference in leverage between banks and non-banks? J. Financ. Quant. Anal. 52, 2677–2702 (2017)

    Article  Google Scholar 

  7. Brunnermeier, M.K.: Deciphering the liquidity and credit crunch 2007–2008. J. Econ. Perspect. 23, 77–100 (2009)

    Article  Google Scholar 

  8. Brunnermeier, M.K., Gorton, G., Krishnamurthy, A.: Risk topography. NBER Macroecon. Annu. 26, 149–176 (2012)

    Article  Google Scholar 

  9. Bush, N., Hambly, B.M., Haworth, H., Jin, L., Reisinger, C.: Stochastic evolution equations in portfolio credit modelling. SIAM J. Financ. Math. 2, 627–664 (2011)

    Article  MathSciNet  MATH  Google Scholar 

  10. Capponi, A., Bo, L.: Systemic risk in interbanking networks. SIAM J. Financ. Math. 6, 386–424 (2015)

    Article  MathSciNet  MATH  Google Scholar 

  11. Carmona, R., Fouque, J.-P., Sun, L.-H.: Mean field games and systemic risk. Commun. Math. Sci. 13, 911–933 (2015)

    Article  MathSciNet  MATH  Google Scholar 

  12. Cho, S., Kim, P., Park, H.: Two-sided estimates on Dirichlet heat kernels for time-dependent parabolic operators with singular drifts in \(C^{1,\alpha }\)-domains. J. Differ. Equ. 252, 1101–1145 (2012)

    Article  MATH  Google Scholar 

  13. Cochrane, J.H.: Toward a run-free financial system. In: Baily, M.N., Taylor, J.B. (eds.) Across the Great Divide, pp. 197–249. Hoover Press, Stanford (2014)

    Google Scholar 

  14. Cont, R., Moussa, A., Santos, E.B.: Network structure and systemic risk in banking systems. In: Fouque, J.-P., Langsam, J. (eds.) Handbook of Systemic Risk, pp. 327–368. Cambridge University Press, Cambridge (2013)

    Chapter  Google Scholar 

  15. Cont, R., Wagalath, L.: Running for the exit: distressed selling and endogenous correlation in financial markets. Math. Finance 23, 718–741 (2013)

    Article  MathSciNet  MATH  Google Scholar 

  16. Cont, R., Wagalath, L.: Fire sale forensics: measuring endogenous risk. Math. Finance 26, 835–866 (2016)

    Article  MathSciNet  MATH  Google Scholar 

  17. Cvitanić, J., Ma, J., Zhang, J.: The law of large numbers for self-exciting correlated defaults. Stoch. Process. Appl. 122, 2781–2810 (2012)

    Article  MathSciNet  MATH  Google Scholar 

  18. Danielsson, J., Shin, H.S., Zigrand, J.-P.: Endogenous extreme events and the dual role of prices. Annu. Rev. Econ. 4, 111–129 (2012)

    Article  Google Scholar 

  19. Dees, S., Henry, J., Martin, R.: STAMPE: Stress-test analytics for macroprudential purposes in the Euro area. ECB (2017). Available online at: https://www.ecb.europa.eu/pub/pdf/other/stampe201702.en.pdf

  20. Delarue, F., Inglis, J., Rubenthaler, S., Tanré, E.: Global solvability of a networked integrate-and-fire model of McKean–Vlasov type. Ann. Appl. Probab. 25, 2096–2133 (2015)

    Article  MathSciNet  MATH  Google Scholar 

  21. Delarue, F., Inglis, J., Rubenthaler, S., Tanré, E.: Particle systems with a singular mean-field self-excitation. Application to neuronal networks. Stoch. Process. Appl. 125, 2451–2492 (2015)

    Article  MathSciNet  MATH  Google Scholar 

  22. Duffie, D.: Financial regulatory reform after the crisis: an assessment. Manag. Sci. 64, 4835–4857 (2018)

    Article  Google Scholar 

  23. Duffie, D., Eckner, A., Horel, G., Saita, L.: Frailty correlated default. J. Finance 64, 2089–2123 (2009)

    Article  Google Scholar 

  24. Eisenberg, L., Noe, T.H.: Systemic risk in financial systems. Manag. Sci. 47, 236–249 (2001)

    Article  MATH  Google Scholar 

  25. Fouque, J.-P., Ichiba, T.: Stability in a model of interbank lending. SIAM J. Financ. Math. 4, 784–803 (2013)

    Article  MathSciNet  MATH  Google Scholar 

  26. Fouque, J.-P., Sun, L.-H.: Systemic risk illustrated. In: Fouque, J.-P., Langsam, J. (eds.) Handbook on Systemic Risk, pp. 444–452. Cambridge University Press, Cambridge (2013)

    Chapter  Google Scholar 

  27. Garnier, J., Papanicolaou, G., Yang, T.-W.: Large deviations for a mean field model of systemic risk. SIAM J. Financ. Math. 4, 151–184 (2013)

    Article  MathSciNet  MATH  Google Scholar 

  28. Garnier, J., Papanicolaou, G., Yang, T.-W.: A risk analysis for a system stabilized by a central agent. Risk Decis. Anal. 6, 97–120 (2017)

    Article  MATH  Google Scholar 

  29. Giesecke, K., Schwenkler, G., Sirignano, J.: Inference for large financial systems. Math. Finance (2019, forthcoming). Available online at: https://doi.org/10.2139/ssrn.3012751

  30. Giesecke, K., Spiliopoulos, K., Sowers, R.: Default clustering in large portfolios: typical events. Ann. Appl. Probab. 23, 348–385 (2013)

    Article  MathSciNet  MATH  Google Scholar 

  31. Giesecke, K., Spiliopoulos, K., Sowers, R., Sirignano, J.: Large portfolio asymptotics for loss from default. Math. Finance 25, 77–114 (2015)

    Article  MathSciNet  MATH  Google Scholar 

  32. Glasserman, P., Young, H.P.: Contagion in financial networks. J. Econ. Lit. 54, 779–831 (2016)

    Article  Google Scholar 

  33. Gornall, W., Strebulaev, I.A.: Financing as a supply chain: the capital structure of banks and borrowers. J. Financ. Econ. 129, 510–530 (2018)

    Article  Google Scholar 

  34. Gyrya, P., Saloff-Coste, L.: Neumann and Dirichlet Heat Kernels in Inner Uniform Domains. Astérisque, vol. 336. Société Mathématique de France, Paris (2011)

    MATH  Google Scholar 

  35. Hambly, B.M., Ledger, S.: A stochastic McKean–Vlasov equation for absorbing diffusions on the half-line. Ann. Appl. Probab. 27, 2698–2752 (2017)

    Article  MathSciNet  MATH  Google Scholar 

  36. Hambly, B.M., Ledger, S., Søjmark, A.: A McKean–Vlasov equation with positive feedback and blow-ups. Ann. Appl. Probab. (2019, forthcoming). Available online at: arXiv:1801.07703

  37. Hellwig, M.F.: Systemic risk in the financial sector: an analysis of the subprime-mortgage financial crisis. Economist 157, 129–207 (2009)

    Article  Google Scholar 

  38. Hurd, T.: Contagion! Systemic Risk in Financial Networks. Springer Briefs in Quantitative Finance. Springer, Berlin (2016)

    Book  MATH  Google Scholar 

  39. Inglis, J., Talay, D.: Mean-field limit of a stochastic particle system smoothly interacting through threshold hitting-times and applications to neural networks with dendritic component. SIAM J. Math. Anal. 47, 3884–3916 (2015)

    Article  MathSciNet  MATH  Google Scholar 

  40. Jakubowski, A.: The almost sure Skorokhod representation for subsequences in nonmetric spaces. Theory Probab. Appl. 42, 167–175 (1998)

    Article  MathSciNet  Google Scholar 

  41. Kallianpur, G., Xiong, J.: Stochastic Differential Equations in Infinite Dimensional Spaces. Lecture Notes, Monograph Series. Institute of Mathematical Statistics, Bethesda (1995)

    MATH  Google Scholar 

  42. Karatzas, I., Shreve, S.E.: Brownian Motion and Stochastic Calculus. Springer, Berlin (1991)

    MATH  Google Scholar 

  43. Lando, D., Nielsen, M.: Correlation in corporate defaults: contagion or conditional independence? J. Financ. Intermed. 19, 355–372 (2010)

    Article  Google Scholar 

  44. Ledger, S.: Sharp regularity near an absorbing boundary for solutions to second order SPDEs in a half-line with constant coefficients. Stoch. Partial Differ. Equ., Anal. Computat. 2, 1–26 (2014)

    MathSciNet  MATH  Google Scholar 

  45. Ledger, S.: Skorokhod’s M1 topology for distribution-valued processes. Electron. Commun. Probab. 21, 1–11 (2016)

    Article  MathSciNet  MATH  Google Scholar 

  46. Moreno-Bote, R., Parga, N.: Response of integrate-and-fire neurons to noisy inputs filtered by synapses with arbitrary time scales: firing rate and correlations. Neural Comput. 22, 1528–1572 (2010)

    Article  MathSciNet  MATH  Google Scholar 

  47. Nadtochiy, S., Shkolnikov, M.: Particle systems with singular interaction through hitting times: application in systemic risk modeling. Ann. Appl. Probab. 29, 89–129 (2019)

    Article  MathSciNet  MATH  Google Scholar 

  48. Ostojic, S., Brunel, N., Hakim, V.: Synchronization properties of networks of electrically coupled neurons in the presence of noise and heterogeneities. J. Comput. Neurosci. 26, 369–392 (2009)

    Article  MathSciNet  Google Scholar 

  49. Pedersen, L.H.: When everyone runs for the exit. Int. J. Cent. Bank. 5, 177–199 (2009)

    Google Scholar 

  50. Schwartz, L.: Radon Measures on Arbitrary Topological Spaces and Cylindrical Measures. Tata Institute of Fundamental Research. Oxford University Press, London (1973)

    MATH  Google Scholar 

  51. Sowers, R., Spiliopoulos, K.: Default clustering in large pools: large deviations. SIAM J. Financ. Math. 6, 86–116 (2015)

    Article  MathSciNet  MATH  Google Scholar 

  52. Spiliopoulos, K.: Systemic risk and default clustering for large financial systems. In: Friz, P., et al. (eds.) Large Deviations and Asymptotic Methods in Finance, pp. 529–557. Springer, Berlin (2015)

    Chapter  Google Scholar 

  53. Whitt, W.: Stochastic-Process Limits: An Introduction to Stochastic-Process Limits and Their Application to Queues. Springer, Berlin (2002)

    Book  MATH  Google Scholar 

Download references

Acknowledgements

We thank two anonymous referees and the Associate Editor for their careful reading and suggestions for improvements, which has helped sharpen the presentation of this paper.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Andreas Søjmark.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

A. Søjmark was supported in part by a grant from the EPSRC [EP/L015811/1].

Appendix

Appendix

1.1 A.1 Some technical lemmas

Lemma A.1

Suppose\(\nu \)is a measure-valued process satisfying the conditions in Definition2.3and let\(g_{s}(x)=g(s,x,\nu _{s},L_{s})\), where\(g\)is any of\(b\), \(\mathfrak{b}\), \(\sigma \)or\(\sigma ^{2}\), with\(L_{s}=1-\nu _{s}(0,\infty )\). Define the error term

$$ \mathcal{E}_{t,\varepsilon }^{g}(x):=\left \langle \nu _{t},g_{t}( \cdot )G_{\varepsilon }(x,\cdot )\right \rangle -g_{t}(x)(\mathcal{T} _{\varepsilon }\nu _{t})(x). $$

Then we have

$$ \mathbb{E}\bigg[\int _{0}^{T}\Vert \mathcal{E}_{t,\varepsilon }^{g} \Vert _{L^{2}(\mathbb{R})}^{2}\,dt\bigg]\longrightarrow 0 \qquad \textit{as } \varepsilon \downarrow 0. $$

Proof

This follows by a simple modification of Lemma 8.1 in [35]. The only thing to note is that the crude bound \(|g_{t}(x)-g_{t}(y)| \leq 2\Vert g_{t}\Vert _{\infty }\) no longer holds, as \(g_{t}\) need not be bounded. However, the arguments from [35] still apply if we instead rely on \(|g_{t}(x)-g_{t}(y)|\leq \Vert \partial _{x}g_{t} \Vert _{\infty }|x-y|\). □

Lemma A.2

Suppose\(\bar{\nu }\)satisfies the whole-space analogues of conditions (iii), (iv) in Definition 2.3. Let\(g_{s}=g(s,\cdot ,\nu _{s},L_{s})\), where\(g\)is any of\(b\), \(\mathfrak{b}\), \(\sigma \)or\(\sigma ^{2}\), and define the error terms

$$\begin{aligned} \bar{\mathcal{E}}_{t,\varepsilon }^{g}(x) & :=\left \langle \bar{ \nu }_{t},g_{t}(\cdot )p_{\varepsilon }(x,\cdot )\right \rangle -g _{t}(x)\partial _{x}(\bar{\mathcal{T}}_{\varepsilon }\bar{\nu }_{t})(x)+ \partial _{x}g_{t}(x)\bar{\mathcal{H}}_{t,\varepsilon }^{g}(x), \\ \bar{\mathcal{H}}_{t,\varepsilon }^{g}(x) & :=\left \langle \bar{ \nu }_{t},(x-\cdot )\partial _{x}p_{\varepsilon }(x-\cdot )\right \rangle . \end{aligned}$$

Then for\(k=1,2\), we have

$$ \mathbb{E}\bigg[\int _{0}^{T}\Vert \bar{\mathcal{E}}_{t,\varepsilon } ^{g}\Vert _{L^{2}(\mathbb{R})}^{2k}\,dt\bigg] \longrightarrow 0,\quad \mathbb{E}\bigg[\int _{0}^{T}\Vert x\bar{\mathcal{E}}_{t,\varepsilon } ^{g}\Vert _{L^{2}(\mathbb{R})}^{2}\,dt\bigg]\longrightarrow 0, \qquad \textit{as }\varepsilon \downarrow 0. $$

Proof

This is immediate from the obvious modifications of [35, Lemma 8.2]. □

Lemma A.3

Fix any positive integer\(k \geq 1\). Let\((\nu _{t}^{{N}})\)be as defined in (2.1) and let\((\tilde{\nu }_{t})\)be any measure-valued process satisfying Definition2.3. Then for every\(a>0\), we have

$$ \mathbb{E}\bigg[\int _{0}^{T}\langle \tilde{\nu }_{t},x^{k}\mathbf{1} _{[\lambda ,\infty )}(x)\rangle\,dt\bigg]=o(\lambda ^{k}e^{-a\lambda }) \qquad \textit{as }\lambda \rightarrow \infty , $$
(A.1)

and likewise, it holds uniformly in\(N\geq 1\)and\(t\in [0,T]\)that

$$ \mathbb{E}[\langle \nu _{t}^{{N}},x^{k}\mathbf{1}_{[\lambda ,\infty )}(x) \rangle ] =o(\lambda ^{k}e^{-a\lambda }) \qquad \textit{as }\lambda \rightarrow \infty . $$
(A.2)

Proof

Fix an arbitrary \(a>0\) and let \(\varepsilon >0\) be given. By Corollary 3.4, there exists \(\lambda _{0}\geq 0\) such that

$$ \mathbb{E}[\nu _{t}^{{N}}(\lambda ,\infty )]\leq \varepsilon e^{-a \lambda }, \qquad \forall \lambda \geq \lambda _{0}, $$
(A.3)

uniformly in \(N\geq 1\) and \(t\in [0,T\text{]}\). Given any \(\lambda \geq \lambda _{0}\), we let \(\{ s_{0},s_{1},\ldots \} \) denote the partition of \([\lambda , \infty ) \) with \(s_{0} = \lambda \) and \(s_{i} - s_{i-1} = 1/a\). By (A.3) and monotone convergence, we get

$$\begin{aligned} \mathbb{E}[\langle \nu _{t}^{{N}},x\mathbf{1}_{[\lambda ,\infty )} \rangle ] & \leq \sum _{i=1}^{\infty }s_{i}\mathbb{E}[\nu _{t}^{{N}}(s _{i-1},\infty )]\leq \varepsilon \sum _{i=1}^{\infty }s_{i}e^{-as_{i-1}}= \varepsilon e\sum _{i=1}^{\infty }s_{i}e^{-as_{i}}. \end{aligned}$$

Now \(x\mapsto xe^{-ax}\) is decreasing for \(x\geq 1/a\); so taking \(\lambda _{0}^{\prime }:=\max \{ 1/a,\lambda _{0}\} \) and noting that \(s_{0}=\lambda \), it holds for all \(\lambda >\lambda _{0}^{\prime }\) that

$$\begin{aligned} \mathbb{E}[\langle \nu _{t}^{{N}},x\mathbf{1}_{[\lambda ,\infty )} \rangle ] & \leq \varepsilon e\int _{\lambda }^{\infty }xe^{-ax}\,dx= \varepsilon e\bigg(\frac{\lambda }{a}+\frac{1}{a^{2}}\bigg)e^{-a \lambda }. \end{aligned}$$

This proves (A.2) for \(k=1\), and the work for \(k\geq 2\) is analogous. The claim in (A.1) follows similarly by relying on the exponential tail property from Definition 2.3. □

Lemma A.4

Let\(\nu \)be a limit point from Theorem 2.4and\(\tilde{\nu }\)any measure-valued process satisfying Definition2.3. Then almost surely (as\(\lambda \rightarrow \infty \)),

$$ \lim _{\varepsilon \downarrow 0}\frac{\nu _{t}(0,\varepsilon )}{\varepsilon }=0, \qquad \nu _{t}(\lambda ,\infty )=O(e^{-\lambda }), \qquad \int _{0}^{T}\tilde{\nu }_{t}(\lambda ,\infty )\,dt=O(e^{-\lambda }). $$

Proof

For the first claim, recall from Proposition 4.3 that there exist \(\delta \in (0,1]\) and \(\beta >0\) such that \(\mathbb{E}[\nu _{t}(0,\varepsilon )]=t^{-\frac{ \delta }{2}}O(\varepsilon ^{1+\beta })\) as \(\varepsilon \rightarrow 0\). Using Markov’s inequality, we thus deduce that for any \(\theta >0\) and \(n\) sufficiently large,

$$ \mathbb{P}[n^{2/\beta }\nu _{t}(0,n^{-2/\beta })>\theta ]\leq \theta ^{-1}n ^{2/\beta }\mathbb{E}[\nu _{t}(0,n^{-2/\beta }]\leq Ct^{-\frac{ \delta }{2}}\theta ^{-1}n^{-2}. $$

Hence the Borel–Cantelli lemma gives \(\limsup _{n}n^{2/\beta }\nu _{t}(0,n ^{-2/\beta })=0\) almost surely. Now, given \(\varepsilon >0\), we have \((n+1)^{-2/\beta }<\varepsilon \leq n ^{-2/\beta }\) for some \(n\geq 1\); so we deduce that

$$ \limsup _{\varepsilon \downarrow 0}\frac{\nu _{t}(0,\varepsilon )}{ \varepsilon }\leq \limsup _{n\to \infty }\frac{\nu _{t}(0,n^{-2/\beta })}{(n+1)^{-2/ \beta }}\leq \limsup _{n\to \infty }\frac{\nu _{t}(0,n^{-2/\beta })}{n ^{-2/\beta }}\bigg(\frac{n+1}{n}\bigg)^{\frac{2}{\beta }}. $$

Since the latter is zero with probability 1, this proves the first claim. The two remaining results follow by analogous considerations for the tail probabilities, using the exponential tail properties from Proposition 4.3 and Definition 2.3. □

1.2 A.2 Proofs of Propositions 3.5 and 3.6

As an immediate consequence of Corollary 6.2, we have the following result.

Lemma A.5

Let\(\varLambda _{t}^{i,{N}}=|X_{t}^{i}|+{ \sum _{j=1}^{{N}}}a_{j}^{{N}}|X _{t}^{j}|\)as in (6.1). Then

$$ \mathbb{P}\bigg[\sup _{t\leq T}\varLambda _{t}^{i,{N}}\geq \lambda \bigg]=o(1) \qquad \textit{as }\lambda \rightarrow \infty , $$

uniformly in\(N\geq 1\)and\(i=1,\ldots ,N\).

Based on this lemma, we can adapt the arguments from [35, Sect. 4] to prove Propositions 3.5 and 3.6.

Proof of Proposition 3.5

Arguing as in the proof of [35, Proposition 4.6], for all \(a>0\) and \(\theta :=\frac{1}{2}(1-r)\), we have

$$ \mathbb{P}[L_{t+h}^{{N}}-L_{t}^{{N}}< \delta ,L_{t}^{{N}}< r]\leq \mathbb{P}[L_{t+h}^{{N}}-L_{t}^{{N}}< \delta ,\;\nu _{t}^{{N}}(0,a)> \theta ]+o(e^{-a}). $$

Let \(E:=\{ L_{t+h}^{{N}}-L_{t}^{{N}}<\delta ,\;\nu _{t}^{{N}}(0,a)> \theta \}\) and define a random index set ℐ by

$$ \mathcal{I}:=\{1\leq i\leq N:X_{t}^{i}< a,\;t< \tau _{i},\;a_{i}^{{N}}> \theta /2N\}. $$

Note that the particles with index in \(\{ 1\leq i\leq N:a_{i}^{{N}} \leq { \frac{1}{2}}\theta N^{-1}\} \) can contribute at most \({ \frac{1}{2}}\theta \) towards \(\nu _{t}^{{N}}(0,a)\). Recalling also that \(a_{i}^{{N}}\leq m/N\) for some \(m>0\), see (2.2), it follows that on the event \(E\), we must have \(|\mathcal{I}|>N\theta /2m\). Hence

$$ \mathbb{P}[E]\leq \sum _{\mathcal{I}_{0}:\left |\mathcal{I}_{0}\right |> \frac{\theta N}{2m}}\mathbb{P}[E\,|\,\mathcal{I}=\mathcal{I}_{0}] \mathbb{P}[\mathcal{I}=\mathcal{I}_{0}]. $$
(A.4)

Furthermore, since \(L_{t+h}^{{N}}-L_{t}^{{N}}<\delta \) on \(E\) while \(a_{i}^{{N}}>\theta /2N\) for \(i\in \mathcal{I}\), we deduce that

$$ \mathbb{P}[E\,|\,\mathcal{I}=\mathcal{I}_{0}]\leq \mathbb{P}\big[ \bigl|\bigl\{ i\in \mathcal{I}_{0}:\inf _{u\leq h}X_{t+u}^{i}\leq 0 \bigr\} \bigr|< 2\delta N/\theta \big| \mathcal{I}=\mathcal{I}_{0} \big]. $$
(A.5)

In order to estimate (A.5), we let \(Z_{t}^{i}=\varUpsilon _{t}(X_{t}^{i})\) as in Lemma 6.3 and recall that \((Z_{t}^{i})\) then satisfies

$$ dZ_{t}^{i}=\hat{b}_{t}^{i}\,dt+dB_{t}^{i}, \qquad |\hat{b}_{t}^{i}|\leq c(1+\varLambda _{t}^{i,{N}}), \qquad |Z_{t}^{i}|\leq c'(1+|X_{t}^{i}|). $$
(A.6)

Using that \(Z_{t}^{i}\leq 0\) if and only if \(X_{t}^{i}\leq 0\), the estimate in (A.5) implies

$$ \mathbb{P}[E\,|\,\mathcal{I}=\mathcal{I}_{0}] \leq \mathbb{P}\big[ \bigl|\bigl\{ i\in \mathcal{I}:\inf _{u\leq h}Z_{t+u}^{i}\leq 0\bigr\} \bigr|< 2\delta N/\theta \big|\mathcal{I}=\mathcal{I}_{0}\big]. $$

Moreover, it is immediate from (A.6) that

$$ Z_{t+u}^{i}\leq Z_{t}^{i}+c_{1}h(1+\tilde{\varLambda }_{h}^{i,{N}})+ \tilde{B}_{u}^{i}, \qquad \tilde{B}_{u}^{i}:=B_{t+u}^{i}-B_{t}^{i}, \qquad \tilde{\varLambda }_{h}^{i,{N}}:=\sup _{u\leq h}\varLambda _{t+u}^{i,{N}}. $$

In particular, if \(\inf _{u\leq h}\tilde{B}_{u}^{i}\leq -Z_{t}^{i}-ch(1+ \tilde{\varLambda }_{h}^{i,{N}})\), then \(\inf _{u\leq h}Z_{t+u}^{i}\leq 0\) and hence

$$ \mathbb{P}[E\,|\,\mathcal{I}=\mathcal{I}_{0}]\leq \mathbb{P}\bigg[ \bigl|\bigl\{ i\in \mathcal{I}:\inf _{u\leq h}\tilde{B}_{u}^{i}\leq -Z _{t}^{i}-ch(1+\tilde{\varLambda }_{h}^{i,{N}})\bigr\} \bigr|< { \frac{2 \delta }{\theta }} N \bigg|\mathcal{I}=\mathcal{I}_{0}\bigg]. $$

Using the bound \(|Z_{t}^{i}| \leq c' |X_{t}^{i}|\) from (A.6), the definition of ℐ implies that there is a constant \(\tilde{c}>0\) such that

$$ \mathbb{P}[E\,|\,\mathcal{I}=\mathcal{I}_{0}]\leq \mathbb{P}\bigg[ \bigl|\bigl\{ i\in \mathcal{I}:\inf _{u\leq h}\tilde{B}_{u}^{i}\leq - \tilde{c}a-\tilde{c}(1+\tilde{\varLambda }_{h}^{i,{N}})\bigr\} \bigr|< { \frac{2\delta }{\theta }} N \bigg| \mathcal{I}=\mathcal{I}_{0} \bigg]. $$

Recalling the fact that

$$ \tilde{B}_{u}^{i}=\int _{t}^{t+u}\rho _{s}\,dW_{s}^{0}+\int _{t}^{t+u}\sqrt{1- \rho _{s}^{2}}\,dW_{s}^{i}=:I_{u}+J_{u}^{i}, $$

we split the above probability on \(\{\sup _{u\leq h}|I_{u}|<\tilde{c}a \}\cap \{\sup _{u\leq h}\varLambda _{t+u}^{i,{N}}<a\}\) and its complement. In this way, we get

$$\begin{aligned} \mathbb{P}[E \,|\,\mathcal{I}=\mathcal{I}_{0}] &\leq \mathbb{P}\bigg[ \bigl|\bigl\{ i\in \mathcal{I}:\inf _{u\leq h}J_{u}^{i}\leq -3 \tilde{c}a-\tilde{c}\bigr\} \bigr|< { \frac{2\delta }{\theta }} N \bigg|\mathcal{I}=\mathcal{I}_{0}\bigg] \\ & \phantom{=:}+\mathbb{P}\bigg[\,\sup _{u\leq h}\left |I_{u}\right | \geq \tilde{c}a\bigg]+\mathbb{P}\bigg[\,\sup _{u\leq h}\varLambda _{t+u} ^{i,{N}}\geq a\bigg]. \end{aligned}$$

From Lemma A.5, we know that the last term is \(o(1)\) as \(a\rightarrow \infty \) uniformly in \(N\geq 1\). Similarly, \((I_{u})\) is a martingale, and hence the second term is also \(o(1)\) as \(a\rightarrow \infty \) by Doob’s maximal inequality. Concerning the first term, we can introduce a time-change to make each \(J^{i}\) an independent Brownian motion. Recalling (A.4), the proof can then be completed via a law of large numbers argument as in [35, Proposition 4.6] by choosing the free parameter \(a\) as a suitable function of \(\delta \). □

Proof of Proposition 3.6

Arguing as in [35, Proposition 4.7], it suffices to show that

$$ \lim _{\delta \rightarrow 0}\lim _{N\rightarrow \infty }\mathbb{P}[L _{t+\delta }^{{N}}-L_{t}^{{N}}\geq \eta ,\;\nu _{t}^{{N}}(0,\varepsilon )< \eta /2]=0. $$

Let \(E:=\{ L_{t+\delta }^{{N}}-L_{t}^{{N}}\geq \eta ,\;\nu _{t}^{{N}}(0, \varepsilon )<\eta /2\} \) and define the random index set

$$ \mathcal{I}:=\{ 1\leq i\leq N:X_{t}^{i}\geq \varepsilon \text{ or } t \geq \tau _{i}\} . $$

By (2.2), there exists \(m\) such that \(a_{i}^{{N}}\leq m/N\); so noting that \(\{ \nu _{t}^{{N}}(0,\varepsilon )<\eta /2\} \) is contained in \(\{ |\mathcal{I}|\geq \frac{N}{m}(1-\frac{ \eta }{2})\} \), we have

$$ \mathbb{P}[E]\leq \sum _{\mathcal{I}_{0}:|\mathcal{I}_{0}|\geq \frac{N}{m}(1-\frac{ \eta }{2})}\mathbb{P}[E\,|\,\mathcal{I}=\mathcal{I}_{0}]\mathbb{P}[ \mathcal{I}=\mathcal{I}_{0}]. $$

The conditional probabilities can be estimated by

$$\begin{aligned} \mathbb{P}[E \,|\,\mathcal{I}=\mathcal{I}_{0}] &\leq \mathbb{P}\bigg[ \bigl|\bigl\{ i\in \mathcal{I}_{0}:\inf _{s\in [t,t+\delta ]}X_{s} ^{i}\leq 0,\quad X_{t}^{i}\geq \varepsilon \bigr\} \bigr|\geq \frac{N \eta }{2m}\bigg|\mathcal{I}=\mathcal{I}_{0}\bigg] \\ & \leq \mathbb{P}\bigg[\bigl|\bigl\{ i\in \mathcal{I}_{0}: \inf _{s\in [t,t+\delta ]}(X_{s}^{i}-X_{t}^{i})\leq -\varepsilon \bigr\} \bigr|\geq \frac{N\eta }{2m}\bigg|\mathcal{I}=\mathcal{I}_{0} \bigg]. \end{aligned}$$
(A.7)

Using the scale transform \(\varUpsilon \) from Lemma 6.3, we let \(U_{s}^{i}:=\varUpsilon _{t+s}(X_{t+s}^{i})-\varUpsilon _{t+s}(X_{t}^{i})\) and note that as in Lemma 6.3,

$$\begin{aligned} dU_{s}^{i}= u_{s}^{i}\,ds+\rho _{t+s}\,dW_{t+s}^{0}+(1-\rho _{t+s}^{2})^{ \frac{1}{2}}\,dW_{t+s}^{i}=:u_{s}^{i}\,ds+dI_{s}+dJ_{s}^{i}, \end{aligned}$$
(A.8)

where the drift satisfies \(|u_{s}^{i}|\leq c_{1}(1+\varLambda _{t+s}^{i, {N}})\). By the construction of \(\varUpsilon \) and the boundedness from below of \(1/\sigma \), say \(1/\sigma \geq c_{2}>0\), it follows from (A.7) that

$$ \mathbb{P}[E \,|\,\mathcal{I}=\mathcal{I}_{0}]\leq \mathbb{P}\bigg[ \bigl|\bigl\{ i\in \mathcal{I}_{0}:\inf _{s\leq \delta }U_{s}^{i} \leq -c_{2}\varepsilon \bigr\} \bigr|\geq \frac{N\eta }{2m}\bigg| \mathcal{I}=\mathcal{I}_{0}\bigg]. $$

Given \(\delta >0\), fix a constant \(a=a(\delta )>0\) to be specified later. Using the decomposition (A.8) of \(U_{s}^{i}\) and the growth estimate for its drift \(u_{s}^{i}\), we see that on the event

$$ \bigg\{ \sup _{s\leq \delta }\varLambda _{t+s}^{i,{N}}< a/\delta \bigg\} \cap \bigg\{ \sup _{s\leq \delta }|I_{s}|< a\bigg\} , $$
(A.9)

if \(i\in \mathcal{I}_{0}\) is such that \(\inf _{s\leq \delta }U_{s}^{i} \leq -c_{2}\varepsilon \), then

$$\begin{aligned} \inf _{s\leq \delta }J_{s}^{i} & \leq -c_{2}\varepsilon +\delta c_{1} \bigg(1+\sup _{s\leq \delta }\varLambda _{t+s}^{i,{N}}\bigg)+ \sup _{s\leq \delta }\left |I_{s}\right | \\ & \leq -c_{2}\varepsilon +c_{1}\delta +(1+c_{1})a\:\leq -c_{3}( \varepsilon -\delta -a). \end{aligned}$$

Consequently, splitting up the desired probability on the event (A.9) and its complement, we get

$$\begin{aligned} \mathbb{P}[E \,|\,\mathcal{I}=\mathcal{I}_{0}] & \leq \mathbb{P} \bigg[\bigl|\bigl\{ i\in \mathcal{I}_{0}:\inf _{s\leq \delta }J_{s} ^{i}\leq -c_{3}(\varepsilon -\delta -a)\bigr\} \bigr|\geq \frac{N \eta }{2m}\bigg|\mathcal{I}=\mathcal{I}_{0}\bigg] \\ & \phantom{=:}+\mathbb{P}\bigg[\,\sup _{s\leq \delta }\varLambda _{t+s} ^{i,{N}}\geq a/\delta \bigg]+\mathbb{P}\bigg[\,\sup _{s\leq \delta } \left |I_{s}\right |\geq a\bigg]. \end{aligned}$$

By Doob’s maximal inequality, the last term is bounded by \(\delta a ^{-2}\); so choosing \(a=a(\delta )\) such that \(\delta a^{-2} \rightarrow 0\) as \(\delta \rightarrow 0\), we can ensure that it vanishes as \(\delta \rightarrow 0\). Moreover, requiring also that \(\delta /a \rightarrow 0\) as \(\delta \rightarrow 0\), it follows from Lemma A.5 that the second term vanishes (uniformly in \(N\)) as \(\delta \rightarrow 0\). These requirements are satisfied by \(a=a(\delta ):=\delta ^{1/2}\log \log (1/\delta )\), and with this choice for \(a\), the proof can now be finished by the same arguments as in [35, Proposition 4.7]. □

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hambly, B., Søjmark, A. An SPDE model for systemic risk with endogenous contagion. Finance Stoch 23, 535–594 (2019). https://doi.org/10.1007/s00780-019-00396-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s00780-019-00396-1

Keywords

Mathematics Subject Classification (2010)

JEL Classification

Navigation