# The price of complexity in financial networks

Edited by Jose A. Scheinkman, Columbia University, New York, NY, and approved July 6, 2016 (received for review November 21, 2015)

## Significance

Estimating systemic risk in networks of financial institutions represents, today, a major challenge in both science and financial policy making. This work shows how the increasing complexity of the network of contracts among institutions comes with the price of increasing inaccuracy in the estimation of systemic risk. The paper offers a quantitative method to estimate systemic risk and its accuracy.

## Abstract

Financial institutions form multilayer networks by engaging in contracts with each other and by holding exposures to common assets. As a result, the default probability of one institution depends on the default probability of all of the other institutions in the network. Here, we show how small errors on the knowledge of the network of contracts can lead to large errors in the probability of systemic defaults. From the point of view of financial regulators, our findings show that the complexity of financial networks may decrease the ability to mitigate systemic risk, and thus it may increase the social cost of financial crises.

### Sign up for PNAS alerts.

Get alerts for new articles, or get an alert when an article is cited.

Several years after the beginning of the so-called Great Recession, regulators warn that we still do not have a satisfactory framework to deal with too-big-to-fail institutions and with systemic events of distress in the financial system (1, 2). The topic is of general societal and scientific interest because assisting financial institutions in recent years has come with a social cost equivalent to a few percent points of gross domestic product for Organization for Economic Cooperation and Development countries (3). In particular, poor estimates of systemic risk are socially costly because regulators and bank managers end up keeping aside either insufficient or redundant buffers (4–7). One of the difficulties is that financial institutions are connected in multilayer networks both directly, via contracts among each other (loans, bonds, repurchasing agreement, derivatives, etc.), and indirectly, via exposures to common assets (8–15). The default probability of one institution depends therefore on the default probability of all of the other institutions in the network. The resulting complexity of the financial system is a potential source of information asymmetries, collective moral hazard, and increased systemic risk (16–20) and, hence, requires deeper understanding. In particular, the determination of the probability of systemic events has remained an open problem so far (21–25). Here, we show that, in a network of financial contracts, the probability of systemic default can be very sensitive to errors on information about contracts as well as on information about the complexity of the network structure. Moreover, the sensitivity depends on the network structure due to a multiplicative interplay of errors along chains of lending. Under certain conditions, network structures with more numerous and longer chains of lending lead to a stronger amplification of errors. Although there may be some intuition for this effect in the literature, it has not been explained analytically or quantified before. Our work aids an understanding of how to contain errors in the estimation of the probability of systemic events.

## Results

We introduce a general model of a network of credit contracts among

*n*financial institutions (hereafter, “banks”). Contracts are over the counter (OTC), i.e., they are not mediated through a central counterparty. They are also collateralized, i.e., banks have to post collateral to receive a loan. In addition to making contracts with each other, banks also hold external assets, i.e., securities that are not issued by the banks in the system. These external assets are the only source of stochasticity in the model. An illustration of the case of three banks connected through credit contracts is sketched in Fig. 1 (see*Methods*for the details of the model).Fig. 1.

### Errors on Contract Characteristics.

We consider a first scenario, named here as “errors on contract characteristics,” in which it is known which assets and which counterparties each bank is exposed to, but there can be an error, for instance, in the recovery rate

