We develop a method to study coupled dynamics of gauge-invariant variables, constructed out of metric and gauge field fluctuations on the background of a AdS5 Reissner-Nordström black brane. Using this method, we compute the numerical spectrum of quasinormal modes associated with fluctuations of spin 0, 1 and 2, non-perturbatively in μ/T . We also analytically compute the spectrum of hydrodynamic excitations in the small chemical potential limit. Then, by studying the spectral curve at complex momenta in every spin channel, we numerically find points at which hydrodynamic and non-hydrodynamic poles collide. We discuss the relation between such collision points and the convergence radius of the hydrodynamic derivative expansion. Specifically in the spin 0 channel, we find that within the range $$ 1.1\underset{\sim }{<}\mu /T\underset{\sim }{<}2 $$
1.1
<
∼
μ
/
T
<
∼
2
, the radius of convergence of the hydrodynamic sound mode is set by the absolute value of the complex momentum corresponding to the point at which the sound pole collides with the hydrodynamic diffusion pole. It shows that in holographic systems at finite chemical potential, the convergence of the hydrodynamic derivative expansion in the mentioned range is fully controlled by hydrodynamic informa- tion. As the last result, we explicitly show that the relevant information about quantum chaos in our system can be extracted from the pole-skipping points of energy density re- sponse function. We find a threshold value for μ/T , lower than which the pole-skipping points can be computed perturbatively in a derivative expansion.