We consider the distribution of the axionic dark matter and the profiles of the axionically induced electric field on the background of a regular magnetic monopole. The background static spherically symmetric spacetime is considered to be regular at the center of the object and has no event horizons, the nonminimal interaction of strong gravitational and monopole type magnetic fields being the source, which provides this regularity. We assume that the pseudoscalar (axion) field also is nonminimally coupled to the electromagnetic and gravitational fields, and we search for solutions to the extended equations of nonminimal axion electrostatics, which satisfy two requirements: first, the axionically induced radial electric field vanishes at the center; second, the pseudoscalar field is regular at the center of the formed axionic dyon. We have found the set of parameters, describing the nonminimal axion-photon coupling, for which both requirements are satisfied. Properties of the two-parameter family of regular solutions to the equations of nonminimal axion electrostatics are analyzed numerically. The regular axionically induced electric field is shown to have the structure, which is typical for a simple electric layer.

We suggest the phenomenological model of emergence of the dynamic aether as a result of decay of the SU(N) symmetric field configuration containing the multiplet of vector fields. The scenario of the transition to the dynamic aether, which is characterized by one unit timelike vector field associated with the aether velocity, is based on the idea of spontaneous color polarization analogous to the spontaneous electric polarization in ferroelectric materials. The mechanism of spontaneous color polarization is described in the framework of anisotropic cosmological model of the Bianchi-I type; it involves into consideration the idea of critical behavior of the eigenvalues of the tensor of color polarization in the course of the Universe accelerated expansion. The interim stage of transition from the color aether to the canonic dynamic aether takes the finite period of time, the duration of which is predetermined by the phenomenologically introduced critical value of the expansion scalar.

We consider the magnetic monopole in the axionic dark matter environment (axionic dyon) in the framework of the Reissner - Nordstr\"om - de Sitter model. Our aim is to study the distribution of the pseudoscalar (axion) and electric fields near the so-called folds, which are characterized by the profiles with the central minimum, the barrier on the left, and the maximum on the right of this minimum. The electric field in the fold-like zones is shown to change the sign twice, i.e., the electric structure of the near zone of the axionic dyon contains the domain similar to a double electric layer. We have shown that the described fold-like structures in the profile of the gravitational potential, and in the profiles of the electric and axion fields can exist, when the value of the dyon mass belongs to the interval enclosed between two critical masses, which depend on the cosmological constant.

We consider the structure of halos of the axionic dark matter, which surround massive relativistic objects with static spherically symmetric gravitational field and monopole-type magneto-electric fields. We work with the model of pseudoscalar field with the extended periodic potential, which depends on additional arguments proportional to the moduli of the Killing vectors; in our approach they play the roles of model guiding functions. The covariant model of the axion field with this modified potential is equipped with the extended formalism of the Killing vector fields, which is established in analogy with the formalism of the Einstein-aether theory, based on the introduction of a unit timelike dynamic vector field. We study the equilibrium state of the axion field, for which the extended potential and its derivative vanish, and illustrate the established formalism by the analysis of two-level axionic dark matter profiles, for which the stage delimiters relate to the critical values of the modulus of the timelike Killing vector field.

In the framework of Einstein-Maxwell-axion theory we consider a static field configuration in the outer zone of a magnetic star. We assume that this field configuration is formed by the following interacting quartet: a spherically symmetric gravitational field of the Reissner-Nordstr\"om type, pseudoscalar field, associated with an axionic dark matter, strong intrinsic magnetic field, and axionically induced electric field. Based on the analysis of the solutions to the master equations of the model, we show that the structure of the formed field configuration reveals the following triple effect. First, if the strong magnetic field of the star has the dipole component, the axionic halo around the magnetic star with spherically symmetric gravitational field is no longer spheroidal, and the distortion is induced by the axion-photon coupling; second, due to the distortion of the axionic halo the magnetic field is no longer pure dipolar, and the induced quadruple, octupole, etc., components appear thus redistributing the magnetic energy between the components with more steep radial profiles than the dipolar one; third, the interaction of the axionic and magnetic fields produces the electric field, which inherits their multipole structure. We attract attention to possible applications of the predicted field restructuring to the theory of axionic rotating magnetars.

In the framework of axionic extension of the Einstein-aether theory we establish the model, which describes a stiff regulation of the behavior of axionic dark matter by the dynamic aether. The aether realizes this procedure via the modified Higgs potential, designed for modeling of nonlinear self-interaction of pseudoscalar (axion) field; the modification of this potential is that its minima are not fixed, and their positions and depths depend now on the square of the covariant derivative of the aether velocity four-vector. Exact solutions to the master equations, modified correspondingly, are obtained in the framework of homogeneous isotropic cosmological model. The effective equation of state for axionic dark matter is of the stiff type. Homogeneous perturbations of the pseudoscalar (axion) field, of the Hubble function and of the scale factor are shown to fade out with cosmological time, there are no growing modes, the model of stiff regulation is stable.