*R*(i.e., the fraction of the face value of the loan that can be recovered after the default of a counterparty). We first study how the systemic default probability ${P}^{\text{sys}}$ varies as a function of parameter errors in the simplest case of two symmetric banks. We explore all pairs of parameter values in a range consistent with empirical evidence on interbank markets (8). In Fig. 2, each pair of curves with the same color represents the maximum and minimum value of the default probability as a function of the deviation in each parameter around a given point. Fig. 2,*Left*refers to the case in which the actual value (i.e., with zero error) of the default probability is 1. For instance, with an error on recovery rate*R*of 20%, one may think that the default probability is 0.4 when it is actually 1. Conversely, Fig. 2,*Right*refers to the case in which the actual value of the default probability is 0. For instance, with an error in the expected asset return*μ*of 20%, one may think that the default probability is 1 when it is actually 0. In both cases, the gap between the possible estimates increases rapidly with the error on the parameters, e.g., an error in parameters of 10% can lead to an error of 100% in the probability. Note that this result does not imply that the error in the default probability is always large compared with errors in the parameters. However, it shows that there exist cases in which the deviation can be very large. Next, to illustrate how the effects of errors on contract characteristics depend also on the structure of the underlying network of contracts, we focus on errors on a single parameter, namely the recovery rate*R*, and we consider three basic architectures with three nodes: a star, a chain, and a ring. Fig. 3 shows the sensitivity $\partial {P}^{\text{sys}}/\partial R$ of the default probability on the recovery rate*R*as a function of the ratio between interbank leverage*β*and the maximal loss $\mathit{\u03f5}\sigma $ that a bank can withstand on external assets. As we can see, the sensitivity is highest for the ring network architecture, followed by the chain and the star architectures. The intuition behind this result is that the systemic default probability in an interbank network depends on the multiplicative interplay among the parameters that matter for the default thresholds. In particular, the multiplication involves the banks located along chains of lending. As a result, more-numerous and longer chains lead to stronger amplification of the errors on the parameters. Indeed, as proved in*Methods*, $\partial {P}^{\text{sys}}/\partial R$ is a polynomial in the interbank leverage*β*, where the leading power depends on the presence of chains or cycles in the network. This result, which has not been reported so far in the literature, illustrates concretely the impact of financial complexity on the determination of systemic risk.Fig. 2.

Fig. 3.

### Errors in the Structure of the Contract Network.

Finally, we would like to investigate how the systemic default probability depends on the complexity of the network of contracts. To this end, we consider a second scenario, named here as “errors in network structure,” in which the information regarding how many contracts a bank has and with which counterparties may be incorrect. More precisely, we are interested in measuring the error in the systemic default probability when the arrangement of contracts (i.e., who trades with whom) is not known and the number of possible arrangements increases, subject to the constraint of a maximum number of possible contracts, referred to as “density cap” (

*Methods*). In Fig. 4, the green area corresponds to a lower expected return on external assets*μ*, and the blue area corresponds to a higher one. Both areas represent the range of possible values of the default probability resulting from all possible network configurations compatible with a given density cap. The tip of each arrow indicates the value of default probability for three specific network architectures, i.e., the star, the circle, and the complete graph. For instance, in the case of lower*μ*(green area), the complete graph yields values of default probability larger than the star architecture. In this context, a simple way to capture the complexity of the network is to count the number of possible network configurations for a given density cap. In Fig. 4, the red curve represents the entropy $\mathrm{\Sigma}\left({C}_{l}\right)$ of the space of possible network configurations (*Methods*) for increasing values of the density cap.Fig. 4.

## Discussion

The complexity of the financial system is a source of potential information asymmetries and collective moral hazard (16, 20). Indeed, if the default of some financial institutions has an impact on the system that is possibly very large but difficult to compute exactly, those institutions are more likely to enjoy a bailout with public funds; they tend therefore to count on being rescued in case of downturn and take more risk than they would otherwise. In other words, such moral hazard leads to a more fragile financial system and to an implicit public subsidy. Therefore, the estimation of the probability of individual and systemic events in a network of contracts is crucial to improve the stability of the financial system. However, the interdependence among asset values and probability of default of all institutions poses conceptual and computational challenges, and little progress has been made in this direction so far. Many previous works on systemic risk build on the approach

*à la*Eisenberg and Noe (8, 26−28) in which a clearing vector of payment and a recovery rate on defaulting banks’ assets is determined endogenously. However, the reason why the recovery rate can be determined endogenously in those works is that the valuation is carried out ex post (i.e., at the maturity of the contracts) and that, in case of default, there is a successful and immediate asset liquidation of the external assets with full recovery. In contrast, asset liquidation implies legal settlements that, in practice, take several months or years, whereas, in the short run (e.g., weeks), the recovery rate is likely to be significantly smaller than in the Eisenberg−Noe approach and to be sensitive to regulators’ interventions. For these reasons, it is important to study the impact of errors on the recovery rate. Moreover, and even more importantly, here we are interested in the valuation that can be carried out*ex ante*, i.e., before the maturity, given the information available. Therefore, we assume that banks lend to each other against a collateral and that the recovery rate on the loans to defaulted banks is smaller than 1 and exogenous (25). Finally, similar to refs. 9, 18, and 19, we intentionally leave aside, at this stage, the question of which configurations and parameters would arise from banks’ dynamic choice, because our method encompasses all possible configurations, including those arising from the uncoordinated individual investment strategies.On the one hand, the ability of a bank to make contracts with any other bank in the system increases its ability to diversify the risk. On the other hand, the resulting complexity comes with the price that “everybody knows less.” Indeed, either a higher complexity of the individual contracts or a higher complexity of the structure of contracts implies that market participants and regulators know less precisely the probability of individual and systemic default. Although there are individual incentives to be part of a complex financial network, this work shows quantitatively the existence of so-called “negative externalities,” which eventually translate into potential social costs.

