Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

Economics

5 Nonlinear Interactions

Authors
Affiliations
New York University
National Graduate Institute for Policy Studies (GRIPS)

Much of what makes network analysis interesting is how ramifications of choices flow across networks. In general, decisions made at a given node ii affect responses of neighboring nodes and, through them, neighbors of neighboring nodes, and so on. As these consequences flow across the network, they in turn affect choices at ii. This is a tail-chasing scenario, which can be unraveled through fixed point theory.

In some network settings, such as the input-output model in Section 2.1, interactions are linear and fixed point problems reduce a system of linear equations. In other settings, interactions are inherently nonlinear and, as a result, we need more sophisticated fixed point theory.

This chapter is dedicated to the study of networks with nonlinear interactions. We begin with relevant fixed point theory and then apply it to a sequence of problems that arise in analysis of economic networks, including production models with supply constraints and financial networks.

5.1Fixed Point Theory

Let SS be any set. Recall from Section 6.1.1.6 that, given a self-map GG on SS, a point x∈Sx \in S is called a fixed point of GG if Gx=xGx = x. (A self-map on SS is a function GG from SS to itself. When working with self-maps it is common to abbreviate G(x)G(x) to GxGx.) In this chapter, we will say that GG is globally stable on SS if GG has a unique fixed point x∗∈Sx^* \in S and Gkx→x∗G^k x \to x^* as k→∞k \to \infty for all x∈Sx \in S. In other words, under this property, the fixed point is not only unique but also globally attracting under iteration of GG.

We have already discussed fixed points, indirectly or directly, in multiple contexts:

In each case, when we introduced these fixed point problems, we immediately needed to consider questions of existence and uniqueness of fixed points. Now we address these same issues more systematically in an abstract setting.

5.1.1Contraction Mappings

In Chapter 2 we studied solutions of the system x=Ax+bx = Ax + b that are fixed points of the affine map Fx=Ax+bFx = Ax + b on Rn\RR^n studied in Example 5.1.1. The Neumann series lemma is, in essence, a statement about existence and uniqueness of fixed points of this map. Here we investigate another fixed point theorem, due to Stefan Banach (1892–1945), that can be thought of as extending the Neumann series lemma to nonlinear systems.

5.1.1.1Contractions

Let SS be a nonempty subset of Rn\RR^n. A self-map FF on SS is called contracting or a contraction of modulus λ\lambda on SS if there exists a λ<1\lambda < 1 and a norm ∥⋅∥\| \cdot \| on Rn\RR^n such that

∥Fu−Fv∥⩽λ∥u−v∥for allu,v∈S.\| Fu - Fv \| \leq \lambda \| u - v \| \quad \text{for all} \quad u, v \in S.

The next example uses a similar idea but is based on a different norm.

Consider again Example 5.1.1. For the affine map Fx=Ax+bFx = Ax + b, the condition ∥A∥<1\| A \| < 1 used to obtain contraction is stronger than the condition r(A)<1r(A) < 1 used to obtain a unique fixed point in the Neumann series lemma (see Exercise 2.3.11). Furthermore, the Neumann series lemma provides a geometric series representation of the fixed point. On the other hand, as we now show, the contraction property can be used to obtain unique fixed points when the map in question is not affine.

5.1.1.2Banach’s Theorem

The fundamental importance of contractions stems from the following theorem.

We complete a proof of Theorem 5.1.1 in stages.

5.1.1.3Eventual Contractions

Let SS be a nonempty subset of Rn\RR^n. A self-map FF on SS is called eventually contracting if there exists a k∈Nk \in \NN such that FkF^k is a contraction on SS. Significantly, most of the conclusions of Banach’s theorem carry over to the case where FF is eventually contracting.

There is a close connection between Theorem 5.1.2 and the Neumann series lemma (NSL). If S=RnS=\RR^n and Fx=Ax+bF x = A x + b with r(A)<1r(A) < 1, then the NSL implies a unique fixed point. We can also obtain this result from Theorem 5.1.2, since, for any k∈Nk \in \NN,

∥Fkx−Fky∥=∥Akx−Aky∥=∥Ak(x−y)∥⩽∥Ak∥∥x−y∥.\| F^k x - F^k y \| = \| A^k x - A^k y \| = \| A^k (x - y) \| \leq \|A^k \| \| x - y \|.