The state of a static spherically symmetric relativistic axionically active multi-component plasma in the gravitational, magnetic and electric fields of an axionic dyon is studied in the framework of the Einstein - Maxwell - Boltzmann - axion theory. We assume that the equations of axion electrodynamics, the covariant relativistic kinetic equations, and the equation for the axion field with modified Higgs-type potential are nonlinearly coupled; the gravitational field in the dyon exterior is assumed to be fixed and to be of the Reissner-Nordstr\"om type. We introduce the extended Lorentz force, which acts on the particles in the axionically active plasma, and analyze the consequences of this generalization. The analysis of exact solutions, obtained in the framework of this model for the relativistic Boltzmann electron-ion and electron-positron plasmas, as well as, for degenerated zero-temperature electron gas, shows that the phenomena of polarization and stratification can appear in plasma, attracting attention to the axionic analog of the known Pannekoek-Rosseland effect.

We study a new exactly solvable model of coupling of the Dark Energy and Dark Matter, in the framework of which the kernel of non-gravitational interaction is presented by the integral Volterra-type operator well-known in the classical theory of fading memory. Exact solutions of this isotropic homogeneous cosmological model were classified with respect to the sign of the discriminant of the cubic characteristic polynomial associated with the key equation of the model. Energy-density scalars of the Dark Energy and Dark Matter, the Hubble function and acceleration parameter are presented explicitly; the scale factor is found in quadratures. Asymptotic analysis of the exact solutions has shown that the Big Rip, Little Rip, Pseudo Rip regimes can be realized with the specific choice of guiding parameters of the model. We show that the Coincidence problem can be solved if we consider the memory effect associated with the interactions in the Dark Sector of the Universe.

The extended Einstein-Maxwell-aether-axion model describes internal interactions inside the system, which contains gravitational, electromagnetic fields, the dynamic unit vector field describing the velocity of an aether, and the pseudoscalar field associated with the axionic dark matter. The specific feature of this model is that the axion field controls the dynamics of the aether through the guiding functions incorporated into the Jacobson's constitutive tensor. Depending on the state of the axion field these guiding functions can control and switch on or switch off the influence of acceleration, shear, vorticity and expansion of the aether flow on the state of physical system as a whole. We obtain new exact solutions, which possess the pp-wave symmetry, and indicate them by the term pp-wave aether modes in contrast to the pure pp-waves, which can not propagate in this field conglomerate. These exact solutions describe a specific dynamic state of the pseudoscalar field, which corresponds to one of the minima of the axion potential, and switches off the influence of shear and expansion of the aether flow; the model does not impose restrictions on the Jacobson's coupling constants and on the axion mass. Properties of these new exact solutions are discussed.

The SU(N)-symmetric generalization of the model of the electromagnetically active dynamic aether is formulated. This generalization is based on the introduction of a Yang-Mills gauge field instead of the Maxwell field, and of a SU(N)-multiplet of vector fields instead of the standard single vector field. In the framework of the second order version of the effective field theory this generalization includes three constitutive tensors, which are the SU(N) extensions of the tensors appeared in the Einstein-Maxwell-aether theory; we reconstructed the full-format set of these constitutive tensors. The total self-consistent system of master equations for the gauge, vector and gravitational fields is obtained by the variation procedure. The general model of the SU(N)-symmetric dynamic aether is reduced to the extended Einstein-Yang-Mills-aether model by the ansatz about spontaneous color polarization of the vector fields. In fact, this ansatz requires the vector fields, which form the SU(N) multiplet, to become parallel in the group (color) space due to a phase transition, and a new selected direction in the group space to appear, thus converting it into the anisotropic color space.

In the framework of the Einstein-Maxwell-aether theory we study the birefringence effect, which can occur in the pp-wave symmetric dynamic aether. The dynamic aether is considered to be latently birefringent quasi-medium, which displays this hidden property if and only if the aether motion is non-uniform, i.e., when the aether flow is characterized by the non-vanishing expansion, shear, vorticity or acceleration. In accordance with the dynamo-optical scheme of description of the interaction between electromagnetic waves and the dynamic aether, we shall model the susceptibility tensors by the terms linear in the covariant derivative of the aether velocity four-vector. When the pp-wave modes appear in the dynamic aether, we deal with a gravitationally induced degeneracy removal with respect to hidden susceptibility parameters. As a consequence, the phase velocities of electromagnetic waves possessing orthogonal polarizations do not coincide, thus displaying the birefringence effect. Two electromagnetic field configurations are studied in detail: longitudinal and transversal with respect to the aether pp-wave front. For both cases the solutions are found, which reveal anomalies in the electromagnetic response on the action of the pp-wave aether mode.

Cosmic Dark Fluid is considered as a non-stationary medium, in which electromagnetic waves propagate, and magneto-electric field structures emerge and evolve. A medium - type representation of the Dark Fluid allows us to involve into analysis the concepts and mathematical formalism elaborated in the framework of classical covariant electrodynamics of continua, and to distinguish dark analogs of well-known medium-effects, such as optical activity, pyro-electricity, piezo-magnetism, electro- and magneto-striction and dynamo-optical activity. The Dark Fluid is assumed to be formed by a duet of a Dark Matter (a pseudoscalar axionic constituent) and Dark Energy (a scalar element); respectively, we distinguish electrodynamic effects induced by these two constituents of the Dark Fluid. The review contains discussions of ten models, which describe electrodynamic effects induced by Dark Matter and/or Dark Energy. The models are accompanied by examples of exact solutions to the master equations, correspondingly extended; applications are considered for cosmology and space-times with spherical and pp-wave symmetries. In these applications we focused the attention on three main electromagnetic phenomena induced by the Dark Fluid: first, emergence of Longitudinal Magneto-Electric Clusters; second, generation of anomalous electromagnetic responses; third, formation of Dark Epochs in the Universe history.

