Currently, the two greatest mysteries in physics seem to be dark matter and dark energy. Almost all models of dark matter resort to largely-unspecified never-discovered subatomic particles or introduce new physical laws  - , and references therein. One exception is papers showing the possibility that dark matter (or at least a part of it) can consist of hydrogen atoms of the 2nd flavor, whose existence is proven by the analysis of atomic experiments and can also explain the latest observations of 21 cm radio line from the early Universe  .
As for dark energy, there are lots of field and particle candidates for it (in addition again to suggestions to introduce new gravitational degrees of freedom, of which massive gravity is one of the examples), but these candidates have never been discovered so far  - , and references therein. This is even more intriguing because the commonly accepted view is that the Universe is currently in the dark energy dominance era (estimated to start about 5 billion years ago)—the era where yet unknown dark energy dominates over the gravitation and is responsible for the observed acceleration of the Universe expansion.
In the present paper, we introduce a scenario providing a possible explanation of the entire history of the Universe expansion: both the era of the decelerating expansion and the current era of the accelerated expansion. In this model, there is no need to resort to dark energy or to new gravitational degrees of freedom (mass gravity etc.).
We present this scenario as the following sequence of logical steps.
1) We consider a “gas” of a large number of gravitating neutral nonrelativistic particles of mass m. The corresponding Hamiltonian is:
where ri and pi are the radius-vector and the momentum of particle number i, respectively; G is the gravitational constant.
At any instant of time, the gas has a subsystem of relatively isolated pairs of particles: a pair of particles is considered relatively isolated if their separation is much smaller than their distance to other particles. The subsystem is open: in the course of time, some pairs could stop qualifying as the subsystem members (because they can no more be considered as relatively isolated), while some other pairs could become relatively isolated and qualify as new members of the subsystem. Here the word “subsystem” does not mean a localized part of the entire system: rather it means a subset of the particle within the ensemble—the subset of pairs (not located in one particular volume) that are relatively isolated.
2) Pairs of the positive total energy (E > 0) lose the energy by the gravitational radiation and reach the state of E = 0. For clarity: we consider non-relativistic particles and by their total energy we mean the difference between the relativistic total energy and the rest mass multiplied by the square of the speed of light.
3) To the motion of pairs of E = 0 applies the virial theorem—contrary to many textbooks stating that the virial theorem applies only to bound states. As shown by Pollard , the boundness of a state is not the necessary condition for the applicability of the virial theorem. He derived a more general condition that is both necessary and sufficient. He specifically showed that the virial theorem applies to pairs of E = 0 interacting by Newton’s gravity.
4) Therefore, the subsystem of pairs of E = 0 is characterized by the following values of the average potential ( ) and kinetic ( ) energies: .
5) Now we apply Dirac’s Generalized Hamiltonian Dynamics (GHD) to the above subsystem of pairs. Dirac designed the GHD for applications to quantum field theory   . The central point of Dirac’s GHD is the addition of a linear combination of constraints to the classical Hamilton function. A brief description of Dirac’s GHD and of its further development is presented in Appendix A.
Oks and Uzer  demonstrated how to apply Dirac’s GHD to microscopic systems of particles, such as e.g., atomic and molecular systems. Specifically, as constraints they used classical integrals of the motion. By doing so, they showed, e.g., that hydrogenic atoms/ions, being described by Dirac’s GHD, exhibit a sequence of classical non-radiating states, i.e., states where the system does not emit the electromagnetic radiation—it is radiationally-stable. The set of energies of these classical non-radiating states turned out to exactly coincide with the corresponding quantum stationary states.
The most important feature of classical non-radiating stable states is the following. In such a state, the radius-vector of the atomic electron r(t) = r0 and the generalized momentum P(t) = P0, where r0 and P0 are some vector constants. Thus, the electron velocity v = 0, but its generalized momentum P0generally differs from zero. This should not be shocking: the generalized momentum is a more complex physical quantity than the velocity. For instance, it is well-known that for a charge e of mass m in an electromagnetic field, characterized by a vector-potential A, it is also possible to have v = P/m – eA/(mc) = 0 while P = eA/c ≠ 0 , the formula after Equation (16.9).
6) The way Oks & Uzer  applied Dirac’s GHD to pairs of charges interacting by the Coulomb law can be similarly applied to pairs of neutral particles interacting by Newton’s gravity. The result is the infinite number of classical non-radiating states of the energies
, , . (2)
Here k is the harmonic number of the primary frequency of the particle moving along an elliptical trajectory. It is convenient to denote q = 1/k, so that
, . (3)
It is obvious that q takes the infinite amount of rational numbers in the segment [0, 1] and that .
7) Thus, within Dirac’s GHD, the pairs of gravitating neutral particles of energies E > 0, lose the energy by the gravitational radiation, reach the state of E = 0, and in this state the gravitational radiation vanishes. The subsystem of pairs of E = 0 is radiationally stable.
8) While the subsystem of pairs of E = 0 has , the rest of the gas has (the subscript “rest” stands for the “rest of gas”). The two subsystems interact and can exchange members. Therefore, the entire system tends to thermalize (due to the exchange of members of the two subsystems), so that for the entire system we have the average kinetic energy .
9) The average kinetic energy of the rest of the gas will continue to decrease in the course of time so that the entire gas will have the average kinetic energy . This means that in the course of time, practically the entire gas would become thermodynamically equivalent to a set of non-radiating pairs of E = 0.
10) Here we come to the following central point. The “freezing” of the pairs into classical states of r = const is equivalent to a partial inhibition of the gravitation. Indeed, without it, the gravitational radiation and the loss of the total energy would continue with 100% effectiveness and the pairs would relatively rapidly collapse (being considered classically). However, the partial inhibition of the gravitation slows down this process.
11) Thus, as the share of the subsystem of pairs (in the entire system) increases in the course of time, the gravitational interaction within the entire system effectively decreases and tends to zero.
12) For this scenario, the particles of the gas should have a practically infinite lifetime and should have zero or very little interaction with the rest of the matter. Therefore, the most probable candidate for these particles is non-relativistic neutrinos. The current technology allows detecting only relativistic neutrinos—because slow-moving neutrinos have very low probabilities of interactions with the rest of the matter. However, nothing prohibits neutrinos from moving at non-relativistic velocities.
13) According to the recent study by Shamir , there is a disparity between galaxies spinning clockwise and counterclockwise so that the Universe as the whole rotates. Thus, there is a centrifugal force. After the gravitation would get partially inhibited at the above stage 10, the centrifugal force would start dominating over the gravitation and the particles of the rotating gas would expand at an increasing speed.
We considered a “gas” of a large number of gravitating neutral particles. We presented a step-by-step scenario, one of the central points of which was the application of Dirac’s GHD to pairs of these particles. Another central point was that the virial theorem is valid for applies to pairs of zero total energy interacting by Newton’s gravity—contrary to many textbooks stating that the virial theorem applies only to bound states. We showed that as a result, the gravitational interaction within the entire system effectively decreases.
The following point should be emphasized. Just the Universe rotation alone cannot explain why in the first 9 billion years or so, the expansion of the Universe was decelerating. The entire situation can be explained by the combination of the inhibition of the gravitational interaction (described above) complemented by the rotation of the Universe. In other words, it took time for the gas to reach the above stage 10 (the “freezing” of the gravitation). During this time, the gravitation was still dominating the opposing centrifugal force caused by the rotation of the Universe, so that the Universe expansion was slowing down. Then, after the gas reached the above stage 10, the centrifugal force, caused by the Universe rotation, started dominating the gravitation, so that the Universe expansion started to accelerate.
Thus, in this model, there is no need to resort to dark energy or to new gravitational degrees of freedom (mass gravity, etc.) for explaining the entire history of the Universe expansion, including the observed acceleration of the expansion in the last 5 billion years or so.
Appendix A. Brief Description of the Application of Dirac’s Generalized Hamiltonian Dynamics (GHD) to Microscopic Systems of Particles
The central point of Dirac’s GHD is the introduction of a classical generalized Hamilton function
, , (A.1)
where is the conventional classical Hamilton function and the 2nd term is a linear combination of constraints (Dirac    ). Here and below, the summation over twice repeated suffix is understood. The coefficients can be determined from so-called consistency conditions
where [f, g] is the usual Poisson bracket. (To distinguish equations that hold after the variation of dynamical variables from equations that do not hold, Dirac used the sign “=” for the former and the sign “≈” for the latter.)
 further developed Dirac’s GHD by using classical integrals of the motion as constraints  . Specifically, for hydrogenic atoms/ions they used components of the angular momentum and of the Runge-Lenz vector as the constraints. This turned out to be sufficient for proving the existence of classical non-radiating states.
