the problem of a single particle interacting with its environment is encountered frequently in condensed matter physics , and in many variations .
such scenarios appear in the form of genuine defects in materials , like magnetic impurities in metals and the associated ubiquitous kondo problem.@xcite impurity problems are also often used as a stepping stone to gain insight into the behavior of strongly correlated systems , such as studying a single hole doped into a mott insulator helps understanding some properties of high - temperature superconductors.@xcite    in recent years , impurity problems have also started to gain some attention in the context of ultracold atomic gases . amongst others ,
the fermionic polaron problem has been studied thoroughly in connection with strongly imbalanced fermi gases close to a feshbach resonance,@xcite where a single `` spin - down '' atom interacting with a fermi sea of `` spin - up '' atoms appears as a polaronic quasiparticle and eventually forms a two - particle bound - state with one of the majority atoms if the interactions are strong enough.@xcite more generally , such polaronic problems appear when we consider the nature of the threshold of spectral functions of gapped excitations above many - body ground states.@xcite    here we are interested in a seemingly similar problem , where a mobile impurity is coupled to a continuum of gapless excitations at a quantum critical point , instead of gapless particle - hole excitations at a fermi surface .
in particular , we want to study the situation where bosons in an optical lattice are tuned to the superfluid - mott insulator transition,@xcite and a single atom of a different species is coupled to the bosons via a density - density interaction .
this scenario is especially interesting due to the recent experimental implementation of quantum - gas microscopes,@xcite which allow tracking the motion of the impurity through the bosonic bath .
indeed , a recent combined experimental and theoretical work has addressed a very similar situation in an effectively one - dimensional system.@xcite we are going to focus on the other experimentally relevant , two dimensional situation instead . from our perspective , this case also happens to be the most interesting scenario , because in @xmath0 spatial dimensions the critical behavior at the sf - mi transition is described by a non - trivial fixed point without well - defined quasiparticle excitations.@xcite recently there has been a renewed interest in this particular transition as experiments succeeded in observing the so - called higgs amplitude mode.@xcite    the main motivation for our work is the prospect of studying aspects of quantum critical transport in a clean , well - defined experimental setting . calculating transport coefficients of systems close to a quantum phase transition remains a challenging problem in theoretical condensed matter physics .
the main reason is that all traditional methods for calculating transport properties rely on a quasiparticle picture of the low - energy excitations above the ground - state , which breaks down at a quantum critical point .
indeed , strongly coupled conformal field theories ( cft ) that describe physical properties of physical systems close to quantum critical points typically have no well defined quasiparticle excitations .
some progress has been made recently using the so called ads / cft correspondence , which circumvents this problem by mapping the strongly coupled cft to a dual , weakly coupled gravitational theory , where transport coefficients can be computed reliably.@xcite    in this paper our main interest is to calculate the diffusion constant @xmath1 of the impurity as a function of temperature in the quantum critical regime , shown in fig .
[ fig : phase ] . this diffusion constant can be measured directly in experiments with quantum - gas microscopes by evaluating the mean square displacement @xmath2 of the impurity atom from its initial position after a time @xmath3 @xmath4 where the brackets denote an average over many experimental realizations and the time @xmath5 has to be larger than the typical collision time @xmath6 , below which the impurity propagates ballistically .
several theoretical works have addressed similar questions using the holographic correspondence in the context of a heavy quark moving through a quark gluon plasma.@xcite here we will use a conventional method based on the boltzmann equation , however , which turns out to provide an adequate description at low temperatures .
the quantum criticality of the bulk critical point has a single characteristic time which controls transport and relaxational processes close to equilibrium @xmath7.@xcite we will show here that the situation with impurity motion is not as universal , and the impurity time scale is determined by an interplay between several couplings which are formally irrelevant at the critical point .
one of the these couplings is the inverse effective mass , @xmath8 , of the mobile impurity ; there is the associated energy scale @xmath9 , where @xmath10 the velocity of the bulk bosonic excitations at the critical point .
the value of @xmath11 controls the quantum wave packet spread of the impurity at short times .
the other `` irrelevant '' parameters couple the impurity to the bulk excitations , and associated energy scales determine the diffusion of the impurity at long times even at the lowest temperatures .
one of our interesting results is that the diffusion constant of the impurity does _ not _ depend on @xmath11 at low temperatures , however .
the rest of the paper is organized as follows . in section [ sec
: ft ] we derive an effective low energy theory for a mobile impurity coupled to bosons at the superfluid  mott - insulator ( sf - mi ) transition and we show that the presence of three - body interactions changes the theory significantly at the particle - hole symmetric point . in section [ sec:2p ]
we study perturbative self - energy corrections to the impurity propagator and calculate its diffusion constant as a function of temperature , neglecting three - body interactions .
we repeat these calculations including three - body interactions in section [ sec:3p ] and show that this changes the results qualitatively at low temperatures .
finally we perform a renormalization group analysis of the effective field theory in sec .
[ sec : rg ] .
we start by deriving the effective low - energy theory for a bose - hubbard model at the sf - mi critical point in @xmath0 spatial dimensions , following ref .  .
a convenient starting point is to consider particle- and hole- excitations on top of a mott insulating state , which are described by the lagrangian @xmath12 here the fields @xmath13 and @xmath14 represent particle- and hole excitations on top of a state with an integer filling of bosons .
the respective energy gaps for these excitations are denoted by @xmath15 and @xmath16 , their effective masses by @xmath17 and @xmath18 .
the last term in eq .   creates or annihilates particles and holes in pairs , as required by the conservation of the total number of bosons , and we do not explicitly show higher order terms in the fields @xmath19 and @xmath20 that are allowed by symmetry .
ultimately we are interested in the experimentally relevant situation where the transition happens at a constant density , i.e.  at the tip of the mott - lobe .
this happens only if the gaps @xmath15 and @xmath16 are equal , such that particles and holes condense at the same time when the gap closes .
for the moment we will consider the general case , however , and parametrize the gaps as @xmath21 now we add a single , mobile impurity atom to the bosons , the dynamics of which is described by the free particle lagrangian @xmath22 where the field @xmath23 represents the impurity and the chemical potential @xmath24 has to be adjusted such that the density of impurities is zero .
this is important only at finite temperatures , however , where we need to ensure that no artificial , thermally excited impurities exist . at zero temperature
we can safely set @xmath25 .    within a microscopic bose - hubbard like model
the coupling between the impurity and the bosons is an onsite density - density interaction , which is a pure two - body interaction in the simplest case .
since the local boson density is given by the difference between particle- and hole- densities ( up to the constant average number of bosons that we neglect in the following , as it only renormalizes the impurity s chemical potential ) , the interaction term in our effective field theory takes the form @xmath26 where the subscript @xmath27 of the coupling constant @xmath28 indicates that this term derives from a two - body interaction .
furthermore we include a term that descends from a three - body interaction @xmath29 because in the case of particle - hole symmetry ( @xmath30 ) this interaction generates couplings that are _ more _ relevant in an rg sense than those deriving from two - body interactions only , as will be shown below . to avoid confusion we note that a three - body interaction in the bose - hubbard model of the form @xmath31 , where @xmath32 denotes the number of bosons on site @xmath33 and @xmath34 is the impurity creation operator , also generates a term of the form in addition to
. for brevity we will always refer to the term in as the three - body interaction , however .
eventually our system is described by the total lagrangian @xmath35 in order to arrive at the more familiar low energy theory for the sf - mi critical point we define the fields @xmath36 and integrate out the field @xmath37 . after a rescaling of @xmath38 and defining @xmath39 as well as @xmath40
we arrive at our final result @xmath41 the terms shown are the leading order terms within a gradient expansion . for the case of particle - hole symmetry ( @xmath42 , i.e.  @xmath43 ) the dynamical critical exponent of the bosons is @xmath44 and the critical point is described by the well known @xmath45 cft in @xmath46 dimensions .
even more interesting , the direct interaction term between the bosons and the impurity @xmath47 vanishes if only two - body interactions are present .
this somewhat counterintuitive result ( in the sense that the particle - hole symmetric interaction term vanishes if the @xmath48-theory is particle - hole symmetric ) can be understood from the fact that the two - body interaction term @xmath49 is invariant under a particle - hole transformation ( @xmath50 ) _ only _ in combination with @xmath51 .
this is not true for the three - body interaction term @xmath52 , however , and that is why the parametrically smaller @xmath53 generates a term that is more relevant than the one generated by two - body interactions at the particle - hole symmetric point .
the leading order terms that are consistent with the symmetry
@xmath54 = \mathcal{l}[\psi^ * ; -\delta ,- u_2,u_3 ] \label{syms}\ ] ] are the ones shown in eq .  .
also note that the two - body interaction @xmath28 couples the impurity density @xmath55 to the time - like component of the bosonic current - density @xmath56 .
we stress again that we focus solely on the particle - hole symmetric case ( @xmath30 ) in the rest of the paper .
a tree - level scaling analysis of the lagrangian using the exact critical exponents of the bosonic cft shows that both interactions @xmath28 and @xmath53 are irrelevant in @xmath0 spatial dimensions at the p / h - symmetric point .
indeed , performing a scaling with @xmath44 gives the scaling dimensions @xmath57 = ( d-1+\eta)/2 $ ] and @xmath58 = d/2 $ ] for the fields , as well as @xmath59 = d+1 - 1/\nu$ ] and @xmath60 = d$ ] for the composite bosonic density- and current operators per definition . here
@xmath61 and @xmath62 denote the anomalous dimension and the correlation length exponent of the bosonic cft , which take the values @xmath63 and @xmath64 at the quantum critical point in @xmath0 spatial dimensions.@xcite from these scaling dimensions we get @xmath65 = 1-d$ ] and @xmath66=1/\nu - d$ ] , thus both couplings are irrelevant in @xmath0 .
also note that the curvature of the impurity dispersion is an irrelevant perturbation , as the scaling dimension of the impurity mass @xmath11 is @xmath67=1 $ ] .
we will perform a more elaborate rg analysis in sec .
[ sec : rg ] .
in this section we analyze the field theory at the particle - hole symmetric point ( @xmath30 ) , assuming that three body interactions are absent ( i.e.  @xmath68 ) .
the lagrangian thus takes the simplified form @xmath69 note that we dropped the chemical potential since it is not important at zero temperature , as argued above .      as a first step towards an understanding of the impurity s properties we calculate its self - energy perturbatively in the interaction strength @xmath28 at zero temperature in the following .
since @xmath28 is irrelevant , we expect the perturbative calculation to give qualitatively correct results .
the terms up to second order are shown in fig .
the first order tadpole contribution turns out to be identically zero .
note that the boson loop in the second order diagram corresponds to the time - like component of the bosonic current - current correlation function of the free bosonic theory due to the particular form of the interaction vertex in . at @xmath70 the form of the bosonic current - current correlator of the interacting @xmath48-theory at the critical point in @xmath0
is fixed by conformal invariance and is given by @xmath71 here @xmath72 is a numerical factor that is renormalized by the boson interaction @xmath73 and takes the value @xmath74 for the free bosonic theory .
the functional form of @xmath75 is independent of @xmath73 , however .
for this reason our results for the impurity self - energy below are perturbative in @xmath28 , but valid to arbitrary order in the boson interaction @xmath73 .    at zero temperature
the imaginary part of the retarded self - energy up to second order in @xmath28 is given by @xmath76 where @xmath77 denotes the impurity dispersion .
at @xmath78 the integral can be evaluated exactly and we obtain @xmath79 \notag \\ & = & - \frac{u_2 ^ 2 \kappa}{3 \pi c^4 } \ , \omega^3 + \mathcal{o}(\omega^4)\end{aligned}\ ] ]    for finite momenta we can derive approximate analytic expressions , which take the form @xmath80 with @xmath81 .
the qualitative change of the self - energy at @xmath82 can be understood from kinematic considerations . the only lifetime limiting process for the impurity at zero temperature happens to be the excitation of a bosonic `` phonon '' mode . due to energy- and momentum conservation this process
is prohibited as long as the group velocity of the impurity @xmath83 is smaller than the `` sound '' velocity @xmath10 of the bosons .
for this reason the impurity spectral function has a sharp delta - function peak for all momenta @xmath84 , which is reflected by the fact that the imaginary part of the self energy scales as @xmath85 close to the quasiparticle pole . on the other hand , for @xmath86 the impurity can scatter strongly by exciting the bosons .
consequently the imaginary part takes a non - zero on - shell value , implying that there is no well defined impurity excitation anymore .
the appearance of this kinematic constraint can be readily seen from eq .  , where the imaginary part of the current correlator is non - zero only if @xmath87 .
it is important to note that this threshold behavior for creating bosonic excitations at @xmath88 holds to all orders in perturbation theory in the impurity - boson coupling @xmath28 .
higher order diagrams with more boson propagators describe processes where the impurity excites multiple boson modes , which in turn have an even higher energetic threshold . from a theoretical perspective
it is an interesting question , if the impurity s properties above this kinematic threshold at @xmath82 are described by a rg fixed point .
this is not the case , however , and we will come back to this problem in sec .  [
sec : rg ] .
an important transport coefficient that can be directly measured in experiments is the diffusion constant @xmath1 of the impurity .
here we calculate @xmath89 as a function of temperature @xmath90 , impurity mass @xmath11 and the mott gap @xmath91 using the boltzmann equation , which gives qualitatively correct results as long as the impurity has a well defined quasiparticle peak .
this condition is satisfied in the temperature regime @xmath92 ( we set @xmath93 from now on ) , where the kinematic threshold discussed in the previous section does nt matter .
let us quickly estimate the experimentally relevant temperature regime before proceeding with the boltzmann calculation .
typical temperatures in experiments@xcite in units of the hubbard interaction @xmath94 are on the order of @xmath95 . using the theoretical results for the hopping amplitude @xmath96 as well as the sound velocity @xmath10 at the critical point of the 2d bose - hubbard model@xcite , which take the values @xmath97 and @xmath98 ,
together with the band mass @xmath99 of the impurity , we obtain @xmath100 , which is typically well below the kinematic threshold .
note that @xmath101 is a relatively large energy scale .
we are interested in a situation without external forces but with an initial density gradient , thus the boltzmann equation takes the form @xmath102 \ , \label{boltzequ}\ ] ] with @xmath103 as the wigner distribution function of the impurity at momentum @xmath104 , position @xmath105 and time @xmath3 .
the collision integral has the standard form @xmath106 = - \sum_{\k ' } \big ( w_{\k,\k ' } f_\k - w_{\k',\k } \ ,
f_{\k ' } \big ) \ , \label{collint}\ ] ] where the transition rates @xmath107 can be calculated using fermi s golden rule and are given by @xmath108   \ . \label{wkk}\end{aligned}\ ] ] here @xmath109 denotes the time - like component of the retarded bosonic current - current correlation function , @xmath110 is the bose - einstein distribution function and @xmath111 denotes the unit step function .
details of the derivation are given in appendix [ app : boltz ] .
in typical experimental situations the impurity is initially localized , released at some time @xmath112 and then is allowed to propagate through the bosonic bath . in principle
the boltzmann equation allows to calculate the full time evolution of the impurity distribution function . at short times collisions
do nt play a role and the impurity propagates ballistically @xmath113 where @xmath6 denotes the typical time between collisions .
if we start at @xmath112 with the particles localized in a gaussian wavefunction @xmath114 then the short time evolution of the wigner distribution function according to eq .
is free particle behavior @xmath115    at long times @xmath5 collisions dominate and the impurity propagates diffusively .
the wigner distribution then takes the form @xmath116 with a maxwell - boltzmann distribution of the momenta . indeed , in the long time limit , where the wigner distribution function of the impurity is close to equilibrium
, we can derive a diffusion equation from the boltzmann equation for arbitrary linear collision integrals of the form .
again , details can be found in appendix  [ app : boltz ] .
the corresponding impurity diffusion constant @xmath1 takes the simple form @xmath117 with a transport relaxation time @xmath6 which is given by @xmath118 f^\text{boltz}_\k \ , \label{reltime}\ ] ] where @xmath119 and @xmath120 is the equilibrium boltzmann distribution function . using the explicit form of the transition rates @xmath107 in eq .
we can cast the expression for the diffusion constant @xmath89 in a scaling form @xmath121 the corresponding scaling function @xmath122 can be calculated numerically and is shown in fig .  [
fig : phij3d ] .
here we used the one loop result for the bosonic current correlator @xmath123 at finite temperature , where the dispersion relation of the bosonic modes @xmath124 acquires a temperature dependent mass term @xmath125 which follows from a @xmath126 expansion in the limit @xmath127 of the bosonic cft.@xcite note that this temperature dependent gap is an artifact of the @xmath127 limit and it is not the scope of this paper to go beyond this limitation .    at a fixed temperature @xmath90 , the diffusion constant @xmath1 takes its minimal value in the quantum critical regime directly above the quantum critical point , where the mott gap vanishes ( @xmath128 ) .
this is consistent with the naive expectation that the phase space for scattering processes between the impurity and the bosonic modes is large , if the gap of the bosonic modes is small .
the scaling function in the quantum critical regime for @xmath128 is shown in more detail in fig .
[ fig : phij ] .
note that @xmath122 approaches a finite value @xmath129 at zero temperature and scales linearly for small @xmath130 , as @xmath131 interestingly , the fact that the scaling function takes a constant value at @xmath70 and @xmath128 implies that the diffusion constant at the critical point does not depend on the impurity mass @xmath11 at very low temperatures .
indeed , the curvature of the impurity dispersion , _ i.e. _  the inverse mass @xmath132 , is an irrelevant perturbation , and here there is no conspiracy between the two irrelevant couplings @xmath132 and @xmath28 that potentially could change the @xmath90 dependence of @xmath1 from that expected from its @xmath133 dependence . only at higher temperatures the mass starts to play a role , leading to a crossover from a @xmath134 behavior at very low temperatures to a @xmath135 dependence at temperatures @xmath136 .
note that this crossover temperature is around the typical temperatures in experiments , where @xmath137 , thus we expect that the @xmath135 behavior at larger temperatures should be easily accessible .    in order to compare our results directly to experiments it would be natural to express the diffusion constant in units of the hopping amplitude @xmath96 ( or more precisely @xmath138 , if the lattice constant @xmath139 and planck s constant are not set to unity ) .
assuming that the hubbard interaction @xmath28 between the impurity and the bosons equals the hubbard-@xmath94 in the bulk , the diffusion constant at a typical temperature @xmath140 above the critical point takes the value @xmath141 .
the corresponding transport relaxation time @xmath6 takes the value @xmath142 , showing that it is experimentally promising to reach the diffusive regime at times @xmath5 .
we note here that in certain parameter regimes a direct comparison of experimental results with our calculations might be difficult due to the harmonic confinement in experiments .
this confinement effectively leads to an averaging of the diffusion constant when the impurity explores different regions of the trap . in most cases
this is not a problem , however . only in the mott regime at temperatures well below the mott gap , where the diffusion constant is large and the impurity propagates almost ballistically , it explores larger regions of the trap and averaging
could become important .
such confinement problems can be circumvented by increasing the temperature slightly , since @xmath1 and @xmath6 are strongly temperature dependent and decrease by an order of magnitude if the temperature is raised from @xmath143 to @xmath144 .
on a formal level , the fact that the scaling function at the critical point @xmath145 is constant can be traced back to the behavior of the current correlator at small frequencies , which scales as @xmath146 for @xmath147 .
it is important to stress that this behavior , which is related to the fact that @xmath148 is an odd function of frequency and expected to be an analytic at @xmath149 , is likely valid in general and not related to our one - loop approximation of the current correlator .
indeed it has been argued that the current - current correlation function takes a hydrodynamic form for @xmath150 @xmath151 the imaginary part of which obviously scales as @xmath152 for small frequencies .
this has also been confirmed by calculations using the ads / cft correspondence.@xcite    in the mott regime away from the critical point , where the excitation gap @xmath91 is finite , the scaling function diverges exponentially as the temperature approaches zero @xmath153 this behavior can be understood from the fact that the impurity does not scatter with bosonic excitations at temperatures well below the mott gap and thus propagates ballistically , which translates to a diverging diffusion constant .
we also note that the mobility @xmath154 of the impurity is independent of its mass at the critical point at low enough temperatures .
the mobility is defined as @xmath155 , where @xmath156 is the impurity s terminal drift velocity in response to an external force @xmath157 .
it is related to the diffusion constant @xmath1 via the einstein relation @xmath158 and thus has the same dependence on the impurity mass as the diffusion constant .      in the superfluid phase
the impurity diffusion constant @xmath1 can be calculated in the same way as in sec .
[ sec : boltz2 ] , using a @xmath126 expansion in the limit @xmath127 for the bosonic current correlation function @xmath123 .
the only difference is that eq .   for the boson mass term
has to be replaced by@xcite @xmath159 in contrast to the mott regime , where the energy scale is set by the mott gap @xmath91 , the energy scale @xmath160 in the superfluid phase is set by the helicity modulus at zero temperature @xmath161.@xcite    in the @xmath126 expansion , the helicity modulus is non - zero only at @xmath70 , and the o(@xmath162 ) symmetry is fully restored at any non - zero temperature .
however , for the @xmath163 case of interest here , the helicity modulus is non - zero for a range of low @xmath90 : the helicity modulus @xmath164 is proportional to the superfluid density @xmath165 , and close to the quantum critical point , @xmath166 is a universal function of @xmath167 and @xmath90 .
nevertheless , we expect the @xmath126 expansion to provide a reasonable description of the impurity atom motion , because the latter couples only to quantities which are invariant under the u(1 ) global symmetry of the superfluid . furthermore , the gap in eq .
( [ mbsuperfluid ] ) is much smaller than @xmath90 , and so is washed out by thermal excitations . in the following
we restrict our discussion to the low temperature regime , where the superfluid density is essentially constant and independent of temperature .
again , we can define a scaling function @xmath122 for the impurity diffusion constant , which takes the same form as eq .  .
the scaling function in the superfluid regime is shown in fig .
[ fig : phijsf ] , together with the results at the critical point and in the mott regime .
note that the temperature dependence of the diffusion constant in the superfluid phase is qualitatively the same as at the critical point , with the only difference that @xmath1 takes numerically slightly smaller values in the superfluid phase than at the critical point .
this can be understood from that fact the impurity couples efficiently to gapless phonons in the superfluid .
we discussed in sec .  [ sec : ft ] that three - body interactions are more relevant than two - body interactions at the particle - hole symmetric quantum critical point .
indeed , the behavior of the impurity s transport coefficients will be dominated by three - body interactions at low enough temperatures , albeit the three - body interaction being parametrically smaller than the two - body interaction . in this section
we thus repeat the analysis of sec .
[ sec:2p ] for the case of three - body interactions @xmath53 . for
the sake of simplicity we set the two - body interaction to zero in this section and start our discussion from the lagrangian in eq .
with @xmath30 and @xmath168 . in reality
both interactions will be present , however , which will lead to a crossover in the temperature dependence of the diffusion constant .
the associated crossover scale depends on the relative strength of the two interactions , which will be discussed briefly at the end of this section .
again , we start by performing a perturbative analysis of the impurity s self - energy at zero temperature , including terms up to second order in @xmath53 .
the two corresponding self - energy diagrams are analogous to the ones shown in fig .
note also that mixed terms @xmath169 vanish by symmetry .
we neglect the first order tadpole contribution as it only renormalizes the chemical potential , which we set to zero in order to ensure a vanishing density of impurities . the second order contribution to the self - energy
can be computed by replacing the current - current correlator @xmath75 in eq .   with the bosonic polarization function @xmath170 .
again we want to obtain a result for the impurity self - energy which is non - peturbative in the boson coupling @xmath73 , thus we use the scaling form @xmath171 for the bosonic polarization function at zero temperature , which follows from the scaling dimensions of the fields discussed in sec .  [
sec : ft ] .
the corresponding integral to eq .
for the imaginary part of the impurity self - energy can be evaluated exactly at @xmath78 . in the limit of small frequencies
@xmath172 we obtain @xmath173 using the exact value for the correlation length exponent@xcite @xmath174 , the imaginary part of the self - energy scales as @xmath175 at low frequencies , i.e.  the impurity spectral function again has a sharp delta - function peak .
note that this is not true in mean - field theory , where @xmath176 and the the imaginary part of the self - energy scales as @xmath152 at small frequencies , which would result in an algebraic singularity instead of a delta - function peak in the spectral function .
at finite momenta we expect a similar behavior of @xmath177 as described in sec .
[ sec : pert2 ] for the case of two - body interactions , where the self - energy acquires a finite imaginary part beyond the kinematic threshold at @xmath178 .
we calculate this imaginary part explicitly in sec .
[ sec : rg ] .      for the computation of the impurity s diffusion
constant the same arguments hold as in sec .  [
sec : boltz2 ] and we restrict our discussion here to the critical point where the mott gap vanishes ( @xmath128 ) .
the transition rates @xmath107 have the same form as in eq .  , with @xmath179 replaced by @xmath180 . using the scaling form of the bosonic polarization function at finite temperature @xmath181
we can express the impurity diffusion constant in the scaling form @xmath182 unfortunately we can not reliably calculate the scaling function @xmath183 in the equation above , because neither the mean - field approximation with @xmath184 , nor the leading order @xmath126 approximation where @xmath185 are good approximations for the polarization function at finite temperature .
in fact , from eq .
[ bospolariz ] one can see that @xmath186 is close to being constant at zero temperature , because the correlation length exponent almost takes the value @xmath187 at the critical point .
however , it is reasonable to assume that the scaling function @xmath188 in eq .
will again take a constant value at @xmath189 , because @xmath190 is an odd function of frequency and expected to scale linearly @xmath191 for small frequencies , similar to the case discussed in sec .
[ sec : boltz2 ] .
based on this assumption the diffusion constant scales as @xmath192 at low temperatures .
if both , two - and three - body interactions are present , there will be a crossover in the temperature dependence of the diffusion constant .
summing up the two inverse relaxation times using matthiessen s rule , we obtain @xmath193 with a constant of order one .
we thus expect a crossover from a @xmath194 to a @xmath195 behavior at temperatures @xmath196 , provided that the mass - dependent crossovers happen at a higher temperature .
the simple tree level - scaling at the end of sec .  [ sec : ft ] showed that both @xmath28 and @xmath53 are irrelevant couplings at the particle - hole symmetric point ( @xmath30 ) , with @xmath53 being less irrelevant than @xmath28 . in this section
we perform a more elaborate zero - temperature rg analysis of the three - body interaction @xmath53 .
our aim is to understand how a finite imaginary part of the self - energy is generated beyond the kinematic threshold discussed in sec .
[ sec : pert2 ] , thus we are going to analyze the flow of the impurity propagator at a fixed external momentum .
the two - body interaction @xmath28 will be set to zero throughout this section .
we are using a functional rg formulation to derive our flow equations.@xcite in order to keep the analysis transparent we use simple parametrizations for the propagators and vertices , however .
since the coupling between the impurity and the bosons is a four point coupling and the frg flow - equations have a one loop structure , we employ the standard trick and consider the flow in the symmetry broken phase in order to generate a non - trivial frequency- and momentum - dependence for the impurity propagator while maintaining a simple parametrization of the interaction vertices.@xcite this approach is similar to the background field method used in field - theoretic rg approaches .
our starting point is to choose a simple truncation of the effective action @xmath197 $ ] at momentum scale @xmath198 in the symmetry broken phase @xmath199 \notag \\   & & + \frac{g_\lambda \phi_\lambda}{2 } \sum_{q , k } \big ( \psi^*_q \psi_{q - k } \psi_k + \text{c.c . }
\big ) \notag \\   & & + \frac{g_\lambda}{4 } \sum_{q , k , p } \psi^*_k \psi^*_{q - k } \psi_{q - p } \psi_p   + \sum_q b^*_q \big ( -b_{\l , q}^{-1 }   \big ) b_q   \notag \\ & & +   \sum_{q , k , p }   u_\lambda(\omega_q ) \ ,
\psi^*_{k } b^*_{q - k } b_{q - p } \psi_{p } \notag \\   & & + \phi_\lambda \sum_{q , k }   u_\lambda(\omega_q ) \left ( b^*_q b_{q - k } \psi_k + \text{c.c . }
\right )   \ . \label{gammaeff}\end{aligned}\ ] ] here @xmath200 is the expectation value of the bosonic field @xmath38 , the subscript @xmath201 denotes scale dependent quantities , we use the shorthand notation @xmath202 for sums over bosonic matsubara frequencies and momenta , and we dropped the subscript on the three - body coupling @xmath203 for notational brevity . note that even though @xmath200 is zero at the critical point , the properly rescaled field expectation value @xmath204 will be non - zero and is related to the anomalous dimension of the bosons .    the effective action in eq .
is a straightforward generalization of eq .   for @xmath30 and @xmath168 to the symmetry broken phase .
the important renormalizations appear in the propagators , which we parametrize as ( we set the velocity of bosonic excitations to unity ( @xmath205 ) from now on ) @xmath206 where @xmath207 and @xmath208 are wave - function renormalizations and @xmath209 denotes the sign function .
note two important modifications of the impurity propagator @xmath210 .
first , we expanded the dispersion relation around a fixed momentum @xmath211 and its corresponding on - shell frequency @xmath212 in order to study the flow of the propagator at a fixed external momentum .
the group velocity @xmath213 is allowed to flow , which in turn captures renormalizations of the impurity s effective mass .
second , we include a damping term @xmath214 , which will flow to non - zero values above the kinematic threshold . as will become clear later , generating a finite imaginary self - energy @xmath214 is the only possible way to cure the mass - shell singularity that arises from the large phase space for scattering processes between the ( almost ) linearly dispersing impurity mode and the linearly dispersing critical bosonic modes above the kinematic threshold .    also note that we did nt include a chemical potential in the impurity propagator . during all our calculations we enforce a vanishing density of impurities by neglecting the contribution of the impurity propagator s pole ( or branch cut , if @xmath215 ) in all diagrams .
this procedure is tantamount to adjusting a chemical potential such that the density of impurities is zero .
apart from the impurity propagator also the three - body coupling @xmath216 in will acquire an important non - analytic frequency dependence during the flow above the kinematic threshold , which can be parametrized as @xmath217    flow equations for all scale dependent quantities can now be derived in the usual way from the exact flow equation for the effective action@xcite @xmath218 = \frac{1}{2 } \text{tr } \left [ \frac{1}{\gamma^{(2)}_\lambda[\phi ] + r_\lambda } \ , \partial_\lambda r_\lambda \right ] \label{wequ}\ ] ] here @xmath219 $ ] is shorthand for the second functional derivative of the effective action with respect to the fields , i.e.  @xmath220)_{kq } = \delta^2 \gamma /(\delta \phi_{k,\alpha } \delta \phi_{q,\beta})$ ] has a @xmath221 matrix structure coming from the field index @xmath222 with @xmath223 , besides its usual frequency / momentum - space matrix structure .
the trace in eq .   is taken with respect to all momentum , frequency and field indices .
@xmath224 denotes the regulator function @xmath225 and we choose an optimized litim cutoff for both , the bosons as well as the impurity@xcite @xmath226 as customary , we neglect derivatives of @xmath227 and @xmath228 in the scale derivative of the regulator @xmath229 .    the explicit form of our flow equations for all scale - dependent parameters is shown in appendix [ app : a ] .
note that the flow equations for the four boson parameters , @xmath230 , @xmath231 , @xmath232 and @xmath233 are decoupled from the impurity flow equations .
this is because we only consider the case of a vanishing impurity density where all diagrams with internal impurity loops evaluate to zero and thus do not give a contribution to the flow equations .
for this reason the properties of the bosons at the critical point are described by the usual wilson - fisher fixed point with critical exponents in the universality class of the 3d xy - model .    without discussing the standard flow equations for the bosons much further
, we now turn to the impurity flow equations , shown in fig .
[ figflow ] in a very schematic way .
the equations for the impurity propagator and the three - point vertex are given by    @xmath234   \ , \label{flowequb } \\ \partial_\l ( u_\l \phi_\l ) & = & - u_\l g_\l \phi_\l \sum_p \big [ 2 \ , \mathcal{g}_p^3 + 3 \ , \mathcal{g}_p^2 \mathcal{f}_p + 6 \ , \mathcal{g}_p \mathcal{f}_p^2+\mathcal{f}_p^3   \big ]   \ , \dot{r}^\psi_{\l , p } \notag \\ & & - 2 \ , u_\l^2 \phi_\l \sum_p \big [ ( \mathcal{g}_p-\mathcal{f}_p)^2 \mathcal{b}_{p } \,\dot{r}^\psi_{\l , p } + ( \mathcal{g}_{p}-\mathcal{f}_{p } ) \mathcal{b}_{p}^2 \ , \dot{r}^b_{\l , p } \big ] \notag \\ & & - 3 \ , u_\l^2 g_\l \phi_\l^3   \sum_p \big[2 ( \mathcal{g}_p-\mathcal{f}_p)^3 \mathcal{b}_{p } \ , \dot{r}^\psi_{\l , p } + ( \mathcal{g}_{p}-\mathcal{f}_{p})^2 \mathcal{b}_{p}^2 \ , \dot{r}^b_{\l , p } \big ] \notag \\ & & - 2 \ , u_\l^3 \phi_\l^3   \sum_p \big[(\mathcal{g}_p-\mathcal{f}_p)^2 \mathcal{b}_{p}^2 \ , \dot{r}^\psi_{\l , p } + 2 ( \mathcal{g}_{p}-\mathcal{f}_{p } ) \mathcal{b}_{p}^3 \ , \dot{r}^b_{\l , p } \big ]   \ , \label{flowequu}\end{aligned}\ ] ]    where @xmath235 is the scale derivative of the regulator and @xmath236 , @xmath237 and @xmath238 denote the regularized normal and anomalous boson green s functions as well as the regularized impurity green s function , respectively .
@xmath239    from the explicit form of the flow equations in app .
[ app : a ] it is clear that the three - body interaction @xmath240 is always irrelevant . this is mainly due to the screening of the three - body interaction by bosonic particle - hole excitations , giving rise to a contribution to the @xmath241-function of the interaction @xmath240 which is linear in @xmath240 and negative , independent of the impurity s properties .    even though @xmath240 always flows to zero , a qualitative change in the rg - flow occurs if the group velocity @xmath242 of the impurity is above the kinematic threshold for creating particle - hole excitations .
in perturbation theory this happens for @xmath243 , as discussed in sec .
[ sec : pert2 ] . in our rg procedure
the threshold velocity is slightly renormalized to @xmath244 in units with @xmath205 , where @xmath245 , @xmath246 and @xmath247 denote the fixed - point values of the rescaled variables ( see appendix [ app : a ] ) . indeed , for @xmath248 the rg generates a finite damping term @xmath214 for the impurity .
typical rg - flows of the impurity parameters in the regime @xmath248 are shown in fig .
[ fig : paramflow ] .
we do nt find a rg fixed - point that controls the flow of @xmath214 above the kinematic threshold .
instead , @xmath214 flows to a finite , non - universal value while @xmath240 flows to zero . in fig .
[ fig : gamma ] we show the value of @xmath214 at the end of the flow ( i.e.   at @xmath249 ) as a function of the initial velocity @xmath250 , showing the sharp onset of @xmath214 for @xmath248 .
formally , generating a flow for @xmath214 is intimately linked to the appearance of a non - integrable mass - shell singularity above the kinematic threshold .
in this paper we have analyzed the problem of a mobile impurity atom interacting with bosons at the superfluid - mott insulator critical point in @xmath0 spatial dimensions .
the effective low energy theory describing this systems is qualitatively different if three - body interactions are present .
we calculated the diffusion constant @xmath1 of the impurity as a function of temperature and showed that @xmath1 does not depend on the impurity mass at low temperatures .
furthermore we showed that the zero - temperature properties of the impurity change if its group velocity exceeds the kinematic threshold for creating particle - hole excitations .
this behavior is not controlled by a rg fixed - point , however .
instead the imaginary part of the impurity s self - energy flows to a non - universal value while the interaction flows to zero .
we gratefully acknowledge interesting discussions with i.  bloch , m.  endres , t.  fukuhara , c.  gross , s.  hild and in particular p.  strack .
this research was supported by the national science foundation under grant dmr-1103860 , and by u.s .
army research office award w911nf-12 - 1 - 0227 .
mp is supported by the erwin schrdinger fellowship j 3077-n16 of the austrian science fund ( fwf ) .
in this appendix we derive the collision integral for the boltzmann equation in sec .  [
sec : boltz2 ] as well as the diffusion equation , which describes the propagation of the impurity at long times . our collision integral will be perturbative in @xmath28 , but we want to make sure that effects to all order in the boson coupling @xmath73 are accounted for , similar to the perturbative calculation in sec .  [
sec : pert2 ] . this can be done without invoking the keldysh formalism by considering an effective model where the impurity is coupled to a bath of phonons . as will be shown in the following
, the spectral density of the phonons needs to be chosen such that it matches the spectral function of the bosonic current - current correlation function .
the collision integral for the boltzmann equation can then be derived straightforwardly from the effective phonon model using fermi s golden rule .
we start from the effective action describing the impurity coupled to a bath of harmonic oscillators @xmath251 where @xmath252 denotes the frequency of the phonon mode @xmath253 , @xmath254 are coupling constants and we use the shorthand notation @xmath255 . in order to match this effective action to our original action we integrate out the phonons and obtain @xmath256 where we have defined the spectral density of the phonons via @xmath257 on the other hand ,
if we formally integrate out the bosons in our original action to second order in @xmath28 , we get the same result as in eq .
[ seffj ] with @xmath258 now that we have matched the spectral density of the effective phonon bath to our original model we can straightforwardly employ fermi s golden rule and obtain the collision integral in eq .  .
close to equilibrium we can solve the boltzmann equation by first employing a fourier transform with respect to the spatial coordinate @xmath105 @xmath259 \ , \label{boltzequq}\ ] ] and making the ansatz @xmath260 here @xmath261 and @xmath262 describes the deviation of the momentum distribution from equilibrium .
this is basically a leading order expansion in spherical harmonics where the angular dependence of @xmath263 on the momentum @xmath104 only enters through the @xmath264 term .
note that the impurity density distribution is given by @xmath265 using the ansatz we can derive a diffusion equation from the boltzmann equation for arbitrary linear collision integrals of the form as follows .
taking an angular average of the boltzmann equation we obtain @xmath266 \notag \\ & = & i[f^0_k(\q ) ]   \ , \equiv \ , 0 \ . \label{equb1}\end{aligned}\ ] ] moreover , taking an angular average over the boltzmann equation multiplied by @xmath267 we get @xmath268 \ , \label{equb2}\ ] ] where we have used @xmath269 to replace @xmath270 with @xmath271 .
the rhs of the above equation can be evaluated straightforwardly for a collision integral of the form @xmath272 = \label{equb3 } \\ & & \ \ \ \   - \frac{q^2 k^2}{2 } \sum_{\k ' } \big [ w_{\k,\k ' } \tilde{f}_k - w_{\k',\k } \frac{k'}{k } \cos \theta_{\k,\k ' } \ , \tilde{f}_{k ' } \big ] \ . \notag\end{aligned}\ ] ] again , using eq .   to replace @xmath270 with @xmath271 in eq .  ,
summing over @xmath104 and neglecting the @xmath273 term , we finally arrive at the diffusion equation @xmath274 with the diffusion constant specified in eqs .   and .
the explicit expression of the scaling function for the diffusion constant in eq . is    @xmath275 e^{-z } \ , \label{scalingfunc}\end{aligned}\ ] ]    where @xmath276 is the scaling function of the current correlator @xmath277 and @xmath278 .
all of the results derived above also apply to the case of three - body interactions if the current correlator @xmath75 is replaced by the density correlator @xmath170 as well as @xmath28 by @xmath53 .
it has been shown in sec .
[ sec : rg ] that the flow equations for the boson parameters decouple from the impurity flow equations . indeed , using the rescaled variables @xmath279 and defining the anomalous dimension of the bosons via @xmath280 the flow equations for the boson parameters
are given by ( @xmath281 ) @xmath282 in order to bring the flow equations into an analytically tractable form we expanded them in @xmath283 .
this basically amounts to an expansion in a small anomalous dimension , which works well in @xmath0 spatial dimensions . at the wilson - fisher fixed point
the rescaled variables take the values @xmath284 note that the value of the anomalous dimension at the critical point @xmath285 agrees pretty well with the expected value for a 3d xy - transition , where @xmath286.@xcite for the impurity flow equations we define the rescaled variables and the anomalous dimension of the impurity as @xmath287 and the flow equations take the form @xmath288 \ , \\ \partial_s \tilde{u}_s ' & = & - \eta _ * \tilde{u}_s ' - \frac{\tilde{u}_s ' \tilde{g}_*}{16 \pi } \bigg [ \frac{3}{(1+\tilde{r}_*)^{5/2}}-\frac{115   \ , \tilde{g } _ * \tilde{\phi}_*^{2}}{4 ( 1+\tilde{r}_*)^{7/2}}\bigg]-\text{re } \bigg [ \frac{(\tilde{u}_s'+i \tilde{u}_s'')^2}{2 \pi } \bigg(\frac{\tilde{v}_s^2+i \tilde{\gamma}_s ( \alpha - i \tilde{\gamma}_s)}{\tilde{v}_s^2 \alpha^3 \big ( ( \alpha - i \tilde{\gamma}_s)^2-\tilde{v}_s^2\big)^{1/2 } } -\frac{i \tilde{\gamma}_s}{\tilde{v}_s^2 \alpha^3 } \bigg ) \bigg ] \notag \\ & & + \frac{3   \tilde{g } _ * \tilde{\phi}_*^{2}}{2 \pi } \text{re}\bigg [ ( \tilde{u}_s'+i \tilde{u}_s'')^2
\bigg(\frac{3 i \tilde{\gamma}_s ( \alpha - i \tilde{\gamma}_s)^3 - 3 \tilde{v}_s+\tilde{v}_s^2 ( 4 \alpha^2 - 9 i \tilde{\gamma}_s \alpha-6 \tilde{\gamma}_s^2)}{4 \tilde{v}_s^2 \alpha^5 \big ( ( \alpha - i \tilde{\gamma}_s)^2-\tilde{v}_s^2\big)^{3/2 } } - \frac{3 i \tilde{\gamma}_s}{4 \tilde{v}_s^2 \alpha^5 } \bigg)\bigg ]   \ , \end{aligned}\ ] ] @xmath289-\text{im } \bigg [ \frac{(\tilde{u}_s'+i \tilde{u}_s'')^2}{2 \pi } \bigg(\frac{\tilde{v}_s^2+i \tilde{\gamma}_s ( \alpha - i \tilde{\gamma}_s)}{\tilde{v}_s^2 \alpha^3 \big ( ( \alpha - i \tilde{\gamma}_s)^2-\tilde{v}_s^2\big)^{1/2 } } -\frac{i \tilde{\gamma}_s}{\tilde{v}_s^2 \alpha^3 } \bigg ) \bigg ] \notag \\ & & + \frac{3   \tilde{g } _ * \tilde{\phi}_*^{2}}{2 \pi } \text{im}\bigg [ ( \tilde{u}_s'+i \tilde{u}_s'')^2
\bigg(\frac{3 i \tilde{\gamma}_s ( \alpha - i \tilde{\gamma}_s)^3 - 3 \tilde{v}_s+\tilde{v}_s^2 ( 4 \alpha^2 - 9 i \tilde{\gamma}_s \alpha-6 \tilde{\gamma}_s^2)}{4 \tilde{v}_s^2 \alpha^5 \big ( ( \alpha - i \tilde{\gamma}_s)^2-\tilde{v}_s^2\big)^{3/2 } } - \frac{3 i \tilde{\gamma}_s}{4 \tilde{v}_s^2 \alpha^5 } \bigg)\bigg ]   \ , \\
\partial_s \tilde{v}_s & = & -\eta_b \tilde{v}_s+\frac { \tilde{\phi}_*^{2}}{2 \pi } \text{re } \bigg[(\tilde{u}_s'+\tilde{u}_s'')^2 \bigg ( \frac{\tilde{v}_s ( \alpha - i \tilde{\gamma}_s)^3+i \tilde{\gamma}_s \tilde{v}_s^3}{\tilde{v}_s^2 \alpha^3 \big ( ( \alpha - i \tilde{\gamma}_s)^2-\tilde{v}_s^2\big)^{3/2}}-\frac{1}{v \alpha^3 } \bigg)\bigg ] \ , \\
\partial_s \tilde{\gamma}_s & = & ( 1-\eta_b ) \tilde{\gamma}_s - \frac{\tilde{u}_s''}{8 \pi } \ , \frac{2 ( 1+\tilde{r}_*)-3 \ , \tilde{g } _ * \tilde{\phi}_*^{2}}{2 ( 1+\tilde{r}_*)^{5/2}}+ \frac{\tilde{\phi}_*^{2}}{4 \pi } \text{im } \bigg [ ( \tilde{u}_s'+i \tilde{u}_s'')^2\bigg ( \frac{\tilde{v}_s^2+i \tilde{\gamma}_s ( \alpha - i \tilde{\gamma}_s)}{\tilde{v}_s^2 \alpha^3 \big ( ( \alpha - i \tilde{\gamma}_s)^2-\tilde{v}_s^2\big)^{1/2 } } -\frac{i \tilde{\gamma}_s}{\tilde{v}_s^2 \alpha^3 }   \bigg ) \bigg ] \ , \end{aligned}\ ] ] where we have defined @xmath290 and fixed all boson parameters at their respective critical values .
again we have expanded boson loops in @xmath283 and moreover we have taken @xmath291 for the sake of simplicity , which amounts to neglecting the small curvature of the impurity dispersion in .
this does not change the results qualitatively .
we also mention that the term @xmath292 cancelled exactly with the @xmath293 contribution in the flow equation for the three - body interaction @xmath216 .
the particular structure with real- and imaginary - parts in the impurity flow equations originates from the non - analytic term @xmath294 in the boson propagator . after setting the external boson frequencies to zero ,
all diagrams with internal impurity lines depend on the external impurity frequency @xmath295 via the damping term @xmath296 .
we can extract the contributions of analytic and non - analytic terms by taking limits @xmath297 from above and below and adding or subtracting the resulting expressions , which gives rise to the particular form of the impurity flow equations above .
c.  p.  herzog , a.  karch , p.  kovtun , c.  kozcaz and l.  g.  yaffe , jhep * 0607 * , 013 ( 2006 ) . j.  casalderrey - solana and d.  teaney , phys .
d * 74 * , 085012 ( 2006 ) .
s.  s.  gubser , phys .
d * 74 * , 126005 ( 2006 ) .
j.  de boer , v.  e.  hubeny , m.  rangamani and m.  shigemori , jhep * 0907 * , 094 ( 2009 ) .
d.  t.  son and d.  teaney , jhep * 0907 * , 021 ( 2009 )
. s.  caron - huot , p.  m.  chesler and d.  teaney , phys .
d * 84 * , 026012 ( 2011 ) .
d.  tong and k.  wong , arxiv:1210.1580 [ hep - th ] .


the global financial crisis that started in 2007 has stimulated an extensive literature on numerous credit - related themes such as risk assessment , financial contagion , regulatory indebtedness or liquidity ratios , misuse of derivatives ( see for example  @xcite ) .
as the crisis unfolded , the understanding of the credit exposure at an individual bank level , as well as at an aggregate systemic standpoint became increasingly important .
a given bank was perceived safer or riskier based not only on its balance sheet figures or growth estimates , but also on the financial strength of the insurance companies who would step in in the case of a debtor s default , the financial distress of the peers with whom the bank had more interbank exposure , liquidity dry - ups which could trigger massive asset fire - sales , or even the assigned probability to a lender of last resort type of bailout from the central government  @xcite .
insurance companies played a key function in the network : they had the commitment to enter the scene if a given asset failed to make promptly payments .
however once it became clear that insurance companies would be not be able to compensate for every risky asset they had insured , it also became apparent that a financial distress event could propagate across the entire financial system @xcite .
banks would have to start recognizing massive asset losses , and if a bank defaulted , other interconnected institutions could follow , potentially triggering a bankruptcy cascade @xcite .
people fearing for its savings could line up and withdraw their deposits causing a so - called bank run , as was for example observed in the united states , united kingdom , iceland , spain , or brazil @xcite .
this context is pertinent to the present work , as nowadays much of the discussion continues to be how to define and regulate higher capital requirements , while at the same time weighting the trade - offs involved , namely higher costs , or lending shifting to the shadow banking ( @xcite ) .
this paper contributes to this academic discussion by using a network - based framework to examine the structure and inner dynamics behind the links between financial institutions in different countries .
the application of network theory to the financial markets @xcite gave place to a new set of methodologies to study interaction in interconnected agents , institutions or financial products ( e.g. portfolios , stock indices , derivatives ) .
however its literature is fairly recent , and to our knowledge none has studied the banking system structure using actual ( non - simulated ) data from banks financial statements across several countries . to the extent that actual interbank data ( e.g. loans , repos , swaps )
is not publicly available , we argue that relevant implications can still be derived upon the common behavior of certain key balance sheet items .
as we have learned from the crisis , two banks can be exposed to each other s risk , even if there is no financial transaction between them : if a bank fails , market conditions ( liquidity , stock prices , mergers or acquisitions expectations ) are all likely to experience cascade consequences . in general , financial institutions ,
be them commercial banks , investment banks , development banks , broker - dealers , or credit unions , fall under the supervision of each country s central bank .
these institutions , henceforth referred to as banks , must usually submit their quarterly or monthly financial statements to the corresponding central bank regulator .
this article takes advantage of that rich , homogeneous and regular available information , and uses balance sheet data from 2005 to 2010 from five emerging countries : argentina , brazil , mexico , south africa , and taiwan .
we then construct banking functional networks from the leverage dependence across financial institutions , and analyze the structure of this networks .
we find configurations that tend to group into large clusters at relatively low correlation levels . a modular structure characterized by one large bank community
, some small ones and isolated banks is also found in all countries .
in addition , we propose and simulate a model of corporate and interbank loans that generates rich and diverse balance sheet growth .
and when we construct simulated banking functional networks according to their leverage dependence , we find results that are consistent with the empirical networks . despite the absence of actual interbank contracts , these findings translate into relevant policy implications in terms of contagion and concentration , as well as
suggest potential avenues for future research .
the paper is divided as follows .
section ii describes the data and the methodology , section iii describes the main results of the paper , section iv simulates a model of interbank loans and compares with empirical data , and section v concludes .
we study the interaction between financial institutions using balance sheet data from 2005 to 2010 from five developing countries : argentina , brazil , mexico , south africa , and taiwan .
we retrieve total assets and total liabilities for each financial institution sampled at a monthly or quarterly basis , whichever is more frequent in its reporting .
let @xmath0 and @xmath1 be the total assets and liabilities , respectively , for bank @xmath2 at time @xmath3 . for each country
, the number of banks @xmath4 , fluctuates in time , and thus @xmath5 is a stochastic birth and death process . throughout the observed period 2005 - 2010
the existing banks at 2005 may fall bankrupt , merge with others , or simply survive until 2010 . on the other side
, new banks can appear in between , say in 2006 , thus modifying @xmath5 as well .
as the purpose of this work is to study the interaction between banks , we shall remove from the dataset those banks with incomplete observations .
for example , if a bank stops or starts reporting data in 2007 , its observations are removed from that country .
let @xmath6 be the banks that satisfy the previous condition , i.e. the effective number of banks studied .
the following table shows both the number of banks in 2005 and 2011 ( @xmath7 and @xmath8 ) , as well as the number of banks present throughout the time investigated ( @xmath6 ) .
table 1 also includes the number of banks births ( deaths ) , @xmath9 ( @xmath10 ) , during the period 2005 - 2011 .    .
financial system size of the five countries studied in the period 2005 - 2011 .
@xmath11one bank from brazil and one from south africa were excluded because of missing data in some time periods . [
cols="^,^,^,^,^,^ " , ]     for instance in brazil 49 banks were created and 28 disappeared between 2005 and 2011 .
only two of these new banks have gone out of the banking system before 2011 . except for brazil ,
the rest of the countries satisfy the following equation : @xmath12 we divide the countries into two groups according to their size @xmath6 , one composed by argentina and brazil , and the other composed by mexico , south africa and taiwan .
the first group represents a large size banking system with @xmath13 and the second a smaller banking system with @xmath14 . also notice the different dynamics between brazil ( also mexico ) and taiwan : whereas in the former country the number of banks increased over time , in the latter decreased .
by contrast , argentina and south africa exhibit very small fluctuations in @xmath5 for the period studied .
the construction of a network of associations out of balance sheet variables requires selecting an interaction metric that meets two general properties .
first , it should be robust to spurious correlations so that we avoid as much as possible to use a non - stationary variable that may increase or decrease simply following a global trend . and second
, it should be relevant to describe an economic process .
in other words , if two banks are associated with a given variable , it should mean that they are subject to certain equivalent market or economic conditions , and therefore may share the same risk .
we argue that balance sheet aggregates such as assets are inappropriate to define a network , and instead propose to use the leverage ratio .
consider for instance total assets , which are regarded as a proxy of bank size , and are subject to a large number of global or macroeconomic circumstances . under normal market conditions ,
the assets of a financial institution are expected to grow with time , as a result of general growth of the economy , inflation , balance sheet effects from expansionary monetary policy , population ( customer base ) growth , increased access to banking services ( banking penetration ) , higher value of the securities held by the bank ( stocks or debt ) .
it seems that almost any pair of banks , regardless whether they are contractually linked or not , will yield a positive assets correlation , and thus will incorrectly suggest an overlinked network .
we instead propose using the degree of indebtedness , henceforth referred to as the leverage ratio , as the metric on which to base the interactions .
we define the leverage as follows @xmath15    where liabilities and equity refer to the balance sheet variables of a given bank which must satisfy the usual accounting identity : @xmath16 , and as measured at book value .
therefore , the leverage ratio of bank @xmath2 , @xmath17 , indicates what fraction of equity and debt this bank is using to finance its assets .
consider for instance panel ( a ) and panel ( b ) of fig .
1 , which shows an example where the leverage increases and decreases , respectively . the former may happen , for instance , through a credit quality shock , such that the value of the existing claims decreases and the bank increases its loan loss provision .
such loan write - offs translate into a reduction in total equity , thus increasing the ratio .
the leverage can also rise if , for example , the bank simply issues debt to finance a potential acquisition or expand its credit portfolio .
what regards to panel ( b ) , a bank may deleverage by prepaying debt with its own cash , where both total assets and liabilities decrease .
it is also possible for the shareholders to conduct a capital increase , with proceeds either paying down debt or cushioning its assets .
these examples illustrate only a few of the many different scenarios why a leverage ratio may fluctuate within a given financial institution .    finally , to address the stationary concern fig .
2 depicts the median ( results remain unchanged with the average ) leverage ratio for each country as a function of time . contrary to what happens with liabilities and assets ( fig .
a.1 ) , the leverage ratio does not exhibit any strong time dependence .
each country exhibits its own leverage range , with more or less fluctuations over time .
taiwan , for example , is characterized by exhibiting a large median leverage ratio ( @xmath18 ) while argentina and brazil exhibit a lower ratio ( @xmath19 ) .
leverage levels across countries are expected to differ due to different bank capital regulation .
each country may impose its own restrictions on hedging operations , ability to buy or sell foreign currency , place debt in the national or international markets ( and at which currency ) , banks may be subject to a particular government program ( e.g. subsidized mortgage or consumer credit programs with predetermined interest rates ) , restrictions on dividend distribution
. the central bank may also have in place a monetary and/or inflation target rate , which in turn modifies the reference rate at which banks lend each other and hence the rates at which households or corporates access to credit .
countries also differ in their accounting consolidation practices or taxation .
overall these factors limit the construction of leverage - based networks at the country level , as opposed to aggregated inter - country networks .
in particular , we construct country - level networks of leverage dependence as follows .    1 .
compute the leverage time series for all banks from 2005 to 2010 in a given country . after accounting for mergers and acquisitions , banks with incomplete data points
are excluded ( either because they appeared later or disappeared sooner than the 2005 - 2010 range ) .
2 .   compute the leverage correlation matrix between all pairs of banks .
3 .   establish links between banks ( nodes ) whenever the correlation exceeds a given threshold @xmath20 .
the structure of the network is studied for different values of @xmath20 .
finally , notice that the leverage ratio also presents the important advantage of not depending on the currency , thus allowing us to use the same variable across countries .
furthermore , any difference in regulation or local practice of an improved indebtedness ratio ( e.g. capital ratio after tier 1 / 2 adjustments ) is neutral within a country ( for it applies to all ) and not relevant across countries , since this paper does not compare absolute leverage levels , but rather the interaction that can be derived from common leverage behaviors .
we next illustrate an example using data from the argentine banking system . given two time series , corresponding to the leverage evolution of two banks , @xmath21 and @xmath22 with @xmath23 ,
we compute the pearson correlation coefficient . where this coefficient takes value @xmath24 [ -1 , 1 ] , and is equal to 1 when the two time series are identical ( @xmath25 for all @xmath3 ) .
the correlation matrix is composed by the @xmath26 pairs of banks correlation coefficients ( interactions ) .
figure 3 ( a ) shows the leverage ratio evolution for six ( from a total of 75 ) different argentine banks .
three colors are used to illustrate the different behaviors observed throughout the sample .
for example the middle ( bottom ) two leverage curves in blue ( green ) correspond to two banks that exhibit an increasing ( decreasing ) leverage evolution .
another type of leverage evolution is shown in red , with a concave - shaped leverage behavior .
the dynamics of the argentine banking system is found to be very rich from an indebtedness perspective . as we shall see later , these rich dynamics lead to a complex banking network .
it also follows that banks c and d have a similar leverage increasing pattern .
the same is found between banks e and f with practically a scale factor of difference between them , and to banks a and b , although with a different type of evolution .
each of these pairs of banks exhibit a large correlation coefficient between its corresponding leverage time series , e.g. @xmath27 , and @xmath28 for the pair of banks a - b , c - d , and e - f respectively .
however , if we compute the leverage correlation between bank c ( or d ) and e ( or f ) this value will be negative , thus stating that the two banks exhibit an opposite behavior . in this last case
the four correlation coefficients range between -0.35 ( for c - e ) and -0.65 ( d - e ) .
the remaining possible interactions between the six banks fluctuate around zero .
the complete characterization of interdependence between the 75 argentine banks is captured by the correlation matrix in panel b , with each column and row representing a bank .
a special algorithm is used to order the banks for a better visualization , such that light yellow `` spots '' correspond to group of banks that share a similar behavior among them ( large @xmath29 coefficient ) , while negative correlations are represented in dark red color .    although the leverage dependence does not imply an interbank transaction ( e.g. as would be a repo agreement ) , a strong and recurring synchronization is still relevant in the understanding of concentration and systemic risk .
it is reasonable to think that a highly leveraged network where many banks are simultaneously leveraged does not imply the same contagion risk as does a network where only a small set of banks are moderately leveraged .
we can also think of the degree of negative and positive correlations as another determinant of how robust the banking system is . in this sense , identifying and studying which banks form clusters that tend to leverage ( deleverage ) all together can also provide useful policy implications at a time where stricter banking regulation is at the core of the debate ( @xcite ) .
interestingly , from a risk - sharing point of view , it would be optimal for the financial institutions not to share identical leverage dynamics , such that a shock affecting the indebtedness of a given bank is not be followed by its peers . for a banking system to remain healthy and resilient to financial contagion ,
its dynamics needs to be heterogeneous .
that is , greater heterogeneity helps to isolate individual shocks by setting a less interconnected financial system .
further , this can be understood as the system s endogenous solution to reduce its risk from financial shocks .
notice that the economic implications behind the positive and negative dependence scenarios are quite different .
a positive correlation between two banks implies that they both increase and decrease their debt - to - equity ratio simultaneously . in this case
, we would expect both banks to share ( or target ) similar growing .
alternatively , banks balance sheets could be affected by exogenous policies along the theories of balance sheet effects ( @xcite ) , bank lending channels ( @xcite ) , and bank runs ( @xcite ) .    on the other hand
, a negative correlation implies that when one bank increases its leverage , the opposite is true for the other bank , and vice versa .
a potential explanation for this case is related to the literature of fly - to - safety , i.e. in periods of financial turmoil households withdraw their deposits from a small bank and transfer them to a bigger ( perceived as safer ) institution
. thus liabilities decrease in the small bank , but increase in the other bank , resulting in a negative leverage correlation between them . as documented in @xcite , this effect was observed in brazil in late 2008 , where depositors ran from smaller banks to larger banks , the latter ones perceived as too - big - to - fail institutions with implicit guarantees from the central bank .
an alternative mechanism comes from what the literature calls the depositors discipline : depositors may punish banks with poor performance either by demanding higher interest rates or by withdrawing their deposits .
@xcite find evidence in argentina , chile , and mexico of the depositors discipline mechanism , in particular that bank deposits growth falls as risk exposure increases .
also @xcite and @xcite document that depositors favor big banks , although it is hard to differentiate whether depositors base their decisions on too - big - to - fail sentiments or bank fundamentals . also ,
@xcite report that during the 2001 argentine convertibility crisis deposit withdrawals were more pronounced in banks with higher risk taking , thus providing further evidence of opposite changes in balance sheet accounts across banks .    .
, scaledwidth=100.0% ]    we now characterize the complete argentine network .
the upper panel of fig .
4 depicts the graphical representation of the argentine financial network with @xmath30 links , or with an average degree @xmath31 equal to @xmath32 .
recall that each node represents a bank , and the @xmath30 highest banks pair correlations are connected by an undirected edge .
put alternatively , an edge ( or link ) between two nodes ( or vertices ) is added whenever its leverage pair correlation exceeds a certain threshold value @xmath20 .
the same banks showed in fig .
3 are highlighted in the network representation of fig .
4 shows clusters of banks , where a cluster is defined as a community of interconnected ( linked ) banks . and each color represents a different leverage ( temporal ) dependence . in particular , the network can be characterized by four homogeneous groups : a large cluster of banks that essentially have an increasing leverage evolution ( e.g. banks c and d ) ; three small groups with various leverage behaviors ( decreasing , concave - shaped , and a mixed of both ) ; and finally a large set of isolated nodes ( 41 from 75 ) , represented with white circles .
interestingly , we also find that the modular structure observed in the argentine banking system is similar to other countries . despite having different size , regulation , or even average indebtedness level ,
the leverage dynamics are consistent across countries .
bottom panel of fig .
4 shows the brazilian banking network using the same methodology as well as the same @xmath33 used in argentina .
it can be observed that brazil s network is similar to argentina s , namely both are characterized by a large community of banks , another set of small clusters of banks , and many isolated banks .    as the correlation threshold decreases ( and the number of links @xmath30 increases ) , both the isolated and the small - cluster nodes start merging into one large interconnected cluster .
this merger - type of behavior between different groups can be observed in fig .
5 , where the fraction of nodes that belongs to the largest cluster ( blue nodes as per fig .
4 ) is shown as a function of the correlation threshold . as we lower the threshold @xmath20 , the largest cluster starts absorbing smaller communities .
big jumps in this graph indicate that two large clusters have merged .
for example , in brazil s case a small change in the correlation threshold around 0.6 causes the two largest clusters ( each one representing approximately 40% of the nodes ) to be either merged or divided . in argentina
this similar and sudden concentration occurs at a greater correlation threshold ( @xmath34 ) .
we also find that in both countries 90% of the banks belong to a unique cluster at a relatively small correlation threshold @xmath35 , thus suggesting that the financial networks may become suddenly too concentrated .
finally , we replicate the analysis to the smaller banking systems of mexico , south africa and taiwan .
( a ) shows the leverage - based networks of these countries for an equal number of links per node , i.e. for @xmath36 .
similar network properties are observed across these countries . in particular , although taiwan s banking system is more volatile ( large fluctuations in @xmath5 ) , the banking network structure is very similar to the south african one , and both have almost the same total number of banks ( @xmath37 ) .
this result is strengthen by the similar relationship between the size of the largest cluster and the correlation threshold shown in panel b.    the mexican network , on the contrary , seems to exhibit certain differences with respect to the other two . for @xmath38 the largest cluster does not have significant weight ( right bottom panel ) , as it only accounts for 30% of the total nodes in the network .
as the correlation threshold decreases , a second cluster grows until is finally merged with the largest one ( left bottom panel ) .
thus we observe that the consolidation of the mexican banking network is characterized by the merger of two large clusters , whereas in the other countries there is only one large cluster which then absorbs either individual nodes or very small clusters .        in summary , our proposed method to examine the banking system from a leverage perspective yields networks composed by a large group of banks that are densely interconnected ( blue cluster in fig . 4 and 6 )
, a few smaller clusters , and many isolated banks . interestingly , this topology is homogeneous across countries , as depicted in fig . 4 and 6 . as
the interaction criterion becomes more and more lax ( smaller @xmath20 , in absolute value ) , banks start re - grouping into a unique financial cluster .
furthermore , the network structure is not dominated by one hub , as it occurs with many nature networks ( @xcite ) . on the contrary , these banking networks exhibit a modular structure where each module is densely connected .
only a few banks function as _ bridges _ , or connectors between different modules .
the existing structure also suggests a somewhat fragile property : a small modification to these connector banks can have sensible effects to the configuration of clusters , as well as to the isolation or propagation of a shock across the network . on the one hand , highly interconnected banks ( hubs ) are not desirable , since a failure in a hub can be spread out the network , potentially causing an aggregate system breakdown .
but on the other hand , some degree of interconnectedness is positive , as it helps to reduce the exposure of an affected node by sharing and loading off risk with its peers .
therefore we argue that the existence of dense network structures may be desirable to some extent in order to increase risk sharing and reduce hub exposure . in this sense , in order to reduce risk and increase shock resilience , regulators should recognize a basic trade - off : as connectivity and clustering increases , individual ( idiosyncratic ) risk may decrease due to better risk sharing , but global or systemic risk could increase thanks to faster channels of propagation .
the degree of optimal bank connectivity , and how can it be reached , are still active questions in the literature .
in this section we intend to contrast the empirical results from section [ sec : results ] with simulations from a model of corporate and interbank loans .
we first start by describing the set up of the model .
we propose that the corporate sector requests loans from the banking system according to a given random process , and links between banks arise due to insufficient liquidity to meet these loans .
this lending effect generates rich balance sheet and leverage dynamics that yield networks that are overall consistent with the empirical results .
our model builds on the interbank model proposed in @xcite . for simplicity
we assume perfect information in the economy , although as explained by @xcite imperfect information can result in an additional source of financial risk .
further , although we think of agents in the real sector as the implicit source of interbank links , the model does not generate balance sheet dynamics from the agents side of the economy .
see @xcite for an example of an agent - based interbank network model .
see also @xcite for a review of ace ( agent - based computational economics ) modeling in economics .
the model is structured as follows .    1
.   the banking system is composed by @xmath4 banks .
assets are defined as the sum of liquidity @xmath39 , iliquid assets @xmath40 , corporate loans @xmath41 , and interbank loans @xmath42 .
liabilities are defined as the sum of deposits @xmath43 and interbank debt @xmath44 . moreover ,
any bank @xmath45 at time @xmath3 must satisfy : @xmath46 @xmath47 2 .
let assets @xmath48 at time @xmath49 be drawn from a uniform distribution , e.g. from 5,000 to 50,000 .
then for each bank set a conservative equity - to - assets ratio selected at random from the interval 0.10 to 0.35 . and
finally , let liquidity @xmath50 at @xmath51 be a fixed share of assets , thus featuring that a given bank can be solvent to provide loans or take additional debt , but still illiquid to do so .
these steps so far are sufficient to characterize the banking system at the initial period .
3 .   next assume the corporate sector requests loans from the banking system according to a poisson process of rate @xmath52 , i.e. the length of time between loan requests is exponentially distributed .
let @xmath53 be the fixed amount of the loan , and @xmath54 be the interest rate of these loans .
each loan is supplied by only one bank @xmath45 , where @xmath45 is chosen at random with equal probability among the @xmath4 banks .
however , if bank @xmath45 is selected and @xmath55 , it lacks enough liquid assets to fulfill the entire loan . to cover the difference
, bank @xmath45 enters into the interbank market and borrows at rate _ r@xmath56 _ ( with @xmath57 ) from bank @xmath58 , where @xmath58 is also selected randomly across the @xmath59 remaining banks .
4 .   similarly ,
if @xmath60 then bank @xmath58 can not supply the interbank loan .
if this occurs then a new bank is chosen from the remaining @xmath61 banks .
as soon as @xmath53 is supplied to the corporate sector , liquidity @xmath62 as well as household deposits @xmath63 increase in a subset of banks @xmath64 by the amount @xmath65 , where @xmath66 are random weights such that @xmath67 , and @xmath68 .
the intuition is that the corporate sector distributes the loan in the economy , for example in the form of wages , and households in turn deposit these proceeds randomly across a small number of banks .
notice that this generates balance sheet ( assets ) growth , which will further enrich the leverage dynamics .
although the @xmath30 banks are drawn randomly each time a loan takes place , @xmath45 s weights are fixed but re - scaled , i.e. @xmath69 @xmath70 .
the rationale is that some banks tend to have a higher ( or lower ) share of deposits that is stable over time .
in addition , we propose that in each period there is a ( low ) probability of a negative shock to a random bank @xmath45 .
when the shock event is triggered , deposits and liquidity decrease in fixed amount , e.g. @xmath71 . 7 .
finally , the corporate and interbank loans are paid back after @xmath72 periods .
therefore at time @xmath73 bank equity @xmath74 will increase by the net amount of interest earned between the loan granted to the corporate sector and , if any , the funds borrowed from the interbank market .
contrary to the loan being initially distributed across the banking system , the funds for repayment come exogenously from the real sector . in other words , we implicitly assume the corporate sector is always productive and therefore is able to generate proceeds elsewhere .
notice that the complete set of bank interactions resulting from the model can be described in the form of a directed network . in this graph nodes
are given by the @xmath4 banks , and the directed links between nodes are given by the direction of the interbank loan .
therefore each point in time can be fully characterized by an adjacency matrix @xmath75 , where @xmath75 is a @xmath76 matrix and in general will not be symmetric .
an entry @xmath77 ( i.e. edge from @xmath45 to @xmath58 ) implies bank @xmath58 has borrowed from bank @xmath45 .
however , we will not use this information for analyzing the data , we will only apply the methodology described before for the simulated leverage time series to study the generated functional networks .
we now simulate the model setting @xmath78 to be comparable to the actual size of the argentine and brazil network .
this allows , without loss of generality , to more easily contrast the results with those in section 3 .
in fact , results remain qualitatively similar if we modify parameters such as @xmath4 , @xmath53 , @xmath79 , @xmath80 , etc .
for instance , consider panel ( a ) and ( b ) of fig .
[ fig : simulacion1 ] which show the simulated dynamics of average assets and leverage over time , respectively . after some transient period ( not shown )
the average leverage become stationary ( panel ( b ) ) .
moreover , the stationary value is 6 which is very similar to the value observed in argentina and brazil ( see fig .
2 ) . the average assets grows over time ( panel ( a ) ) , obtaining a growth of 4.5 for the time period studied .
this last value is compatible to the brazil assets growth ( fig .
now , that the average assets and leverage are similar to the empirical one , we study the individual bank behavior applying the network methodology described before for the individual simulated leverage time series .
panel ( c ) in fig .
7 shows the graphical representation of the simulated leverage network for a correlation threshold @xmath81 .
the topology is very similar to the network in fig .
in particular , the network is characterized by a large interconnected cluster , together with smaller communities of banks , and many isolated banks .
furthermore , not only banks tend to merge with the largest cluster as @xmath20 decreases , but also it presents the discontinuity property that leads to a unique large cluster ( fig . 7 ( d ) )
the proportion of nodes that belongs to the largest cluster jumps from 0.4 to over 0.8 at around @xmath35 .
the result shown in panel ( d ) is very similar to that obtained in fig .
5 .     and
5,000 periods . _
( a ) banks assets growth . ( b )
average leverage over time .
( c ) representation of simulated network .
( d ) proportion of banks that belongs to the largest cluster as a function of the correlation threshold.,scaledwidth=105.0% ]    therefore the interbank model generates modular structures that resemble the empirical networks .
but can we say more about the mechanisms that drive the links between banks ?
in particular , we want to know whether the model can shed light on the dynamics behind the leverage dependence . in order to answer this
we propose to study two banks , and in particular , the two most correlated banks . to avoid drawing conclusions from one simulation , we instead simulate the model 40 times , and for each of them we identify the two banks with the largest leverage dependence
. then compute the total leverage and assets growth , that is , the growth between the last period @xmath82 and @xmath51 .
we compute this ratio for both banks in each of the 40 replications .
the histogram of these ratios is depicted in panels ( a ) and ( b ) of fig .
interestingly , we find that the two banks exhibit large growth in both variables . in the model this takes place when banks receive greater frequency of loan requests . and at the same time , since these banks tend to have sufficient liquidity , they are also more likely to lend to other banks .
the model therefore suggests that banks with more aggressive growth strategies may be more interconnected , and hence share additional risks between them .
we then translate this measure to the analog of the empirical samples . in panels ( c ) and ( d ) of fig .
8 we compute the leverage and the total assets growth for each bank in argentina and brazil , respectively .
where banks are ordered for better visualization .
recall we are comparing the simulations against argentina and brazil because they share a similar size @xmath78 .
red and green dots identify the two most correlated banks in each country .
notice that both banks indeed exhibit a higher leverage growth .
a higher assets growth is also found in argentina , and to a less extent in brazil s network .
after the recent and still ongoing financial crisis an important bulk of research has turned its eyes into new ways to model financial risk and contagion between financial institutions , and how to improve regulation so as to reduce risk at both the individual and aggregate level . at the time of writing ,
banking supervision codes are being modified in order to incorporate tighter capital and liquidity requirements such that financial institutions are much better prepared to absorb losses in future credit - related events .
this is particularly important since a crisis starting at a given place can very easily amplify or spread across the system through various channels : counterparty risk due to interbank lending , generalized fall in asset prices ( fire sales to boost capital / margin requirements ) , liquidity hoarding or funding constraints ( higher margins , run on repo / money market ) , self - fulfilling prophecies ( bank runs ) , non - government bailout , or underestimation of risk in , and exposure to , derivate securities .
this paper studies the banking system within five different countries , and then compares the results across them .
we construct networks using bank - level balance sheet data , where links between each pair of banks are assigned upon the leverage ratio dependence .
although this interaction does not imply an interbank financial transaction , our findings do illustrate that recurring leverage dynamics ( e.g. cluster of banks that simultaneously leverage or deleverage ) is relevant in the understanding of risk concentration .
interestingly , we find that in all countries ( argentina , brazil , mexico , south africa , and taiwan ) the network topology is characterized by a large cluster of interconnected banks , together with other small communities and isolated banks . as the correlation threshold gets lower , the proportion of banks that belongs to the largest cluster increases almost dramatically : in both argentina and brazil , for instance , 90% of the banks group into a unique cluster at a 0.5 threshold .
we also find that these networks are not dominated by a hub ( extremely interconnected node ) , but by a modular structure where communities are connected by bridge banks .
this novel finding can be thought of as an optimal strategy to diversify risk within a group of banks that share similar characteristics . in other words ,
this subset restricts the spread of a distress cascade when one bank fails , as well as limits ex - post government assistance to a reduced set of banks .    to test the robustness of our results we present a model of corporate and interbank loans .
the model generates rich balance sheet and lending dynamics , which certainly do not intend to capture all the intricacies between the financial system and the rest of the economy .
our objective is much smaller , and is to generate lending dynamics that can shed light on the complexity of interbank relationships .
interestingly , in doing so we obtain leverage networks that closely follow the empirical ones .
our simulations also suggest that banks with recent rapid or aggressive growth strategies are more likely to be interconnected , and hence share more risks .
understanding risk and leverage using a network framework can also be of interest to policymakers and regulators , who for instance can build upon these tools to locate heavily interconnected banks and evaluate its importance to the aggregate system . when credit events occur , the indebtedness ratio is often the most straightforward and commonly used proxy to assess sustainability .
policymakers should enrich this analysis by incorporating additional indicators of direct exposure such as interbank loans or repo agreements , which can also be used to construct more detailed networks .
finally , understanding how individual or common shocks trigger financial crises is not only necessary to design countercyclical or containment policies , but also to design policies that prevent or limit the surge of such crises in the first place .
as suggested in @xcite , @xcite and @xcite , institutions that are more likely to destabilize the system could be subject to stricter supervision , higher capital requirements or lower interbank exposure .
this situation would certainly reduce moral hazard incentives for excessive risk - taking at institutions that are too large to be risked to fail .
the g-20 , the financial stability oversight committee and several countries on their own are already in the process of drafting special supervision guidelines for these so - called  systematically important financial institutions " or sifis ( @xcite and @xcite ) .
we are grateful to staff members from the corresponding central or reserve banks for providing us the data .
we thank d. heymann , j.c .
hallak , and two anonymous referees for helpful comments .
32 blanchard , o. j. `` the crisis : basic mechanisms and appropriate policies '' , _ imf working papers _ ,
april , 1- 22 ( 2009 ) .
duffie , darrell.innovations in credit risk transfer : implications for financial stability._bis working paper _ 255 ( 2008 ) .
haldane , r.m .
systemic risk in banking ecosystems .
_ nature _ 469:351 - 355 ( 2011 ) .
buiter , willem , `` lessons from the 2007 financial crisis '' , background paper submitted to the uk treasury select committee , december 11 ( 2007 ) .
gorton , g.b . and
metrick , a. .
`` securitized banking and the run on repo '' .
yale icf working paper no .
09 - 14 ( 2010 ) .
brunnermeier , m.k .
`` deciphering the liquidity and credit crunch 2007 - 2008 '' .
_ journal of economic perspectives _
, 23 , 1 , 77 - 100 ( 2009 ) .
klodt h , lehment h , eds .
`` the crisis and beyond .
lessons from the 2007 financial crisis '' .
kiel institute for the world economy , e - books , ( 2009 ) .
fender , ingo , and jacob gyntelberg .
`` overview : global financial crisis spurs unprecedented policy actions '' . _ bis quarterly review _ , december ( 2008 ) .
battiston , s. , gatti , d.d , gallegati , m. , greenwald , b.c.n . , and stiglitz , j.e .
`` liaisons dangereuses : increasing connectivity , risk sharing and systemic risk '' .
_ journal of economic dynamics and control _ , 36 , 1121 - 1141 ( 2012 ) .
shin , hyun song . `` reflections on northern rock : the bank run that heralded the global financial crisis '' .
_ journal of economic perspectives _ , 23(1 ) : 101 - 19 ( 2009 ) .
basel committee on banking supervision .
basel ii : international convergence of capital measurement and capital standards : a revised framework ( 2004 ) .
basel committee on banking supervision .
basel iii : a global regulatory framework for more resilient banks and banking systems ( 2010 ) .
hannoun , h. the basel iii capital framework : a decisive breakthrough .
bank for international settlements , ( 2010 ) .
s. heise and r. kuhn .
`` derivatives and credit contagion in interconnected networks '' .
_ the european physical journal b - condensed matter and complex systems _ , 85 , number 4 , 115 ( 2012 ) .
bastos e santos , e. and cont , r .. `` the brazilian interbank network structure and systemic risk '' .
_ banco central do brasil _ , working paper no . 219 ( 2010 )
nier , e. , yang , j. , yorulmazer , t. , and alentorn , a. `` network models and financial stability '' . _ journal of economic dynamics and control _ , 31(6 ) : 2033 - 2060 ( 2007 ) .
iori , g. , jafarey , s. , and padilla , f. `` systemic risk on the interbank market '' .
_ journal of economic behaviour and organisation _ , 61(4 ) : 525 - 542 ( 2006 ) .
kiyotaki , n. and moore , j .. `` credit cycles '' . _
journal of political economy _
, 105 , 2 , 211 - 248 ( 1997 ) .
krugman , p. `` balance sheets , the transfer problem , and financial crisis '' . _ international tax and public finance _ , 6 , 459 - 472 ( 1999 ) . l.f .
cspedes , r. chang , a. velasco .
`` balance sheets and exchange rate policy '' . _ american economic review _ * 94 * , 1183 - 1193 ( 2004 ) .
f.s . mishkin .
`` symposium on the monetary transmission mechanism '' . _ the journal of economic perspectives _ * 9 * , 4 , 3 - 10 ( 1995 ) .
kishan , t.p .
`` bank size , bank capital , and the bank lending channel '' .
_ journal of money , credit , and banking _ * 32 * , 1 , 121 - 141 ( 2000 ) .
s. morris , h.s .
shin.``financial regulation in a system context '' .
_ brookings papers on economic activity _ , 229 - 274 ( 2008 ) .
r. oliveira , r.f .
schiozer , l.a.b .
`` too big to fail perception by depositors : an empirical investigation '' .
_ banco central do brasil _ , working paper no . 233 ( 2011 ) .
martinez peria , s.l .
`` do depositors punish banks for bad behavior ?
market discipline , deposit insurance , and banking crises '' .
_ the journal of finance _ * 56 * , 3 , 1029 - 1051 ( 2001 ) .
a.m. maechler , k.m .
`` dynamic depositor discipline in us banks '' .
_ journal of banking and finance _ * 30 * 7 , 1871 - 1898 ( 2006 )
. m. imai .
`` market discipline and deposit insurance reform in japan '' .
_ journal of banking and finance _
, 30 , 3433 - 3452 ( 2006 ) .
`` using market information in prudential bank supervision : a review of the u.s . empirical evidence '' . _ journal of money , credit and banking _ * 30 * , 3 , 273 - 305 ( 1998 ) .
e. levy - yeyati , m.s .
martinez peria , s.l .
`` market discipline under systemic risk : evidence from bank runs in emerging economies '' . _ the world bank _ , policy research working paper series 3400 ( 2004 ) .
s. galiani , d. heymann , m. tommasi.``great expectations and hard times : the argentine convertibility plan '' .
_ economa _ * 3 * , 2 , 109 - 147 ( 2003 ) .
calvo , a. izquierdo , e. talvi .
`` sudden stops , the real exchange rate , and fiscal sustainability : argentina s lessons '' .
nber working paper no .
9828 ( 2003 ) .
l. damato , t. burdisso , v. cohen .
`` the argentine banking and exchange rate crisis of 2001 : can we learn something from financial crises ? '' .
banco central de la repblica argentina _ , working paper ( 2002 ) .
m. newman , a.l .
barabasi , and d.j . watts .
`` the structure and dynamics of networks '' .
princeton studies in complexity , ( 2006 ) .
g. caldarelli .
`` scale - free networks : complex webs in nature and technology '' .
oxford finance , ( 2007 ) .
bank for international settlements .
global systemically important banks : assessment methodology and the additional loss absorbency requirement .
( 2011 ) bank of england , the role of macroprudential policy .
discussion paper , november ( 2009 ) .
bank of england , instruments of macroprudential policy .
discussion paper , december ( 2011 ) .
m. raberto , f. rapallo , and e. scalas .
`` semi - markov graph dynamics '' .
plos one 6.8 : e23370 , ( 2011 ) .
`` banking panics : the role of the first - come , first - served rule and information externalities '' . _ journal of political economy _ , 107 ( 5 ) , 946968 , ( 1999 ) .
m. lengnick , s. krug , and h.w .
`` money creation and financial instability : an agent - based credit network approach '' .
economics : the open - access , open - assessment e - journal , 7 ( 2013 - 32 ) : 1 - 44 , ( 2013 ) .
b. lebaron , and l. tesfatsion .
`` modeling macroeconomics as open - ended dynamic systems of interaction agents '' . _ american economic review _ , 98(2 ) : 246 - 50 , ( 2006 ) .
gai , p. , haldane , a. , kapadia , s. `` complexity , concentration and contagion '' .
_ journal of monetary economics _ , 58(5 ) , 453 - 470 ( 2011 ) .
arinaminpathy , n. , kapadia , s. , may , r. m .. `` size and complexity in model financial systems '' . _
proceedings of the national academy of sciences _ , ( 2012 ) .
upper , c , `` simulation methods to assess the danger of contagion in interbank markets '' .
_ journal of financial stability _ , 7(3 ) , 111 - 125 ( 2011 ) .
bair , s.c .
`` we must resolve to end too big to fail '' .
proceedings , federal reserve bank of chicago , may , 9 - 15 ( 2011 ) .
elliot , d. j. and litan , r. e. `` identifying and regulating systemically important financial institutions : the risks of under and over identification and regulation '' .
washington : the brookings institution , ( 2011 ) .