We extend the Einstein-aether theory to take into account the interaction between a pseudoscalar field, which describes the axionic dark matter, and a time-like dynamic unit vector field, which characterizes the velocity of the aether motion. The Lagrangian of the Einstein-aether-axion theory includes cross-terms based on the axion field and its gradient four-vector, on the covariant derivative of the aether velocity four-vector, and on the Riemann tensor and its convolutions. We follow the principles of the Effective Field theory, and include into the Lagrangian of interactions all possible terms up to the second order in the covariant derivative. Interpretation of new couplings is given in terms of irreducible parts of the covariant derivative of the aether velocity, namely, the acceleration four-vector, the shear and vorticity tensors, and the expansion scalar. A spatially isotropic and homogeneous cosmological model with dynamic unit vector field and axionic dark matter is considered as an application of the established theory; new exact solutions are discussed, which describe models with a Big Rip, Pseudo Rip and de Sitter-type asymptotic behavior.

We consider a pp-wave symmetric model in the framework of the Einstein-Maxwell-aether-axion theory. Exact solutions to the equations of axion electrodynamics are obtained for the model, in which pseudoscalar, electric and magnetic fields were constant before the arrival of a gravitational pp-wave. We show that dynamo-optical interactions, i.e. couplings of the electromagnetic field to a dynamic unit vector field, attributed to the velocity of a cosmic substratum (aether, vacuum, dark fluid...), provide the response of axionically active electrodynamic system to display anomalous behavior.

We discuss the problem of identification of coupling constants, which describe interactions between photons and space-time curvature, using exact regular solutions to the extended equations of the nonminimal Einstein-Maxwell theory. We argue the idea that three nonminimal coupling constants in this theory can be reduced to the single guiding parameter, which plays the role of nonminimal radius. We base our consideration on two examples of exact solutions obtained earlier in our works: the first of them describes a nonminimal spherically symmetric object (star or black hole) with regular radial electric field; the second example represents a nonminimal Dirac-type object (monopole or black hole) with regular metric. We demonstrate that one of the inflexion points of the regular metric function identifies a specific nonminimal radius, thus marking the domain of dominance of nonminimal interactions.

We extend the Einstein-aether theory to include the Maxwell field in a nontrivial manner by taking into account its interaction with the time-like unit vector field characterizing the aether. We also include a generic matter term. We present a model with a Lagrangian that includes cross-terms linear and quadratic in the Maxwell tensor, linear and quadratic in the covariant derivative of the aether velocity four-vector, linear in its second covariant derivative and in the Riemann tensor. We decompose these terms with respect to the irreducible parts of the covariant derivative of the aether velocity, namely, the acceleration four-vector, the shear and vorticity tensors, and the expansion scalar. Furthermore, we discuss the influence of an aether non-uniform motion on the polarization and magnetization of the matter in such an aether environment, as well as on its dielectric and magnetic properties. The total self-consistent system of equations for the electromagnetic and the gravitational fields, and the dynamic equations for the unit vector aether field are obtained. Possible applications of this system are discussed. Based on the principles of effective field theories, we display in an appendix all the terms up to fourth order in derivative operators that can be considered in a Lagrangian that includes the metric, the electromagnetic and the aether fields.

We establish a new model of coupling between a cosmic dark fluid and electrodynamic systems, based on an analogy with effects of electric and magnetic striction, piezo-electricity and piezo-magnetism, pyro-electricity and pyro-magnetism, which appear in classical electrodynamics of continuous media. Extended master equations for electromagnetic and gravitational fields are derived using Lagrange formalism. A cosmological application of the model is considered, and it is shown that a striction-type interaction between the dark energy (the main constituent of the dark fluid) and electrodynamic system provides the universe history to include the so-called unlighted epochs, during which electromagnetic waves can not propagate and thus can not scan the universe interior.

We extend the Einstein-Maxwell-axion theory including into the Lagrangian cross-terms of the dynamo-optical type, which are quadratic in the Maxwell tensor, linear in the covariant derivative of the macroscopic velocity four-vector, and linear in the pseudoscalar (axion) field or its gradient four-vector. We classify the new terms with respect to irreducible elements of the covariant derivative of the macroscopic velocity four-vector of the electromagnetically active medium: the expansion scalar, acceleration four-vector, shear and vorticity tensors. Master equations of the extended axion electrodynamics are used for the description of the response of an axionically active electrodynamic system, induced by a pp-wave gravitational background. We show that this response has a critical character, i.e., the electric and magnetic fields, dynamo-optically coupled to the axions, grow anomalously under the influence of the external pp-wave gravitational field.

We study the problem of axion-photon coupling in the magnetic field influenced by gravitational radiation. We focus on exact solutions to the equations for axion electrodynamics in the pp-wave gravitational background for two models with initially constant magnetic field. The first model describes the response of an initially constant magnetic field in a gravitational-wave vacuum with unit refraction index; the second model is characterized by a non-unit refraction index prescribed to the presence of ordinary and/or dark matter. We show that both models demonstrate anomalous behavior of the electromagnetic field generated by the axion-photon coupling in the presence of magnetic field, evolving in the gravitational wave background. The role of axionic dark matter in the formation of the anomalous response of this electrodynamic system is discussed.

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.

