We establish a new self-consistent system of equations accounting for a nonminimal coupling of the cooperative gravitational, electromagnetic and pseudoscalar (axion) fields in a multi-component relativistic plasma. The axionic extension of the nonminimal Einstein-Maxwell-Vlasov theory is based on two consistent procedures. First, we use the Lagrange formalism to obtain nonminimal equations for the gravitational, electromagnetic and pseudoscalar fields with the additional sources generated in plasma. Second, we use the Vlasov version of the relativistic kinetic theory of the plasma, guided by the cooperative macroscopic electromagnetic, gravitational and axionic fields, to describe adequately the response of the plasma on the variations of these fields. In order to show the self-consistency of this approach we check directly the compatibility conditions for the master equations for the cooperative fields. Using these compatibility conditions we reconstruct the ponderomotive force, which acts on the plasma particles, and discuss the necessary conditions for existence of the distribution function of the equilibrium type.

We consider cosmological applications of a new self-consistent system of equations, accounting for a nonminimal coupling of the gravitational, electromagnetic and pseudoscalar (axion) fields in a relativistic plasma. We focus on dispersion relations for electromagnetic perturbations in an initially isotropic ultrarelativistic plasma coupled to the gravitational and axion fields in the framework of isotropic homogeneous cosmological model of the de Sitter type. We classify the longitudinal and transversal electromagnetic modes in an axionically active plasma and distinguish between waves (damping, instable or running), and nonharmonic perturbations (damping or instable). We show that for the special choice of the guiding model parameters the transversal electromagnetic waves in the axionically active plasma, nonminimally coupled to gravity, can propagate with the phase velocity less than speed of light in vacuum, thus displaying a possibility for a new type of resonant particle-wave interactions.

A three-parameter toy-model, which describes a non-minimal coupling of gravity field with electromagnetic field of a relativistic two-component electrically neutral plasma, is discussed. Resonance interactions between particles and transversal waves in plasma are shown to take place due to the curvature coupling effect.

The spherically symmetric perturbations in the spatially flat Friedman models are considered. It is assumed that the Friedmannian density and pressure are related through a linear equation of state. The perturbation is joined smoothly with an unperturbed Friedmann's background at the sound horizon of perturbation. Such junction is in accordance with the "birth" of a local perturbation as a result of the redistribution of matter. The solution of the Einstein's equations is obtained in linear approximation on a Friedmann's background near the the sound horizon of perturbation.