For obtaining further results, Oks & Uzer  employed Plank’s hypothesis (later incorporated by Einstein into his photoelectric law  ), which can be formulated as follows. In the process of the radiation, the change of the energy between the initial and final states DE should be proportional to the average frequency of the radiation in this process. As a result, Oks & Uzer obtained the expression for the infinite set of classical non-radiating states Ek (where ) coinciding with the corresponding set of quantum stationary states.
Later on, Dirac’s GHD was successfully applied to other microscopic systems and to finding classical non-radiating states in these systems. Example is pairs of particles interacting by a modified Coulomb potential  and a spherical harmonic oscillator .
 de Martino, I., Chakrabarty, S.S., Cesare, V., Gallo, A., Ostorero, L. and Diaferio, A. (2020) Dark Matters on the Scale of Galaxies. Universe, 6, Article No. 107.
 Oks, E. (2020) Alternative Kind of Hydrogen Atoms as a Possible Explanation of the Latest Puzzling Observation of the 21 cm Radio Line from the Early Universe. Research in Astronomy and Astrophysics, 20, Article No. 109.
 Brax, P. (2017) What Makes the Universe Accelerate? A Review on What Dark Energy Could Be and How to Test It. Reports on Progress in Physics, 81, Article No. 016902. https://doi.org/10.1088/1361-6633/aa8e64
 Oks, E. and Uzer, T. (2002) Application of Dirac’s Generalized Hamiltonian Dynamics to Atomic and Molecular Systems. Journal of Physics B: Atomic, Molecular and Optical Physics, 35, 165-173.
 Shamir, L. (2020) Patterns of Galaxy Spin Directions in SDSS and Pan-STARRS Show Parity Violation and Multipoles. Astrophysics and Space Science, 365, Article No. 136. https://doi.org/10.1007/s10509-020-03850-1