Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Distress Propagation in Complex Networks: The Case of Non-Linear DebtRank

  • Marco Bardoscia ,

    mb@lims.ac.uk

    Affiliations Department of Banking and Finance, University of Zürich, Zürich, Switzerland, London Institute for Mathematical Sciences, London, United Kingdom

  • Fabio Caccioli,

    Affiliations Department of Computer Science, University College London, London, United Kingdom, Systemic Risk Centre, London School of Economics and Political Sciences, London, United Kingdom

  • Juan Ignacio Perotti,

    Affiliation IMT: Institute for Advanced Studies, Lucca, Italy

  • Gianna Vivaldo,

    Affiliation IMT: Institute for Advanced Studies, Lucca, Italy

  • Guido Caldarelli

    Affiliations IMT: Institute for Advanced Studies, Lucca, Italy, London Institute for Mathematical Sciences, London, United Kingdom, CNR-ISC: Institute for Complex Systems, Rome, Italy

Abstract

We consider a dynamical model of distress propagation on complex networks, which we apply to the study of financial contagion in networks of banks connected to each other by direct exposures. The model that we consider is an extension of the DebtRank algorithm, recently introduced in the literature. The mechanics of distress propagation is very simple: When a bank suffers a loss, distress propagates to its creditors, who in turn suffer losses, and so on. The original DebtRank assumes that losses are propagated linearly between connected banks. Here we relax this assumption and introduce a one-parameter family of non-linear propagation functions. As a case study, we apply this algorithm to a data-set of 183 European banks, and we study how the stability of the system depends on the non-linearity parameter under different stress-test scenarios. We find that the system is characterized by a transition between a regime where small shocks can be amplified and a regime where shocks do not propagate, and that the overall stability of the system increases between 2008 and 2013.

Introduction

Complex networks [13] have proved useful to describe systems characterised by pair-wise interactions. Properties of dynamical processes on networks can be strongly affected by the underlying topology [4]. Examples include spread of news [5], rumours [6], diseases [7], financial distress [8], random walkers travelling the graph [911], and avalanches [12, 13].

In these cases stylized models, despite their apparent simplicity, can give meaningful indications on the large scale dynamics of the system [7], also helping to shed light on the importance of the network topology [14]. For example, models of epidemic contagion (such as SIS or SIR [15]) display dramatically different behaviors depending weather they take place on regular lattices or on complex networks. Similarly, also the spread of distress [1619] in financial networks is deeply dependent on the pattern of connections among financial institutions. In particular, it is not clear yet if a single topology can be considered robust with respect to different types of shocks [20] or not [21, 22].

Financial institutions are strongly interconnected in a variety of ways (e.g. ownership relationships [23, 24], common asset holdings [2527], trading of derivatives [28], possible arbitrage opportunities to exploit [29]) through which distress can propagate and lead to amplification phenomena, such as default cascades. Here we focus on a single layer of interconnectedness, namely that associated with interbank loans. To cope with fluctuations of liquidity, banks constantly lend money to each other, at different maturities. Hence, lenders are subject to counterparty risk, i.e. the risk that their borrowers could default and therefore not be able to fulfill their obligations. This, in turn, could lead to the default of lenders, resulting in a further wave of distress.

In the literature on financial contagion, a bank is represented by its balance sheet, consisting of assets with a positive economic value (such as loans, derivatives, stocks, bonds, real estate) and of liabilities with a negative economic value (such as customers’ deposits, debits). The balance sheet identity for bank i defines its equity as the difference between its total assets Ai and its total liabilities Li: Ei = AiLi. A bank with a negative equity would not be able to pay back its debtors, even assuming that it could sell all of its assets. Therefore, usually a negative equity is considered a good proxy for the default of a bank. An interbank loan extended by bank i to bank j is an asset for bank i and a liability for bank j. Hence, the relationship between a lender and a borrower is pair-wise in nature and a convenient way to represent it is by means of a directed weighted network [30, 31] in which edges of weight correspond to a loan of amount from bank i to bank j. We call all the other assets and liabilities external and we denote them with and respectively (see Fig 1).

thumbnail
Fig 1. Interbank network and balance sheet.

