The physical mechanism responsible for driving accretion flows in astrophysical accretion disks is commonly thought to be related to the development of plasma instabilities and turbulence. A key question is therefore the determination of consistent equilibrium configurations for accretion-disk plasmas and investigation of their stability properties. In the case of collisionless plasmas kinetic theory provides the appropriate theoretical framework. This paper presents a kinetic description of low-frequency and long-wavelength axisymmetric electromagnetic perturbations in non-relativistic, strongly-magnetized and gravitationally-bound axisymmetric accretion-disk plasmas in the collisionless regime. The analysis, carried out within the framework of the Vlasov-Maxwell description, relies on stationary kinetic solutions of the Vlasov equation which allow for the simultaneous treatment of non-uniform fluid fields, stationary accretion flows and temperature anisotropies. It is demonstrated that, for such solutions, no axisymmetric unstable perturbations can exist occurring on characteristic time and space scales which are long compared with the Larmor gyration time and radius. Hence, these stationary configurations are actually stable against axisymmetric kinetic instabilities of this type. As a fundamental consequence, this rules out the possibility of having the axisymmetric magneto-rotational or thermal instabilities to arise in these systems.A fundamental issue in the physics of accretion disks (ADs) concerns the stability of equilibrium or quasi-stationary configurations occurring in AD plasmas. The observed transport phenomena giving rise to the accretion flow are commonly ascribed to the existence of instabilities and the subsequent development of fluid or MHD turbulence [1][2][3][4]. In principle, these can include both MHD phenomena (such as drift instabilities driven by gradients of the fluid fields) and kinetic ones (due to velocity-space anisotropies, including, for example, trapped-particle modes, cyclotron and Alfven waves, etc.). Possible candidates for the angular momentum transport mechanism are usually identified either with the magneto-rotational instability (MRI) [5][6][7] or the thermal instability (TMI) [8][9][10][11], caused by unfavorable gradients of rotation/shear and temperature respectively. The validity of the above identifications needs to be checked in this case, because they usually rely on incomplete physical descriptions, which ignore the microscopic (kinetic) plasma behavior. In fact, "stand-alone" fluid and MHD approaches which are not explicitly based on kinetic theory and/or do not start from consistent kinetic equilibria, may become inadequate or inapplicable for collisionless or weakly-collisional plasmas. Apart from possible gyrokinetic and finite Larmor-radius effects (which are typically not included for MRI and TMI), this concerns consistent treatment of the kinetic constraints which must be imposed on the fluid fields (see related discussion in Refs. [12,13]). This concerns,...