The fidelity susceptibility measures sensitivity of eigenstates to a change of an external parameter. It has been fruitfully used to pin down quantum phase transitions when applied to ground states (with extensions to thermal states). Here we propose to use the fidelity susceptibility as a useful dimensionless measure for complex quantum systems. We find analytically the fidelity susceptibility distributions for Gaussian orthogonal and unitary universality classes for arbitrary system size. The results are verified by a comparison with numerical data.The discovery of many body localization (MBL) phenomenon resulting in non-ergodicity of the dynamics in many body systems [1] restored also the interest in purely ergodic phenomena modeled by Gaussian random ensembles (GRE) [2] and in possible measures to characterize them. The gap ratio between adjacent level spacings [3] was introduced precisely for that purpose as it does not involve the so called unfolding [4] necessary for meaningful studies of level spacing distributions and yet often leading to spurious results [5]. Still, the level spacing distribution belongs to the most popular statistical measures used for single particle quantum chaos studies [6][7][8][9] and also in the transition to MBL [10-13]. A particular place among different measures was taken by those characterizing level dynamics for a Hamiltonian H(λ) dependent on some parameter λ. In Pechukas-Yukawa formulation [14,15] energy levels are positions of a fictitious gas particles, derivatives with respect to the fictitious time λ are velocities (level slopes), the second derivatives describe curvatures of the levels (accelerations). Simons and Altschuler [16] put forward a proposition that the variance of velocities distribution is an important parameter characterizing universality of level dynamics. This led to predictions for distributions of avoided crossings [17] and, importantly, curvature distributions postulated first on the basis of numerical data for GRE [18] and then derived analytically via supersymmetric method by von Oppen [19,20] (for alternative techniques see [21,22]). Curvature distributions were recently addressed in MBL studies [23,24].Apart from quantum chaos studies in the eighties and nineties of the last millennium, another "level dynamics" tool has been introduced in the quantum information area, i.e. the fidelity [25]. It compares two close (possibly mixed) quantum states. If these states are dependent on a parameter λ it is customary to introduce a fidelity susceptibility χ. For sufficiently small λ, in a finite system, one hasFidelity susceptibility is directly related to the quantum Fisher information (QFI), G, being directly proportional to the Bures distance between density matrices at slightly differing values of λ [26,27], with G(λ) = 4χ. Fidelity susceptibility emerged as a useful tool to study quantum phase transitions as at the transition point the ground state changes rapidly leading to the enhancement of χ [27][28][29][30][31][32][33][34]. All of these studies w...