Researchers began to question the theory of gravity after the advent of the theory of general relativity (GR). Weyl (1919) and Eddington (1923) considered modifications to the theory by including higher-order invariants in the action  -  . In the 1960s, the more complex gravitational action with many advantages appeared. GR was not renormalizable at that time and cannot be conventionally quantized. In 1962, Utiyama and De Witt showed that renormalization at one loop demands the Einstein Hilbert action be supplemented by high-order curvature  . Then, Stelle showed that higher-order actions are renormalizable  . There is higher-order curvature action for the effective low-energy gravitational action, when quantum corrections or string theory are considered    .
GR correction is not an easy task. First, there are many naive GR corrections, which are unrealistic  -  . The best-known example is an alternative to GR, i.e. scalar-tensor theory. There are also many methods of gravity correction  -  . The typical examples are Dvali-Gabadadze-Porrati gravity  , brane- world gravity  , vector-scalar tensor theory  and Einstein-Aether theory  . However, there are many different gravity theories. These theories are the summary of the Einstein-Hilbert action,
where, G is the gravitational constant, g is the determinant of the metric and R is the Ricci scalar (). With the general function of R, the following is obtained.
There are two motivations for GR correction: (1) adding higher-order gravitational action in high-energy physics, and (2) applying the GR correction to cosmology and astrophysics.
However, there are two problems. The first problem is why specifically actions and not more general ones, which include other higher-order invariants, such as. The answer is twofold. First, actions are sufficiently general to encapsulate some of the basic characteristics of higher-order gravity, and at the same time they are simple enough to be easily handle. For instance, viewing as a series expansion, i.e.,
When the coefficients have certain values, the action shows some interesting phenomenology. In sum, there are some advantages of theory in gravity correction. It can help the high-order gravitational theory to avoid the fatal Ostrogradski instability   .
The second problem is that it is related to a possible loophole. First, how can high-energy corrections of the gravitational action have nothing to do with late- time cosmological phenomenology? Would not effective field theory considerations require that the coefficients in Equation (3) be such, as to make any modifications to the standard Einstein-Hilbert term important only near the Planck scale?  
The observed late-time acceleration of the Universe poses one challenge to theoretical physics. In principle, this phenomenon may be the result of unknown physical processes. For instance, it involves either the correction of gravitational theory or the existence of new fields in high-energy physics. Although the latter one is most commonly used, the correction of gravitational theory is an attractive and complementary approach to explain this phenomenon, known as gravity  . Some researchers added the weight of Ricci scalar in the Einstein-Hilbert Lagrangian for the GR correction  -  .
In  , a gravitational theory of a scalar field ϕ with non-minimal derivative coupling to curvature is considered. The results show that a cosmological model with non-minimal derivative coupling is able to explain in a unique manner both a quasi-de Sitter phase and an exit from it without any fine-tuned potential. In  , the authors approached the problem of testing dark energy and alternative gravity models to general relativity by cosmography. The results show that degeneration among parameters can be removed by accurate data analysis of large data samples and also present the examples.
Several different forms for have been suggested in the literatures  -  . These different -gravity theories have also been discussed in the stability conditions  -  , inflationary epoch  , compatibility with solar-system tests and galactic data  -  and the late-time cosmological evolution  -  . Additional constraints to theories may be caused by imposing the so-called energy conditions    , for example, the phantom field potentials  , expansion history of the Universe      , as well as evolution of the deceleration parameter and their confrontation with supernovae observations   . If gravity is considered as a step toward a more complicated theory, that which generalization would be more straightforward will depend on the chosen representation (see also Sotiriou et al., 2008 for a discussion)  -  .
In this paper, we introduce a new action, the new action effective amount of inclusion of Quintessence and Phantom can solve the problem. Thus, this model is a more general model of dark energy. We obtained the general sound speed in the evolutions of the Universe, and give the exact expressions for the exact energy- momentum tensor, pressure, energy, and different gravity theories.
The rest of this paper is organized as follows. Section 2 investigates a general action of gravity and Phantom and Quintessence fields coupled to scalar curvature. Section 3 studies the exact energy-momentum tensor and sound speed of the new general single field action of gravity. Sect.4 shows different cosmological evolutions with single field inflation. Section 5 studies general gravity theory with general matter action and its applications. Section 6 gives general gravity theory with many general Phantom and Quintessence fields. Section 7 presents summary and conclusions.
2. General Action of f(R) Gravity and Phantom and Quintessence Fields Coupled to the Scalar Curvature
The first idea is to combine the actions for Phantom and Quintessence fields into one action, by adding a parameter α into the general action of gravity. Thus, we have a new action as follows:
(i) When, Equation (4) is just a general action of gravity and Phantom field;
(ii) When, Equation (4) is a general action of gravity and Quintessence field.
Further, the scalar curvature R can be coupled to the Phantom and Quintessence fields. Thus, the second idea is to add scalar curvature R into Equation (4)
where is related to the scalar curvature. Because general characteristics of the scalar curvature are coupled to Phantom and Quintessence fields, the scalar and curvature are naturally coupled. Therefore, Equation (5) is a new action with gravity, Phantom and Quintessence fields and scalar curvature R. Compared with our previous researches, this action provides the possibility to explore the energy-momentum tensor and sound speed of the single field action of gravity. A general function of scalar curvature (i.e.) is considered for Phantom and Quintessence fields. It is very natural to consider coupling scalar curvature R to matter field, e.g., Phantom and Quintessence fields. Because Equation (5) satisfies the invariance of the general coordinate transformation (specially coupling of R to scalar field is the constant Einstein-Hilbert action, which has all invariant properties of Einstein’s RG. Its coupling strength can be fitted by coupling parameter from physical experiments), it means that the whole laws remain effective in the whole spacetime and are important in field theories. Thus, Equation (5) is consistent.
The variance of Equation (5) is
For the deducing details see Appendix A.
To study Equation (6), we discuss generally as follows 
We rewrite Equation (9) as
Substituting Equation (11) into Equation (8), we have
The deducing details of Equation (12) are in Appendix B.
Substituting Equation (12) into Equation (7), we obtain
Substituting Equation (13) into Equation (6), we have
Considering a general partial integration, we have
The deducing details of Equation (15) are in Appendix C. Thus, we have
Substituting Equation (16) into Equation (14), we have
Using Equation (17), we deduce Euler-Lagrange Equations
where and. we can rewrite Equation (18) as
Therefore, a general action of gravity and Phantom and Quintessence fields is generally presented. We generalize Equation (4) to a general action of gravity and Phantom and Quintessence fields coupled to scalar curvature by making variance of the general Lagrangian. We further deduce Euler-La- grange Equations of different fields, matter tensor, effective matter tensor, etc.
3. Exact Energy-Momentum Tensor and Sound Speed of the New General Single Field Action of f(R) Gravity
From Equation (5) and Equations (20)-(23), we can obtain the exact energy- momentum tensor.
To simplify the exact energy-momentum tensor, we can rewrite Equation (24) as
where. Similar to  , we can generally define
Substituting Equation (26) into Equation (25), we obtain a general energy- momentum tensor
where is different from the past expression.
In general, using, and Friedman- Robertson-Walker metric
Due to, we have. Further, using Equations (29) and (32), we obtain a new general density related to as follows
Furthermore, we can obtain the general sound speed
Using Equations (30) and (31), we obtain
Thus, we can rewrite Equations (33) and (34) as
where. The deducing details are in Appendix D. And then we have
In conclusion, we obtain the general energy-momentum tensor, pressure, density and speed sound of the new general single field action of gravity and Phantom and Quintessence fields coupled to scalar curvature.
4. Different Cosmological Evolutions with Single Field Inflation
The most hopeful models for the different evolutions of the Universe are the cosmological models of the initial evolution and subsequent development. They are supported by the most comprehensive and accurate explanations based on the current scientific evidences and observations. According to the observations on the current Universe, the dark matter accounts for 24% of the mass-energy density of the observable Universe, the dark energy amounts for nearly 72% and the ordinary matter only accounts for about 4%.
The equation of state (EOS) is a powerful way to describe the matter and the evolutions of the Universe. In cosmology, the EOS of a perfect fluid is characterized by a dimensionless number that is equal to the ratio of its pressure to its energy density. It is closely related to the thermodynamic EOS and ideal gas law.
Therefore, with EOSs of matter, dark energy and dark matter, the new general action can be used to explain the different cosmological evolutions: (I) Big Rip Universe; (II) De Sitter Universe; (III) Harmonic Universe.
In the case of cosmological inflation, using Equations (28), (30) and (5), we deduce the EOS
Substituting into Equation (35), it follows that
To discuss the different evolutions of the Universe, we use the Friedman equations 
and then we have
Now accelerating expansion requires smallness of the variation of the Hubble parameter, as defined by the parameter 
Using Equation (33), we rewrite Equation (35) as
Equation (42) can describe the different evolution characteristics of the Universe:
(i) Big Rip Universe, i.e., the case where, the inflation and the accelerating expansion of the Universe can be characterized by the EOS of dark energy.
Generally, the expansion of the Universe is accelerating when.
When EOS for Phantom energy is, the Big Rip will occur. According to Equation (42), the total pressure of our Universe is negative, and then the relation between pressure P and is deduced as follows
(ii) De Sitter Universe, i.e. the case where, the relation between pressure P and is
The de Sitter Universe is a solution to Einstein’s field equations of General Relativity, which is named after Willem de Sitter. When one considers the Universe as spatially flat and neglects ordinary matter, the dynamics of the Universe would be dominated by the cosmological constant, corresponding to dark energy. If the current acceleration of our Universe is due to a cosmological constant, then the universe would continue to expand. All the matter and radiation will be diluted. Eventually, there will be almost nothing left except the cosmological constant, and the Universe will become a de Sitter Universe.
(iii) Harmonic Universe, i.e., the case where, the relation between pressure P and is
The EOS for ordinary non-relativistic matter is, which means that it is diluted as. This is natural for ordinary non-relativistic matter.
The EOS of radiation and matter is in the very early Universe, and
then the Universe is diluted as. In the expanding universe, the energy density decreases more quickly than the volume expansion.
Substituting Equation (26) into Equation (42), we can further deduce a concrete expression for Equation (42)
which shows their concrete evolution details. For example, we focus on the case of different potentials as follows:
(a) When, the EOS is
(b) When, the EOS is
(c) When, the EOS is
Thus, we deduce the EOS, i.e., Equations (35), (36), (42), (47)-(49), of gravity and Phantom and Quintessence fields coupled to scalar curvature. We investigate different cosmological evolutions with single field inflation including the Big Rip Universe, De Sitter Universe and Harmonic Universe. We further study the cases of different potentials for different EOSs.
5. General f(R) Gravity Theory with General Matter Action and Its Applications
We begin with a general action of gravity plus general matter action
where is the determinant of the metric, is the Ricci scalar and is a general matter action. Varying this action with respect to the metric, we obtain the field equation  - 
where a prime denotes differentiation with respect to and. To use an approach to the null energy condition (NEC) and strong energy condition (SEC) similar to that in general relativity (GR) context, we note that Equation (51) can be rewritten as
Multiplying Equation (53) by, we have
Then, we get
In addition, we can have
Using Equations (56) and (57), we deduce
Combining Equation (56) with Equation (58), we have
The deducing details are in Appendix E.
Equations (56) and (59) are consistent. This is because by multiplying Equation (59) with we have
The deducing details are in Appendix F. Equation (60) is just Equation (56), so Equations (56) and (59) are consistent.
In addition, we have
When, the virtue of Einstein’s Equation (57) implies
where is an effective Equation (59).
Using Equation (61), we similarly have
For the homogeneous and isotropic Friedman-Lemaitre-Robertson-Walker (FLRW) metric with scale factor and for a perfect fluid, and  . Substituting Equation (61) into Equation (62), we have
where and. (For the deducing details see Appendix G). Therefore, we obtain an effective SEC
The evolution equation for the expansion of a null geodesic congruence is defined by a vector field, which has the same form as the Raychaudhuri equation, with in the last term. In this case, the condition for the convergence (geodesic focusing) of hyper-surface orthogonal congruences of null geodesics along with Einsteins’s equation implies
where. The deducing details are in Appendix H. Therefore, we obtain an effective NEC
Using Equation (59), we concretely have
Substituting Equations (68) and (69) into Equations (65) and (67), we deduce inequalities
For the deducing details see Appendix I.
Thus, we have
For, Equations (70) and (71) give
where. Therefore, Equations (68) and (69) can return to the well-known forms of the SEC and NEC in GR      . Therefore, the above investigations are consistent.
In sum, we investigate a general gravity with general matter action and obtain the different field equations, general matter tensor and effective matter tensor. Further, we obtain an effective SEC and an effective NEC. When approaches to R, the effective SEC and the effective NEC approaches, respectively, to the usual SEC and the usual NEC. Thus, these studies are consistent.
6. General f(R) Gravity Theory with General Scalar Fields
To investigate more general cases and extend the applications of the new action proposed in Section 2, the single scalar field is changed into a more general form .
We now generalize Equation (4) to a general action of gravity, and many general scalar fields are coupled to scalar curvature as follows.
When and, Equation (74) is a generic action of gravity and Quintessence fields as follows.
Equation (76) is a generic Lagrangian of gravity and Quintessence fields.
When and, Equation (74) is a generic action of gravity and Phantom fields as follows
Equation (78) is a generic Lagrangian of gravity and Phantom fields.
Using Equation (74) or Equation (77), we can calculate and obtain all corresponding results similar to all investigations above Equation (74).
7. Summary and Conclusions
This paper introduces the development process of gravity theory, which is based on the Quintessence, Phantom and theories  . We generally present a general action of gravity, Phantom and Quintessence fields and generalizes the action to a new general action of gravity, Phantom and Quintessence fields coupled to scalar curvature. Further, we deduce Euler- Lagrange Equations of different fields, and give matter tensor and effective matter tensor.
Then, this paper obtains the general energy-momentum tensor, pressure, density and speed sound of the new general single field action of gravity, Phantom and Quintessence fields coupled to scalar curvature and so on.
Further, this paper deduces the equations of states of gravity, Phantom and Quintessence fields coupled to scalar curvature. We also investigate different cosmological evolutions with single field inflation including the Big Rip Universe, De Sitter Universe and Harmonic Universe. Besides, we study the cases of different potentials for different equations of states.
In addition, this paper investigates a general gravity theory with general matter action and obtains the different field equations, general matter tensor and effective matter tensor. We further obtain an effective Strong Energy Condition and an effective Null Energy Condition. We prove that when f(R) approaches to R, the effective Strong Energy Condition and the effective Null Energy Condition approach, respectively, to the usual Strong Energy Condition and the usual Null Energy Condition. Thus, these investigations are consistent.
The Hawking-Penrose singularity theorems invoke the weak and strong energy conditions, whereas the proof of the second law of black hole thermodynamics requires the null energy condition. More recently, several researchers used the classical energy conditions of GR to investigate cosmological issues.
In the cosmology, these theories provide an alternative way to explain the cosmic speed-up. The freedom in building different functional forms of causes the problem of how to constrain these many possible gravities from theoretical or observational aspects. Recently, this possibility has been explored by testing the cosmological viability of some specific forms of gravities. Specially, all calculations with appendixes in this paper are helpful for beginners in this field.
Finally, we generalize Equation (4) to a general action of gravity and many scalar fields, and obtain a general action of gravity, Quintessence field and Phantom field. Then more applications can be done. For levels and multi-field coupling, we will carry out further study  .
We would like to thank Prof. Rong-Gen Cai for the helpful discussions and comments.
The Conflict of Interests
The authors declare no conflict of interest regarding the publication of this paper.
The work is supported by National Natural Science Foundation of China (Grant Nos. 11275017 and 11173028).