More generally, our results show that higher interdependence on the credit market among banks decreases the ability to make estimates on default probabilities and hence to correctly price debt instruments. One possible approach to this problem could be to increase the complexity of the regulation in order to match the complexity of the financial market. However, here we have shown that there are intrinsic limitations in the accuracy of estimating the probability of default, which implies that increasing the complexity of regulation may not be an appropriate way to address market complexity (29). Hence, a potential tradeoff emerges between financial stability and market complexity.

## Methods

We consider a financial network with OTC contracts among

*n*banks, including secured credit contracts. We distinguish between contracts within the banking system itself (“interbank”) and contracts of banks on securities outside the banking system (“external”).### Model Timing.

The timing of the model is as follows. At time 1, banks raise funds and make investments in external and interbank assets. At time 2, the values of the external assets are shocked and updated. Although the shock distribution is known at time 1, shocks are only observed at time 2. At time 2, the interbank contracts mature, and their value is also updated depending on the shocks that have occurred. Fig. 1 sketches the structure of our model. For each bank

*i*, the main quantities are the following.### Assets and Liabilities.

Assets and liabilities of

*i*on the external markets are denoted as ${a}_{i}^{E}$ and ${\ell}_{i}^{E}$. Assets and liabilities of*i*on the interbank credit market are denoted as ${a}_{i}^{B}$ and ${\ell}_{i}^{B}$. Total liabilities are denoted as ${\ell}_{i}$. At time 1, each bank*i*allocates its external assets in a portfolio of securities on the external market, ${E}_{ik}$, denoting the fraction of*i*’s external assets invested at time 1 in the security*k*. The unitary value of the external security*k*is ${x}_{k}^{E}$. Without loss of generality, at time 1, ${x}_{k}^{E}\left(1\right)=1$ for all*k*, and ${x}_{K}^{E}\left(2\right)$ is a random variable drawn from a given distribution. At time 2, then, the external assets of bank*i*is a sum of random variables, ${a}_{i}^{E}\left(2\right)={a}_{i}^{E}\left(1\right){\displaystyle {\sum}_{k}{E}_{ik}{x}_{K}^{E}\left(2\right)}$. For our purposes, it is sufficient to assume that we can express the external assets of bank*i*as follows: ${a}_{i}^{E}\left(2\right)={a}_{i}^{E}\left(1\right)\left(1+\mu +\sigma \hspace{0.17em}{u}_{i}\right)$, where ${u}_{i}$ is a random variable drawn from a given distribution with mean zero and finite variance, the parameter ${\mu}_{i}$ is the expected return of the portfolio, and ${\sigma}_{i}$ is a scaling factor controlling the effective magnitude of the shocks. We also assume to know the joint probability distribution $p\left({u}_{1},\dots ,{u}_{n}\right)$.At time 1, each bank

*i*allocates its interbank assets among the other banks, ${B}_{ij}$, denoting the fraction of*i*’s interbank assets invested at time 1 in the liability of bank*j*. These investments are secured via collateralization, i.e., bank*j*posts as collateral for the loan an asset that bank*i*will collect in case bank*j*defaults. In spirit, this approach is similar to ref. 19, although, for the purpose of our study, we exclude rehypothecation, i.e., assets used as collateral are kept aside and cannot be reused. The value of the collateral is assumed to stay constant during the two periods, and it can be ignored in the default condition described in*Default Condition*.The unitary value of the interbank liability of bank