We analyze late-time evolution of the Universe in the framework of the self-consistent model, in which the dark matter is influenced by the Archimedean-type force proportional to the four-gradient of the dark energy pressure. The dark energy is considered as a fluid with the equation of state of the relaxation type, which takes into account a retardation of the dark energy response to the Universe accelerated expansion. The dark matter is guided by the Archimedean-type force, which redistributes the total energy of the dark fluid between two its constituents, dark energy and dark matter, in the course of the Universe accelerated expansion. We focus on the constraints for the dark energy relaxation time parameter, for the dark energy equation of state parameter, and for the Archimedean-type coupling constants, which guarantee the Big Rip avoidance. In particular, we show that the Archimedean-type coupling protects the Universe from the Big Rip scenario with asymptotically infinite negative dark energy pressure, and that the Little Rip is the fate of the Universe with the Archimedean-type interaction inside the dark fluid.

We consider mathematical aspects of the axion electrodynamics in application to the problem of evolution of geomagnetic and terrestrial electric fields, which are coupled by relic axions born in the early Universe and (hypothetically) forming now the cold dark matter. We find axionic analogs of the Debye potentials, well-known in the standard Faraday - Maxwell electrodynamics, and discuss exact solutions to the equations of the axion electrodynamics describing the state of axionically coupled electric and magnetic fields in a spherical resonator Earth-Ionosphere. We focus on the properties of the specific electric and magnetic oscillations, which appeared as a result of the axion-photon coupling in the dark matter environment. We indicate such electric and magnetic field configurations as longitudinal electro-magnetic clusters.

We establish a new self-consistent Einstein-Maxwell-axion model based on the Lagrangian, which is linear in the pseudoscalar (axion) field and its four-gradient and includes the four-vector of macroscopic velocity of the axion system as a whole. We consider extended equations of the axion electrodynamics, modified gravity field equations, and discuss nonstationary effects in the phenomenon of optical activity induced by axions.

We establish an extended version of the Einstein - Maxwell - axion model by introducing into the Lagrangian cross-terms, which contain the gradient four-vector of the pseudoscalar (axion) field in convolution with the Maxwell tensor. The gradient model of the axion-photon coupling is applied to cosmology: we analyze the Bianchi-I type Universe with an initial magnetic field, electric field induced by the axion-photon interaction, cosmological constant and dark matter, which is described in terms of the pseudoscalar (axion) field. Analytical, qualitative and numerical results are presented in detail for two distinguished epochs: first, for the early Universe with magnetic field domination; second, for the stage of late-time accelerated expansion.

During the evolution of the universe there are at least two epochs during which electromagnetic waves cannot scan the universe's internal structure neither bring information to outside observers. The first epoch is when photons are in local thermodynamic equilibrium with other particles, and the second is when photon scattering by charged particles is strong. One can call these two periods of cosmological time as standard unlighted epochs. After the last scattering surface, photons become relic photons and turn into a source of information about the universe. Unlighted cosmic epochs can also appear when one considers non-minimal theories, i.e., theories in which the electromagnetic field is coupled in an intricate way with the cosmological gravitational field. By considering a cosmological model where the dark sector, i.e., the dark energy and dark matter, self-interacts via an Archimedean-type force, and taking into account a non-minimal coupling theory for the electromagnetic field, we discuss the appearance of unlighted epochs. In the framework of our non-minimal theory, a three-parameter non-minimal Einstein-Maxwell model, the curvature coupling can be formulated in terms of an effective refraction index n(t). Then, taking advantage of a well-known classical analogy, namely, in a medium with n^2<0 electromagnetic waves do not propagate and their group velocity, i.e., energy transfer velocity, has zero value at the boundary of the corresponding zone, one can search for the unlighted epochs arising in the interacting dark fluid cosmological model. We study here, both analytically and numerically, cosmological models admitting unlighted epochs.

Pseudoscalar-photon interactions were proposed in the study of the relations among equivalence principles. The interaction of pseudoscalar axion with gluons was proposed as a way to solve the strong CP problem. Subsequent proposal of axion as a dark matter candidate has been a focus of search. Motivation from superstring theories add to its importance. After a brief introduction and historical review, we present (i) the current status of our optical experiment using high-finesse Fabry-Perot resonant cavity - Q & A experiment - to detect pseudoscalar-photon interactions, (ii) the constraints on pseudoscalar-photon interactions from astrophysical and cosmological observations on cosmic polarization rotation, and (iii) theoretical models of non-minimal interactions of gravitational, electromagnetic and pseudoscalar (axion) fields, and their relevance to cosmology.

