In the present paper we propose a straightforward way to quantization of the gravitational field when the cosmological constant is non-vanishing within the Einstein equations.
In Section 2 we propose an energy-momentum interpretation of the negative Einstein tensor, when it is shifted from the l.h.s. to the r.h.s., so that the total energy-momentum including matter and vacuum (or dark) contribution is identically null.
In Section 3 we represent the metric tensor onto the tetrad of its ortho-normalized eigenvectors and propose a correspondence rule according to which quantization is obtained replacing the eigenvectors with the creation and annihilation operators for the field. Some care is required in evaluating the product of the complex operators in order to provide a Hermitian result, the metric tensor being real. More, it is shown how the space-time interval expectation values are discretized.
In Section 4 we apply the proposed quantization method to physical solutions such as Schwarzschild-De Sitter, Robertson-Walker-De Sitter and Kerr-De Sitter metrics.
In Section 5 some conclusions and perspectives are suggested.
2. Energy-Momentum Interpretation of the Equations of the Gravitational Field
Let us consider the Einstein field equations of the gravitational field, including the cosmological constant   in presence of observable matter and non-gravitational interaction fields:
the universal constant being the Einstein gravitational constant equal to:
with , Newton gravitational constant.
As usual is the metric tensor of signature , , is the Ricci tensor, where:
of trace , being the contravariant components of the metric tensor (inverse matrix):
the symmetric connection coefficients (Christoffel symbols) in a space-time torsionless manifold. Moreover is the energy-momentum tensor of the observable matter and non-gravitational fields, while is the cosmological constant.
On introducing the notation , where:
which we may interpret as the energy-momentum tensor of the gravitational field, we can write the field Equation (1) as an energy-momentum balance equation.
More, it results to be convenient to split the gravitational contribution in two parts, the former due to ordinary attractive gravitation (geometric), and the latter due to the repulsive gravitation (cosmological or dark)    :
Of course, in absence of observable matter and non-gravitational fields the energy-momentum components of the full gravitational field vanish.
As we will see soon, the representation (7) of the field equations has the non-trivial advantage of facilitating quantization of the gravitational field.
We emphasize that, as a consequence of the equivalence principle, the description of a physical gravitational field within a pseudo-Euclidean (i.e., Minkowskian) space-time is equivalent to the geometry of a suitable non-Euclidean (i.e., Riemannian) space-time. Therefore the Einstein equations may be interpreted either geometrically as a set of partial differential equations governing metric, connection and curvature of space-time ( ), or as energetic conditions ( ) governing the energy-momentum of a physical gravitational field living within a flat space-time. Adopting the latter interpretation we have split the energy-momentum tensor of the free physical gravitational field (which is null because of the Einstein equations and then results to be non-localizable) into two non-vanishing localizable contributions ( ) and ( ). Quantization will be performed on the Hamiltonian of the cosmological part, so providing quantization also of the opposite Riemannian part. It is remarkable that even if the null Hamiltonian of the whole gravitational field cannot be quantized being non-localizable, each one of its half opposite contributions can be quantized separately.
3. Approaching the Quantization of the Gravitational Field
As a first step, we will be concerned with the gravitation in absence of external non-gravitational fields, i.e., when:
Then Equation (7) becomes simply:
In order to open a way towards a method for gravity quantization resembling the well known procedures for Abelian (electromagnetic) and non-Abelian Yang-Mills (weak and strong interaction) fields, we represent the metric tensor on the tetrad of its ortho-normalized eigenvectors :
Then the components of the energy-momentum tensor result:
Let us now introduce the notation:
where V is any constant volume in the physical space and the reduced Planck constant. We emphasize that the eigenvalues are positive, the signature of the metric tensor being involved in .
It follows in (16):
In order to quantization we assume the following correspondence relation:
where are complex operators, and replace the non-quantized relation (19) with the correspondent quantum equation:
symmetrization being required to preserve hermiticity of the energy-momentum operator . The relation between the creation and annihilation operators with the co-ordinates and momenta operators:
provides the energy-momentum tensor for a set of harmonic oscillators:
Definitions (22) ensure, as usual in quantum field theory, the commutation relation:
The the operators result to be the creation and respectively annihilation operators for the quantized gravitational field. The inidces are to be interpreted as labels of the oscillation modes. Finally we obtain:
The Hamiltonian operator follows by integration on the region of volume V of the time-time component :
Or expanding the sum:
The previous results provide soon, thanks to (11) also the information, holding in empty space-time:
and in presence of matter:
the r.h.s. terms been known.
Eventually we point out that, after quantization, the metric tensor itself becomes an operator, resulting:
Therefore the interval becomes an operator too:
The expectation values of discretized interval are given by:
4. Applications to Physical Solutions
The previous results are easily applied to the known physical solutions to Einstein equations in presence of cosmological constant.
1) Schwarzschild-De Sitter metric (spherical co-ordinates )
In correspondence to the diagonal Schwarzschild-De Sitter solution for which it results:
the characteristic frequencies provided by (18) become:
Then the proper frequencies of the 4 characteristic modes depend on the mass M of the oscillating body, the volume V into which it is confined and on the cosmological constant.
2) Robertson-Walker-De Sitter metric (spherical coordinates )
In correspondence to the diagonal Robertson-Walker-De Sitter solution we have:
and the characteristic frequencies provided by (18) become:
Here, the proper frequencies of oscillation of the observed region of the universe depend on the evolution function , the volume of the region and, of course on the cosmological constant.
3) Kerr-De Sitter metric (untwisted co-ordinates )
In untwisted co-ordinates   the Kerr metric reduces to a Schwarzschild-De Sitter like metric and the results are the formally similar.
The frequencies are:
Now, the proper frequencies of oscillation of the observed region of the universe depend on the volume V of the region and on the cosmological constant, the mass being hidden within the resized ray coordinate y.
The proposed approach could offer a suggestion also for quantization of multidimensional gravity in space-time manifolds  , and  or higher power theories of gravitation .
We presented a very simple method to approach the hard problem of quantization of the gravitational field, suggesting a way based on the presence of a non-vanishing cosmological constant, a condition which is today universally recognized. Quantization of the cosmological contribution owed to vacuum energy-momentum provides immediately also the quantization of the geometric gravity which is opposite in sign respect to the vacuum term, possibly incremented by other sources of matter and non-gravitational fields.
Our approach could be applied, in further investigations, even to generalized theories of gravitation.
I sincerely thank the referee for suggesting hints to improve the paper.
 Ellis, G.F.R. (2003) A Historical Review of How the Cosmological Constant Has Fared in General Relativity and Cosmology. Chaos, Solitons & Fractals, 16, 505-512.
 Farnes, J.S. (2018) A Unifying Theory of Dark Energy and Dark Matter: Negative Masses and Matter Creation within a Modified ΛCDM Framework. Astronomy & Astrophysics, 620, A92.
 Arbey, A. and Mahmoudi, F. (2021) Dark Matter and the Early Universe: A Review. Journal Reference: Progress in Particle and Nuclear Physics, 119, Article ID: 103865.
 Strumia, A. (2020) Fundamental Fields Included into Geometry in Sixteen Space-time Dimensions. In: Theory and Practice of Mathematics and Computer Science, Vol. 4. Book Publisher International.