As r(A)<1r(A) < 1, we can choose kk such that ∥Ak∥<1\| A^k \| < 1 (see Section 2.3.3). Hence FF is eventually contracting and Theorem 5.1.2 applies.

As mentioned above, contractions and eventual contractions have much wider scope than the NSL, since they can also be applied in nonlinear settings. At the same time, the NSL is preferred when its conditions hold, since it also gives inverse and power series representations of the fixed point.

5.1.1.4A Condition for Eventual Contractions

The result below provides a useful test for the eventual contraction property. (In the statement, the absolute value of a vector is defined pointwise, as in Section 6.1.2.2.)

5.1.2Shortest Paths Revisited

Consider again the shortest path problem introduced in Section 3.1. One modification that sometimes appears in applications is the addition of discounting during travel between vertices. For example, if the vertices are international ports and travel takes place by sea, then port-to-port travel time is measured in weeks or even months. It is natural to apply time discounting to future costs associated with that travel, to implement the idea that paying a given dollar amount in the future is preferable to paying it now.

Recall that, without discounting, the Bellman equation for the shortest path problem takes the form q(x)=min⁡y∈O(x){c(x,y)+q(y)}q(x) = \min_{y \in \oO(x)} \{ c(x, y) + q(y) \} for all x∈Vx \in V, where cc is the cost function VV is the set of vertices and qq is a candidate for the cost-to-go function. We showed that the minimum cost-to-go function q∗q^* satisfies the Bellman equation and is the unique fixed point of the Bellman operator.

The Bellman equation neatly divides the problem into current costs, embedded in the term c(x,y)c(x, y), and future costs embedded in q(y)q(y). To add discounting, we need only discount q(y)q(y). We do this by multiplying it by a discount factor β∈(0,1)\beta \in (0, 1). The Bellman equation is then q(x)=min⁡y∈O(x){c(x,y)+βq(y)}q(x) = \min_{y \in \oO(x)} \{ c(x, y) + \beta q(y) \} for all x∈Vx \in V and the Bellman operator is

Tq(x)=min⁡y∈O(x){c(x,y)+βq(y)}(x∈V).Tq(x) = \min_{y \in \oO(x)} \{ c(x, y) + \beta q(y) \} \qquad (x \in V).

In Section 3.1.2, without discounting, we had to work hard to show that the Bellman operator has a unique fixed point in UU, the set of all q ⁣:V→R+q \colon V \to \RR_+ with q(d)=0q(d)=0. With discounting the proof is easier, since we can leverage the Banach contraction mapping theorem.

In what follows, we identify the vertices in VV with integers 1,…,n1, \ldots, n, where dd is identified with nn. We then understand UU as all nonnegative vectors qq in Rn\RR^n with q(n)=0q(n) = 0. (We continue to write q(x)q(x) for the xx-th element of the vector qq, but now xx is in [n]\natset{n}.)

Now let ∥⋅∥∞\| \cdot \|_\infty be the supremum norm on Rn\RR^n (see Section 2.3.1). We claim that TT is a contraction on UU of modulus β\beta. To see that this is so, fix p,q∈Up, q \in U and observe that, pointwise,

Tq=T(p+q−p)⩽T(p+∥q−p∥∞1)⩽Tp+β∥q−p∥∞1,Tq = T(p + q - p) \leq T(p + \|q - p\|_\infty \1) \leq Tp + \beta \|q - p\|_\infty \1,

where the first inequality is by the order-preserving of TT and the second follows from Exercise 5.1.9. Hence

Tq−Tp⩽β∥q−p∥∞1.Tq - Tp \leq \beta \|q - p\|_\infty \1.

Reversing the roles of pp and qq gives the reverse inequality. Hence

∣Tq(x)−Tp(x)∣⩽β∥q−p∥∞|Tq (x) - Tp(x) | \leq \beta \|q - p\|_\infty

for all x∈[n]x \in \natset{n}. Taking the maximum on the left hand side yields ∥Tq−Tp∥∞⩽β∥q−p∥∞\|Tq - Tp \|_\infty \leq \beta \|q - p\|_\infty, which shows that TT is a contraction of modulus β\beta. Hence Banach’s theorem applies and a unique fixed point exists.

5.1.3Supply Constraints

While the input-output model from Section 2.1.2 has many useful applications, its linear structure can be a liability. One natural objection to linearity is supply constraints: if sector jj doubles its orders from sector ii, we cannot always expect that sector ii will quickly meet this jump in demand.