We establish a new self-consistent model in order to explain from a unified viewpoint two key features of the cosmological evolution: the inflation in the early Universe and the late-time accelerated expansion. The key element of this new model is the Archimedean-type coupling of the dark matter with dark energy, which form the so-called cosmic dark fluid. We suppose that dark matter particles immersed into the dark energy reservoir are affected by the force proportional to the four-gradient of the dark energy pressure. The Archimedean-type coupling is shown to play a role of effective energy-momentum redistributor between the dark matter and the dark energy components of the dark fluid, thus providing the Universe's evolution to be a quasiperiodic and/or multistage process. In the first part of the work we discuss a theoretical base and new exact solutions of the model master equations. Special attention is focused on the exact solutions for which the scale factor is presented by the anti-Gaussian function: these solutions describe the late-time acceleration and are characterized by a nonsingular behavior in the early Universe. The second part contains qualitative and numerical analysis of the master equations; we focus there on the solutions describing a multi-inflationary Universe.

In this (second) part of the work we present the results of numerical and qualitative analysis, based on a new model of the Archimedean-type interaction between dark matter and dark energy. The Archimedean-type force is linear in the four-gradient of the dark energy pressure and plays a role of self-regulator of the energy redistribution in a cosmic dark fluid. Because of the Archimedean-type interaction the cosmological evolution is shown to have a multistage character. Depending on the choice of the values of the model guiding parameters,the Universe's expansion is shown to be perpetually accelerated, periodic or quasiperiodic with finite number of deceleration/acceleration epochs. We distinguished the models, which can be definitely characterized by the inflation in the early Universe, by the late-time accelerated expansion and nonsingular behavior in intermediate epochs, and classified them with respect to a number of transition points. Transition points appear, when the acceleration parameter changes the sign, providing the natural partition of the Universe's history into epochs of accelerated and decelerated expansion. The strategy and results of numerical calculations are advocated by the qualitative analysis of the instantaneous phase portraits of the dynamic system associated with the key equation for the dark energy pressure evolution.

We discuss new exact solutions of a three-parameter nonminimal Einstein-Maxwell model. The solutions describe static spherically symmetric objects with and without center, supported by an electric field nonminimally coupled to gravity. We focus on a unique one-parameter model, which admits an exact solution for a traversable electrically charged wormhole connecting two universes, one asymptotically flat the other asymptotically de Sitter ones. The relation between the asymptotic mass and charge of the wormhole and its throat radius is analyzed. The wormhole solution found is thus a nonminimal realization of the idea of Wheeler about charge without charge and shows that, if the world is somehow nonminimal in the coupling of gravity to electromagnetism, then wormhole appearance, or perhaps construction, is possible.

We establish a new self-consistent system of equations accounting for a non-minimal interaction of gravitational, electromagnetic and axion fields. The procedure is based on a non-minimal extension of the standard Einstein-Maxwell-axion action. The general properties of a ten-parameter family of non-minimal linear models are discussed. We apply this theory to the models with pp-wave symmetry and consider propagation of electromagnetic waves non-minimally coupled to the gravitational and axion fields. We focus on exact solutions of electrodynamic equations, which describe quasi-minimal and non-minimal optical activity induced by the axion field. We also discuss empirical constraints on coupling parameters from astrophysical birefringence and polarization rotation observations.

We establish a new model, which takes into account a dynamic (inertial) self-interaction of gravitating systems. The model is formulated by introduction of a new function depending on the square of the covariant derivative of the velocity four-vector of the system as a whole into the Lagrangian. This term is meant for description of both self-action of the system irregularly moving in the gravitational field, and back-reaction of the motion irregularities on the gravity field. We discuss one example of exact solution to the extended master equations in the framework of cosmological model of the FLRW type with vanishing cosmological constant. It is shown that accelerated expansion of the Universe can be driven by traditional matter with positive pressure (e.g., dust, ultrarelativistic fluid) due to the back-reaction of the gravity field induced by irregular motion of the system as a whole; this back-reaction is shown to be characterized by the negative effective pressure.

We discuss a class of exact solutions of a three-parameter non-minimally extended Einstein-Maxwell model, which are attributed to non-minimal magnetic monopoles of the Dirac type. We focus on the investigation of the gravitational field of Dirac monopoles for those models, for which the singularity at the central point is hidden inside of an event horizon independently on the mass and charge of the object. We obtained the relationships between the non-minimal coupling constants, for which this requirement is satisfied. As explicit examples, we consider in detail two one-parameter models: first, non-minimally extended Reissner-Nordstr\"om model (for the magnetically charged monopole), second, the Drummond-Hathrell model.

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.

We establish a new non-minimal Einstein-Yang-Mills-dilaton model, for which the Lagrangian is linear in the curvature and contains eight arbitrary functions of the scalar (dilaton) field. The self-consistent system of equations for the non-minimally coupled gauge, scalar and gravitational fields is derived. As an example of an application we discuss the model with pp-wave symmetry. Two exact explicit regular solutions of the whole system of master equations, belonging to the class of pp-wave solutions, are presented.