Sketch of a portion of the interbank network and stylised representation of the balance sheet of a bank, with interbank assets and liabilities highlighted. The difference between assets and liabilities is the equity. A negative equity is usually considered a good proxy for the default of a bank.

https://doi.org/10.1371/journal.pone.0163825.g001

The study of the interbank network has attracted considerable attention, also for its practical importance. Two widely recognized algorithms to quantify losses due to financial contagion are the Furfine algorithm [12, 32] and DebtRank [3339]. The former is essentially a threshold model according to which a bank propagates distress to its creditors only after its default. In contrast, DebtRank was introduced precisely to account for shock propagations occurring also in absence of default. To this end, relative losses in the equity of a borrower translate into the same relative devaluation of interbank assets of the corresponding lender. Those two mechanisms represent two extremes. On one hand, the Furfine algorithm is likely to underestimate the build-up of systemic risk. On the other hand, in DebtRank even tiny variations in the equity (as those deriving from daily market fluctuations) have a sizeable impact on the value of interbank assets. As a realistic scenario is likely to lie in-between those two extremes, in this paper we propose a model that interpolates between them and use it to perform stress tests to the European banking system. We will refer to the introduced model as non-linear DebtRank.

The paper is organized as follows: In the section Results we specify the model used and present a detailed characterization of its behavior within the context of a case study. In the section Discussion we discuss the main implications of our results, also from a policy-making perspective, and point out some limitations of our approach. We refer the reader interested in the details about the data used and a derivation of the algorithm to the section Materials and Methods.

Results

We perform stress tests on N = 183 publicly traded European banks using data from their balance sheets for the years from 2008 to 2013 (see the section Materials and Methods: Data for a detailed description of data). Since data on bilateral exposures are not publicly available, we employ a reconstruction technique to infer plausible values [40, 41] and sample for each year 100 instances of interbank networks for given values of connectivity p, defined as the number of reconstructed edges divided by the number of possible edges (N(N − 1)) (see the section Materials and Methods: Data for more details about the reconstruction of data).

A stress test consists in applying an initial exogenous shock to the system and to measure its response in terms of the resulting equity losses. From the point of view of risk management, the relevant quantity is the relative equity loss of bank i at time t: (1) The corresponding quantity at the aggregate level is the total relative equity loss: (2) In the context of financial contagion, the initial exogenous shock amounts to a relative devaluation of external assets, which corresponds to setting the initial condition hi(1) (see the section Materials and Methods: Model dynamics). If banks hold no interbank asset they will not incur any additional equity loss and their final equity losses will be equal only to the devaluation of external assets. In contrast, let us suppose that bank i extends a loan to bank j and therefore that it holds the corresponding interbank asset . If bank j suffers an equity loss, its probability of default will increase and thus the probability that it will be able to fully pay back its loan will decrease. Bank i will account for the possibility of not being fully paid back by reducing the value of the interbank asset in its balance sheet (in financial jargon the interbank asset will be marked-to-market). However, as the value of the interbank asset decreases, also the equity of bank i will decrease by the same amount. Lenders of bank i will now re-evaluate their interbank assets towards bank i, using the same mechanism used by bank i to re-evaluate its interbank asset towards bank j. As a consequence, the process of re-evaluation of interbank assets and equities proceeds recursively from borrowers to lenders, until convergence.