*j*to other banks is ${x}_{j}^{B}$. Without loss of generality, at time 1, ${x}_{j}^{B}\left(1\right)=1$ for all*j*. The liabilities of bank*j*are constant in value from the perspective of bank*i*, i.e., the debt agreed upon in the contract at time 1. However, from the point of view of counterparties of*j*, ${x}_{j}^{B}\left(2\right)=1$ if bank*j*honors its obligation, and ${x}_{j}^{B}\left(2\right)={R}_{j}$ otherwise, where ${R}_{j}$ is the recovery rate, i.e., the fraction of the interbank asset that is covered by the collateral and that the lender can recover after the default of*j*. Accordingly, at time 2, the interbank assets of bank*i*are ${a}_{i}^{B}\left(2\right)={a}_{i}^{B}\left(1\right){\displaystyle {\sum}_{j}{B}_{ij}{x}_{j}^{B}\left(2\right)}$. Notice that this approach differs from previous models (26, 28) that are based on the determination ex post of the clearing vector of payments.### Default Condition.

The standard balance sheet identity in financial accounting states that equity of bank

*i*, ${e}_{i}$, is the difference between assets and liabilities. Hence, ${e}_{i}\left(2\right)={a}_{i}^{E}\left(2\right)+{a}_{i}^{B}\left(2\right)-{\ell}_{i}={a}_{i}^{E}\left(1\right)\left(1+\mu +\sigma \hspace{0.17em}{u}_{i}\right)+{a}_{i}^{B}\left(1\right){\displaystyle {\sum}_{j}{B}_{ij}{x}_{j}^{B}\left(2\right)}-{\ell}_{i}$. It is also standard to assume that the default of bank*i*occurs when equity becomes negative, i.e., if ${e}_{i}\left(2\right)<0$. In the following, we are interested in the probability of default of individual banks. Notice that we assume ${e}_{i}\left(1\right)>0$, and thus ${e}_{i}\left(2\right)<0$ iff $\left[{e}_{i}\left(2\right)\right]/\left[{e}_{i}\left(1\right)\right]<0$. It is then convenient to write the default condition as ${\epsilon}_{i}\left(1+\mu +\sigma \hspace{0.17em}{u}_{i}\right)+{\beta}_{i}{\displaystyle {\sum}_{j}{B}_{ij}{x}_{j}^{B}\left(2\right)}-{\lambda}_{i}<0$, where the parameter ${\epsilon}_{i}=\left[{a}_{i}^{E}\left(1\right)\right]/\left[{e}_{i}\left(1\right)\right]$ measures the leverage with respect to external assets, i.e., the magnitude, per unit of initial equity of bank*i*, of the investments of bank*i*in external assets (30). Similarly, the parameter ${\beta}_{i}=\left[{a}_{i}^{B}\left(1\right)\right]/\left[{e}_{i}\left(1\right)\right]$ measures the interbank leverage, i.e., the magnitude, per unit of initial equity, of*i*’s investments in interbank assets (30), and the parameter ${\lambda}_{i}=\left[{l}_{i}\left(1\right)\right]/\left[{e}_{i}\left(1\right)\right]$ measures the magnitude, per unit of initial equity, of*i*’s total liabilities. Let us define a default indicator ${\chi}_{i}$, with ${\chi}_{i}=1$ in case of default of bank*i*and ${\chi}_{i}=0$ otherwise. Because the only variable that is exogenously stochastic is the shock ${u}_{i}$ on each bank’s external assets, we finally write the default condition as follows:$${u}_{i}<{\theta}_{i}\equiv \frac{1}{{\epsilon}_{i}\sigma}\left({\lambda}_{i}-1-{\mu}_{i}-{\beta}_{i}{\displaystyle \sum _{j}{B}_{ij}{x}_{j}^{B}\left({\chi}_{j}\right)}\right),$$

[1]

where ${\theta}_{i}$ denotes the default threshold. Notice that we have dropped the time in the notation, and we have emphasized in the formula that the value of the interbank liability of a counterparty

*j*, ${x}_{j}^{B}$, depends on the default indicator of*j*, ${\chi}_{j}$, to recall that it is ${x}_{j}^{B}\left({\chi}_{j}=0\right)=1$ and ${x}_{j}^{B}\left({\chi}_{j}=1\right)={R}_{j}$. Thus, depending on the magnitude and sign of the shocks ${u}_{i}$ that hit all banks, some of them can default on their obligations, possibly causing other banks to default. We can now express the default indicators ${\chi}_{i}$ of all banks as a system of equations$$\forall i\hspace{1em}{\chi}_{i}=\mathrm{\Theta}\left[{u}_{i}-{\theta}_{i}\left({\chi}_{1},\dots ,{\chi}_{n}\right)\right],$$