We formulate a self-consistent non-minimal five-parameter Einstein-Yang-Mills-Higgs (EYMH) model and analyse it in terms of effective (associated, color and color-acoustic) metrics. We use a formalism of constitutive tensors in order to reformulate master equations for the gauge, scalar and gravitational fields and reconstruct in the algebraic manner the so-called associated metrics for the Yang-Mills field. Using WKB-approximation we find color metrics for the Yang-Mills field and color-acoustic metric for the Higgs field in the framework of five-parameter EYMH model. Based on explicit representation of these effective metrics for the EYMH system with uniaxial symmetry, we consider cosmological applications for Bianchi-I, FLRW and de Sitter models. We focus on the analysis of the obtained expressions for velocities of propagation of longitudinal and transversal color and color-acoustic waves in a (quasi)vacuum interacting with curvature; we show that curvature coupling results in time variations of these velocities. We show, that the effective metrics can be regular or can possess singularities depending on the choice of the parameters of non-minimal coupling in the cosmological models under discussion. We consider a physical interpretation of such singularities in terms of phase velocities of color and color-acoustic waves, using the terms ``wave stopping'' and ``trapped surface''.

Using a Lagrangian formalism, a three-parameter non-minimal Einstein-Maxwell theory is established. The three parameters, $q_1$, $q_2$ and $q_3$, characterize the cross-terms in the Lagrangian, between the Maxwell field and terms linear in the Ricci scalar, Ricci tensor, and Riemann tensor, respectively. Static spherically symmetric equations are set up, and the three parameters are interrelated and chosen so that effectively the system reduces to a one parameter only, $q$. Specific black hole and other type of one-parameter solutions are studied. First, as a preparation, the Reissner-Nordstr\"om solution, with $q_1=q_2=q_3=0$, is displayed. Then, we seek for solutions in which the electric field is regular everywhere as well as asymptotically Coulombian, and the metric potentials are regular at the center as well as asymptotically flat. In this context, the one-parameter model with $q_1\equiv -q$, $q_2=2q$, $q_3=-q$, called the Gauss-Bonnet model, is analyzed in detail. The study is done through the solution of the Abel equation (the key equation), and the dynamical system associated with the model. There is extra focus on an exact solution of the model and its critical properties. Finally, an exactly integrable one-parameter model, with $q_1\equiv -q$, $q_2=q$, $q_3=0$, is considered also in detail. A special sub-model, in which the Fibonacci number appears naturally, of this one-parameter model is shown, and the corresponding exact solution is presented. Interestingly enough, it is a soliton of the theory, the Fibonacci soliton, without horizons and with a mild conical singularity at the center.

Non-minimal interactions in the pp-wave Einstein - Yang - Mills - Higgs (EYMH) model are shown to give rise to color cross-effects analogous to the magneto-electricity in the Maxwell theory. In order to illustrate the significance of these color cross-effects, we reconstruct the effective (associated, color and color-acoustic) metrics for the pp-wave non-minimal seven-parameter EYMH model with parallel gauge and scalar background fields. Then these metrics are used as hints for obtaining explicit exact solutions of the non-minimally extended Yang-Mills and Higgs equations for the test fields propagating in the vacuum interacting with curvature. The influence of the non-minimal coupling on the test particle motion is interpreted in terms of the so-called trapped surfaces, introduced in the Analog Gravity theory.

We establish a nonminimal Einstein-Yang-Mills-Higgs model, which contains six coupling parameters. First three parameters relate to the nonminimal coupling of non-Abelian gauge field and gravity field, two parameters describe the so-called derivative nonminimal coupling of scalar multiplet with gravity field, and the sixth parameter introduces the standard coupling of scalar field with Ricci scalar. The formulated six-parameter nonminimal Einstein-Yang-Mills-Higgs model is applied to cosmology. We show that there exists a unique exact cosmological solution of the de Sitter type for a special choice of the coupling parameters. The nonminimally extended Yang-Mills and Higgs equations are satisfied for arbitrary gauge and scalar fields, when the coupling parameters are specifically related to the curvature constant of the isotropic spacetime. Basing on this special exact solution we discuss the problem of a hidden anisotropy of the Yang-Mills field, and give an explicit example, when the nonminimal coupling effectively screens the anisotropy induced by the Yang-Mills field and thus restores the isotropy of the model.

Based on the analogy with non-minimal SU(2) symmetric Wu-Yang monopole with regular metric, the solution describing a non-minimal U(1) symmetric Dirac monopole is obtained. In order to take into account the curvature coupling of gravitational and electromagnetic fields, we reconstruct the effective metrics of two types: the so-called associated and optical metrics. The optical metrics display explicitly that the effect of birefringence induced by curvature takes place in the vicinity of the non-minimal Dirac monopole; these optical metrics are studied analytically and numerically.

We discuss a non-minimal Einstein-Yang-Mills-Higgs model with uniaxial anisotropy in the group space associated with the Higgs field. We apply this theory to the problem of propagation of color and color-acoustic waves in the gravitational background related to the non-minimal regular Wu-Yang monopole.

The self-consistent system of master equations describing the kinetics of a relativistic non-Abelian plasma, influenced by curvature interactions, is formulated. Non-minimal (curvature induced) coupling is shown to modify all the subsystems of the model: the gauge field equations, the gravity field equations and the kinetic equation for colored particles.

A self-consistent extended Einstein-Maxwell model for relativistic non-stationary polarizable-magnetizable anisotropic media is presented. Based on the analogy with relativistic extended irreversible (transient) thermodynamics, the extended constitutive equations for the electrodynamics of continua are formulated phenomenologically, the convective derivatives of the first, second, etc. orders being taken into account. The master equations for the gravity field contain a modified effective (symmetric) stress-energy tensor of the electromagnetic field in a material medium, the use of this tensor being motivated both by historical analogies and direct variational procedure. By way of example we consider the exact solution of the extended Einstein-Maxwell model, describing the isotropic cosmological model with hidden non-vanishing electromagnetic field, electric polarization and magnetization.