In this sector we investigate the impact of supply constraints on equilibrium. These constraints introduce nonlinear relationships between nodes that affect equilibria and make analysis more challenging.

We recall from Section 2.1.2 that did_i is final demand for good ii, xix_i is total sales of sector ii, zijz_{ij} is inter-industry sales from sector ii to sector jj, and aij=zij/xja_{ij} = z_{ij}/x_j is dollar value of inputs from ii per dollar output from jj.

Departing from our previous formulation of equilibrium in the input-output model, suppose that, in the short run, the total output value of sector ii is constrained by positive constant xˉ(i)\bar x(i). Holding prices fixed (in the short run), this means that sector ii has a capacity constraint in terms of unit output. For the purposes of our model, the vector of capacity constraints xˉ:=(xˉ(i))i=1n\bar x := (\bar x(i))_{i=1}^n can be any vector in R+n\RR^n_+.

For each sector ii, we modify the accounting identity (2.1) to

xi=min⁡{∑j=1noij+di,  xˉ(i)},x_i = \min \left\{ \sum_{j=1}^n o_{ij} + d_i, \; \bar x(i) \right\},

where oijo_{ij} is the value of orders from sector ii made by sector jj. Thus, if the capacity constraint in sector ii is not binding, then output is the sum of orders from other sectors and orders from final consumers. If xˉ(i)\bar x(i) is less than this sum, however, then sector ii produces to capacity xˉ(i)\bar x(i).

An equilibrium for this model is one where all orders are met, subject to capacity constraints. The fact that orders are met means that oij=zij=aijxjo_{ij} = z_{ij} = a_{ij} x_j. Substituting this equality into (5.6) and rewriting as a vector equality, (5.6) can equivalently be formulated as

x=GxwhereGx:=(Ax+d)∧xˉ.x = G x \quad \text{where} \quad G x := (A x + d) \wedge \bar x.

The following exercise is key to solving the fixed point problem (5.7).

We are now ready to prove existence of a unique fixed point under the assumption that every sector has positive value added.

5.1.4Fixed Points and Monotonicity

Banach’s fixed point theorem and its extensions are foundations of many central results in pure and applied mathematics. For our purposes, however, we need to search a little further, since not all mappings generated by network models have the contraction property. In this section, we investigate two fixed point results that drop contractivity in favor of monotonicity.

5.1.4.1Existence

Without contractivity, one needs to work harder to obtain even existence of fixed points, let alone uniqueness and convergence. This is especially true if the map in question fails to be continuous. If, however, the map is order preserving, then existence can often be obtained via some variation on the Knaster–Tarski fixed point theorem.

Here we present a version of this existence result that is optimized to our setting, while avoiding unnecessary excursions into order theory. In stating the theorem, we recall that a closed order interval in Rn\RR^n is a set of the form

[a,b]:={x∈Rn:a⩽x⩽b}[a, b] := \setntn{x \in \RR^n}{a \leq x \leq b}

where aa and bb are vectors in Rn\RR^n. Also, we call (xk)⊂Rn(x_k) \subset \RR^n increasing (resp., decreasing) if xk⩽xk+1x_k \leq x_{k+1} (resp., xk⩾xk+1x_k \geq x_{k+1}) for all kk.

For a self-map GG on S⊂RnS \subset \RR^n, we say that x∗x^* is a least fixed point (resp., greatest fixed point) of GG on SS if x∗x^* is a fixed point of GG in SS and x∗⩽xx^* \leq x (resp., x⩽x∗x \leq x^*) for every fixed point xx of GG in SS. Finally, we say that GG is

Here xk↑xx_k \uparrow x means that (xk)(x_k) is increasing and xk→xx_k \to x. The definition of xk↓xx_k \downarrow x is analogous. In the next theorem, S:=[a,b]S := [a, b] is a closed order interval in Rn\RR^n and GG is a self-map on SS.

5.1.4.2Du’s Theorem

Theorem 5.1.5 is useful because of its relatively weak assumptions. At the same time, it fails to deliver uniqueness. Hence its conclusions are considerably weaker than the results we obtained from contractivity assumptions in Section 5.1.1.

In order to recover uniqueness without imposing contractivity, we now consider order-preserving maps that have additional shape properties. In doing so, we use the definition of concave and convex functions in Section 6.1.5.2.