[2]

where $\mathrm{\Theta}$ denotes the step function or Heaviside function (i.e., equal to 1 if the argument is positive, and zero otherwise). A solution of the system, denoted as ${\chi}^{*}$, depends on the vector of shocks

*u*and on the initial condition ${\chi}^{0}$, which represents the initial belief of the banks in other banks’ defaults. The existence and uniqueness of the solution is discussed in*Existence and Uniqueness of Solutions*.The determination of the fixed point of the map above can become computationally cumbersome if we want to sample, at a fine resolution, the shock space of an arbitrary number of banks. However, in the default condition of the counterparties

*j*of bank*i*, we can simplify the computation by replacing the value of the second-order counterparties’ credit obligations (i.e., the obligations of neighbors of order 2 in the contract network) with their expected value. In other words, we replace the stochastic variable ${x}_{k}^{B}$ with its expected value $E\left[{x}_{k}\right]=R{P}_{k}+\left(1-{P}_{k}\right)$, where ${P}_{k}$ is the default probability of bank*k*. Indeed, expected values are commonly used in banking practice to estimate the future value of assets. Notice also that this approximation does not remove the effects of correlations across shocks on banks.### Default Probability.

The probability ${P}_{i}$ of default at time 2 of bank

*i*is then simply the integral over the shock space of the default indicator$$\forall i\hspace{1em}{P}_{i}={\displaystyle \int {\chi}_{i}^{*}\left(u,{\chi}^{0}\right)p\left(u\right)du},$$

[3]

where $p\left(u\right)$ denotes the joint density function of the shocks and accounts for possible correlations across shocks. Finally the probability of systemic default is

$${P}^{\text{sys}}={\displaystyle \int {\chi}^{\text{sys}}\left[{\chi}^{*}\left(u,{\chi}^{0}\right)\right]p\left(u\right)du},$$

[4]

where ${\chi}^{\text{sys}}$ is the systemic default indicator. The choice of the systemic default identification can vary, and there is no consensus on what should be defined as a systemic event. For the sake of simplicity, here we consider the extreme but intuitive case of all banks defaulting, i.e., ${\chi}^{\text{sys}}={\mathrm{\Pi}}_{i}{\chi}_{i}^{*}$.

### Existence and Uniqueness of Solutions.

In our model, banks’ default conditions are described by a system of nonlinear equations, $\chi =\mathrm{\Theta}\left(u,\chi \right)$, for which no closed-form solution exists. However, finding all of the solutions of such a system is equivalent to studying, for any point

*u*of the shock space and for any initial condition on*χ*, the deterministic map of the finite set ${\left\{0,\hspace{0.17em}1\right\}}^{n}$ in itself: $\mathrm{\Theta}:{\left\{0,\hspace{0.17em}1\right\}}^{n}\to {\left\{0,\hspace{0.17em}1\right\}}^{n}$. Because the map is deterministic and the set is finite, for any given initial condition, the map must have either a unique fixed point or a limit cycle; otherwise the set would have to be nonfinite (31). If banks only make regular credit contracts with each other, then there is a unique solution. Indeed, the default state of a bank can only be affected adversely by the default of some counterparty. Hence, the default state of each bank is a nondecreasing function of the default states of the other banks. As a result, the map cannot enter a cycle involving more than one state. Otherwise, the default state of at least one bank would have to revert from default to no default, which is not possible, by assumption. Hence, the fixed point is unique for any given initial condition on*χ*.### Contract Errors.

By contract errors, we mean that the estimates of the contract characteristics available to the regulator who computes the systemic default probability deviate from the actual value. We define the relative error on a given parameter