Extended Einstein-Maxwell model and its application to the problem of evolution of magnetized Bianchi-I Universe are considered. The evolution of medium magnetization is governed by a relaxation type extended constitutive equation. The series of exact solutions to the extended master equations is obtained and discussed. The anisotropic expansion of the Bianchi-I Universe is shown to become non-monotonic (accelerated/decelerated) in both principal directions (along the magnetic field and orthogonal to it). A specific type of expansion, the so-called evolution with hidden magnetic field, is shown to appear when the magnetization effectively screens the magnetic field and the latter disappears from the equations for gravitational field.

We discuss exact solutions of three-parameter non-minimal Einstein-Yang-Mills model, which describe the wormholes of a new type. These wormholes are considered to be supported by SU(2)-symmetric Yang-Mills field, non-minimally coupled to gravity, the Wu-Yang ansatz for the gauge field being used. We distinguish between regular solutions, describing traversable non-minimal Wu-Yang wormholes, and black wormholes possessing one or two event horizons. The relation between the asymptotic mass of the regular traversable Wu-Yang wormhole and its throat radius is analysed.

We discuss new exact spherically symmetric static solutions to non-minimally extended Einstein-Yang-Mills equations. The obtained solution to the Yang-Mills subsystem is interpreted as a non-minimal Wu-Yang monopole solution. We focus on the analysis of two classes of the exact solutions to the gravitational field equations. Solutions of the first class belong to the Reissner-Nordstr{\"o}m type, i.e., they are characterized by horizons and by the singularity at the point of origin. The solutions of the second class are regular ones. The horizons and singularities of a new type, the non-minimal ones, are indicated.

A self-consistent non-minimal non-Abelian Einstein-Yang-Mills model, containing three phenomenological coupling constants, is formulated. The ansatz of a vanishing Yang-Mills induction is considered as a particular case of the self-duality requirement for the gauge field. Such an ansatz is shown to allow obtaining an exact solution of the self-consistent set of equations when the space-time has a constant curvature. An example describing a pure magnetic gauge field in the de Sitter cosmological model is discussed in detail.

We establish a new self-consistent system of equations for the gravitational and electromagnetic fields. The procedure is based on a non-minimal non-linear extension of the standard Einstein-Hilbert-Maxwell action. General properties of a three-parameter family of non-minimal linear models are discussed. In addition, we show explicitly, that a static spherically symmetric charged object can be described by a non-minimal model, second order in the derivatives of the metric, when the susceptibility tensor is proportional to the double-dual Riemann tensor

Motivated by the structure of one-loop vacuum polarization effects in curved spacetime we discuss a non-minimal extension of the Einstein-Maxwell equations. This formalism is applied to Bianchi I models with magnetic field. We obtain several exact solutions of the non-minimal system including those which describe an isotropization process. We show that there are inflationary solutions in which the cosmological constant is determined by the non-minimal coupling parameters. Furthermore, we find an isotropic de Sitter solution characterized by a "screening" of the magnetic field as a consequence of the non-minimal coupling.

Using relativistic kinetic theory, we study spherically symmetric, static equilibrium configurations of a collisionless Maxwell-Boltzmann gas with non-standard self-interactions, modelled by an effective one--particle force. The resulting set of equilibrium conditions represents a generalization of the classical Tolman-Oppenheimer-Volkov equations. We specify these conditions for two types of Lorentz--like forces: one coupled to the 4-acceleration and the 4--velocity and the other one coupled to the Riemann tensor. We investigate the weak field limits in each case and show that they lead to various Newtonian type configurations that are different from the usual isothermal sphere characterizing the conventional Newtonian Maxwell--Boltzmann gas. These configurations could provide a plausible phenomenological and theoretical description of galactic dark matter halo structures. We show how the self--interaction may act phenomenologically as an effective cosmological constant and discuss possible connections with Modified Newtonian Dynamics(MOND).

A new exactly solvable model for the evolution of relativistic kinetic system interacting with an internal stochastic reservoir under the influence of a gravitational background expansion is established. This model of self-interaction is based on the relativistic kinetic equation for the distribution function defined in the extended phase space. The supplementary degree of freedom is described by the scalar stochastic variable (Langevin source), which is considered to be the constructive element of the effective one-particle force. The expansion of the Universe is shown to be accelerated for the suitable choice of the non-minimal self-interaction force.

An exact solution to the dynamic equations for a massive boson traveling in a pp-wave gravitational background under the influence of the force induced by curvature, is presented. We focus on the effect of anomalous polarization-curvature interaction and consider models in which the coupling constant of such an interaction is treated to be either a deterministic quantity or a random variable.

A curvature self-interaction of the cosmic gas is shown to mimic a cosmological constant or other forms of dark energy, such as a rolling tachyon condensate or a Chaplygin gas. Any given Hubble rate and deceleration parameter can be traced back to the action of an effective curvature force on the gas particles. This force self-consistently reacts back on the cosmological dynamics. The links between an imperfect fluid description, a kinetic description with effective antifriction forces, and curvature forces, which represent a non-minimal coupling of gravity to matter, are established.

Specific internal self-interactions in gaseous cosmic fluids are shown to give rise to effective negative pressures which may violate the strong energy condition. On this basis we discuss the transition from an initial de Sitter phase to a subsequent Friedmann-Lema\^{\i}tre-Robertson-Walker period as a non-equilibrium configuration of an ultrarelativistic gas. An accelerated expansion of the present universe is obtained as the consequence of a negative internal friction force which is self-consistently exerted on the non-relativistic cold dark matter (CDM) type microscopic constituents of the cosmic gas.

Accelerated expansion of the Universe may result from an anti-frictional force that is self-consistently exerted on cold dark matter (CDM). Cosmic anti-friction is shown to give rise to an effective negative pressure of the cosmic medium. While other models introduce a component of dark energy besides ``standard'' CDM, we resort to a phenomenological one-component model of CDM with internal self-interactions. We demonstrate how the dynamics of the LambdaCDM model may be recovered as a special case of cosmic anti-friction.

Introducing an effective refraction index of an isotropic cosmic medium, we investigate the cosmological fluid dynamics which is consistent with a conformal, timelike symmetry of a corresponding ``optical'' metric. We demonstrate that this kind of symmetry is compatible with the existence of a negative viscous pressure and, consequently, with cosmological entropy production. We establish an exactly solvable model according to which the viscous pressure is a consequence of a self-interacting one-particle force which is self-consistently exerted on the microscopic particles of a relativistic gas. Furthermore, we show that a sufficiently high decay rate of the refraction index of an ultrarelativistic cosmic medium results in an inflationary expansion of the universe.

A covariant criterion for the Cherenkov radiation emission in the field of a non-linear gravitational wave is considered in the framework of exact integrable models of particle dynamics and electromagnetic wave propagation. It is shown that vacuum interacting with curvature can give rise to Cherenkov radiation. The conically shaped spatial distribution of radiation is derived and its basic properties are discussed.

We explain an accelerated expansion of the present universe, suggested from observations of supernovae of type Ia at high redshift, by introducing an anti-frictional force that is self-consistently exerted on the particles of the cosmic substratum. Cosmic anti-friction, which is intimately related to ``particle production'', is shown to give rise to an effective negative pressure of the cosmic medium. While other explanations for an accelerated expansion (cosmological constant, quintessence) introduce a component of dark energy besides ``standard'' cold dark matter (CDM) we resort to a phenomenological one-component model of CDM with internal self-interactions. We demonstrate how the dynamics of the LambdaCDM model may be recovered as a special case of cosmic anti-friction. We discuss the connection with two-component models and obtain an attractor behavior for the ratio of the energy densities of both components which provides a possible phenomenological solution to the coincidence problem.

We describe the ``deflationary'' evolution from an initial de Sitter phase to a subsequent Friedmann-Lema\^{\i}tre-Robertson-Walker (FLRW) period as a specific non-equilibrium configuration of a self-interacting gas. The transition dynamics corresponds to a conformal, timelike symmetry of an ``optical'' metric, characterized by a refraction index of the cosmic medium which continously decreases from a very large initial value to unity in the FLRW phase.

We establish a general thermodynamic scheme for cosmic fluids with internal self-interactions and discuss equilibrium and non-equilibrium aspects of such systems in connection with (generalized) symmetry properties of the cosmological dynamics. As an example we construct an exactly solvable gas dynamical model of a ``deflationary'' transition from an initial de Sitter phase to a subsequent Friedmann-Lema\^{\i}tre-Robertson-Walker (FLRW) period. We demonstrate that this dynamics represents a manifestation of a conformal symmetry of an ``optical'' metric, characterized by a specific effective refraction index of the cosmic medium.

We show that a relativistic gas may be at ``global'' equilibrium in the expanding universe for any equation of state $0 < p \leq \rho /3$, provided that the gas particles move under the influence of a self-interacting, effective one-particle force in between elastic binary collisions. In the force-free limit we recover the equilibrium conditions for ultrarelativistic matter which imply the existence of a conformal timelike Killing vector.

We show that a de Sitter spacetime is a solution of Einstein's field equations with the energy momentum tensor of a self-interacting, classical Maxwell-Boltzmann gas in collisional equilibrium. The self-interaction is described by a four-force which is quadratic in the (spatially projected) particle four-momenta. This force does not preserve the particle number and gives rise to an exponential increase in the comoving entropy of the universe while the temperature of the latter remains constant. These properties of a gas universe are related to the existence of a ``projector-conformal'' timelike Killing vector representing a symmetry which is ``in between'' the symmetries characterized by a Killing vector and those characterized by a conformal Killing vector.

The particles of a classical relativistic gas are supposed to move under the influence of a quasilinear (in the particle four-momenta), self-interacting force inbetween elastic, binary collisions. This force which is completely fixed by the equilibrium conditions of the gas, gives rise to an effective viscous pressure on the fluid phenomenological level. Earlier results concerning the possibility of accelerated expansion of the universe due to cosmological particle production are reinterpreted. A phenomenon such as power law inflation may be traced back to specific self-interacting forces keeping the particles of a gas universe in states of generalized equilibrium.