A proof of Theorem 5.1.6 was obtained in an more abstract setting in Du (1990). Interested readers can consult that article for a proof.

To illustrate how these results can be applied, consider the constrained production model without assuming positive value added, so that r(A)<1r(A) < 1 is not enforced. In Exercise 5.1.12 we obtained existence. With Theorem 5.1.6 in hand, we can also show uniqueness whenever d≫0d \gg 0 and xˉ≫0\bar x \gg 0.

Indeed, we have already seen that GG is a self-map on S:=[0,xˉ]S := [0, \bar x] and, when this last condition holds, we have G0=d∧xˉ≫0G 0 = d \wedge \bar x \gg 0. Hence the conclusions of Theorem 5.1.6 will hold if we can establish that GG is concave.

Here is a small extension of Du’s theorem that will prove useful soon:

5.2Financial Networks

Given the long history of crises in financial markets around the globe, economists and business analysts have developed many tools for assessing the credit-worthiness of banks and other financial institutions. After the major financial crises of 2007-2008, originating in the subprime market in the US and the sudden collapse of Lehmann Brothers, it became clear that the financial health of individual institutions cannot be assessed in isolation. Rather, it is essential to analyze solvency and credit-worthiness in terms of the entire network of claims and liabilities within a highly interconnected financial system. In this section, we review financial crises and apply network analysis to study how they evolve.

5.2.1Contagion

Some financial crises have obvious causes external to the banking sector. A prominent example is the hyperinflation that occurred in Weimar Germany around 1921–1923, which was driven by mass printing of bank notes to meet war reparations imposed under the Treaty of Versailles. Here the monetary authority played the central role, while the actions of private banks were more passive.

Other crises seem to form within the financial sector itself, driven by interactions between banks, hedge funds, and asset markets. In many cases, the crisis follows a boom, where asset prices rise and economic growth is strong. Typically, the seeds of the crisis are laid during this boom phase, when banks extend loans and firms raise capital on the basis of progressively more speculative business plans. At some point it becomes clear to investors that these businesses will fail to meet expectations, leading to a rush for the exit.

The last phase of this cycle is painful for the financial sector, since rapidly falling asset values force banks and other financial institutions to generate short-term capital by liquidating long-term loans, typically with large losses, as well as selling assets in the face of falling prices, hoarding cash and refusing to roll over or extend short term loans to other institutions in the financial sector. The financial crisis of 2007–2008 provides a textbook example of these dynamics.

One key aspect of the financial crisis of 2007–2008, as well as other similar crises, is contagion, which refers to the way that financial stress spreads across a network of financial institutions. If one institution becomes stressed, that stress will often spread to investors or counterparties to which this institution is indebted. The result of this process is not easy to predict, since, like equilibrium in the production networks studied in Chapter 2, there is a tail chasing problem: stress spreads from institution AA to institutions BB, CC and DD, which may in turn increase stress on AA, and so on.

In this section we study financial contagion, beginning with a now-standard model of default cascades.

5.2.2Default Cascades

Default cascades are a form of financial contagion, in which default by a node in a network leads to default by some of its counterparties, which then spreads across the network. Below we present a model of default cascades and analyze its equilibria.

5.2.2.1Network Valuation

Consider a financial network G=(V,E,w)\gG = (V, E, w), where V=[n]V = \natset{n} is a list of nn financial institutions called banks, with an edge (i,j)∈E(i, j) \in E indicating that jj has extended credit to ii. The size of that loan is w(i,j)w(i, j). Thus, an edge points in the direction of a liability, as in Figure 1.17: edge (i,j)(i, j) indicates a liability for ii and an asset for jj. As in Section 1.4.1, the set of all direct predecessors of i∈Vi \in V will be written as I(i)\iI(i), while the set of all direct successors will be denoted O(i)\oO(i).

Banks in the network have both internal and external liabilities, as well as internal and external assets. Internal (i.e., interbank) liabilities and assets are given by the weight function ww, in the sense that w(i,j)w(i,j) is a liability for ii, equal to the size of its loan from jj, and also an asset for jj. Positive weights indicate the presence of counterparty risk: when jj holds an asset of book value w(i,j)w(i, j), whether or not the loan is repaid in full depends on the stress placed on bank ii and rules that govern repayment in the event of insolvency.

We use the following notation for the primitives of the model: xi:=∑j∈O(i)w(i,j)x_i := \sum_{j \in \oO(i)} w(i, j) is total interbank liabilities of bank ii,

