Stability of large complex systems with heterogeneous relaxation dynamics

We study the probability of stability of a large complex system of size N within the framework of a generalized May model, which assumes a linear dynamics of each population size n _i (with respect to its equilibrium value): $\frac{\mathrm{d}{n}_{i}}{\mathrm{d}t}=-{a}_{i}{n}_{i}-\sqrt{T}{\sum }_{j}{J}_{ij}{n}_{j}$ . The a _i > 0's are the intrinsic decay rates, J _ij is a real symmetric (N × N) Gaussian random matrix and $\sqrt{T}$ measures the strength of pairwise interaction between different species. Our goal is to study how inhomogeneities in the intrinsic damping rates a _i affect the stability of this dynamical system. As the interaction strength T increases, the system undergoes a phase transition from a stable phase to an unstable phase at a critical value T = T _c. We reinterpret the probability of stability in terms of the hitting time of the level b = 0 of an associated Dyson Brownian motion (DBM), starting at the initial position a _i and evolving in 'time' T. In the large N → ∞ limit, using this DBM picture, we are able to completely characterize T _c for arbitrary density μ(a) of the a _i 's. For a specific flat configuration ${a}_{i}=1+\sigma \frac{i-1}{N}$ , we obtain an explicit parametric solution for the limiting (as N → ∞) spectral density for arbitrary T and σ. For finite but large N, we also compute the large deviation properties of the probability of stability on the stable side T < T _c using a Coulomb gas representation.