As we show in the section Materials and Methods: Model dynamics, the dynamic equation for the relative equity loss that includes the aforementioned recursive re-evaluation of interbank assets is the following: (3) where pD(t) is the probability of default of bank j at time t and is the interbank leverage matrix. In order to completely define the iterative map Eq (3), we need to establish a relationship between and hj(t). In the case of Furfine’s algorithm, the probability of default is equal to one only if the equity is smaller than or equal to zero, and it is equal to zero otherwise, while in the linear DebtRank . Here we take a practical approach by introducing the following functional form: (4) in which α is a free parameter. Such approach, albeit “phenomenological” in the sense that Eq (4) is intended to provide an effective description of how shocks propagate, has several merits. It depends on a single parameter, which has a clear financial meaning: it is the inverse of the typical relative equity loss after which banks start to propagate distress to their creditors. Thus, 1/α can be interpreted as a soft threshold: losses smaller that 1/α have negligible impact on the probability of the default of banks. The value α should be therefore tuned by practitioners and regulators so that the probabilities of default can be calibrated (in principle also diversely) for each bank. Estimating it is beyond the scope of this paper. However, we will perform a sensitivity analysis by exploring the behaviour of the model as a function of α. Moreover, Eq (4) leads to probabilities of defaults that are convex with respect to the equity (see Fig 2), which is their expected behaviour. In fact small fluctuations in the equity of borrowers reasonably trigger tiny variations in the probability of the default of lenders, while when borrowers experience larger losses the marginal effect on lenders can be dramatic. Finally, Eq (4) easily allows to interpolate between two of the mostly widely used contagion models: the linear DebtRank, which is recovered for α = 0, and the Furfine algorithm, which is recovered for α → ∞. In Fig 2 we plot for these cases the probabilities of default as a function of the relative equity loss. The importance of the interbank leverage matrix in the context of distress propagation has been already highlighted by [36, 39, 42], pointing out that the stability of the system is determined by its largest eigenvalue. Also in our case, using the general framework discussed in [43] it is easy to show that the system is stable if λmax eα < 1, i.e. if α > log λmax, where λmax is the largest eigenvalue of the interbank leverage matrix Λ.

thumbnail
Fig 2. Comparison of different algorithms.

Probability of default as a function of the relative equity loss hj(t) for different algorithms. The non-linear DebtRank interpolates between the Furfine algorithm and the non-linear DebtRank.

https://doi.org/10.1371/journal.pone.0163825.g002

Our stress test consists in assuming a devaluation of external assets by a factor xshock for a fraction of banks equal to pshock. All presented results are averaged both over the network samples and over the set of initially shocked banks (10 realizations of such set for each network in the sample). Therefore, independently for each realization, we proceed to shock pshockN banks. Each shocked bank suffers an initial loss in the external assets: (see the section Materials and Methods: Model dynamics).

A first analysis is focused on the most critical year: 2008. In Fig 3 we show the effects of the initial shock, as its propagation through the network unravels in time. More in details, Fig 3 shows S(t), the fraction of stressed banks, D(t), the fraction of defaulted banks, and H(t), the total relative equity loss experienced by the system. Stressed banks are those which have experienced equity losses, but have not defaulted yet. Defaulted banks at time t are those for which hi(t) = 1. From left to right, plots correspond to α = 0, α = 1, and α = 2.

thumbnail
Fig 3. Unraveling of shocks propagation over time.

Fraction of stressed banks S(t) (blue line), fraction of defaulted banks D(t) (red line), and H(t) (violet line), total relative equity loss experienced by the system as a function of the time t over which shocks propagate. Banks experience a shock in the external assets, which suffer a relative loss equal to xshock = 0.5%. All points are averaged over a sample of 100 reconstructed networks with connectivity p = 0.05 and compatible with 2008 balance sheets, and over 10 realisations of the shock in which each bank is shocked with probability pshock = 0.05. Error bars span three standard errors. α = 0 in panel A and the algorithm reduces to the linear DebtRank, while α = 1 in panel B, and α = 2 in panel C. We see that the dynamics unravels within a few time steps in the panels A and C, while it takes considerably more time steps in panel B.

https://doi.org/10.1371/journal.pone.0163825.g003

The qualitative behaviour of both stressed and defaulted banks is shared by all panels. Stressed banks sharply increase in the first time steps and decrease afterwards, as defaults start to occur. This is consistent with the fact that stress propagates even in absence of defaults. However, a clear dependence from α emerges. The most striking feature is that the time scale over which the system reaches its steady state is a non-monotonous function of α. In fact, in both panels A and C convergence is reached before the first 20 time steps, while in panel B the dynamics is slower. This phenomenon can be intuitively understood as follows: α is related to a soft threshold in the value of the relative equity loss that a bank needs to attain before it can propagate a shock. Such threshold is zero in the linear case and approaches one in the strongly non-linear regime (α ≫ 1). In the first case shocks are quickly propagated through the network, while in the second case shocks are easily dampened. In the intermediate regime the build-up of the stress happens gradually. Nevertheless, the total relative equity loss can still reach values comparable to those of the linear case (see panels A and B).

