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 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.

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.

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.