The theory of solitons in general relativity was first elaborated by G. N. Shikin in 1995. He formulated the requirements to be fulfilled by solitons . His research work allowed an intensive study on soliton in general relativity by many authors. In a series of papers, exact plane-symmetric solutions to the spinor and gravitational field equations have been obtained. The role of gravitational field and the nonlinear terms in the formation of the field configurations with limited total energy, spin and charge have been thoroughly investigated. Let us emphasize that the total charge and spin of the self-consistent system of spinor and gravitational field equations are unlimited. The divergence of the total charge and spin is related to the non-consideration of the torsion and the properties of the static plane-symmetric utilized   . For an excellent review of the interacting scalar and spinor fields in plane-symmetric metric refer to   . Note that in    the obtained solutions are exact regular with localized energy density and finite total energy. Nevertheless, the total charge and spin are not bounded. The unlimited problem of the charge and the spin is resolved in a series of interesting articles  - . The gravitational field is given by a spherically symmetric metric.
The aim of the paper was to study the role of the interaction of nonlinear spinor, scalar and gravitational fields in the formation of configurations with localized energy density and limited energy, spin and charge of the spinor field.
The paper is organized as follows. Section 2 deals with general relativistic equations. The lagrangian of the self-consistent interaction spinor, scalar and gravitational fields has been defined. The gravitational field in our study is given by a spherically symmetric metric. Section 3 addresses the main results. The fundamental field equations have been solved and their localisation properties have been examined. Section 4 deals with the discussion of the main results. Thus, we studied in detail the spinor, scalar and gravitational interaction equations by choosing the concrete form of the arbitrary function under the form . Finally, concluding remarks and future work are outlined in Section 5.
2. General Relativistic Equations
The lagrangian of the self-consistent system of interaction between spinor, scalar and gravitational fields may be written under the following form :
; ; and correspond respectively to gravitational field lagrangian, free spinor field lagrangian, free scalar field lagrangian and interaction lagrangian. They are defined as follows:
R denotes the scalar curvature; is Einstein’s gravitational constant; G is Newton’s universal gravitational constant; c is the velocity of light; is an arbitrary function of the invariant .
In the sequel, in ordor to simplify the expressions, we shall consider:
The grvitational field in our case is given by a spherical symmetric space-time via the metric which is defined under the following form:
For simplicity reason, the speed of light has been taken to be unity (c = 1). The metric functions , and are time and angular coordinates and independent. They are functions of spatial variable alone which is defined as in  in the form , where r stands for the radial component of the spherical symmetric metric. These metric functions obey the harmonic coordinates conditions:
Variation of (1) with respect to the spinor field and its conjugate , we establish the nonlinear spinor field equations as follows:
Similarly, varying of (1) with respect to the scalar field we obtain the following scalar field equation:
Then, the general form of Einstein’s field equation is:
In virtue of (8) and (9), expression (13) becomes:
where prime ( ) denotes differentiation with respect to , is the energy momentum tensor of the spinor, scalar fields and its interaction, is Einstein’s tensor, is Ricci’s tensor and is Kronecker’s symbol which is 0 if and 1 if .
The metric energy-momentum tensor of the interaction of the spinor and scalar fields field can be written as follows:
Taking into account (10) and (11), we rewrite under the form:
Taking into account (19) and (21), let us try to rewrite explicitly the non null components of the metric energy-momentum tensor . In this optic, we obtain:
In expressions (10), (11) and (19), is the derivative covariant of spinor. It is connected to the spinor affine connection matrices  :
In the above equations, are Dirac’s matrices in curved space-time. In order to define , let us use the equalities:
where is the metric of Minkowski and are tetradic 4-vectors.
With the relation (25), we have:
For the matrices in flat space-time, we take :
At present, let us define the affine connection matrices of the spinor. To this end, the general form of is:
In (27), are Christoffel’s symbols. From expression (27), we have the components of the affine connection matrices of the spinor in the following way:
Equation (12) has solution:
According to Einstein’s convention of sommation, we get:
Then, using expressions (24), (29) and (30), we can rewrite Equations (10) and (11) as follows:
Further setting with , from (31), we obtain the following set of equations:
where function is:
The functions , , and are connected by the relation:
The following section deals with the main results.
3. Main Results
Summing the set of Equations (33)-(36), we infer that the invariant function:
satisfies a firt order differential equation as follows:
Equation (40) has solution:
Expression (41), reflects the natural link between the nonlinear spinor field of elementary particles and their own gravitational field.
By Considering the spinor field equation in the form (31) and the conjugate one in the form (32), from (23), the momentum may be rewritten in the following form:
In the following paragraph, we shall solve Einstein’s field equations in order to determine the general expressions of the metric functions , and and then the link which exits between them.
In this perspective, in virtue of , substraction of Einstein’s Equations (14) and (16) leads to the following equation:
The previous equation can be transformed into a Liouville equation type to produce the following solutions (refer to , page 30):
A and D are integration constants. T is a function. It is defined under the form:
where h and are integration constants.
From (44), (45) and (9), one finds the explicit form of the metric function as follows:
Then, the metric functions , and are connected as follows:
Let us note that Einstein’s Equation (15) is a first integral of Equations (14) and (16). It is also a first order differential equation. By substituting, (48) and (42) into (15), we obtain:
Substituting and , into (49), we have:
Equation (50) has the general analytical solutions under the form:
Setting a concrete form of the function , from (7) we can find . Knowing , we can determine , from (37). Thus, from (51), the invariant function can be found. Then, if is known, we can determine from (41). Finally, the metric functions and can be determined from expression (48).
Considering the concrete expression of the invariant function , we can establish the regularity properties of the solutions obtained. Studying the distribution of the energy per unit invariant volume , we can establish the localization properties of the solutions.
In the sequel, we shall get a concrete form of the functions by solving Equations (33)-(36) in a more compact form if we pass to the functions , with . To this end, introducing the new functions into (33)-(36), we obtain the following set of equations:
Let us re-express Equations (52)-(55) to the function of argument , i.e. , . In these conditions, we have the following set of equations for the functions :
Let us pass from the set of first-order differential Equation (57)-(60) to a set of second-order differential equations. By doing so, differentiating Equation (57) in S and introducing into the result the expressions of the function and , we obtain a second-order differential equation for the function as follows:
Similarly differentiating Equation (60) and introducing into the result the expression of and the expression of its derivative, we obtain the second-order differential equation for the function :
Summing (63)-(64) and setting , we obtain the following second-order differential equations of the function :
Under the condition with , the solution of Equation (65) is:
Substracting Equations (57) and (60) and taking into account (66), we get:
From expressions (66) and (68), we define the functions and as follows:
The same operating on Equations (58) and (59) leads to the following expressions:
It follows that, the concrete form of the functions is:
The exact functions are regular and bounded.
In the sequel, we shall use the concrete form of the functions in order to determine the total charge and total spin.
The general form of the spinor current vector is:
From (78), we deduce the components of the spinor current vector as follows:
In this study, the configuration is assumed static. Thus, only the component is non null. With this assumption, we obtain: ; ; . In addition to this, we have .
The component determines the charge density of the spinor field and scalar field in interaction:
The charge density is continuous localized function.
The total charge of the interaction system of spinor and scalar fields is:
In the integral (84), denotes the center of the fields configurations.
The spin tensor of the spinor field reads:
The analytical expression defines the spatial density of the spin vector. In virtue of (85), we have:
It follows from (86) that:
The chronometricaly invariant spatial density of the spinor vector is:
The projection of the spin vector on the axis is given by the following expression:
According to expressions (84) and (93), the total charge Q and the total spin are finite.
In the following section, we shall examine the Einstein’s, spinor and scalar fields equations by choosing a concrete form of the arbitrary function .
In this section, we shall analyze the general results obtained in the previous section for concrete form of the function . The function is chosen under the form:
THus, the expression of the function becomes:
By substituting (95) into (51), without losing the generality we can consider massless spinor and scalar fields, according to the theory of Heisenberg , we obtain:
The function is a continuous and limited function.
Introducing the relations (94) and (96) into (22), we define the energy density of the interaction of spinor and scalar fields:
where the function is defined by the following expression:
The distribution of the energy density per unit invariant volume is given by the expression:
From (98) and (100), the energy density and the distribution of the energy density per unit invariant volume are continuous and localized functions. Moreover, the total energy is limited. The solutions of Equation (31) are soliton-like when the arbitrary function is chosen under the form .
Let us emphasize that when the interaction lagangian is null i.e the invariant function is null in the relation (1), the obtained solutions are exact and regular. Nevertheless, the energy density is not localized. The lagrangian interaction is very important in order to obtain the regular solutions with localized energy density.
In order to clarify the role of the nonlinear terms of the interaction of spinor, scalar and gravitational fields in the nonlinear fields equations in the formation of regular localized soliton-like solutions, we obtained the soslutions to the fields equations in the case where . Let us note that the obtained solutions are exact regular but their energy density is not bounded. In the present case soliton-like solutions are absent.
5. Concluding Remarks
Taking into account the proper gravitational field of elementary particles, we have obtained and examined the exact general solutions of interaction scalar, spinor and gravitational field equations. These solutions describe the configurations of scalar and spinor fields with a localized energy density , finite total energy E. In addition, the total charge and spin are limited.
The solitons-like solutions with limited total spin and charge exist in the case where the arbitrary function is under the general form .
 Adomou, A. and Shikin, G.N. (1998) Exact Self-Consistent Plane-Symmetric Solutions of the Equations of Two Interacting Fields (Spinor Field And Scalar Field). Russian Physics Journal, 41, 672–678.
 Adomou, A., Edou, J. and Massou, S. (2019) plane Symmetric Solutions to the Nonlinear Spinor Field Equations in General Relativity Theory. Journal of Modern Physics, 10, 1222-1234.
 Adomou, A. and Shikin, G.N. (1997) Intearcting Spinor and Scalar Fields: Exact Self-Consistent Plane-Symmetric Solutions. 15th International Conference of General Relativity and gravitation, 16-21 December 1997, Pune.
 Adanhoumè, A., Adomou, A., Codo, F.P. and Hounkonnou, M.N. (2012) Nonlinear Spinor Field Equations in Gravitaional Theory: Spherical Symmetric Soliton-Like Solutions. Journal of Modern Physics, 3, 935-942.
 Adomou, A., Edou, J. and Massou, S. (2019) Soliton-Like Spherical Symmetric Solutions of the Nonlinear Spinor Field Equations Depending on the Invariant Function IP = P2 in the General Relativity Theory. Journal of Applied Mathematics and Physics, 7, 2818-2835.
 Massou, S., Adomou, A. and Edou, J. (2019) Soliton-Like Spherical Symmetric Solutions of the Nonlinear Spinor Field Equations in General Relativity. International Journal of Applied Mathematics and Theoretical Physics, 5, 118-128.
 Adomou, A., Edou, J., Hontinfinde S.I.V. and Massou, S. (2020) Exact Soliton-Like Spherical Symmetric Solutions of Heisenberg-Ivanenko Nonlinear Spinor Field Equation in Gravitational Theory. Journal of Applied Mathematics and Physics, 8, 1236-1254.
 Edou, J., Adomou, A., Hontinfinde, S.I.V. and Massou, S. (2020) Analytical Soliton-Like Solutions to Nonlinear DIRAC Equations of Spinor Field in Spherical Symmetric Metric. Journal of High Energy Physics, Gravitation and Cosmology, 6, 623-639.