*p*, the ratio ${\mathrm{\Delta}}_{p}=\left(\tilde{p}-{p}^{*}\right)/{p}^{*}$, i.e., the relative deviation of the regulator’s estimation, $\tilde{p}$, from the actual value ${p}^{*}$ of the variable*p*. For a given level of deviation ${\mathrm{\Delta}}_{p}^{*}$, we consider a realistic number giving a ${\mathrm{\Delta}}_{p}$ precision of $5\%$ around the original value ${p}^{*}$: ${\mathrm{\Delta}}_{p}\in \left[{p}^{*}-{\mathrm{\Delta}}_{p}^{*},{p}^{*}+{\mathrm{\Delta}}_{p}^{*}\right]$. For instance, a deviation of 10% on the interbank leverage*β*(i.e., ${\mathrm{\Delta}}_{\beta}^{*}=0.1$) will yield ${\mathrm{\Delta}}_{\beta}\in \left\{-10\%,-9\%,\dots ,0,\dots ,10\%\right\}$. For each value of ${\mathrm{\Delta}}_{\beta}$, we compute the corresponding probability of systemic default using $\tilde{\beta}={\beta}^{*}\left(1+{\mathrm{\Delta}}_{\beta}\right)$, and we record the minimum and the maximum values. Additionally, we consider the effect of combined errors where more than one parameter is subject to changes at the same time. Take as an example the leverage parameters*ε*and*β*. As $\tilde{\epsilon}$ and $\tilde{\beta}$ can each take 21 different values, we have ${21}^{2}$ combinations. In the spirit of what has been done previously, we compute all probabilities of systemic default using all possible combinations of $\tilde{\epsilon}$ and $\tilde{\beta}$ given a level of deviation ${\mathrm{\Delta}}_{\left(\epsilon ,\beta \right)}^{*}$.### Analytical Example.

To understand how the uncertainty of specific parameters can be amplified by the network structure, consider the following three basic architectures with three nodes: a star, a chain, and a ring. Here we compute the probability ${P}^{\text{sys}}$ of the event in which all banks default under the assumption that the shocks hitting the banks are independent and drawn from the same uniform distribution in the space $\left[-1,\hspace{0.17em}1\right]$. Denote by ${\theta}_{i}^{+}$ the specific value of the threshold ${\theta}_{i}$ in Eq.

**1**in the case that all of the counterparties of a given bank*i*default. Conversely, denote by ${\theta}_{i}^{-}$ the case of no defaulting counterparties. In particular, if a bank has no counterparties (e.g., a leaf node in a tree network structure), then it has ${\theta}_{i}={\theta}_{i}^{-}$. The computation yields ${P}^{\text{sys},\text{star}}=\left(1/{2}^{3}\right)\left(1+{\theta}_{1}^{+}\right)\left(1+{\theta}_{2}^{-}\right)\left(1+{\theta}_{3}^{-}\right)$; ${P}^{\text{sys},\text{chain}}=\left(1/{2}^{3}\right)\left(1+{\theta}_{1}^{+}\right)\left(1+{\theta}_{2}^{+}\right)\left(1+{\theta}_{3}^{-}\right)$; and ${P}^{\text{sys},\text{ring}}=\left(1/{2}^{3}\right)\left(1+{\theta}_{1}^{+}\right)\left(1+{\theta}_{2}^{+}\right)\left(1+{\theta}_{3}^{+}\right)$. Note that ${\theta}_{i}^{+}\ge {\theta}_{i}^{-}$ because a bank needs more positive shocks on its external assets to survive when counterparties on the interbank have defaulted. It follows that ${P}^{\text{sys},\text{ring}}\ge {P}^{\text{sys},\text{chain}}\ge {P}^{\text{sys},\text{star}}$. It also follows that the sensitivity of the default probability on the recovery rate*R*depends on the network structure $\partial {P}^{\text{sys},\text{ring}}/\partial R\propto {\left[\beta /\left(\mathit{\u03f5}\hspace{0.17em}\sigma \right)\right]}^{3}$; $\partial {P}^{\text{sys},\text{chain}}/\partial R\propto {\left[\beta /\left(\mathit{\u03f5}\hspace{0.17em}\sigma \right)\right]}^{2}$; and $\partial {P}^{\text{sys},\text{star}}/\partial R\propto \beta /\left(\mathit{\u03f5}\hspace{0.17em}\sigma \right)$. In this example, as long as $\beta /\left(\mathit{\u03f5}\hspace{0.17em}\sigma \right)>1$, which is empirically plausible in many cases, the sensitivity on errors on the recovery rate can be much larger in the ring than in the star.### Network Errors.