Next, we present a comprehensive characterization of the steady state. Fig 4 shows the surface plots of H, the total relative equity loss in the steady state, as a function of α and xshock, for network connectivity p = 0.05, and for different choices of pshock. n order to meaningfully compare the results for different values of pshock, we tune the range of xshock spanned so that the ranges of total shock pshockxshock affecting the system are equal across all the cases considered.

thumbnail
Fig 4. Equity loss in 2008.

Surface plot of H, the total relative equity loss in the steady state, as a function of the size of the shock suffered by external assets of banks xshock and of the parameter α, which tunes the non-linearity of the algorithm. All points are averaged over a sample of 100 reconstructed networks with connectivity p = 0.05, and compatible with 2008 balance sheets, and over 10 realization of the shock in which each bank is shocked with probability pshock = 0.05 (panel A) and pshock = 1 (panel B). Note that the range of the total size of the shock pshockxshock is the same for both panels. As α increases, the propagation of the shock is dampened, resulting in smaller losses, and in two different regimes, whose separation is especially evident for pshock = 1, i.e. when all banks are shocked.

https://doi.org/10.1371/journal.pone.0163825.g004

Let us start from focusing on the behaviour of H as a function of α for fixed values of xshock. From the panel A of Fig 4 we can clearly see that, apart from the numerical fluctuations due to the finite number of realizations, for any given value of xshock, H decreases monotonically with α. Starting from large values of α and moving towards smaller values, H reaches a plateaux in which most of the equity in the system has been lost (H > 0.9). As expected, focusing instead on the behaviour of H as a function of xshock, for fixed values of α, H increases with xshock. Overall we can detect two different regimes, one in which (almost) all the equity in the system has been lost (for small values of α) and one in which (almost) no equity has been lost (for large values of α and small values of xshock). For smaller and smaller values of xshock the crossover between such two regimes is sharper and sharper. Finally, in the limit xshock → 0 the system displays a transition between a stable regime in which no losses occur, and an unstable regime in which also infinitesimal shocks can lead to large total relative losses. The presence of such transition can be better appreciated from panel B of Fig 4, where we present the results in the case in which all banks suffer the same initial shock (i.e. pshock = 1). This can be interpreted as a shock to a risk factor common to all banks, such as a sudden change in interest rates or similar to that experienced during a major macroeconomic downturn. We note here that we have performed simulations for different values of the connectivity parameter p ranging between p = 0.05 to p = 1 (fully connected network). Interestingly, we observed that systems with very different connectivities behave in a similar way. A possible explanation is that, due to the reconstruction technique used (see the section Materials and Methods: Data), p = 0.05 is already enough to connect systemically important banks.

Finally, in Fig 5 we adopt the same setting as in Fig 4 for different years. Overall, we observe that H markedly decreases from 2008 to 2013. It clearly emerges that the system was more prone to amplify shocks in 2008, when a region in parameters space in which H ≃ 1 exists. This is consistent with the intuition that banks in 2008 were more fragile.

thumbnail
Fig 5. Equity loss in 2010 and 2013.

Analogous of Fig 4, but for the years 2010 and 2013. Here the connectivity of the reconstructed networks is p = 0.05 and the fraction of shocked banks on average is pshock = 0.05. Going from 2008 (see Fig 4, panel A) to 2013 H is less and less sensitive to changes in α and it is generally smaller, meaning that banks are more and more robust.

https://doi.org/10.1371/journal.pone.0163825.g005

Discussion

In the present study, a general shock propagation mechanism is applied to an interbank network of 183 publicly traded European banks. With probability pshock each bank is subject to an initial shock consisting in the devaluation of their external assets by a factor xshock. The system reaction to shocks is measured in terms of the total relative equity loss, which takes into account the contribution of each bank to the relative equity loss of the system. The dynamics of shock propagation that we consider (non-linear DebtRank) interpolates through the parameter α between two stress test algorithms: the Furfine algorithm and the linear DebtRank.