Πij:={w(i,j)/xi if xi>00 otherwise\Pi_{ij} := \begin{cases} w(i, j) / x_i & \text{ if } x_i > 0 \\ 0 & \text{ otherwise} \end{cases}

is the matrix of relative interbank liabilities, aia_i is external assets held by bank ii and did_i is external liabilities.

When considering the interbank assets of bank jj, we need to distinguish between the book value ∑i∈I(j)w(i,j)\sum_{i \in \iI(j)} w(i, j) of its claims on other banks and the realized value in the face of partial or complete default by its counterparties I(j)\iI(j). To this end, we introduce a clearing vector p∈R+np \in \RR^n_+, which is a list of proposed payments by each bank in the network. In particular, pip_i is total payments made by bank ii to its counterparties within the banking sector. Under the choice of a particular clearing vector, the actual payments received by bank jj on its internal loan portfolio are ∑i∈VpiΠij\sum_{i \in V} p_i \Pi_{ij}.

The last statement is an assumption about the legal framework for the banking sector. It means that the actual payment piΠijp_i \Pi_{ij} from ii to jj is proportional to the amount that ii owes jj, relative to its total interbank obligations. The idea is that all counterparties in the banking sector have equal seniority, so that residual funds are spread across claimants according to the relative size of the claims.

Let p^j\hat p_j be the amount of funds bank jj makes available to repay all of its debts, both interbank and external. This quantity is

p^j=min⁡{aj+∑i∈VpiΠij,  dj+xj}.\hat p_j = \min \left\{ a_j + \sum_{i \in V} p_i \Pi_{ij}, \; d_j + x_j \right\}.

The right hand term inside the min⁡\min operator is total debts of bank jj. The left hand side is the amount on hand to repay those debts, including external assets and repayments by other banks. The bank repays up to—but not beyond—its ability to pay.

External liabilities are assumed to be senior to interbank liabilities, which means that for bank jj we also have

pj=max⁡{p^j−dj,  0}.p_j = \max\{\hat p_j - d_j, \; 0\}.

Thus, interbank payments by jj are a remainder after external debts are settled. If these debts exceed the bank’s ability to pay, the bank becomes insolvent and pays nothing to internal creditors. This is a form of limited liability.

Combining (5.10) and (5.11) and rearranging slightly yields

pj=max⁡{min⁡{aj−dj+∑i∈VpiΠij,  xj},  0}p_j = \max \left\{ \min \left\{ a_j - d_j + \sum_{i \in V} p_i \Pi_{ij}, \; x_j \right\}, \; 0 \right\}

Now let’s take pp, aa, dd and xx as row vectors in Rn\RR^n and write this collection of equations, indexed of jj, in vector form. With max⁡\max and min⁡\min taken pointwise, and using the symbols ∨\vee and ∧\wedge for max and min, we get

p=((a−d+pΠ)∧x)∨0.p = ((a - d + p\Pi) \wedge x) \vee 0.

A solution to this equation is called an equilibrium clearing vector for the banking system.

5.2.2.2Existence and Uniqueness of Fixed Points

In order to analyze existence and uniqueness of equilibria, we introduce the operator T ⁣:Rn→RnT \colon \RR^n \to \RR^n defined by

Tp=((e+pΠ)∧x)∨0,Tp = ((e + p\Pi) \wedge x) \vee 0,

where e:=a−de := a - d represents net external assets. Evidently p∈R+np \in \RR^n_+ is an equilibrium clearing vector if and only if it is a fixed point of TT.

Using this operator, establishing existence of at least one equilibrium clearing vector is not problematic for this model, regardless of the values of the primitives and configuration of the network:

While existence is automatic in this model, uniqueness is not:

There are several approaches to proving uniqueness of fixed points of TT. Here is one:

This leads us to the following result:

There are several possible assumptions about the structure of the bank network G\gG that imply the conditions of Proposition 5.2.1. For example, it would be enough that G\gG is strongly connected and has at least one bank with zero interbank liabilities. Below we investigate another sufficient condition, related to cyclicality.

5.2.2.3Nonnegative External Equity

In this section we investigate the special case in which

(E1) each bank has nonzero interbank debt, so that Π\Pi is a stochastic matrix, and

(E2) external net assets are nonnegative, in the sense that e=a−d⩾0e = a - d \geq 0.

In view of (E1), we cannot hope to use Proposition 5.2.1, since we always have r(Π)=1r(\Pi)=1 when Π\Pi is stochastic. Nonetheless, we can still obtain global stability under certain restrictions on ee and the topology of the network. Here is a relatively straightforward example, which we will later try to refine:

The condition e≫0e \gg 0 is rather strict. Fortunately, it turns out that we can obtain global stability under significantly weaker conditions. To this end, we will say that node jj in a financial system G\gG, is cash accessible if there exists an i∈Vi \in V such that i→ji \to j and e(i)>0e(i) > 0. In other words, jj is downstream in the liability chain from at least one bank with positive net assets outside of the banking sector.

With the result from Exercise 5.2.6 in hand, the next lemma is easy to establish.

Stronger results are available, with a small amount of effort. In fact, under (E1)–(E2), there is a strong sense in which uniqueness of the fixed point is obtained without any further assumptions, provided that we rule out an ambiguity related to what happens when e=0e=0. That ambiguity is discussed in the next exercise and further references can be found in Section 5.3.

Although Exercise 5.2.7 suggests ambiguity about outcomes when e=0e = 0, it is natural to adopt the convention that the equilibrium clearing vector p∗p^* obeys p∗=0p^*=0 whenever e=0e=0. If the entire banking sector has no zero net assets, then no positive payment sequence can be initiated (without outside capital).

In the next exercise, we say that U⊂VU \subset V is accessible from i∈Vi \in V if there exists a j∈Uj \in U such that jj is accessible from ii.

5.2.3Equity-Cross Holdings

In this section we analyze the model of default cascades constructed by Elliott et al., 2014. The model differs from the one studied in Section 5.2.2 in several ways. One is that financial institutions are linked by share cross-holdings: firm ii owns a fraction cijc_{ij} of firm jj for i,j∈V:=[n]i, j \in V := \natset{n}. This implies that failure of firm jj reduces the market value of firm ii, which in turn reduces the market value of other firms, and so on.

The second—and ultimately more significant—difference is the introduction of failure costs that add significant nonlinearities to the model. Failure costs reinforce the impact of each firm failure, leading to greater shock propagation across the network. This feature ties into the intuitive idea that, when many firms are financially stressed, a failure by one firm can trigger a wave of bankruptcies.

5.2.3.1Book and Market Value

We now describe the features of the model. Let C=(cij)i,j∈VC = (c_{ij})_{i, j \in V} be the matrix of fractional cross-holdings, as mentioned above, with 0⩽cij⩽10 \leq c_{ij} \leq 1 for all i,ji,j.

Assumption 5.2.1 implies that firms are not solely owned by other firms in the network: investors outside the network hold at least some fraction of firm ii for all ii.

The book value of firm ii is given by

bi=ei+∑jcijbj(i∈V).b_i = e_i + \sum_j c_{ij} b_j \qquad (i \in V).

Here the first term ei⩾0e_i \geq 0 is external assets of firm ii and the second represents the value of firm ii’s cross-holdings. In vector form, the last equation becomes b=e+Cbb = e + C b.

There is a widely used argument that cross-holdings artificially inflate the value of firms, in the sense that the sum of book values of firms exceeds ∑iei\sum_i e_i, the sum of underlying equity values. The next exercise illustrates:

Due to this artificial inflation, we distinguish between the book value bib_i of a firm and its market value vˉi\bar v_i, which is defined as ribir_i b_i with ri:=1−∑kckir_i := 1 - \sum_{k} c_{ki}. The value rir_i, which gives the share of firm ii held by outsider investors, is strictly positive for all ii by Assumption 5.2.1. With R:=diag⁡(r1,…,rn)R := \diag(r_1, \ldots, r_n), we can write the vector of market values as vˉ=Rb\bar v = R b. Substituting in (5.16) gives

vˉ:=Ae,whereA:=R(I−C)−1.\bar v := A e , \quad \text{where} \quad A : =R (I - C)^{-1}.

5.2.3.2Failure Costs

So far the model is very straightforward, with the market value of firms being linear in external assets ee. However, since bankruptcy proceedings are expensive, it is reasonable to assume that firm failures are costly. Moreover, when the market value of a firm falls significantly, the firm will experience difficulty raising short term funds, and will often need to cease revenue generating activities and sell illiquid assets well below their potential value.

We now introduce failure costs. As before, vˉi\bar v_i is market value without failure costs, as determined in (5.17), while viv_i will represent market value in the presence of failure costs. Failure costs for firm ii are modeled as a threshold function

f(vi)=β1{vi<θvˉi}(i=1,…,n),f(v_i) = \beta \1\{v_i < \theta \bar v_i\} \qquad (i = 1, \ldots, n),

where θ∈(0,1)\theta \in (0, 1) and β>0\beta > 0 are parameters. Thus, costs are zero when viv_i is large and −β-\beta when they fall below the threshold θvˉi\theta \bar v_i. In particular, a discrete failure cost of −β-\beta is incurred when firm value falls far enough below the no-failure market value vˉi\bar v_i. The larger is θ\theta, more prone firms are to failure.

The book value of firm ii without failure costs was defined in (5.15). The book value of firm ii in the presence of failure costs is defined as

bi=ei+∑jcijbj−f(vi).b_i = e_i + \sum_j c_{ij} b_j - f(v_i).

Written in vector form, with ff applied pointwise to the vector vv, we get b=e+Cb−f(v)b = e + C b - f (v). Solving for bb gives b=(I−C)−1(e−f(v))b = (I-C)^{-1} (e - f (v)). The corresponding market value is

v=Rb=A(e−f(v)).v = R b = A (e - f(v)).

Notice that, when no firms fail, we have vi=vˉiv_i = \bar v_i, as expected.

5.2.3.3Equilibria

Equation (5.18) is a nonlinear equation in nn unknowns. A vector v∈Rnv \in \RR^n solves (5.18) if and only if it is a fixed point of the operator T ⁣:Rn→RnT \colon \RR^n \to \RR^n defined by

Tv=A(e−f(v)).Tv = A (e - f(v)).

In what follows we set d:=A(e−β1)d := A(e - \beta \1) and S:=[d,vˉ]S := [d, \bar v].

Note that v∗=v∗∗v^* = v^{**} is a possibility, in which case TT has a unique fixed point in SS.

Why does (Tkd)k∈N(T^k d)_{k \in \NN} converge up to v∗v^* in a finite number of steps? First, as you saw in the solution to Exercise 5.2.12, the map TT is an order-preserving self-map on [d,vˉ][d, \bar v], so Td∈[d,vˉ]Td \in [d, \bar v]. In particular, d⩽Tdd \leq Td. Iterating on this inequality and using the order-preserving property gives d⩽Td⩽T2d⩽⋯d \leq Td \leq T^2 d \leq \cdots, so (Tkd)(T^k d) is indeed increasing. Moreover, the range of TT is a finite set, corresponding to all vectors of the form

u=A(e−βw),u = A(e - \beta w),

where ww is an nn-vector containing only zeros and ones. Finiteness holds because there are only finitely many binary sequences of length nn.

Similar logic can be applied to prove that (Tkvˉ)k∈N(T^k \bar v)_{k \in \NN} converges down to v∗∗v^{**} in a finite number of steps.

If we set v0=vˉv^0 = \bar v and vk+1=Tvkv^{k+1} = T v^k, we can consider the sequence of valuations (vk)(v^k) as a dynamic process, and the number of currently failing firms mk:=∑i1{vik<θ}m^k := \sum_i \1\{v^k_i < \theta\} can be understood as tracking waves of bankruptcies. Failures of firms in the first wave put stress on otherwise healthy firms that have exposure to the failed firms, which in turn causes further failures and so on.

As discussed in Proposition 5.2.3, the sequence (vk)(v^k), which can also be written as (Tkvˉ)(T^k \bar v), is decreasing pointwise. In other words, the value of each firm is non-increasing. Hence, if vik<θv^k_i < \theta for some kk, then vik+j<θv^{k+j}_i < \theta for all j⩾0j \geq 0.

Figure 5.1 illustrates a growing wave of failures that can arise in a financial network. Firms with lighter colors have better balance sheets. Black firms have failed. The code for generating this figure, along with details on parameters, can be found in the code book.

Waves of bankruptcies in a financial network

Figure 5.1:Waves of bankruptcies in a financial network

5.3Chapter Notes

Shin (2010) gives an excellent overview of systemic risk and the financial crisis of 2007–2008. Battiston et al. (2012) study connectedness in financial networks and introduce a measure of systemic impact called DebtRank. Bardoscia et al. (2015) provide a dynamic theory of instability related to DebtRank. Demange (2018) provides a threat index for contagion in financial networks related to Katz centrality. Bardoscia et al. (2019) analyze risks associated with solvency contagion during crises. Jackson & Pernoud (2019) study investment in risky portfolios by banks in financial networks with debt and equity interdependencies. Jackson & Pernoud (2020) analyze optimal bailouts in financial networks. Amini & Minca (2020) offer an introduction to modeling of clearing systems, analyzing equilibria for network payments and asset prices. Jackson & Pernoud (2021) provide a survey of the relationship between financial networks and systemic risk.

The setting in Section 5.2.2.3 was investigated by Eisenberg & Noe (2001), one of the first papers on a network approach to default cascades. Additional details on stability properties can be found in Stachurski (2022). As already mentioned, Section 5.2.3 is based on Elliott et al., 2014, which also includes an interesting discussion of how the level of integration across a network affects equilibria. Klages-Mundt & Minca (2021) discuss optimal intervention in economic networks via influence maximization methods, using Elliott et al. (2014) as an example. Acemoglu et al. (2021) study how anticipation of future defaults can result in "credit freezes."

A general discussion of contraction maps and related fixed point theory can be found in Goebel & Kirk (1990) and Cheney (2013). For more on fixed point methods for order preserving operators, see, for example, Guo et al. (2004), Zhang (2012), Marinacci & Montrucchio (2019) or Deplano et al. (2020).

Footnotes
  1. Hint: Theorem 5.1.1 is self-improving: it implies this seemingly stronger result. The proof is not trivial but see if you can get it started. You might like to note that FkF^k has a unique fixed point u∗u^* in SS. (Why?) Now consider the fact that ∥Fu∗−u∗∥=∥FFnku∗−u∗∥\| F u^*- u^* \| = \| F F^{nk} u^*- u^*\| for all n∈Nn \in \NN.

  2. Here’s a hint, if you get stuck: show that r(Π)<1r(\Pi) < 1 implies TT is an eventual contraction on [0,x][0, x] using one of the propositions presented in this chapter.

  3. Hint: every stochastic matrix has at least one stationary distribution.

References
  1. Davey, B. A., & Priestley, H. A. (2002). Introduction to lattices and order. Cambridge university press.
  2. Du, Y. (1990). Fixed points of increasing operators in ordered Banach spaces and applications. Applicable Analysis, 38(01–02), 1–20.
  3. Elliott, M., Golub, B., & Jackson, M. O. (2014). Financial networks and contagion. American Economic Review, 104(10), 3115–3153.
  4. Shin, H. S. (2010). Risk and liquidity. Oxford University Press.
  5. Battiston, S., Puliga, M., Kaushik, R., Tasca, P., & Caldarelli, G. (2012). Debtrank: Too central to fail? financial networks, the fed and systemic risk. Scientific Reports, 2(1), 1–6.
  6. Bardoscia, M., Battiston, S., Caccioli, F., & Caldarelli, G. (2015). DebtRank: A microscopic foundation for shock propagation. PloS One, 10(6), e0130406.
  7. Demange, G. (2018). Contagion in financial networks: a threat index. Management Science, 64(2), 955–970.
  8. Bardoscia, M., Barucca, P., Codd, A. B., & Hill, J. (2019). Forward-looking solvency contagion. Journal of Economic Dynamics and Control, 108, 103755.
  9. Jackson, M. O., & Pernoud, A. (2019). Investment Incentives and Regulation in Financial Networks [Techreport]. SSRN 3311839.
  10. Jackson, M. O., & Pernoud, A. (2020). Credit freezes, equilibrium multiplicity, and optimal bailouts in financial networks [Techreport]. arxiv 2012.12861.
  11. Amini, H., & Minca, A. (2020). Clearing financial networks: Impact on equilibrium asset prices and seniority of claims. Tutorials in Operations Research, 154–175.
  12. Jackson, M. O., & Pernoud, A. (2021). Systemic risk in financial networks: A survey. Annual Review of Economics, 13, 171–202.
  13. Eisenberg, L., & Noe, T. H. (2001). Systemic risk in financial systems. Management Science, 47(2), 236–249.
  14. Stachurski, J. (2022). Systemic Risk in Financial Systems: Properties of Equilibria.
  15. Klages-Mundt, A., & Minca, A. (2021). Optimal intervention in economic networks using influence maximization methods. European Journal of Operational Research.