By network errors, we mean that the regulator does not have information on the arrangement of the contracts among banks but only their maximum number. We denote by ${C}_{l}$ the cap on the maximum number of possible contracts in the market (i.e., the network density). By increasing the cap ${C}_{l}$, the total number of possible configurations grows as follows: $\overline{\mathcal{N}}\left({C}_{l}\right)={\displaystyle {\sum}_{l=0}^{{C}_{l}}\left(n\left(n-1\right)/l\right)}$. For a given ${C}_{l}$, we inspect all possible configurations of links arrangement and compute the probability of systemic default. As a measure of the market complexity we take the so-called network entropy, relative to a given density cap (32), i.e., the logarithm of the number of configurations normalized by the number of nodes, $\mathrm{\Sigma}\left({C}_{l}\right)=\left(1/n\right)\hspace{0.17em}log\overline{\mathcal{N}}\left({C}_{l}\right)$.

### Available Code and Data.

Analyses and figures can be reproduced using Python, Matlab and C++ scripts and data files available in the public GitHub repository (https://github.com/troukny/price-of-complexity) following the instructions provided.

## Acknowledgments

S.B. and G.C. acknowledge support from the European Commission’s Global Systems Science (GSS) initiative and the Future and Emerging Technologies programme (FET) through the Grants SIMPOL 610704, DOLFINS 640772, MULTIPLEX 317532, and SoBigData 654024. S.B. acknowledges the Swiss National Fund Professorship Grant PP00P1-144689. T.R. acknowledges the support of Fonds National de la Recherche Scientifique. J.E.S. and S.B. acknowledge the support of the grant on Financial Stability from the Institute for New Economic Thinking (INET).

## References

1

AG Haldane, RM May, Systemic risk in banking ecosystems.

*Nature***469**, 351–355 (2011).2

; Bank of England

*A Framework for Stress Testing the UK Banking System*(Bank England, London, 2013).3

FJ Lambert, K Ueda, P Deb, DF Gray, P Grippa, How big is the implicit subsidy for banks considered too important to fail?

*Global Financial Stability Report*(Int Monetary Fund, Washington, DC), Chap 3. (2014).4

JC Rochet

*Why Are There So Many Banking Crises? the Politics and Policy of Bank Regulation*(Princeton Univ Press, Princeton, 2008).5

AG Haldane,

*Rethinking Financial Networks*(Bank of England, London). Available at www.bankofengland.co.uk/archive/Documents/historicpubs/speeches/2009/speech386.pdf. Accessed August 5, 2016. (2009).6

; International Monetary Fund

*Global Financial Stability Report: Meeting New Challenges to Stability and Building a Safer System*(Int Monetary Fund, Washington, DC, 2010).7

E Nier, J Yang, T Yorulmazer, A Alentorn, Network models and financial stability.

*J Econ Dyn Control***31**, 2033–2060 (2007).8

H Elsinger, A Lehar, M Summer, Risk assessment for banking systems.

*Manage Sci***52**, 1301–1314 (2006).9

N Beale, et al., Individual versus systemic risk and the Regulator’s Dilemma.

*Proc Natl Acad Sci USA***108**, 12647–12652 (2011).10

RM May, N Arinaminpathy, Systemic risk: The dynamics of model banking systems.

*J Roy Soc Interface***7**, 823–838 (2010).11

K Anand, P Gai, M Marsili, Rollover risk, network structure and systemic financial crises.

*J Econ Dyn Control***36**, 1088–1100 (2012).12

D Acemoglu, A Ozdaglar, A Tahbaz-Salehi, Systemic risk and stability in financial networks.

*Am Econ Rev***105**, 564–608 (2015).13

M Elliott, B Golub, MO Jackson, Financial networks and contagion.

*Am Econ Rev***104**, 3115–3153 (2014).14

TA Peltonen, M Scheicher, G Vuillemey, The network structure of the CDS market and its determinants.

*J Financ Stab***13**, 118–133 (2014).15

A Capponi, M Larsson, Price contagion through balance sheet linkages.

*Rev Asset Pric Stud***5**, 227–253 (2015).16

S Battiston, G Caldarelli, C-P Georg, R May, J Stiglitz, Complex derivatives.

*Nat Phys***9**, 123–125 (2013).17

RJ Caballero, A Simsek, Fire sales in a model of complexity.

*J Finance***68**, 2549–2587 (2013).18

N Arinaminpathy, S Kapadia, RM May, Size and complexity in model financial systems.

*Proc Natl Acad Sci USA***109**, 18338–18343 (2012).19

P Gai, A Haldane, S Kapadia, Complexity, concentration and contagion.

*J Monet Econ***58**, 453–470 (2011).20

E Farhi, J Tirole, Collective Moral hazard, maturity mismatch and systemic bailouts.

*Am Econ Rev***102**, 60–93 (2012).21

JE Stiglitz, Risk and global economic architecture: Why full financial integration may be undesirable.

*Am Econ Rev***100**, 388–392 (2010).22

P Gai, S Kapadia, Contagion in financial networks.

*Proc Roy Soc A Math Phys Eng Sci***466**, 2401–2423 (2010).23

S Battiston, D Delli Gatti, M Gallegati, B Greenwals, JE Stiglitz, Liaisons dangereuses: Increasing connectivity, risk sharing, and systemic risk.

*J Econ Dyn Control***36**, 1121–1141 (2012).24

S Battiston, M Puliga, R Kaushik, P Tasca, G Caldarelli, DebtRank: Too central to fail? Financial networks, the FED and systemic risk.

*Sci Rep***2**, 541 (2012).25

R Cont, A Moussa, EB Santos, Network structure and systemic risk in banking systems.

*Handbook of Systemic Risk*, eds J-P Fouque, JA Lansam (Cambridge Univ Press, Cambridge, UK), pp. 327–368 (2013).26

L Eisenberg, TH Noe, Systemic risk in financial systems.

*Manage Sci***47**, 236–249 (2001).27

R Cifuentes, G Ferrucci, HS Shin, Liquidity risk and contagion.

*J Eur Econ Assoc***3**, 556–566 (2005).28

LCG Rogers, LAM Veraart, Failure and rescue in an interbank network.

*Manage Sci***59**, 882–898 (2013).29

AG Haldane,

*Speech: The dog and Frisbee*(Bank of England, London). Available at www.bankofengland.co.uk/archive/Documents/historicpubs/speeches/2012/speech596.pdf. Accessed August 5, 2016. (2012).30

S Battiston, G Caldarelli, M D’Errico, S Gurciullo, Leveraging the network: A stress-test framework based on DebtRank.

*Statistics and Risk Modelling*, in press. (2016).31

R Laubenbacher, B Pareigis, Equivalence relations on finite dynamical systems.

*Adv Appl Math***26**, 237–251 (2001).32

K Anand, G Bianconi, Entropy measures for networks: toward an information theory of complex topologies.

*Phys Rev E Stat Nonlin Soft Matter Phys***80**, 045102 (2009).## Information & Authors

### Information

#### Published in

#### Classifications

#### Copyright

Freely available online through the PNAS open access option.

#### Submission history

**Published online**: August 23, 2016

**Published in issue**: September 6, 2016

#### Keywords

#### Acknowledgments

S.B. and G.C. acknowledge support from the European Commission’s Global Systems Science (GSS) initiative and the Future and Emerging Technologies programme (FET) through the Grants SIMPOL 610704, DOLFINS 640772, MULTIPLEX 317532, and SoBigData 654024. S.B. acknowledges the Swiss National Fund Professorship Grant PP00P1-144689. T.R. acknowledges the support of Fonds National de la Recherche Scientifique. J.E.S. and S.B. acknowledge the support of the grant on Financial Stability from the Institute for New Economic Thinking (INET).

#### Notes

This article is a PNAS Direct Submission.

### Authors

#### Competing Interests

The authors declare no conflict of interest.

## Metrics & Citations

### Metrics

#### Citation statements

#### Altmetrics

### Citations

If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.

#### Cited by

Loading...

## View Options

### View options

#### PDF format

Download this article as a PDF file

DOWNLOAD PDF### Get Access

#### Login options

Check if you have access through your login credentials or your institution to get full access on this article.

Personal login Institutional Login#### Recommend to a librarian

Recommend PNAS to a Librarian#### Purchase options

Purchase this article to get full access to it.