We notice that the propagation of shocks strongly depends on the parameter α. In particular, in all stress scenarios that we have considered we observe a crossover between a regime of large losses (for small α), in which potentially all banks could default, and a regime of small losses (for large α), in which most banks survive the shock. The width of the intermediate region shrinks as the fraction of banks affected by the shock approaches one.

The model also shows that the interbank network was significantly more fragile in 2008, when the financial crisis took place, than in the subsequent years. This observation holds qualitatively for all values of model parameters and connectivities explored and it is also in agreement with other empirical studies [44].

In addition to the properties of the steady state we have also looked into the dynamics of quantities such as the number of stressed and defaulted banks, whose behavior highlights the existence of different time scales, depending on the model parameters. For instance, we observe that the time needed to reach the steady state is a non-monotonic function of α: in certain cases the shock produces a slow drive of the interbank network towards its collapse, while in other cases the crash occurs immediately after the shock.

Clearly, establishing a coherent mapping between probability of defaults and changes in equity opens several possible directions for future research. Obviously, calibrating α, possibly extracting a different value for each bank, would represent a major achievement. Beyond the “phenomenological” approach adopted here, one could try to derive such relationship in the context of standard financial risk management theory. Moreover, here we have limited our analysis to direct exposures due interbank landing. A proper assessment of systemic risk should account for additional types of interconnectedness, such as that associated with overlapping portfolios, exchange of derivates, and ownership structure. Hence, another future extension of the model could be based on a multilayer network that incorporates those effects. Complex interactions across different layers could lead to non-trivial amplification phenomena [45]. Finally, studying the stability of the system, which is becoming increasingly possible due to recent works [43, 46], could allow regulators to develop faster qualitative stress-testing frameworks to continuously monitor systemic risk.

Materials and Methods

Model dynamics

In order to derive Eq (3) we start from the balance sheet identity introduced in the section Introduction in which we distinguish between external and interbank assets and liabilities: (5) where, as customary in the literature on financial contagion, assets depend explicitly on time, while liabilities do not. In fact, while external assets are subject to exogenous fluctuations, the value of interbank assets will change accordingly (marked-to-market), depending on the probability of default of the borrower j. However, the fact that bank i reassesses the value of its interbank claim towards bank j does not change the value of the debt that bank j owes to bank i: hence interbank liabilities (and analogously external liabilities) always keep their face value (i.e. the value initially present in the balance sheet), and therefore do not depend on time. As a consequence, using Eqs (1) and (5) we can readily compute: (6) If the probability of default of bank j is equal to zero, obviously its creditor bank i can expect to fully recover the face value of its loan. In contrast, if the probability of default of bank j is equal to one, bank j has defaulted and its creditor bank i can expect to recover only a fraction of its loan, where Rij ∈ [0, 1]. A natural way to interpolate between these two extreme cases is to identify the value of interbank assets at time t + 1 with their expected value at time t: (7) By making the conservative assumption [18, 47] that creditors do not actually recover any amount of money in case of default of their borrowers, we have: (8) which implies that bank i updates the value of its interbank claims towards bank j such that it is equal to their face value if the probability of default of the borrower bank j is zero and it decreases proportionally to otherwise. By plugging Eq (8) into Eq (6) and using the definition of the interbank leverage matrix [36, 39]: (9) we can immediately compute: (10) In our stress test scenario we will initially shock external assets by a relative amount xshock, i.e. , such that . However, after the initial shock external assets do not change anymore and the evolution of the relative equity losses is entirely due to the re-assessment of the value interbank assets. As a consequence, the first term in the right-hand side of Eq (10) is equal to zero, for t > 1. Finally, when the equity of bank i becomes equal to zero the bank defaults and is not able neither to further propagate shocks nor to sustain any additional losses. Hence, the maximum value attainable by the relative equity losses is one, which leads to Eq (3).

From Eq (3), we see that the results will depend on the relationship that we assume between the relative loss in equity of a bank and its probability of default. In the Furfine algorithm we have that: (11) while in the linear DebtRank: (12) In the non-linear DebtRank we interpolate between this two extreme cases by means of the parameter α (see Eq (4)): for α = 0 we recover the linear DebtRank, while for α → ∞ we recover the Furfine algorithm.

Data

