1. Introduction
The recently proposed classification of Shapovalov wave spacetimes [
1] provides additional mathematical tools for the construction of exact integrable models of gravitational waves, including primordial gravitational waves in Bianchi universes [
2,
3] and other exact models of plane gravitational waves [
4,
5]. Recent advances in gravitational wave astronomy in the detection of gravitational waves and the resulting astrophysical information [
6,
7,
8] have also increased interest in the mathematical aspects of gravitational wave research.
The relevance of these lines of research is also related to the importance of the discovery of primordial gravitational waves for constructing a theory of the early universe and the possible confirmation of the stage of inflation predicted by some theories. The interest in gravitational wave research is also due to possible secondary physical effects in a gravitational wave, such as the formation of black holes [
9,
10], the capture of astrophysical objects by a gravitational wave [
11], gravitational wave lensing, and other effects.
The Shapovalov wave spacetimes allow for the presence of special symmetries of the Hamilton–Jacobi equation, leading to the existence of privileged coordinate systems, where the equations of motion of test particles and the eikonal equation for radiation in the Hamilton–Jacobi formalism allow for the exact integration by the method of separation of variables. Moreover, among the separable variables on which the spacetime metric depends (non-ignored variables), there are wave variables along which the spacetime interval vanishes. The gravitational wave velocity, according to the gravitational wave detection data from neutron star mergers [
12], equals the velocity of light. Therefore, the existence of the possibility of separating null variables indicates the wave nature of such spacetime models, and the use of wave variables in describing gravitational waves is experimentally justified.
Shapovalov wave spaces in the four-dimensional cases include three main types according to the number of commuting killing vectors they allow and, accordingly, in privileged coordinate systems, the number of killing vectors determines the number of non-ignorable variables on which the spacetime metric depends. For privileged coordinate systems, Vladimir Shapovalov obtained the structure of the space metric [
13,
14,
15], which leads to the separation of variables in the Hamilton–Jacobi equation, which allows one to construct exact integrable models of gravitational waves in various theories of gravity in which test particles move along geodesic curves of spacetime.
Currently, there is interest in studying the properties of modified theories of gravity that could provide corrections to Einstein’s theory of gravity in the early stages of the universe (quantum corrections, quadratic theories,
-theories of gravity, theories with a scalar field, etc.) and could provide a theoretical description of the phenomena of “inflation”, “dark matter”, and “dark energy” [
16,
17,
18,
19]. Therefore, the use of additional mathematical tools for the comparative analysis of exact models of gravitational waves in different theories of gravity also enables the selection of the most realistic theories and models.
In this direction, one of the first models of modified theories of this type was the Brans–Dicke scalar-tensor theory of gravity [
20]. Scalar-tensor theories have taken a firm place in theoretical research on gravity and cosmology [
21], determined by attempts to use an additional scalar field to describe various possible scenarios for the dynamics of the universe and other theoretical constructions [
22,
23,
24,
25].
In this work, Shapovalov wave spaces of type I are considered in the scalar-tensor theory of Brans–Dicke gravity. An additional scalar field in the theory complicates the field equations but also provides additional possibilities. The resulting equation for a Klein–Gordon equation-type scalar field also requires further investigation. As Shapovalov [
13,
14] showed, the scalar equation can be integrated by separating the variables in the same privileged coordinate systems as the Hamilton–Jacobi equation. Therefore, the use of Shapovalov wave spaces also provides additional opportunities for the exact solution of the scalar equation [
26,
27,
28,
29]. Recently, new results have also been obtained in the study of the symmetries of the Klein–Gordon–Fock equations. [
30,
31,
32,
33].
Type I Shapovalov metrics are the most general gravitational wave models of Shapovalov spacetimes since they depend on three variables in the privileged coordinate system, including the wave variable. Note that metrics of this type lead to degeneracy when using Einstein’s vacuum equations since the number of non-ignorable variables in the metric decreases as the vacuum field equations are solved, and the space becomes either a Shapovalov wave space of type II with two non-ignorable variables in the metric, or a space of type III with one wave variable in the metric. The preservation of type I metrics is possible if additional fields and matter exist in the models under consideration. The existence of non-degenerate type I Shapovalov models in the Brans–Dicke theory of gravity and the possibility of detecting such gravitational waves in observations would therefore be additional evidence for the scalar-tensor theories of gravity.
The authors investigate exact wave models of spacetime to find out the differences that may arise in these gravitational wave models for Einstein’s theory and for modified theories of gravity. Shapovalov type I wave spaces have internal symmetries that relate them to mathematical models of gravitational waves, but still, these models of gravitational waves do not lead to the appearance of exact solutions of Einstein’s vacuum equations, unlike the type II and III models. Additionally, in the modified Brans–Dicke theory of gravity, exact wave solutions occur for certain spaces, as we show below. This allows us to observe differences in the detection of gravitational waves by existing and future detectors. This is also important for analyzing the stochastic gravitational wave noise detected by existing gravitational wave detectors. In addition, primordial gravitational waves of this type could leave an imprint on the cosmic microwave background. Detecting “traces” of such gravitational waves would argue that the scalar-tensor gravity theories are more realistic than Einstein’s theory. The possibility of the exact integration of the equations of motion of the test particles and the equations for geodesic deviation in these models of gravitational waves allows gravitational wave detectors to more efficiently analyze and select such disturbances in signals with high noise.
2. Gravitational-Wave Spacetimes of Shapovalov
Shapovalov wave spaces allow for the exact integration of the equations of motion of test particles in the Hamilton–Jacobi formalism
with the separation of non-ignorable wave variables (the space metric depends on these variables) along which the spacetime interval vanishes [
1]. Here,
S is the action function of the test particle,
m is the mass of the particle, and
c is the speed of light. Let us also choose a system of units in which the speed of light is equal to one.
Type I Shapovalov spacetimes admit a killing vector contained in the so-called “complete set” of vectors and second-rank killing tensors that determine integrals of motion linearly and quadratically in the momenta. In a privileged coordinate system where the Equation (
1) allows for the complete separation of variables, the type I Shapovalov spacetime metric depends on three variables, including the null wave variable denoted by
. Another null variable
is ignored (cyclic) and is not included in the metric.
Consider a model with a type I Shapovalov wave space metric, which can be written in a privileged coordinate system in the following form:
where
The scalar curvature
R takes the following form:
where the top prime means the ordinary derivative with respect to the variable on which the function depends.
The nonzero components of the Riemann curvature tensor have the following form:
Moreover, all nonzero components of Weyl’s conformal curvature tensor are proportional to the scalar curvature R, and when the scalar curvature vanishes, they also vanish, leading to a conformally flat spacetime.
Note that the ability to exactly integrate the Hamilton–Jacobi equations for Shapovalov wave spaces makes it possible to obtain the trajectories of test particles, as well as to find exact solutions for the geodesic deviation equation and the exact form of the tidal accelerations in these spaces. These possibilities allow one to determine all of the physical effects in gravitational waves.
3. Brans–Dicke Scalar-Tensor Theory of Gravity
The Lagrange function for the Brans–Dicke theory of gravity (BDT) can be written in the following form [
21]:
where
is a scalar field,
is a constant parameter of the theory, and
is a cosmological constant.
The field equations of the Brans–Dicke scalar-tensor gravity theory with cosmological constant
can be written as follows:
Here, the comma denotes the ordinary partial derivative, and the semicolon denotes the covariant derivative.
To obtain exact gravitational wave solutions, consider the vacuum model
. It is convenient to pass to the form of the scalar field
; then, the BDT equations will take the following more compact form:
Note that the higher derivatives of the scalar field
can be found in Equation (
10) and substituted into the scalar Equation (
11) to reduce the order of this equation.
The equations of the scalar-tensor theory of gravity of Brans–Dicke for the considered metric lead to the following form of field equations:
Equation (
19) arises as a consequence of scalar Equation (
9) when the higher derivatives of the scalar field
from the other field equations of the theory are substituted into it.
3.1. Scalar Field of the Form
and
Assume that the scalar field depends on all variables, including the null variable , on which the spacetime metric does not depend.
In this case, the system of BDT equations is quite cumbersome, but it allows one to find the conditions for the compatibility of field equations, which leads to relations of the following form:
The requirement
imposes strict conditions on the dependence of the scalar field on the variables
and
, so the case of
in (
20) leads to the degeneration of the space into a flat one for any
.
Thus, the scalar field for
can only have the following form:
In this case, the system of BDT equations takes the following form:
Consider Equation (
23). For
, the scalar field is in the form
however, as a result of substitution into other equations of the system and the separation of the coefficients of
, it turns out that
, and the scalar field does not depend on
. Consequently, a solution that satisfies the stated requirements is possible only for
.
The integration of the Equation (
23) and the separation of the coefficients for
for
leads to the scalar field
,
. Substitution into other equations of the system allows one to find all of the unknown functions.
The final solution for the case
is as follows:
where
,
,
,
, and
are constants.
In the resulting solution, we can set
by scaling transformations; then,
The cosmological constant
remains arbitrary, the Riemann curvature tensor
and the Weyl conformal curvature tensor
do not vanish, and for the scalar curvature
R we obtain:
Thus, a nontrivial solution for Shapovalov type I wave spaces in the Brans–Dicke theory is obtained with a scalar field depending on null variables, with a cosmological constant and an exponential conformal factor depending on the wave variable. Note that a vacuum solution of this type does not arise in general relativity [
1].
3.2. Scalar Field of the Form
Consider the case where the scalar field
does not depend on the “ignored” null variable
, on which the metric also does not depend. The system of BDT Equations (
10) and (
11), resolved in terms of higher derivatives, has the following form:
The study of the compatibility of the subsystem of Equations (
34)–(
36) leads to relations that are further used to classify possible solutions:
The compatibility Equation (
43) lead to the following two possible cases:
- (1)
;
- (2)
, .
In the following, we will consider these cases separately.
3.2.1. Scalar Field of the Form and
From the compatibility conditions, we obtain
, and now we can integrate Equations (
34)–(
36). In this case, two cases are possible depending on the values of the constant
:
- (A)
,
. Substitution into Equation (
41) allows us to separate the coefficients at
. As a result, we obtain
, which, in view of (
6) and (
43), means the degenerate case since spacetime becomes flat in this case.
- (B)
,
,
, where the function
must satisfy the following equations:
Let us consider the solution when the scalar field
admits the separation of variables:
From here, we obtain the relations:
Substituting the relations (
49) into Equation (
46) and separating the variables, we obtain
From this, we obtain relations relating the functions
and
of the following form:
where
a,
b, and
c are constants.
Since the functions
enter the metric only as a difference, and the constants
add only a constant term to the scalar field
, then the constants
can be chosen in the simplest way. As a result, we have an intermediate result for the scalar field
:
Having this, excluding the second derivatives from the Equations (
44), (
45), and (
47), we obtain the following relation:
The case
leads to a flat solution, and it is easy to see this, taking into account the relations (
6) and (
41); thus,
The sum of Equations (
44) and (
45) after substituting (
53) and (
55) allows us to separate variables, and we obtain
the solution is the following:
Substituting this solution into Equations (
44) and (
45) leads to a simple condition
; furthermore, the scaling constants
can be converted to the simplest form, and as a result, the final solution has the following form:
For the resulting exact solution (
59) of the equations of the Brans–Dicke theory, the Riemann curvature tensor, the scalar curvature, and the Weyl conformal curvature tensor do not vanish:
3.2.2. Scalar Field of the Form
Let us consider the case where the scalar field
depends only on the non-ignorable wave variable
(
), on which the wave metric itself depends. Equations (
37) and (
41) imply that the cosmological constant
and the scalar curvature
R vanish in this case:
Then, one can immediately separate the variables in the Equation (
41), which allows you to find differential equations for the functions
:
where
a,
b, and
c are constants.
The system of field equations is left with only one Equation (
34) connecting the function
(scale factor of the metric) and the scalar field
:
Equation (
64) obviously admits the solutions for the scalar field
for a given conformal factor
of the metric.
In the case under consideration, the Weil tensor
vanishes, and only two nonzero components of the Riemann curvature tensor remain
:
Thus, this solution is the case of conformally flat space, and the metric and components of the Riemann curvature tensor depend on the wave variable
. The scalar field
and the conformal factor of the metric
remain arbitrary functions of the wave variable
but are related by the Equation (
64).
Note that Equation (
64) admits a de Sitter-type solution. Consider the scale factor of the metric in the following form:
Then, Equation (
64) takes the following form:
Equation (
68) for scalar field
for metric (
2) with scale factor (
67) has three solutions depending on the values of the constant
:
- (1)
- (2)
- (3)
. Then,
where
is a constant parameter.
Thus, we have obtained a set of exact solutions to the equations of the scalar-tensor Brans–Dicke theory in vacuum for the type I Shapovalov wave space. Recall that there are no exact solutions of this type for the case of the vacuum Einstein equations [
1] since Shapovalov spaces of type I degenerate in the case of vacuum Einstein equations.
Therefore, scalar-tensor theories of gravity can provide such exact models of gravitational waves that do not appear in general relativity.
4. Equations of Motion and Trajectories of Test Particles
Let us consider the equation for the motion of test particles in a gravitational field in the Hamilton–Jacobi formalism (the speed of light is chosen to be unity):
The complete integral of the Hamilton–Jacobi Equation (
1) for the Shapovalov spacetimes in the privileged coordinate system can be represented in a separated form:
moreover, for the ignored variable
, on which the metric in the privileged coordinate system does not depend, the function
can be reduced by admissible coordinate transformations to the form
, where
k is a constant.
From Equation (
1), separating the variable
, we obtain:
Finally, separating the variables
and
we obtain:
Thus, for the complete integral of the particle action function
we obtain the following expression:
where
k,
p, and
q are constant parameters.
The equations of the trajectories of the test particles for the considered Shapovalov spaces of type I in the Hamilton–Jacobi formalism take the following form:
where
,
, and
are new independent constant parameters of the test particle motion determined by the initial conditions. Thus, the motion of test particles in the considered spacetime models is determined by a set of constant parameters
k,
p,
q,
,
, and
, given by the initial conditions.
Note that the proper time of the test particle
using the obtained trajectory equations can be written in the following form in terms of the coordinate variables on the trajectory:
Here, we set the additive constant equal to zero by choosing the origin. Thus, proper time, as expected, is a linear combination of the null variables and .
If we substitute (
79) into Equation (
76), then on the trajectories of test particles, we obtain:
Thus, Equations (
79) and (
80), together with Equations (
77) and (
78), define the coordinates of the test particle
on the trajectory as a function of the proper time of the particle
.