The original source for raw data about balance sheets of banks is the Bureau van Dijk Bankscope database, from which we extract data for a subset of 183 among the largest European banks. In particular we use data about the total interbank assets (liabilities) () and we compute external assets (liabilities) as the difference between total assets (liabilities) and interbank assets (liabilities): (). The same data have already been used in [36, 39]. See [36] for all the details about the handling of missing data. Even though European banks do not necessarily constitute an isolated system [48], they are weakly correlated to American and Asian banks [49]. Moreover, banks outside Europe fall under different regulating authorities, and therefore they might be subject to different conditions (e.g. during bailouts). The analysis of Ref. [49] makes it clear that European banks are also connected to financial institutions other than banks. To account for the additional channels of contagion due to the interaction between financial institutions of different types, one could extend our approach to the case of multilayer networks, as suggested in the section Discussion.

As already pointed out, the balance sheet of bank i contains only data about its total interbank lending and borrowing, i.e. the values of and . As a consequence, the full matrix needs to be reconstructed by making some assumptions. Here we proceed as in [39] and make the reasonable assumption, successfully tested in [40, 41], that the probability that bank i extended a loan to bank j is proportional to , the total amount of interbank lending of bank i and to , the total amount of interbank borrowing of bank j. The fitness model [50] allows us to compute the values of the probabilities {pij} that an edge ij exists for a given value of connectivity p. We then use {pij} to build a sample of 100 directed un-weighted networks for each year. For each network in the sample we assign weights using the iterative RAS algorithm [51]. See [39] for a full account of the procedure.

Acknowledgments

All authors thank Marco D’Errico for sharing the data and for fruitful discussions and Stefano Battiston for the numerous exchanges of views. MB, JIP, GV, and GC acknowledge support from FP7-ICT project MULTIPLEX nr. 317532, FP7-ICT project SIMPOL nr. 610704, and Horizon 2020 project DOLFINS nr. 640772. FC acknowledges support of the Economic and Social Research Council (ESRC) in funding the Systemic Risk Centre (ES/K002309/1). GC acknowledges additional support from Horizon 2020 projects SoBIGData nr. 654024 and CoEGSS nr. 676547.

Author Contributions

  1. Conceptualization: MB FC JIP GV GC.
  2. Data curation: JIP GV.
  3. Formal analysis: MB FC JIP GV GC.
  4. Investigation: MB FC JIP GV GC.
  5. Methodology: MB FC JIP GV GC.
  6. Software: JIP GV.
  7. Validation: MB FC JIP GV GC.
  8. Visualization: MB FC JIP GV GC.
  9. Writing – original draft: MB FC JIP GV GC.
  10. Writing – review & editing: MB FC JIP GV GC.

References

  1. 1. Albert R, Barabási AL. Statistical Mechanics of Complex Networks. Reviews of Modern Physics. 2002;74(1):47–97.
  2. 2. Caldarelli G. Scale-Free Networks: Complex Webs in Nature and Technology. Oxford University Press; 2007. https://doi.org/10.1093/acprof:oso/9780199211517.001.0001
  3. 3. Barabási AL. The network takeover. Nature Physics. 2011;8(1):14–16.
  4. 4. Barrat A, Barthélemy M, Vespignani A. Dynamical processes on complex networks. Cambridge University Press; 2008. https://doi.org/10.1017/CBO9780511791383
  5. 5. Leskovec J, Backstrom L, Kleinberg J. Meme-tracking and the Dynamics of the News Cycle. In: Proceedings of the 15th ACM SIGKDD International Conference on Knowledge Discovery and Data Mining; 2009. p. 497–506.
  6. 6. Nekovee M, Moreno Y, Bianconi G, Marsili M. Theory of rumour spreading in complex social networks. Physica A. 2007;374(1):457–470.
  7. 7. Pastor-Satorras R, Castellano C, Van Mieghem P, Vespignani A. Epidemic processes in complex networks. Reviews of Modern Physics. 2015;87(3):925–979.
  8. 8. Schweitzer F, Fagiolo G, Sornette D, Vega-Redondo F, Vespignani A, White DR. Economic Networks: The New Challenges. Science. 2009;325(5939):422–425. pmid:19628858
  9. 9. Motter AAE, Zhou C, Kurths J. Network synchronization, diffusion, and the paradox of heterogeneity. Physical Review E. 2005;71(1):1–10. pmid:15697667
  10. 10. Ehrhardt G, Marsili M, Vega-Redondo F. Diffusion and growth in an evolving network. International Journal of Game Theory. 2006;34:383–397.
  11. 11. Scholtes I, Wider N, Pfitzner R, Garas A, Tessone CJ, Schweitzer F. Causality-driven slow-down and speed-up of diffusion in non-Markovian temporal networks. Nature Communications. 2014;5:5024. pmid:25248462
  12. 12. Watts DJ. A simple model of global cascades on random networks. Proceedings of the National Academy of Sciences. 2002;99(9):5766–5771. pmid:16578874
  13. 13. Lee KM, Yang JS, Kim G, Lee J, Goh KI, Kim Im. Impact of the topology of global macroeconomic network on the spreading of economic crises. PloS One. 2011;6(3):e18443. pmid:21483794
  14. 14. Van den Broeck W, Gioannini C, Gonçalves B, Quaggiotto M, Colizza V, Vespignani A. The GLEaMviz computational tool, a publicly available software to explore realistic epidemic spreading scenarios at the global scale. BMC Infectious Diseases. 2011;11:37. pmid:21288355
  15. 15. Allen LJS. Some discrete-time SI, SIR, and SIS epidemic models. Mathematical Biosciences. 1994;124(1):83–105. pmid:7827425
  16. 16. Iori G, Jafarey S, Padilla FG. Systemic risk on the interbank market. Journal of Economic Behavior & Organization. 2006;61(4):525–542.
  17. 17. Nier E, Yang J, Yorulmazer T, Alentorn A. Network models and financial stability. Journal of Economic Dynamics and Control. 2007;31(6):2033–2060.
  18. 18. Gai P, Kapadia S. Contagion in financial networks. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences. 2010;466(2120):2401–2423.
  19. 19. Battiston S, Delli Gatti D, Gallegati M, Greenwald B, Stiglitz JE. Liaisons dangereuses: Increasing connectivity, risk sharing, and systemic risk. Journal of Economic Dynamics and Control. 2012;36(8):1121–1141.
  20. 20. Tran HAQ, Namatame A. Mitigating Cascading Failure with Adaptive Networking. New Mathematics and Natural Computation. 2015;11(02):151–163.
  21. 21. Roukny T, Bersini H, Pirotte H, Caldarelli G, Battiston S. Default cascades in complex networks: Topology and systemic risk. Scientific Reports. 2013;3:2759. pmid:24067913
  22. 22. Caccioli F, Catanach T, Farmer JD. Heterogeneity, correlations and financial contagion. Advances in Complex Systems. 2012;15(2):1250058.
  23. 23. Elliott M, Golub B, Jackson MO. Financial networks and contagion; 2014. Available from: http://ssrn.com/abstract=2175056.
  24. 24. Vitali S, Glattfelder JB, Battiston S. The Network of Global Corporate Control. PLoS One. 2011 10;6(10):e25995. pmid:22046252
  25. 25. Huang X, Vodenska I, Havlin S, Stanley HE. Cascading Failures in Bi-partite Graphs: Model for Systemic Risk Propagation. Scientific Reports. 2013;3:1219. pmid:23386974
  26. 26. Caccioli F, Shrestha M, Moore C, Farmer JD. Stability analysis of financial contagion due to overlapping portfolios. Journal of Banking & Finance. 2014;46:233–245.
  27. 27. Guo W, Minca A, Wang L. The Topology of Overlapping Portfolio Networks; 2015. Available from: http://ssrn.com/abstract=2619514.
  28. 28. Heise S, Kühn R. Derivatives and credit contagion in interconnected networks. The European Physical Journal B. 2012;85(4).
  29. 29. Bardoscia M, Livan G, Marsili M. Financial instability from local market measures. Journal of Statistical Mechanics: Theory and Experiment. 2012;2012(08):P08017.
  30. 30. De Masi G, Iori G, Caldarelli G. Fitness model for the Italian interbank money market. Physical Review E. 2006;74(6):066112. pmid:17280126
  31. 31. Cont R, Moussa A, Santos EB. Network Structure and Systemic Risk in Banking Systems. In: Fouque JP, J L, editors. Handbook of Systemic Risk. Cambridge University Press; 2010. p. 327–368.
  32. 32. Furfine CH. Interbank Exposures: Quantifying the Risk of Contagion. Journal of Money, Credit & Banking. 2003;35(1):111–129.
  33. 33. Battiston S, Puliga M, Tasca P, Caldarelli G, Kaushik R, Tasca P, et al. DebtRank: too central to fail? Financial networks, the FED and systemic risk. Scientific Reports. 2012;2:541. pmid:22870377
  34. 34. Di Iasio G, Battiston S, Luigi I, Federico P. Capital and contagion in financial networks; 2013. Available from: http://mpra.ub.uni-muenchen.de/52141/.
  35. 35. Thurner S, Poledna S. DebtRank-transparency: controlling systemic risk in financial networks. Scientific Reports. 2013;3:1888. pmid:23712454
  36. 36. Battiston S, Caldarelli G, Gurciullo S. Leveraging the network: a stress-test framework based on DebtRank; 2015. Available from: https://arxiv.org/abs/1503.00621.
  37. 37. Battiston S, D’Errico M, Gurciullo S. DebtRank and the network of leverage. Journal of Alternative Investments. 2016;18(4):68–81.
  38. 38. Battiston S, Caldarelli G, D’Errico M. The financial system as a nexus of interconnected networks. In: Garas A, editor. Interconnected networks. Springer; 2016. p. 195–229.
  39. 39. Bardoscia M, Battiston S, Caccioli F, Caldarelli G. DebtRank: A Microscopic Foundation for Shock Propagation. PloS One. 2015;10(6):e0134888. pmid:26231004
  40. 40. Cimini G, Squartini T, Garlaschelli D, Gabrielli A. Systemic risk analysis in reconstructed economic and financial networks. Scientific Reports. 2016;5:15758. pmid:26507849
  41. 41. Cimini G, Squartini T, Gabrielli A, Garlaschelli D. Estimating topological properties of weighted networks from limited information. Physical Review E. 2015;92:040802. pmid:26565153
  42. 42. Markose S, Giansante S, Shaghaghi AR. ‘Too interconnected to fail’ financial network of US CDS market: Topological fragility and systemic risk. Journal of Economic Behavior and Organization. 2012;83(3):627–646.
  43. 43. Bardoscia M, Battiston S, Caccioli F, Caldarelli G. Pathways towards instability in financial networks; 2016. Available from: http://arxiv.org/abs/1602.05883.
  44. 44. Memmel C, Sachs A. Contagion in the interbank market and its determinants. Journal of Financial Stability. 2013;9(1):46–54.
  45. 45. Caccioli F, Farmer JD, Foti N, Rockmore D. Overlapping portfolios, contagion, and financial stability. Journal of Economic Dynamics and Control. 2015;51:50–63.
  46. 46. Tan JPL. Symmetric and Asymmetric Tendencies in Stable Complex Systems. Scientific Reports. 2016;6:31762. pmid:27545722
  47. 47. Markose, Sheri M. Systemic Risk from Global Financial Derivatives: A Network Analysis of Contagion and Its Mitigation with Super-Spreader Tax; 2012. Available from: https://www.imf.org/external/pubs/ft/wp/2012/wp12282.pdf.
  48. 48. Kubelec C, Sá F. The geographical composition of national external balance sheets: 1980–2005. Bank of England Working Paper; 2010. 384.
  49. 49. Sandoval L. Structure of a global network of financial companies based on transfer entropy. Entropy. 2014;16(8):4443–4482.
  50. 50. Musmeci N, Battiston S, Caldarelli G, Puliga M, Gabrielli A. Bootstrapping topological properties and systemic risk of complex networks using the fitness model. Journal of Statistical Physics. 2013;151(3–4):720–734.
  51. 51. Upper C, Worms A. Estimating Bilateral Exposures in the German Interbank Market: Is there a Danger of Contagion? European Economic Review. 2004;48(4):827–849.